Attenuation of clinical and immunological outcomes during SARS‐CoV‐2 infection by ivermectin

Abstract The devastating pandemic due to SARS‐CoV‐2 and the emergence of antigenic variants that jeopardize the efficacy of current vaccines create an urgent need for a comprehensive understanding of the pathophysiology of COVID‐19, including the contribution of inflammation to disease. It also warrants for the search of immunomodulatory drugs that could improve disease outcome. Here, we show that standard doses of ivermectin (IVM), an anti‐parasitic drug with potential immunomodulatory activities through the cholinergic anti‐inflammatory pathway, prevent clinical deterioration, reduce olfactory deficit, and limit the inflammation of the upper and lower respiratory tracts in SARS‐CoV‐2‐infected hamsters. Whereas it has no effect on viral load in the airways of infected animals, transcriptomic analyses of infected lungs reveal that IVM dampens type I interferon responses and modulates several other inflammatory pathways. In particular, IVM dramatically reduces the Il‐6/Il‐10 ratio in lung tissue and promotes macrophage M2 polarization, which might account for the more favorable clinical presentation of IVM‐treated animals. Altogether, this study supports the use of immunomodulatory drugs such as IVM, to improve the clinical condition of SARS‐CoV‐2‐infected patients.


Introduction
Coronaviruses cause respiratory disease in a wide variety of hosts. During the ongoing pandemic of SARS-CoV-2 causing coronavirus disease 19 , clinical signs other than respiratory symptoms have been linked to infection, frequently associated with neurological symptoms such as anosmia and ageusia. These features have been related to an over-responsiveness of patients' immune system to SARS-CoV-2 (Bhaskar et al, 2020;Han et al, 2020;Qiu et al, 2020). Consequently, there is an urgent need to understand the hallmarks of this over-responsiveness and to find novel therapeutics or repurpose drugs to improve the clinical condition of COVID-19 patients (Batalha et al, 2021).
Ivermectin (IVM), a macrocyclic lactone, is a commercially available anti-parasitic drug which prevents infection by a wide range of endo-and ectoparasites (Sajid et al, 2006;Heidary & Gharebaghi, 2020). IVM is an efficient positive allosteric modulator of the α-7 nicotinic acetylcholine receptor (nAChR) (Krause et al, 1998) and of several ligand-gated ion channels, including the muscle receptor for glutamate (GluCl) in worms (Hibbs & Gouaux, 2011). Furthermore, IVM has been shown to exert an immunomodulatory effect in humans and animals (Sajid et al, 2006;Heidary & Gharebaghi, 2020) under conditions that are known to involve the α-7 nAChR (Pavlov & Tracey, 2012), even though its underlying mechanisms are yet to be established (Laing et al, 2017). A direct or indirect interaction of SARS-CoV-2 with nAChR has also been hypothesized, in particular because of sequence homologies between SARS-CoV-2 spike proteins and nAChR ligands such as snake venom toxins . IVM has been shown to be active beyond its anti-parasitic activity in a wide variety of pathologies, including cancer, allergy, and viral infections (Laing et al, 2017). Recently, IVM has been reported to reduce viral load and improve the clinical status of mice infected by an animal coronavirus, the mouse hepatitis virus (MHV) (Ar evalo et al, 2021). In vitro inhibition of SARS-CoV-2 replication by IVM in Vero/hSLAM cells has also been reported (Caly et al, 2020), albeit at much higher concentrations (50-to 100-fold) than those clinically attainable in humans (150-400 µg/kg) (Guzzo et al, 2002;Bray et al, 2020;Chaccour et al, 2020).
The aim of this study is to investigate the impact of IVM on the pathogenesis of COVID-19, in a SARS-CoV-2 infection model, the golden Syrian hamster. This species is naturally permissive to this virus and the most reliable and affordable animal model for COVID-19 (Chan et al, 2020;Muñoz-Fontela et al, 2020). Moreover, it was recently used to demonstrate the importance of lowering the inflammation with intranasal administration of type I IFN to prevent disease progression (Hoagland et al, 2020). Male and female adult golden Syrian hamsters were intranasally inoculated with 6 × 10 4 PFU of SARS-CoV-2 [BetaCoV/France/IDF00372/2020]. This inoculum size was selected as it invariably causes symptomatic infection in golden Syrian hamster, with a high incidence of anosmia and high viral loads in the upper and lower respiratory tracts within 4 days post-inoculation (dpi) (de Melo et al, 2021). At the time of infection, animals received a single subcutaneous injection of IVM at the anti-parasitic dose of 400 µg/kg, commonly used in human clinical setting, and were monitored over 4 days.
Here, we show that the modulation of the host's inflammatory response using IVM as a repurposed drug strongly diminished the clinical score and severity of the disease (including anosmia) observed in these animals, although it has no impact on viral load. IVM-treated animals presented a strong modulation in several signaling pathways, including a significant reduction of the type I and III interferon response and of the Il-6/Il-10 ratio, along with the presence of M2 macrophages in the lung. These effects were mostly compartmentalized and sex-dependent, and treated infected females exhibited better clinical outcomes.

