A molnupiravir-associated mutational signature in global SARS-CoV-2 genomes

Molnupiravir, an antiviral medication widely used against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), acts by inducing mutations in the virus genome during replication. Most random mutations are likely to be deleterious to the virus and many will be lethal; thus, molnupiravir-induced elevated mutation rates reduce viral load1,2. However, if some patients treated with molnupiravir do not fully clear the SARS-CoV-2 infections, there could be the potential for onward transmission of molnupiravir-mutated viruses. Here we show that SARS-CoV-2 sequencing databases contain extensive evidence of molnupiravir mutagenesis. Using a systematic approach, we find that a specific class of long phylogenetic branches, distinguished by a high proportion of G-to-A and C-to-T mutations, are found almost exclusively in sequences from 2022, after the introduction of molnupiravir treatment, and in countries and age groups with widespread use of the drug. We identify a mutational spectrum, with preferred nucleotide contexts, from viruses in patients known to have been treated with molnupiravir and show that its signature matches that seen in these long branches, in some cases with onward transmission of molnupiravir-derived lineages. Finally, we analyse treatment records to confirm a direct association between these high G-to-A branches and the use of molnupiravir.

Molnupiravir versus placebo in unvaccinated and vaccinated patients with early SARS-CoV-2 infection in the UK (AGILE CST-2): a randomised, placebo-controlled, double-blind, phase 2 trial

BACKGROUND

The antiviral drug molnupiravir was licensed for treating at-risk patients with COVID-19 on the basis of data from unvaccinated adults. We aimed to evaluate the safety and virological efficacy of molnupiravir in vaccinated and unvaccinated individuals with COVID-19.

METHODS

This randomised, placebo-controlled, double-blind, phase 2 trial (AGILE CST-2) was done at five National Institute for Health and Care Research sites in the UK. Eligible participants were adult (aged ≥18 years) outpatients with PCR-confirmed, mild-to-moderate SARS-CoV-2 infection who were within 5 days of symptom onset. Using permuted blocks (block size 2 or 4) and stratifying by site, participants were randomly assigned (1:1) to receive either molnupiravir (orally; 800 mg twice daily for 5 days) plus standard of care or matching placebo plus standard of care. The primary outcome was the time from randomisation to SARS-CoV-2 PCR negativity on nasopharyngeal swabs and was analysed by use of a Bayesian Cox proportional hazards model for estimating the probability of a superior virological response (hazard ratio [HR]>1) for molnupiravir versus placebo. Our primary model used a two-point prior based on equal prior probabilities (50%) that the HR was 1·0 or 1·5. We defined a priori that if the probability of a HR of more than 1 was more than 80% molnupiravir would be recommended for further testing. The primary outcome was analysed in the intention-to-treat population and safety was analysed in the safety population, comprising participants who had received at least one dose of allocated treatment. This trial is registered in ClinicalTrials.gov, NCT04746183, and the ISRCTN registry, ISRCTN27106947, and is ongoing.

FINDINGS

Between Nov 18, 2020, and March 16, 2022, 1723 patients were assessed for eligibility, of whom 180 were randomly assigned to receive either molnupiravir (n=90) or placebo (n=90) and were included in the intention-to-treat analysis. 103 (57%) of 180 participants were female and 77 (43%) were male and 90 (50%) participants had received at least one dose of a COVID-19 vaccine. SARS-CoV-2 infections with the delta (B.1.617.2; 72 [40%] of 180), alpha (B.1.1.7; 37 [21%]), omicron (B.1.1.529; 38 [21%]), and EU1 (B.1.177; 28 [16%]) variants were represented. All 180 participants received at least one dose of treatment and four participants discontinued the study (one in the molnupiravir group and three in the placebo group). Participants in the molnupiravir group had a faster median time from randomisation to negative PCR (8 days [95% CI 8-9]) than participants in the placebo group (11 days [10-11]; HR 1·30, 95% credible interval 0·92-1·71; log-rank p=0·074). The probability of molnupiravir being superior to placebo (HR>1) was 75·4%, which was less than our threshold of 80%. 73 (81%) of 90 participants in the molnupiravir group and 68 (76%) of 90 participants in the placebo group had at least one adverse event by day 29. One participant in the molnupiravir group and three participants in the placebo group had an adverse event of a Common Terminology Criteria for Adverse Events grade 3 or higher severity. No participants died (due to any cause) during the trial.

