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

Direct RNA sequencing of respiratory syncytial virus infected human cells generates a detailed overview of RSV polycistronic mRNA and transcript abundance.

To characterize species of viral mRNA transcripts generated during respiratory syncytial virus (RSV) infection, human fibroblast-like MRC-5 lung cells were infected with subgroup A RSV for 6, 16 and 24 hours. In addition, we characterised the viral transcriptome in infected Calu-3 lung epithelial cells at 48 hours post infection. Total RNA was harvested and polyadenylated mRNA was enriched and sequenced by direct RNA sequencing using an Oxford nanopore device. This platform yielded over 450,000 direct mRNA transcript reads which were mapped to the viral genome and analysed to determine the relative mRNA levels of viral genes using our in-house ORF-centric pipeline. We examined the frequency of polycistronic readthrough mRNAs were generated and assessed the length of the polyadenylated tails for each group of transcripts. We show a general but non-linear decline in gene transcript abundance across the viral genome, as predicted by the model of RSV gene transcription. However, the decline in transcript abundance is not uniform. The polyadenylate tails generated by the viral polymerase are similar in length to those generated by the host polyadenylation machinery and broadly declined in length for most transcripts as the infection progressed. Finally, we observed that the steady state abundance of transcripts with very short polyadenylate tails less than 20 nucleotides is less for N, SH and G transcripts in both cell lines compared to NS1, NS2, P, M, F and M2 which may reflect differences in mRNA stability and/or translation rates within and between the cell lines

Genome evolution of dengue virus serotype 1 under selection by <i>Wolbachia pipientis</i> in <i>Aedes aegypti</i> mosquitoes.

The introgression of antiviral strains of Wolbachia into Aedes aegypti mosquito populations is a public health intervention for the control of dengue. Plausibly, dengue virus (DENV) could evolve to bypass the antiviral effects of Wolbachia and undermine this approach. Here, we established a serial-passage system to investigate the evolution of DENV in Ae. aegypti mosquitoes infected with the wMel strain of Wolbachia. Using this system, we report on virus genetic outcomes after twenty passages of serotype 1 of DENV (DENV-1). An amino acid substitution, E203K, in the DENV-1 envelope protein was more frequently detected in the consensus sequence of virus populations passaged in wMel-infected Ae. aegypti than wild-type counterparts. Positive selection at residue 203 was reproducible; it occurred in passaged virus populations from independent DENV-1-infected patients and also in a second, independent experimental system. In wild-type mosquitoes and human cells, the 203K variant was rapidly replaced by the progenitor sequence. These findings provide proof of concept that wMel-associated selection of virus populations can occur in experimental conditions. Field-based studies are needed to explore whether wMel imparts selective pressure on DENV evolution in locations where wMel is established

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

Acknowledgements: We thank all data contributors, that is, the authors and their originating laboratories responsible for obtaining the specimens, and their submitting laboratories for generating the genetic sequence and metadata and sharing via the GISAID initiative, on which this research is based. We also thank everyone who contributed to the generation of the genomes deposited in the INSDC databases, on which this research is also based. We thank A. Hinrichs and colleagues for access to an UShER mutation-annotated tree built with all available genomic data. We thank NHS England for providing the Blueteq data on treatment records. We thank the UKHSA COVID-19 Therapeutics Programme Team past and present, in particular J. Charlesworth, A. Lackenby, A. Demirjian, M. Chand and C. Brown. We thank J. Bloom, M. Lin, R. Neher, K. Harris and F. Débarre for useful discussions. T.S. was supported by the Wellcome Trust (no. 210918/Z/18/Z) and the Francis Crick Institute, which receives its core funding from Cancer Research UK (no. FC001043), the UK Medical Research Council (MRC) (no. FC001043) and the Wellcome Trust (no. FC001043). This research was funded in whole, or in part, by the Wellcome Trust (nos. 210918/Z/18/Z, FC001043). For the purpose of open access, the authors have applied a CC-BY public copyright licence to any author-accepted manuscript resulting from this article. I.D.-B. is supported by PhD funding from the National Institute for Health and Care Research (NIHR) Health Protection Research Unit in Emerging and Zoonotic Infections at the University of Liverpool in partnership with Public Health England (now UKHSA), in collaboration with the Liverpool School of Tropical Medicine and the University of Oxford (award no. 200907). The views expressed are those of the authors and not necessarily those of the Department of Health and Social Care or NIHR. Neither the funders nor the trial sponsor were involved in study design, data collection, analysis, interpretation or preparation of the manuscript. T.P.P. was funded by the G2P-UK National Virology Consortium, which is funded by the MRC (no. MR/W005611/1). C.R. was supported by a Fondation Botnar Research Award (programme grant no. 6063), the UK Cystic Fibrosis Trust (Innovation Hub Award 001) and funding from the Oxford Martin School.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

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

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

Amplicon and Metagenomic Analysis of Middle East Respiratory Syndrome (MERS) Coronavirus and the Microbiome in Patients with Severe MERS.