COVID-19 clinical outcome is attenuated by ivermectin
In order to study the effects of IVM chemical therapy on clinical outcome, we assessed body weight, clinical score, and olfactory performance daily for 4 days post-infection. The golden hamster model reproduces a moderate-to-severe COVID-19 in humans, with a clinical phase lasting 5-6 days post-infection followed by a complete recovery by 2 weeks post-infection, with no death occurring. The viral titers and viral RNA loads in the airways of infected animals are elevated after 2-4 days post-infection, but they drastically drop by 7 days post-infection (Chan et al, 2020;Sia et al, 2020;de Melo et al, 2021).
The IVM-treated infected animals, of both sexes, showed a decrease in body weight similar to that observed in saline-treated infected hamsters ( Fig EV1A). However, IVM-treated infected animals exhibited a significant reduction in the severity of the clinical score, in a sex-dependent manner: Infected males presented a significantly reduced clinical score, and it fully returned to normal in infected females ( Fig 1A). Remarkably, IVM treatment reduced the olfactory deficit in infected animals (Figs 1B and EV1B-D): 66.7% (12/18) of the saline-treated infected hamsters presented with hyposmia/anosmia, whereas only 22.2% (4/18) of IVM-treated infected hamsters presented with signs of olfactory dysfunction . The olfactory performance was also influenced by sex: 83.3% (10/12) of the saline-treated infected males presented with hyposmia/anosmia, against only 33.3% (4/12) of IVM-treated infected males . In females, while 33.3% (2/6) of saline-treated infected animals presented with hyposmia/anosmia, no olfactory deficit was observed in IVMtreated infected females (0/6; Figs 1B and EV1B-D).
Since males presented a higher incidence of anosmia/hyposmia, we subsequently performed a dose-response curve to test the effect of IVM on the clinical presentation and olfactory functions of infected males: Lower doses of IVM (100 or 200 µg/kg) elicited similar clinical outcomes as the anti-parasitic dose of 400 µg/kg ( Fig  EV2). As expected, no signs of olfactory deficit were observed in mock-infected hamsters (Fig EV1B-D).

Ivermectin treatment does not influence SARS-CoV-2 load in the respiratory tract of infected hamsters
To evaluate the effect of IVM treatment on the viral load in the respiratory tract, we tested the nasal turbinates and lungs of infected hamsters using both classical RT-qPCR ( Fig 1C) and the highly sensitive technique of digital droplet PCR (Suo et al, 2020;Fig 1D, Appendix Fig S1). Surprisingly, the viral RNA load in the respiratory tract remained unaffected by IVM treatment in both samples in both sexes. Furthermore, IVM treatment did not influence viral replication rate, as evaluated by the ratio between structural and nonstructural gene transcription (Fig 1D, Appendix Fig S1). Finally, IVM treatment did not alter infectious viral titers in the lungs (Fig 1E). These results illustrate that no antiviral activity of IVM is detected in vivo at the standard anti-parasitic dose of 400 µg/kg, in contrast to a previous report suggesting that IVM inhibits the replication of SARS-CoV-2 in vitro, albeit used at far higher concentrations (Caly et al, 2020). Therefore, the action of IVM on COVID-19 clinical signs of infected golden hamsters does not result from a decrease in viral replication.
◀ Figure 1. Clinical presentation, olfaction test, viral load and immune profile in the nasal turbinates of SARS-CoV-2-infected hamsters with and without ivermectin treatment.
A Clinical signs in infected hamsters. The clinical score is based on a cumulative 0-4 scale: ruffled fur; slow movements; apathy; and absence of exploration activity.
Symbols indicate the median AE interquartile range. B Olfactory performance in infected hamsters. The olfaction test is based on the hidden (buried) food finding test. Curves represent the percentage of animals that did not find the buried food. Food finding assays were performed at 3 days post-infection (dpi). Data were obtained from three independent experiments for males and two independent experiments for females. C Viral load in the nasal turbinates and in the lungs at 4 dpi. D Ratio between the CPD (copy per droplets, normalized to γ-actin and Hprt reference gene relative expression) of structural [N, nucleocapsid] and non-structural [IP4: RdRp, RNA-dependent RNA polymerase] viral gene expression determined by digital droplet PCR (ddPCR) in the nasal turbinates and in the lungs at 4 dpi. E Infectious viral titer in the lung at 4 dpi expressed as plaque-forming units (PFU)/g of tissue. F Cytokine and chemokine transcripts in the nasal turbinates at 4 dpi in male and female SARS-CoV-2-infected hamsters, treated with saline or with 400 µg/kg ivermectin. Ivermectin therapy modulates local immune responses in infected hamsters' nasal turbinates Anosmia is a typical symptom of COVID-19 in humans, with some sex-dependent differences Qiu et al, 2020;Xydakis et al, 2020). Inflammation in the nasal cavity, following olfactory sensory neurons infection and deciliation, has been shown to be an underlying factor for smell loss during SARS-CoV-2 infection (de Melo et al, 2021), and the chemokine Cxcl10 could be directly implicated due to its neurotoxic potential (Oliviero et al, 2020). We therefore tested a possible modulation by IVM of the local inflammatory response in hamsters and in particularly in the nasal turbinates, the primary target tissue of SARS-CoV-2 infection (de Melo et al, 2021), that could correlate its effect on the olfactory score. To this aim, a panel of cytokines (Il-6, Il-10, Il-1β, Tnf-α, Ifn-β, Ifn-γ, and Ifn-λ) and chemokines (Cxcl10 and Ccl5) were used to assess the impact of IVM treatment on immune responses in the nasal turbinates of SARS-CoV-2-infected hamsters at 4 dpi. Upon treatment with IVM, females presented a significant downregulation of Il-6, Il-10, and Tnf-α, which are key inflammatory mediators of prognostic value in COVID-19 patients (McElvaney et al, 2020a), and of Cxcl10 (Fig 1F), in line with their better olfactory performance observed in the food finding tests (Fig 1B). The differences between sex groups are illustrated by the increase in three pro-inflammatory mediators (Ifn-γ, Ifn-λ, and Ccl5) only in males ( Fig 1F). No difference for the Il-6/Il-10 ratio was observed in the nasal turbinates.

