Genetic consequences of effective and suboptimal dosing with mutagenic drugs in a hamster model of SARS-CoV-2 infection.

Mutagenic antiviral drugs have shown promise against multiple viruses, but concerns have been raised about whether their use might promote the emergence of new and harmful viral variants. Recently, genetic signatures associated with molnupiravir use have been identified in the global SARS-COV-2 population. Here, we examine the consequences of using favipiravir and molnupiravir to treat SARS-CoV-2 infection in a hamster model, comparing viral genome sequence data collected from (1) untreated hamsters, and (2) from hamsters receiving effective and suboptimal doses of treatment. We identify a broadly linear relationship between drug dose and the extent of variation in treated viral populations, with a high proportion of this variation being composed of variants at frequencies of less than 1 per cent, below typical thresholds for variant calling. Treatment with an effective dose of antiviral drug was associated with a gain of between 7 and 10 variants per viral genome relative to drug-free controls: even after a short period of treatment a population founded by a transmitted virus could contain multiple sequence differences to that of the original host. Treatment with a suboptimal dose of drug showed intermediate gains of variants. No dose-dependent signal was identified in the numbers of single-nucleotide variants reaching frequencies in excess of 5 per cent. We did not find evidence to support the emergence of drug resistance or of novel immune phenotypes. Our study suggests that where onward transmission occurs, a short period of treatment with mutagenic drugs may be sufficient to generate a significant increase in the number of viral variants transmitted.

Detection of a historic reservoir of bedaquiline/clofazimine resistance-associated variants in Mycobacterium tuberculosis.

BACKGROUND: Drug resistance in tuberculosis (TB) poses a major ongoing challenge to public health. The recent inclusion of bedaquiline into TB drug regimens has improved treatment outcomes, but this advance is threatened by the emergence of strains of Mycobacterium tuberculosis (Mtb) resistant to bedaquiline. Clinical bedaquiline resistance is most frequently conferred by off-target resistance-associated variants (RAVs) in the mmpR5 gene (Rv0678), the regulator of an efflux pump, which can also confer cross-resistance to clofazimine, another TB drug. METHODS: We compiled a dataset of 3682 Mtb genomes, including 180 carrying variants in mmpR5, and its immediate background (i.e. mmpR5 promoter and adjacent mmpL5 gene), that have been associated to borderline (henceforth intermediate) or confirmed resistance to bedaquiline. We characterised the occurrence of all nonsynonymous mutations in mmpR5 in this dataset and estimated, using time-resolved phylogenetic methods, the age of their emergence. RESULTS: We identified eight cases where RAVs were present in the genomes of strains collected prior to the use of bedaquiline in TB treatment regimes. Phylogenetic reconstruction points to multiple emergence events and circulation of RAVs in mmpR5, some estimated to predate the introduction of bedaquiline. However, epistatic interactions can complicate bedaquiline drug-susceptibility prediction from genetic sequence data. Indeed, in one clade, Ile67fs (a RAV when considered in isolation) was estimated to have emerged prior to the antibiotic era, together with a resistance reverting mmpL5 mutation. CONCLUSIONS: The presence of a pre-existing reservoir of Mtb strains carrying bedaquiline RAVs prior to its clinical use augments the need for rapid drug susceptibility testing and individualised regimen selection to safeguard the use of bedaquiline in TB care and control.

Favipiravir induces HuNoV viral mutagenesis and infectivity loss with clinical improvement in immunocompromised patients.

Chronic human norovirus (HuNoV) infections in immunocompromised patients result in severe disease, yet approved antivirals are lacking. RNA-dependent RNA polymerase (RdRp) inhibitors inducing viral mutagenesis display broad-spectrum in vitro antiviral activity, but clinical efficacy in HuNoV infections is anecdotal and the potential emergence of drug-resistant variants is concerning. Upon favipiravir (and nitazoxanide) treatment of four immunocompromised patients with life-threatening HuNoV infections, viral whole-genome sequencing showed accumulation of favipiravir-induced mutations which coincided with clinical improvement although treatment failed to clear HuNoV. Infection of zebrafish larvae demonstrated drug-associated loss of viral infectivity and favipiravir treatment showed efficacy despite occurrence of RdRp variants potentially causing favipiravir resistance. This indicates that within-host resistance evolution did not reverse loss of viral fitness caused by genome-wide accumulation of sequence changes. This off-label approach supports the use of mutagenic antivirals for treating prolonged RNA viral infections and further informs the debate surrounding their impact on virus evolution.

Haplotype assignment of longitudinal viral deep sequencing data using covariation of variant frequencies.