INTERPRETATION

We found molnupiravir to be well tolerated and, although our predefined threshold was not reached, we observed some evidence that molnupiravir has antiviral activity in vaccinated and unvaccinated individuals infected with a broad range of SARS-CoV-2 variants, although this evidence is not conclusive.

FUNDING

Ridgeback Biotherapeutics, the UK National Institute for Health and Care Research, the Medical Research Council, and the Wellcome Trust.

Analysis of SARS-CoV-2 known and novel subgenomic mRNAs in cell culture, animal model, and clinical samples using LeTRS, a bioinformatic tool to identify unique sequence identifiers

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has a complex strategy for the transcription of viral subgenomic mRNAs (sgmRNAs), which are targets for nucleic acid diagnostics. Each of these sgmRNAs has a unique 5' sequence, the leader-transcriptional regulatory sequence gene junction (leader-TRS junction), that can be identified using sequencing. High-resolution sequencing has been used to investigate the biology of SARS-CoV-2 and the host response in cell culture and animal models and from clinical samples. LeTRS, a bioinformatics tool, was developed to identify leader-TRS junctions and can be used as a proxy to quantify sgmRNAs for understanding virus biology. LeTRS is readily adaptable for other coronaviruses such as Middle East respiratory syndrome coronavirus or a future newly discovered coronavirus. LeTRS was tested on published data sets and novel clinical samples from patients and longitudinal samples from animal models with coronavirus disease 2019. LeTRS identified known leader-TRS junctions and identified putative novel sgmRNAs that were common across different mammalian species. This may be indicative of an evolutionary mechanism where plasticity in transcription generates novel open reading frames, which can then subject to selection pressure. The data indicated multiphasic abundance of sgmRNAs in two different animal models. This recapitulates the relative sgmRNA abundance observed in cells at early points in infection but not at late points. This pattern is reflected in some human nasopharyngeal samples and therefore has implications for transmission models and nucleic acid-based diagnostics. LeTRS provides a quantitative measure of sgmRNA abundance from sequencing data. This can be used to assess the biology of SARS-CoV-2 (or other coronaviruses) in clinical and nonclinical samples, especially to evaluate different variants and medical countermeasures that may influence viral RNA synthesis.

Mutations that adapt SARS-CoV-2 to mink or ferret do not increase fitness in the human airway

SARS-CoV-2 has a broad mammalian species tropism infecting humans, cats, dogs, and farmed mink. Since the start of the 2019 pandemic, several reverse zoonotic outbreaks of SARS-CoV-2 have occurred in mink, one of which reinfected humans and caused a cluster of infections in Denmark. Here we investigate the molecular basis of mink and ferret adaptation and demonstrate the spike mutations Y453F, F486L, and N501T all specifically adapt SARS-CoV-2 to use mustelid ACE2. Furthermore, we risk assess these mutations and conclude mink-adapted viruses are unlikely to pose an increased threat to humans, as Y453F attenuates the virus replication in human cells and all three mink adaptations have minimal antigenic impact. Finally, we show that certain SARS-CoV-2 variants emerging from circulation in humans may naturally have a greater propensity to infect mustelid hosts and therefore these species should continue to be surveyed for reverse zoonotic infections.

Deep splicing plasticity of the human adenovirus type 5 transcriptome drives virus evolution

Viral genomes have high gene densities and complex transcription strategies rendering transcriptome analysis through short-read RNA-seq approaches problematic. Adenovirus transcription and splicing is especially complex. We used long-read direct RNA sequencing to study adenovirus transcription and splicing during infection. This revealed a previously unappreciated complexity of alternative splicing and potential for secondary initiating codon usage. Moreover, we find that most viral transcripts tend to shorten polyadenylation lengths as infection progresses. Development of an open reading frame centric bioinformatics analysis pipeline provided a deeper quantitative and qualitative understanding of adenovirus's genetic potential. Across the viral genome adenovirus makes multiple distinctly spliced transcripts that code for the same protein. Over 11,000 different splicing patterns were recorded across the viral genome, most occurring at low levels. This low-level use of alternative splicing patterns potentially enables the virus to maximise its coding potential over evolutionary timescales.