Lung immunometabolism is affected by SARS-CoV-2 infection and modulated by ivermectin
IVM attenuates lung pathology and inflammation pathways, including cholinergic synapse-related genes In order to further study the mode of action of IVM on clinical signs, we performed at 4 dpi a comparative agnostic transcriptomic approach using RNA-seq in the lower respiratory tract in hamsters treated or not with IVM. In SARS-CoV-2-infected lungs, male and female hamsters exhibited an overall similar pattern, although there were slightly more KEGG pathways modulated by infection in males (Fig 2A). A high number of common dysregulated inflammatory and metabolic pathways in SARS-CoV-2-infected males and females ( Fig  EV3) were similar to those observed in lung cells from human COVID-19 patients, such as "cytokine-cytokine receptor interaction", "TNF signaling pathway", and "insulin resistance" (Dey et al, 2020;Islam et al, 2020) including many genes of the IFN response. The transcriptomic profiles observed in the lungs from IVM-treated hamsters compared to non-treated infected hamsters presented a striking sex difference: Lungs of IVM-treated females presented 1,206 downregulated genes and 1,428 upregulated genes, whereas only 36 downregulated and 51 upregulated genes were detected in the lungs of IVM-treated males. This sex difference is also illustrated by KEGG and GO enrichments representations ( Fig EV4). Several KEGG pathways were significantly regulated in IVMtreated females: "TNF signaling pathway" and "cholinergic synapse", in line with the activation of vagus nerve cholinergic antiinflammatory pathway (CAP) (Pavlov & Tracey, 2012;De Virgiliis & Di Giovanni, 2020), and "platelet activation", related to prevention of thrombosis (Zhang et al, 2020b;Chen & Pan, 2021) (Figs 2A and EV4A). We also observed a modulation of the "Wnt signaling pathway" KEGG pathway in the IVM treatment group, which has been shown to correlate with lung homeostasis and regeneration (Raslan & Yoon, 2020), obesity, type-2 diabetes, and cancer (Melotti et al, 2014;Aamir et al, 2020) (Figs 2A and EV4A). Further modulated pathways in IVM-treated female hamsters were "renin secretion", "vascular smooth muscle contraction", and "adrenergic signaling in cardiomyocytes" which could be related to the prevention of vasoconstriction and therefore attenuation of lung pathology ( Fig EV4A) (Potus et al, 2020;Vaduganathan et al, 2020). Moreover, carbohydrate metabolism and insulin resistance were frequent terms observed in RNA-seq analyses of the lung of female hamsters during SARS-CoV-2 infection, linked to hyperglycemia, metabolic syndrome, and impairments in the immune system, similarly to what is observed in COVID-19 in humans (Pavlov & Tracey, 2012;Gianchandani et al, 2020). In IVM-treated hamsters, specifically in females, the term "insulin resistance" was no longer observed in the KEGG enrichment, and the insulin secretion pathway was downregulated. Several other related pathways linking hyperglycemia, inflammation, and lung pathology were downregulated, such as Rap1, AGE-RAGE, and cGMP-PKG signaling pathways (Figs EV4 and EV5) (Oczypok et al, 2017;Pei et al, 2020;Isidori et al, 2021). In IVM-treated males, significant modulated pathways determined by KEGG and GO enrichments were related to tight junctions, muscular and epithelial cells ( Fig EV4).
IVM exerts a strong regulation of type I and III IFN-related genes and lowers down the Il-6/Il-10 ratio in the lungs of SARS-CoV-2infected hamsters To analyze in more details the modulated genes in IVM-treated animals, we selected representative target genes from these relevant ◀ Figure 2. Transcriptomic profile in the lung of SARS-CoV-2-infected hamsters with and without ivermectin treatment at 4 days post-infection.