Longitudinal deep sequencing of viruses can provide detailed information about intra-host evolutionary dynamics including how viruses interact with and transmit between hosts. Many analyses require haplotype reconstruction, identifying which variants are co-located on the same genomic element. Most current methods to perform this reconstruction are based on a high density of variants and cannot perform this reconstruction for slowly evolving viruses. We present a new approach, HaROLD (HAplotype Reconstruction Of Longitudinal Deep sequencing data), which performs this reconstruction based on identifying co-varying variant frequencies using a probabilistic framework. We illustrate HaROLD on both RNA and DNA viruses with synthetic Illumina paired read data created from mixed human cytomegalovirus (HCMV) and norovirus genomes, and clinical datasets of HCMV and norovirus samples, demonstrating high accuracy, especially when longitudinal samples are available.

Evolution of viral variants in remdesivir-treated and untreated SARS-CoV-2-infected pediatrics patients.

Detailed information on intrahost viral evolution in SARS-CoV-2 with and without treatment is limited. Sequential viral loads and deep sequencing of SARS-CoV-2 from the upper respiratory tract of nine hospitalized children, three of whom were treated with remdesivir, revealed that remdesivir treatment suppressed viral load in one patient but not in a second infected with an identical strain without any evidence of drug resistance found. Reduced levels of subgenomic RNA during treatment of the second patient, suggest an additional effect of remdesivir on viral replication. Haplotype reconstruction uncovered persistent SARS-CoV-2 variant genotypes in four patients. These likely arose from within-host evolution, although superinfection cannot be excluded in one case. Although our dataset is small, observed sample-to-sample heterogeneity in variant frequencies across four of nine patients suggests the presence of discrete viral populations in the lung with incomplete population sampling in diagnostic swabs. Such compartmentalization could compromise the penetration of remdesivir into the lung, limiting the drugs in vivo efficacy, as has been observed in other lung infections.

The combined treatment of Molnupiravir and Favipiravir results in a potentiation of antiviral efficacy in a SARS-CoV-2 hamster infection model.

BACKGROUND: Favipiravir and Molnupiravir, orally available antivirals, have been reported to exert antiviral activity against SARS-CoV-2. First efficacy data have been recently reported in COVID-19 patients. METHODS: We here report on the combined antiviral effect of both drugs in a SARS-CoV-2 Syrian hamster infection model. The infected hamsters were treated twice daily with the vehicle (the control group) or a suboptimal dose of each compound or a combination of both compounds. FINDINGS: When animals were treated with a combination of suboptimal doses of Molnupiravir and Favipiravir at the time of infection, a marked combined potency at endpoint is observed. Infectious virus titers in the lungs of animals treated with the combination are reduced by ∼5 log10 and infectious virus are no longer detected in the lungs of >60% of treated animals. When start of treatment was delayed with one day a reduction of titers in the lungs of 2.4 log10 was achieved. Moreover, treatment of infected animals nearly completely prevented transmission to co-housed untreated sentinels. Both drugs result in an increased mutation frequency of the remaining viral RNA recovered from the lungs of treated animals. In the combo-treated hamsters, an increased frequency of C-to-T mutations in the viral RNA is observed as compared to the single treatment groups which may explain the pronounced antiviral potency of the combination. INTERPRETATION: Our findings may lay the basis for the design of clinical studies to test the efficacy of the combination of Molnupiravir/Favipiravir in the treatment of COVID-19. FUNDING: stated in the acknowledgment.

Systematic Review and Patient-Level Meta-Analysis of SARS-CoV-2 Viral Dynamics to Model Response to Antiviral Therapies.

Severe acute respiratory syndrome-coronavirus 2 (SARS-CoV-2) viral loads change rapidly following symptom onset, so to assess antivirals it is important to understand the natural history and patient factors influencing this. We undertook an individual patient-level meta-analysis of SARS-CoV-2 viral dynamics in humans to describe viral dynamics and estimate the effects of antivirals used to date. This systematic review identified case reports, case series, and clinical trial data from publications between January 1, 2020, and May 31, 2020, following Preferred Reporting Items for Systematic Reviews and Meta-Analyses guidelines. A multivariable Cox proportional hazards (Cox-PH) regression model of time to viral clearance was fitted to respiratory and stool samples. A simplified four parameter nonlinear mixed-effects (NLME) model was fitted to viral load trajectories in all sampling sites and covariate modeling of respiratory viral dynamics was performed to quantify time-dependent drug effects. Patient-level data from 645 individuals (age 1 month to 100 years) with 6,316 viral loads were extracted. Model-based simulations of viral load trajectories in samples from the upper and lower respiratory tract, stool, blood, urine, ocular secretions, and breast milk were generated. Cox-PH modeling showed longer time to viral clearance in older patients, men, and those with more severe disease. Remdesivir was associated with faster viral clearance (adjusted hazard ratio (AHR) = 9.19, P < 0.001), as well as interferon, particularly when combined with ribavirin (AHR = 2.2, P = 0.015; AHR = 6.04, P = 0.006). Combination therapy should be further investigated. A viral dynamic dataset and NLME model for designing and analyzing antiviral trials has been established.