Middle East respiratory syndrome coronavirus (MERS-CoV) is a zoonotic infection that emerged in the Middle East in 2012. Symptoms range from mild to severe and include both respiratory and gastrointestinal illnesses. The virus is mainly present in camel populations with occasional zoonotic spill over into humans. The severity of infection in humans is influenced by numerous factors, and similar to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), underlying health complications can play a major role. Currently, MERS-CoV and SARS-CoV-2 are coincident in the Middle East and thus a rapid way of sequencing MERS-CoV to derive genotype information for molecular epidemiology is needed. Additionally, complicating factors in MERS-CoV infections are coinfections that require clinical management. The ability to rapidly characterize these infections would be advantageous. To rapidly sequence MERS-CoV, an amplicon-based approach was developed and coupled to Oxford Nanopore long read length sequencing. This and a metagenomic approach were evaluated with clinical samples from patients with MERS. The data illustrated that whole-genome or near-whole-genome information on MERS-CoV could be rapidly obtained. This approach provided data on both consensus genomes and the presence of minor variants, including deletion mutants. The metagenomic analysis provided information of the background microbiome. The advantage of this approach is that insertions and deletions can be identified, which are the major drivers of genotype change in coronaviruses. IMPORTANCE Middle East respiratory syndrome coronavirus (MERS-CoV) emerged in late 2012 in Saudi Arabia. The virus is a serious threat to people not only in the Middle East but also in the world and has been detected in over 27 countries. MERS-CoV is spreading in the Middle East and neighboring countries, and approximately 35% of reported patients with this virus have died. This is the most severe coronavirus infection so far described. Saudi Arabia is a destination for many millions of people in the world who visit for religious purposes (Umrah and Hajj), and so it is a very vulnerable area, which imposes unique challenges for effective control of this epidemic. The significance of our study is that clinical samples from patients with MERS were used for rapid in-depth sequencing and metagenomic analysis using long read length sequencing

The P323L substitution in the SARS-CoV-2 polymerase (NSP12) confers a selective advantage during infection.

BACKGROUND: The mutational landscape of SARS-CoV-2 varies at the dominant viral genome sequence and minor genomic variant population. During the COVID-19 pandemic, an early substitution in the genome was the D614G change in the spike protein, associated with an increase in transmissibility. Genomes with D614G are accompanied by a P323L substitution in the viral polymerase (NSP12). However, P323L is not thought to be under strong selective pressure. RESULTS: Investigation of P323L/D614G substitutions in the population shows rapid emergence during the containment phase and early surge phase during the first wave. These substitutions emerge from minor genomic variants which become dominant viral genome sequence. This is investigated in vivo and in vitro using SARS-CoV-2 with P323 and D614 in the dominant genome sequence and L323 and G614 in the minor variant population. During infection, there is rapid selection of L323 into the dominant viral genome sequence but not G614. Reverse genetics is used to create two viruses (either P323 or L323) with the same genetic background. L323 shows greater abundance of viral RNA and proteins and a smaller plaque morphology than P323. CONCLUSIONS: These data suggest that P323L is an important contribution in the emergence of variants with transmission advantages. Sequence analysis of viral populations suggests it may be possible to predict the emergence of a new variant based on tracking the frequency of minor variant genomes. The ability to predict an emerging variant of SARS-CoV-2 in the global landscape may aid in the evaluation of medical countermeasures and non-pharmaceutical interventions

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

Characterisation of SARS-CoV-2 genomic variation in response to molnupiravir treatment in the AGILE Phase IIa clinical trial.

Molnupiravir is an antiviral, currently approved by the UK Medicines and Healthcare products Regulatory Agency (MHRA) for treating at-risk COVID-19 patients, that induces lethal error catastrophe in SARS-CoV-2. How this drug-induced mechanism of action might impact the emergence of resistance mutations is unclear. To investigate this, we used samples from the AGILE Candidate Specific Trial (CST)-2 (clinical trial number NCT04746183). The primary outcomes of AGILE CST-2 were to measure the drug safety and antiviral efficacy of molnupiravir in humans (180 participants randomised 1:1 with placebo). Here, we describe the pre-specified exploratory virological endpoint of CST-2, which was to determine the possible genomic changes in SARS-CoV-2 induced by molnupiravir treatment. We use high-throughput amplicon sequencing and minor variant analysis to characterise viral genomics in each participant whose longitudinal samples (days 1, 3 and 5 post-randomisation) pass the viral genomic quality criteria (n = 59 for molnupiravir and n = 65 for placebo). Over the course of treatment, no specific mutations were associated with molnupiravir treatment. We find that molnupiravir significantly increased the transition:transversion mutation ratio in SARS-CoV-2, consistent with the model of lethal error catastrophe. This study highlights the utility of examining intra-host virus populations to strengthen the prediction, and surveillance, of potential treatment-emergent adaptations