A Heatmaps showing the differentially expressed genes according to the selected KEGG pathways calculated in comparison with mock-infected hamsters. * indicates
Benjamini-Hochberg-adjusted P-value < 0.05 in the comparison between saline and ivermectin within the same sex. Color gradient represents the transcription log 2 fold change comparing infected and mock-infected. Complete analyses are listed in Dataset EV1. B Validation targets in the lung at 4 dpi. Horizontal lines indicate medians. The P value is indicated in bold when significant at a 0.05 threshold. Mann-Whitney test. pathways that were highly regulated in IVM-treated females and we compared by RT-qPCR their respective transcription levels in the lungs of the different groups of animals. Among these genes, several were also significantly modulated between IVM-and saline-treated infected males including Tnfaip3, Sfrp4, Epha2, Gnai1, Hgf, and Fos.
Others presented a similar regulation in males compared to females (from KEGG and GO enrichments) although their modulation was not invariably significant: Casp3, Plcb1, Chrna7, Chrnb4, Adra1d, Grin2d/Nmdar2d, Grid1, Gabrr1, Pik3r5, Igf1r, Wnt11, Wnt3a, Il-2, Il-2ra, Prkg1, Krt4, Creb5, and Ager. The modulation of these targets, together with other genes from relevant inflammatory mediators taken from the literature (Boudewijns et al, 2020;Hoagland et al, 2020) (Il-6, Il-10, Il-1β, Tnf-α, Ifn-β, Ifn-λ, Ifn-γ, Tgf-β, Cxcl10, Ccl5, and Mx2), was confirmed by RT-qPCR (Figs 2B and EV5). IVM limited the expression of Ifn-β (males and females), Ifn-λ (females), and IFN-stimulated gene Mx2 (males and females) in the lung of infected and IVM-treated hamsters compared to infected and saline-treated hamsters. In contrast, saline-treated hamsters presented an increased expression of Ifn-β, Ifn-λ, and Mx2 compared to non-infected hamsters (males and females). This is expected as type I and III IFN signaling pathways have already been shown to correlate with lung pathology severity in SARS-CoV-2-infected hamsters, possibly resulting from a STAT2-dependent response (Boudewijns et al, 2020). In contrast, type I and III IFNs are differently expressed in the nasal turbinates, where Ifn-λ is upregulated in treated males (Fig 1F). This difference between upper and lower airways may be explained by two factors: (i) the specificity of Ifn-λ in the nasal cavity's antiviral response (Okabayashi et al, 2011;Klinkhammer et al, 2018), whereas type I and III IFNs perform redundant actions in the lung (Stanifer et al, 2020), and (ii), following studies with influenza A virus, type I and III IFNs have a sequential activation, where Ifn-λ comes at first to fight infection, and type I IFNs come later, with pro-inflammatory and reinforced antiviral effects (Galani et al, 2017;Stanifer et al, 2020).
Interestingly, some of these genes are not only involved in the previously mentioned pathways but also play a major role in other molecular functions that may be critical in the pathogenesis of SARS-CoV-2 infection and the efficacy of the IVM treatment in hamsters. Among them, the significant overexpression of Il-10 was a common feature of IVM-treated males and females (Fig 2B). This effect may be related to a modulation of the inflammatory response in the lung (downregulation of Tnf-α in males, Tnfaip3 in both sexes, and Il-6 in females) associated with reduced clinical signs. Tnf-α reduction and better clinical presentation are common features caused by the IVM treatment in another model of coronavirus infection, where mice were infected with MHV (Ar evalo et al, 2021). Additionally, differently from the nasal turbinates, the Il-6/Il-10 ratio in the lung of IVM-treated hamsters was significantly lower than in non-treated animals (Fig 2B), which may relate to their comparatively better clinical presentation. Of note, lower plasmatic Il-6/Il-10 ratios were detected in hospitalized COVID-19 patients who did not require intensive care (McElvaney et al, 2020a;McElvaney et al, 2020b).

