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

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

Sequencing data and code - Prince et al 2023

&lt;p&gt;Raw sequencing data, code and intermediate analysis files from &quot;Antiviral activity of molnupiravir precursor NHC against SARS-CoV-2 Variants of Concern (VOCs) and implications for the therapeutic window and resistance&quot; (Prince et al, 2023).&lt;/p&gt; &lt;p&gt;Please see sample_metadata.xlsx for all metadata relating to the files contained in this repository.&lt;/p&gt; &lt;p&gt;Code for data visualisation can be found in: mut_sub_all_pts_serial-pass.R. Specific paths to data will have to be changes to refer to where you have downloaded the data in this repository. Metadata for use with the R script is&nbsp;serial-pass-nimagen-metadata-forR.csv.&lt;/p&gt

'Queer in Microbiology': a Microbiology Society members' endeavour for creating a safe and inclusive environment for LGBTQIA+ microbiologists.

The past decade has seen growing awareness of the challenges faced by LGBTQIA+ scientists, including discrimination in the workplace and the lack of representation. Initiatives such as 500 Queer Scientists, Pride in STEM and the Microbiology Society's LGBTQIA+ events have been instrumental in promoting inclusivity in science, technology, engineering, mathematics and medicine (STEMM). The Microbiology Society and its members have played a pivotal role in these efforts and summarized here are their initiatives towards safer and more inclusive scientific and research environments. Starting with a series of interviews and blog posts about the experiences of LGBTQIA+ microbiologists in research, the Society has promoted the organization of networking and social events and developed guidelines for creating more inclusive scientific conferences. These initiatives have not only improved the representation and visibility of LGBTQIA+ individuals in microbiology, but have also served as a blueprint for similar efforts in other scientific areas. Nevertheless, despite improvements in some areas, full inclusion of LGBTQIA+ scientists is still hindered by societal and institutional policies around the world. Here, we propose novel measures to support and empower LGBTQIA+ microbiological communities within learned societies

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

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

SARS-CoV-2 population dynamics in immunocompetent individuals in a closed transmission chain shows genomic diversity over the course of infection.

BackgroundSARS-CoV-2 remains rapidly evolving, and many biologically important genomic substitutions/indels have characterised novel SARS-CoV-2 lineages, which have emerged during successive global waves of the pandemic. Worldwide genomic sequencing has been able to monitor these waves, track transmission clusters, and examine viral evolution in real time to help inform healthcare policy. One school of thought is that an apparent greater than average divergence in an emerging lineage from contemporary variants may require persistent infection, for example in an immunocompromised host. Due to the nature of the COVID-19 pandemic and sampling, there were few studies that examined the evolutionary trajectory of SARS-CoV-2 in healthy individuals.MethodsWe investigated viral evolutionary trends and participant symptomatology within a cluster of 16 SARS-CoV-2 infected, immunocompetent individuals with no co-morbidities in a closed transmission chain. Longitudinal nasopharyngeal swab sampling allowed characterisation of SARS-CoV-2 intra-host variation over time at both the dominant and minor genomic variant levels through Nimagen-Illumina sequencing.ResultsA change in viral lineage assignment was observed in individual infections; however, there was only one indel and no evidence of recombination over the period of an acute infection. Minor and dominant genomic modifications varied between participants, with some minor genomic modifications increasing in abundance to become the dominant viral sequence during infection.ConclusionsData from this cohort of SARS-CoV-2-infected participants demonstrated that long-term persistent infection in an immunocompromised host was not necessarily a prerequisite for generating a greater than average frequency of amino acid substitutions. Amino acid substitutions at both the dominant and minor genomic sequence level were observed in immunocompetent individuals during infection showing that viral lineage changes can occur generating viral diversity

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