Mixed cytomegalovirus genotypes in HIV-positive mothers show compartmentalization and distinct patterns of transmission to infants.

Cytomegalovirus (CMV) is the commonest cause of congenital infection and particularly so among infants born to HIV-infected women. Studies of congenital CMV infection (cCMVi) pathogenesis are complicated by the presence of multiple infecting maternal CMV strains, especially in HIV-positive women, and the large, recombinant CMV genome. Using newly developed tools to reconstruct CMV haplotypes, we demonstrate anatomic CMV compartmentalization in five HIV-infected mothers and identify the possibility of congenitally transmitted genotypes in three of their infants. A single CMV strain was transmitted in each congenitally infected case, and all were closely related to those that predominate in the cognate maternal cervix. Compared to non-transmitted strains, these congenitally transmitted CMV strains showed statistically significant similarities in 19 genes associated with tissue tropism and immunomodulation. In all infants, incident superinfections with distinct strains from breast milk were captured during follow-up. The results represent potentially important new insights into the virologic determinants of early CMV infection.

Favipiravir at high doses has potent antiviral activity in SARS-CoV-2-infected hamsters, whereas hydroxychloroquine lacks activity.

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) rapidly spread around the globe after its emergence in Wuhan in December 2019. With no specific therapeutic and prophylactic options available, the virus has infected millions of people of which more than half a million succumbed to the viral disease, COVID-19. The urgent need for an effective treatment together with a lack of small animal infection models has led to clinical trials using repurposed drugs without preclinical evidence of their in vivo efficacy. We established an infection model in Syrian hamsters to evaluate the efficacy of small molecules on both infection and transmission. Treatment of SARS-CoV-2-infected hamsters with a low dose of favipiravir or hydroxychloroquine with(out) azithromycin resulted in, respectively, a mild or no reduction in virus levels. However, high doses of favipiravir significantly reduced infectious virus titers in the lungs and markedly improved lung histopathology. Moreover, a high dose of favipiravir decreased virus transmission by direct contact, whereas hydroxychloroquine failed as prophylaxis. Pharmacokinetic modeling of hydroxychloroquine suggested that the total lung exposure to the drug did not cause the failure. Our data on hydroxychloroquine (together with previous reports in macaques and ferrets) thus provide no scientific basis for the use of this drug in COVID-19 patients. In contrast, the results with favipiravir demonstrate that an antiviral drug at nontoxic doses exhibits a marked protective effect against SARS-CoV-2 in a small animal model. Clinical studies are required to assess whether a similar antiviral effect is achievable in humans without toxic effects.

Emergence of genomic diversity and recurrent mutations in SARS-CoV-2.

SARS-CoV-2 is a SARS-like coronavirus of likely zoonotic origin first identified in December 2019 in Wuhan, the capital of China's Hubei province. The virus has since spread globally, resulting in the currently ongoing COVID-19 pandemic. The first whole genome sequence was published on January 5 2020, and thousands of genomes have been sequenced since this date. This resource allows unprecedented insights into the past demography of SARS-CoV-2 but also monitoring of how the virus is adapting to its novel human host, providing information to direct drug and vaccine design. We curated a dataset of 7666 public genome assemblies and analysed the emergence of genomic diversity over time. Our results are in line with previous estimates and point to all sequences sharing a common ancestor towards the end of 2019, supporting this as the period when SARS-CoV-2 jumped into its human host. Due to extensive transmission, the genetic diversity of the virus in several countries recapitulates a large fraction of its worldwide genetic diversity. We identify regions of the SARS-CoV-2 genome that have remained largely invariant to date, and others that have already accumulated diversity. By focusing on mutations which have emerged independently multiple times (homoplasies), we identify 198 filtered recurrent mutations in the SARS-CoV-2 genome. Nearly 80% of the recurrent mutations produced non-synonymous changes at the protein level, suggesting possible ongoing adaptation of SARS-CoV-2. Three sites in Orf1ab in the regions encoding Nsp6, Nsp11, Nsp13, and one in the Spike protein are characterised by a particularly large number of recurrent mutations (>15 events) which may signpost convergent evolution and are of particular interest in the context of adaptation of SARS-CoV-2 to the human host. We additionally provide an interactive user-friendly web-application to query the alignment of the 7666 SARS-CoV-2 genomes.