IVM increases the infiltration of monocytes/macrophages and promote M2 polarization in the lung of SARS-CoV-2-infected hamsters
To assess directly lung pathology in infected animals, we performed histopathological analyses. The lungs of SARS-CoV-2-infected and saline-treated hamsters exhibited substantial lesions with focal areas of edema, congestion, microhemorrhages, associated with mononuclear cells, hyaline membranes, alveolar walls thickening (Fig 3A), in line with previous reports (Chan et al, 2020;Sia et al, 2020). In contrast, the lungs of SARS-CoV-2-infected and IVM-treated hamsters exhibited with reduced degrees of edema and congestion, yet with greater amounts of mononuclear cells in the alveolar spaces ( Fig 3A). SARS-CoV-2 infection caused an increase in the lung weight/body weight ratio; however, this increase was lowered in IVM-treated as compared to saline-treated infected animals ( Fig 3E).
The better clinical presentation observed in IVM-treated animals could be correlated with the activation of the cholinergic antiinflammatory pathway (CAP) under vagus nerve control (Pavlov & Tracey, 2012). Indeed, the interaction of IVM with the α-7 nAChR in macrophages may account for the immune response modulation in IVM-treated hamsters (Gahring et al, 2017;Zhao et al, 2019), eliciting Il-10 production, promoting M2 polarization (dampening inflammation and triggering tissue-repair), and counterbalancing M1 action (pro-inflammatory, associated with tissue damages and microbicidal activity) (Sang et al, 2015). This interpretation is supported by the upregulation of Il-10 observed in the lung of IVM-◀ Figure 3. Identification of macrophages in the lung of SARS-CoV-2-infected hamsters with and without ivermectin treatment and characterization of their transcriptomic profile related to M1/M2 polarization.
A Representative histopathology photomicrographies of lungs according to the different groups: mock_saline, CoV_saline, and CoV_ivermectin. Top panels: whole lung sections. Bottom panels: high magnification. CoV_saline section exhibits important congestion (*), edema associated with few mononuclear cells (white arrowheads). Note the thickening of the alveolar walls. CoV_ivermectin section exhibits important amounts of mononuclear cells (black arrowheads) and less marked signs of congestion or edema. Hematoxylin and eosin. Scale bars = 1 mm (top panels) and 20 μm (bottom panels). B Representative immunofluorescence photomicrographies of neutrophils (Ly-6G), monocytes/macrophages (Iba1), M2 macrophages (Arg1), and SARS-CoV-2 (NP) in the lung. Scale bars = 50 μm. C Quantification of Iba1 + cells, Arg1 + cells, and Iba1 + Arg1 + cells in the lungs. mock_saline n = 3 (males), mock_ivermectin n = 4 (males), CoV_saline n = 9 (6 males and 3 females), and CoV_ivermectin n = 6 (4 males and 2 females). D Heatmaps showing the differentially expressed genes related to the M1/M2 polarization in comparison with mock-infected hamsters. *indicates Benjamini-Hochberg-adjusted P-value < 0.05 in the comparison between saline and ivermectin within the same sex. Color gradient represents the transcription log 2 fold change comparing infected and mock-infected. Complete analyses are listed in Dataset EV1. E Lung weight-to-body weight ratio in the different groups (n = 4/sex/group). treated animals, along with the reduction of key M1 proinflammatory mediators, such as Il-6, Tnf-α, and Cxcl10 (Fig 2B), in similar ways as observed in other viral infections (Sang et al, 2015). Whereas few Iba1 + cells were observed in mock-infected animals, most likely resident interstitial and alveolar macrophages, we observed a large number of Iba1 + myeloid cells in the lungs of IVM-treated animals (Fig 3B). No distinguishable differences were observed in neutrophils population (Ly-6G + ) between saline-treated and IVM-treated animals (Fig 3B). Part of these Iba1 + cells were also Arg1 + cells (Fig 3B), a marker of M2-polarized macrophages. IVM treatment was associated with an increase of both Iba1 + and Arg1 + cells in the lungs of SARS-CoV-2-infected male hamsters (Fig 3C). In contrast, in females, while infection promoted an increase of both Iba1 + and Arg1 + cells (Fig 3C), IVM had no impact on the recruitment of these cells. Additionally, RNA-seq analyses in the lung identified the upregulation of key M2-related genes (Arg1, Cd209a/DC-SIGN, Clec7a/Dectin-1, and Myc/c-Myc) along with classical M1 markers (Cd86 and Cd38; Fig 3D), giving additional support to the M2 polarization tendency caused by the IVM treatment.

Conclusions
Our results demonstrate that IVM improves clinical outcome in SARS-CoV-2-infected animals and is associated with a reduced inflammatory status, but with no impact of SARS-CoV-2 loads in the upper and lower respiratory tracts. Thus, in hamsters, as in humans The results presented herein are consistent with a role of type I and III IFN responses in the pathogenesis of SARS-CoV-2-associated lung disease in hamsters. They show that IVM administration limits IFN response and lung inflammation, even though defects in the type I IFN pathways have been associated with severe COVID-19 (Bastard et al, 2020;Zhang et al, 2020a). This result may suggest that while IFN signaling is crucial to control viral replication and prevent severe disease, in infected hamsters, which only develop a moderate disease, IFN signaling may actually increase tissue damage and associated signs such as anosmia. This is consistent with previous reports in animal models where inhibition of IFN pathways was shown to increase viral replication but with lower lung pathology (Boudewijns et al, 2020;Sun et al, 2020).
Even if the effects observed may, to some extent, share similarities with the impact of dexamethasone (Horby et al, 2020) and tocilizumab (anti-IL-6) (Rossotti et al, 2020) on COVID-19, IVM action is steady and strong in the golden hamsters and is not expected to block the effectors of inflammation but rather to dampen its initiation. Interestingly, the activation of nuclear factor erythroid 2related transcription factor (Nrf2) via low-dose radiotherapy (LDRT) treatment has been proposed to cause a shift from M1 to M2 macrophages and a blockade of NLRP3 inflammasomes (Calabrese et al, 2021) and to be potentially beneficial on the lungs of COVID-19 patients. Along the same line, we noticed via our RNA-seq analyses that IVM treatment increased the gene expression of Nrf2 (or Nfe2l2) in the lungs of female hamsters and slightly in males (Dataset EV1), which gives additional evidence of the broad and upstream activity of IVM during SARS-CoV-2 infection. Further, our data show that these effects are compartmentalized and that the upper and lower respiratory tracts of hamsters respond differently to IVM treatment.
Considerable sex differences are observed in terms of clinical presentation, inflammatory profile, and transcriptomic signatures in the lungs of hamsters, as seen in COVID-19 human patients, where men tend to develop more severe disease than women (Jin et al, 2020; Takahashi  Moreover, the data presented herein are consistent with those reported in human clinical trials with IVM (Kory et al, 2021). In humans, IVM is widely used as anti-helminthic and anti-scabies at therapeutic doses (150-400 µg/kg) (Guzzo et al, 2002) that are in the range of those used in our hamster experiments. Further, several clinical trials on COVID-19 using IVM have been declared, where IVM has been associated with reduction of inflammatory markers and disease severity (preprint: Hill et al, 2021;Kaur et al, 2021). Interestingly, in a long-term care facility where the residents received IVM to control a scabies outbreak, no death or severe COVID-19 was observed (Bernigaud et al, 2021). IVM has already been administered to hospitalized COVID-19 patients, with contrasting outcomes: one study related no efficacy of late IVM administration (8-18 days after symptom onset) in severe COVID-19 patients treated in combination with other drugs (hydroxychloroquine, azithromycin, tocilizumab, steroids) (Camprubí et al, 2020), whereas another study reported lower mortality, especially in severe COVID-19 patients treated with IVM in addition to other treatments (hydroxychloroquine, azithromycin, or both) (Rajter et al, 2020). Importantly, in a study that administered IVM alone within 72 h of symptom onset, the authors noticed an important diminution of anosmia/hyposmia in the IVM group without differences in PCR positivity between IVM and placebo groups .
The presently available data support the view that IVM pharmacokinetics is quite stable across species, generally with slow absorption, broad distribution, low metabolism, and slow excretion, even if some conditions, such as route of administration, formulation, and body condition, may modify these features. In humans, the oral route is the only approved, and the administration of a single dose of 30 mg of IVM (corresponding to 347-594 µg/kg) leads to a C max (maximum concentration) of 84.8 ng/ml, a t max (time to reach the C max ) of 4.3 h and a half-life of 20.1 h (Guzzo et al, 2002). In animal models, the subcutaneous route is most often used, and comparatively, golden hamsters that received a subcutaneous injection of 400 µg/kg of IVM (as in the present study) showed comparable results with a C max of 80.2 ng/ml and a t max of˜4 h (Hanafi et al, 2011). Further, based on veterinary data 8 of 14 EMBO Molecular Medicine 13: e14122 | 2021 ª 2021 The Authors (Lifschitz et al, 2000), simulations using a minimal physiologically based pharmacokinetic (mPBPK) model revealed that the lungs would be exposed to IVM concentrations 2.7× greater than those found in the plasma (Jermain et al, 2020). Yet, this dose did not suffice to achieve the range of antiviral concentrations reported in vitro (Caly et al, 2020). Consequently, considering the results observed in the golden hamster model, IVM may be considered as a therapeutic agent against COVID-19, which would not strongly affect SARS-CoV-2 replication but limit the pathophysiological consequences of the infection in vivo, potentially mediated by type I and III IFN responses and several other related signaling pathways, and a favorable M1/M2 myleoid cells ratio in the lungs. A characteristic modulation of the immune response in the lower airways was observed in IVM-treated hamsters characterized by a transcriptomic profile similar to that observed in humans exhibiting less severe symptoms and a better prognosis (preprint: Masood et al, 2020;McElvaney et al, 2020b). Our data are consistent with the hypothesis that this effect is mediated by the cholinergic anti-inflammatory action of IVM on the vagus nerve reflex Tizabi et al, 2020), that should be addressed experimentally. In particular, the precise contribution of the nAChR in IVM action should be elucidated in comparison with that of other possible IVM targets (Zemkova et al, 2014). Altogether, this study brings the proof of concept that an IVM-based immunomodulatory therapy improves the clinical condition of SARS-CoV-2-infected hamsters, and in clinical trials, it alleviates symptoms of COVID-19 in humans and possibly limits post-COVID-19 syndrome (also known as long COVID) via an antiinflammatory action.

SARS-CoV-2 model and ivermectin treatment of hamsters
Male and female Syrian hamsters (Mesocricetus auratus; RjHan: AURA) of 5-6 weeks of age (average weight 60-80 g) were purchased from Janvier Laboratories and handled under specific pathogen-free conditions. The animals were housed and manipulated in isolators in a Biosafety level-3 facility, with ad libitum access to water and food. Before manipulation, animals underwent an acclimation period of 1 week.
Infected and mock-infected animals were housed in separate isolators, and all hamsters were followed up daily during 4 days at which the body weight and the clinical score were noted. The clinical score was based on a cumulative 0-4 scale: ruffled fur, slow movements, apathy, and absence of exploration activity.
At day 3 post-infection (dpi), animals underwent a food finding test to assess olfaction as previously described (Lazarini et al, 2012;de Melo et al, 2021). Briefly, 24 h before testing, hamsters were fasted and then individually placed into a fresh cage (37 × 29 × 18 cm) with clean standard bedding for 10 min. Subsequently, hamsters were placed in another similar cage for 2 min when about five pieces of cereals were hidden in 1.5 cm bedding in a corner of the test cage. The tested hamsters were then placed in the opposite corner, and the latency to find the food (defined as the time to locate cereals and start digging) was recorded using a chronometer. The test was carried out during a 15-min period. As soon as food was uncovered, hamsters were removed from the cage. One minute later, hamsters performed the same test but with visible chocolate cereals, positioned upon the bedding. The tests were realized in isolators in a Biosafety level-3 facility that were specially equipped for that.
At 4 dpi, animals were euthanized with an excess of anesthetics (ketamine and xylazine) and exsanguination (AVMA, 2020), and samples of nasal turbinates and lungs were collected and immediately frozen at −80°C. Fragments of lungs were also collected and fixed in 10% neutral buffered formalin.

RNA isolation and transcriptional analyses by quantitative PCR from golden hamsters' tissues
Frozen tissues were homogenized with TRIzol (15596026, Invitrogen) in Lysing Matrix D 2-ml tubes (116913100, MP Biomedicals) using the FastPrep-24 ™ system (MP Biomedicals) at the speed of 6.5 m/s during 1 min. Total RNA was extracted using the Direct-zol RNA MicroPrep Kit (R2062, Zymo Research: nasal turbinates) or MiniPrep Kit (R2052, Zymo Research: lung) and reverse-transcribed to first-strand cDNA using the SuperScript ™ IV VILO ™ Master Mix (11766050, Invitrogen). qPCR was performed in a final volume of 10 μl per reaction in 384-well PCR plates using a thermocycler ª 2021 The Authors EMBO Molecular Medicine 13: e14122 | 2021 (QuantStudio 6 Flex, Applied Biosystems). Briefly, 2.5 μl of cDNA (12.5 ng) was added to 7.5 μl of a master mix containing 5 μl of Power SYBR Green Mix (4367659, Applied Biosystems) and 2.5 μl of nuclease-free water with nCoV_IP2 primers (nCoV_IP2-12669Fw: 5 0 -ATGAGCTTAGTCCTGTTG-3 0 ; nCoV_IP2-12759Rv: 5 0 -CTCCCTTTGT TGTGTTGT-3 0 ) at a final concentration of 1 μM (WHO, 2020). The amplification conditions were as follows: 95°C for 10 min, 45 cycles of 95°C for 15 s and 60°C for 1 min; followed by a melt curve, from 60 to 95°C. Viral load quantification of hamster tissues was assessed by linear regression using a standard curve of eight known quantities of plasmids containing the RdRp sequence (ranging from 10 7 to 10 0 copies). The threshold of detection was established as 200 viral copies/µg of RNA. The Golden hamster gene targets were selected for quantifying host inflammatory mediator transcripts in the tissues using the Hprt (hypoxanthine phosphoribosyltransferase), the γactin, and/or the actinB genes as reference (Appendix Table S1). Variations in gene expression were calculated as the n-fold change in expression in the tissues from the infected hamsters compared with the tissues of the uninfected ones using the 2 ÀΔΔCt method (Pfaffl, 2001).

Droplet digital PCR (ddPCR)
Reverse transcription 200 ng of RNA was reverse-transcribed using iScript Advanced cDNA Synthesis kit for RT-qPCR (1702537, Bio-Rad) according to the manufacturer's specifications.
Quantitative PCR for γ-actin and Hprt reference genes Real-time PCR was performed in a CFX96 qPCR machine (Bio-Rad). All samples were measured in duplicate. The 10 μl PCR included 0.8 ng of cDNA, 1× PowerUp PCR master mix (A25742, Applied Biosystems), and 0.5 µM of each primer (Appendix Table S1). The reactions were incubated in a 96-well optical plate at 95°C for 2 min, followed by 40 cycles of 95°C for 15 s and 60°C for 1 min.
Droplet digital PCR ddPCRs were performed on the QX200 Droplet Digital PCR system according to the manufacturer's instructions (Bio-Rad). Briefly, reaction mixture consisted in 10 μl ddPCR Supermix for probe no dUTP (1863023, Bio-Rad), 0.25-1 ng of cDNA, primers and probes for E/ IP4 and N/nsp13 duplex reactions used at concentration listed in Appendix Table S2 in a final volume of 20 μl. PCR amplification was conducted in a iCycler PCR instrument (Bio-Rad) with the following condition: 95°C for 10 min, 40 cycles of 94°C for 30 s with a ramping of 2°/s, 59°C for 1 min with a ramping of 2°/s, followed by 98°C for 5 min with a ramping of 2°/s and a hold at 4°C. After amplification, the 96-well plate was loaded onto the QX200 droplet reader (Bio-Rad) that measures automatically the fluorescence intensity in individual droplets. Generated data were subsequently analyzed with QuantaSoft ™ software (Bio-Rad) based on positive and negative droplet populations. Data are expressed as CPD (copy per droplets) normalized to γ-actin and Hprt reference gene relative expression.

Transcriptomics analysis in golden hamsters' lung
RNA preparation was used to construct strand-specific single-end cDNA libraries according to the manufacturers' instructions (TruSeq Stranded mRNA sample prep kit, Illumina). Illumina NextSeq 500 sequencer was used to sequence libraries. The RNA-seq analysis was performed with the Sequana framework (Cokelaer et al, 2017). We used the RNA-seq pipeline (v0.9.16), which is available online (https://github.com/sequana/sequana_rnaseq). It is built on top of Snakemake 5.8.1 (Köster & Rahmann, 2012). Reads were trimmed from adapters using Cutadapt 2.10 (Martin, 2011) and then mapped to the golden hamster MesAur1.0 genome assembly from Ensembl using STAR 2.7.3a (Dobin et al, 2012). FeatureCounts 2.0.0 (Liao et al, 2014) was used to produce the count matrix, assigning reads to features using annotation MesAur1.0.100 with strand-specificity information. Quality control statistics were summarized using MultiQC 1.8 (Ewels et al, 2016). Statistical analysis on the count matrix was performed to identify differentially regulated genes, comparing infected versus non-infected samples considering all samples and separating by sex. Clustering of transcriptomic profiles was assessed using a principal component analysis (PCA). Differential expression testing was conducted using DESeq2 library 1.24.0 (Love et al, 2014) scripts based on SARTools 1.7.0 (Varet et al, 2016) indicating the significance (Benjamini-Hochberg-adjusted P-values, false discovery rate FDR < 0.05) and the effect size (fold change) for each comparison. Finally, enrichment analysis was performed using modules from Sequana, first by converting golden hamster ensembl ids to gene names and then using human annotations for GO terms and KEGG pathways. The GO enrichment module uses PantherDB (Mi et al, 2019) and QuickGO (Huntley et al, 2014) services; the KEGG pathways enrichment uses gseapy (https://github.com/zqfa ng/GSEApy/), EnrichR (Chen et al, 2013), KEGG (Kanehisa & Goto, 2000), and BioMart services. All programmatic accesses to the online web services were performed via BioServices (Cokelaer et al, 2013).

Histopathology
Lung fragments fixed in 10% neutral buffered formalin were embedded in paraffin. Four-µm-thick sections were cut and stained with hematoxylin and eosin staining. The slides were then scanned using Axioscan Z1 Zeiss slide scanner, using the Zen 2 blue edition software.

Statistics
Statistical analysis was performed using Prism software (GraphPad, version 9.0.0, San Diego, USA), with P < 0.05 considered significant. Quantitative data were compared across groups using log-rank test or two-tailed Mann-Whitney test. Randomization and blinding were not possible due to pre-defined housing conditions (separated isolators between infected and non-infected animals). Ex vivo analysis was blinded (coded samples). All animals were included, and data were provided from 2 replications, except food finding in males, that were replicated 3 times.
Expanded View for this article is available online.

Acknowledgements
The SARS-CoV-2 strain was supplied by the National Reference Centre for

Conflict of interest
The authors declare that they have no conflict of interest.

References
Aamir K, Khan HU, Sethi G, Hossain MA, Arya A (2020) Wnt signaling mediates TLR pathway and promote unrestrained adipogenesis and The paper explained

Problem
The current pandemic of COVID-19 has caused more than 3.5 million deaths and more than 150 million laboratory-confirmed cases worldwide since December 2019 (as of May 2021). COVID-19, caused by SARS-CoV-2, commonly brings about upper airways and pulmonary symptoms and in severe cases can lead to respiratory distress and death. Different therapeutic approaches have been proposed to fight this disease but comprehensive therapeutic studies are still lacking.

Results
We report that ivermectin, used at the standard anti-parasitic dose of 400 µg/kg, protects infected hamsters from developing clinical signs and from losing the sense of smell during SARS-CoV-2 infection. The treated animals exhibited a specific inflammatory response, presenting a reduced type I/III interferon stimulation and a modulation in several intracellular signaling pathways, with an important reduction of the Il-6/Il-10 ratio and promoting M2 polarization of myeloid cells recruited to the lung. These effects are strongly influenced by sex, with treated females exhibiting the best outcome. Surprisingly, ivermectin treatment did not limit viral replication, as treated and nontreated animals presented similar amounts of SARS-CoV-2 in the nasal cavity and in the lungs.

Impact
The results of this study establish that irrespective of viral load, the symptoms and severity of COVID-19 highlight the critical role played by host inflammatory response in COVID-19 severity and highlight that reduced type I/III interferon and Il-6/Il-10 and the presence of M2 macrophages might account for a more favorable clinical presentation, contributing to a better understanding of COVID-19 pathophysiology. Ivermectin might then be considered as promising therapeutic agent against COVID-19 with no impact on SARS-CoV-2 replication but alleviating inflammation and ensuing symptoms.