Tiled-ClickSeq for targeted sequencing of complete coronavirus genomes with simultaneous capture of RNA recombination and minority variants

RpoS activates formation of Salmonella Typhi biofilms and drives persistence in the gall bladder

The development of strategies for targeting the asymptomatic carriage of Salmonella Typhi in chronic typhoid patients has suffered owing to our basic lack of understanding of the molecular mechanisms that enable the formation of S. Typhi biofilms. Traditionally, studies have relied on cholesterol-attached biofilms formed by a closely related serovar, Typhimurium, to mimic multicellular Typhi communities formed on human gallstones. In long-term infections, S. Typhi adopts the biofilm lifestyle to persist in vivo and survive in the carrier state, ultimately leading to the spread of infections via the fecal-oral route of transmission. In the present work, we studied S. Typhi biofilms directly, applied targeted as well as genome-wide genetic approaches to uncover unique biofilm components that do not conform to the CsgD-dependent pathway established in S. Typhimurium. We undertook a genome-wide Tn5 mutation screen in a highly successful parental lineage of S. Typhi, strain H58, in gallstone-mimicking conditions. We generated New Generation Sequencing libraries based on the ClickSeq technology to identify the key regulators, IraP and RpoS, and the matrix components as Sth fimbriae, Vi capsule and lipopolysaccharide. We discovered that the starvation sigma factor, RpoS, was required for the transcriptional activation of matrix-encoding genes in vitro, and for S. Typhi colonization in persistent infections in vivo, using a heterologous fish larval model. An rpoS null mutant failed to colonize the gall bladder in chronic zebrafish infections. Overall, our work uncovered a novel RpoS-driven paradigm for the formation of cholesterol-attached Typhi biofilms, and emphasized the role(s) of stress signaling pathways for adaptation in chronic infections. Our identification of the biofilm regulators in S. Typhi paves the way for the development of drugs against typhoid carriage, which will ultimately control the increased incidence of gall bladder cancer in typhoid carriers.

Characterization of SARS-CoV-2 genetic evolution in vaccinated and non-vaccinated patients from the Kenyan population

Vaccination is a key control measure of COVID-19 by preventing severe effects of disease outcomes, reducing hospitalization rates and death, and increasing herd immunity. However, vaccination can affect the evolution and adaptation of SARS-CoV-2, largely through vaccine-induced immune pressure. Here we investigated the recombination events and single nucleotide polymorphisms (SNPs) on SARS-CoV-2 genome in non-vaccinated and vaccinated patients in the Kenyan population. We identified recombination hotspots in the S, N, and ORF1a/b genes and showed the genetic evolution landscape of SARS-CoV-2 by comparing within-wave and inter-wave recombination events from the beginning of the pandemic (June 2020) to (October 2022) in Kenya. An in-depth analysis of (SNPs) on the S, ORf1a/b, and N genes identified previously unreported mutations. We detected a minority variant in non-vaccinated patients in Kenya, that contained immune escape mutation S255F of the spike gene and showing a differential recombination pattern within the non-vaccinated patients. Detailed analysis of recombination between waves suggested an association between increased population immunity and declining risk of emergence of variants of concern. Overall, this work identified unique mutations in SARS-CoV-2 which could have significant implications for virus evolution, virulence, and immune escape.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

Generation and Functional Analysis of Defective Viral Genomes during SARS-CoV-2 Infection

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from transcriptome sequencing (RNA-seq) data sets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hot spots were identified for DVG recombination, and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single-cell RNA-seq analysis indicated the interferon (IFN) stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the next-generation sequencing (NGS) data set from a published cohort study and observed a significantly higher amount and frequency of DVG in symptomatic patients than those in asymptomatic patients. Finally, we observed exceptionally diverse DVG populations in one immunosuppressive patient up to 140 days after the first positive test of COVID-19, suggesting for the first time an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and into how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection. IMPORTANCE Defective viral genomes (DVGs) are generated ubiquitously in many RNA viruses, including SARS-CoV-2. Their interference activity to full-length viruses and IFN stimulation provide the potential for them to be used in novel antiviral therapies and vaccine development. SARS-CoV-2 DVGs are generated through the recombination of two discontinuous genomic fragments by viral polymerase complex, and this recombination is also one of the major mechanisms for the emergence of new coronaviruses. Focusing on the generation and function of SARS-CoV-2 DVGs, these studies identify new hot spots for nonhomologous recombination and strongly suggest that the secondary structures within viral genomes mediate the recombination. Furthermore, these studies provide the first evidence for IFN stimulation activity of de novo DVGs during natural SARS-CoV-2 infection. These findings set up the foundation for further mechanism studies of SARS-CoV-2 recombination and provide evidence to harness the immunostimulatory potential of DVGs in the development of a vaccine and antivirals for SARS-CoV-2.

ViReMa: a virus recombination mapper of next-generation sequencing data characterizes diverse recombinant viral nucleic acids

BACKGROUND

Genetic recombination is a tremendous source of intrahost diversity in viruses and is critical for their ability to rapidly adapt to new environments or fitness challenges. While viruses are routinely characterized using high-throughput sequencing techniques, characterizing the genetic products of recombination in next-generation sequencing data remains a challenge. Viral recombination events can be highly diverse and variable in nature, including simple duplications and deletions, or more complex events such as copy/snap-back recombination, intervirus or intersegment recombination, and insertions of host nucleic acids. Due to the variable mechanisms driving virus recombination and the different selection pressures acting on the progeny, recombination junctions rarely adhere to simple canonical sites or sequences. Furthermore, numerous different events may be present simultaneously in a viral population, yielding a complex mutational landscape.

FINDINGS

We have previously developed an algorithm called ViReMa (Virus Recombination Mapper) that bootstraps the bowtie short-read aligner to capture and annotate a wide range of recombinant species found within virus populations. Here, we have updated ViReMa to provide an "error density" function designed to accurately detect recombination events in the longer reads now routinely generated by the Illumina platforms and provide output reports for multiple types of recombinant species using standardized formats. We demonstrate the utility and flexibility of ViReMa in different settings to report deletion events in simulated data from Flock House virus, copy-back RNA species in Sendai viruses, short duplication events in HIV, and virus-to-host recombination in an archaeal DNA virus.

Generation and functional analysis of defective viral genomes during SARS-CoV-2 infection.. [PMID: 36172120 PMCID: PMC9516852 DOI: 10.1101/2022.09.22.509123] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Defective viral genomes (DVGs) have been identified in many RNA viruses as a major factor influencing antiviral immune response and viral pathogenesis. However, the generation and function of DVGs in SARS-CoV-2 infection are less known. In this study, we elucidated DVG generation in SARS-CoV-2 and its relationship with host antiviral immune response. We observed DVGs ubiquitously from RNA-seq datasets of in vitro infections and autopsy lung tissues of COVID-19 patients. Four genomic hotspots were identified for DVG recombination and RNA secondary structures were suggested to mediate DVG formation. Functionally, bulk and single cell RNA-seq analysis indicated the IFN stimulation of SARS-CoV-2 DVGs. We further applied our criteria to the NGS dataset from a published cohort study and observed significantly higher DVG amount and frequency in symptomatic patients than that in asymptomatic patients. Finally, we observed unusually high DVG frequency in one immunosuppressive patient up to 140 days after admitted to hospital due to COVID-19, first-time suggesting an association between DVGs and persistent viral infections in SARS-CoV-2. Together, our findings strongly suggest a critical role of DVGs in modulating host IFN responses and symptom development, calling for further inquiry into the mechanisms of DVG generation and how DVGs modulate host responses and infection outcome during SARS-CoV-2 infection.

QTQTN motif upstream of the furin-cleavage site plays a key role in SARS-CoV-2 infection and pathogenesis

The furin cleavage site (FCS), an unusual feature in the SARS-CoV-2 spike protein, has been spotlighted as a factor key to facilitating infection and pathogenesis by increasing spike processing. Similarly, the QTQTN motif directly upstream of the FCS is also an unusual feature for group 2B coronaviruses (CoVs). The QTQTN deletion has consistently been observed in in vitro cultured virus stocks and some clinical isolates. To determine whether the QTQTN motif is critical to SARS-CoV-2 replication and pathogenesis, we generated a mutant deleting the QTQTN motif (ΔQTQTN). Here, we report that the QTQTN deletion attenuates viral replication in respiratory cells in vitro and attenuates disease in vivo. The deletion results in a shortened, more rigid peptide loop that contains the FCS and is less accessible to host proteases, such as TMPRSS2. Thus, the deletion reduced the efficiency of spike processing and attenuates SARS-CoV-2 infection. Importantly, the QTQTN motif also contains residues that are glycosylated, and disruption of its glycosylation also attenuates virus replication in a TMPRSS2-dependent manner. Together, our results reveal that three aspects of the S1/S2 cleavage site-the FCS, loop length, and glycosylation-are required for efficient SARS-CoV-2 replication and pathogenesis.

ViReMaShiny: an interactive application for analysis of viral recombination data

MOTIVATION

Recombination is an essential driver of virus evolution and adaption, giving rise to new chimeric viruses, structural variants, sub-genomic RNAs and defective RNAs. Next-generation sequencing (NGS) of virus samples, either from experimental or clinical settings, has revealed a complex distribution of recombination events that contributes to intrahost diversity. We and others have previously developed alignment tools to discover and map these diverse recombination events in NGS data. However, there is no standard for data visualization to contextualize events of interest, and downstream analysis often requires bespoke coding.

RESULTS

We present ViReMaShiny, a web-based application built using the R Shiny framework to allow interactive exploration and point-and-click visualization of viral recombination data provided in BED format generated by computational pipelines such as ViReMa (Viral-Recombination-Mapper).

AVAILABILITY AND IMPLEMENTATION

The application is hosted at https://routhlab.shinyapps.io/ViReMaShiny/ with associated documentation at https://jayeung12.github.io/. Code is available at https://github.com/routhlab/ViReMaShiny.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203-205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral 'RG' motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2's continued adaptation to human infection.

DVGfinder: A Metasearch Tool for Identifying Defective Viral Genomes in RNA-Seq Data

The generation of different types of defective viral genomes (DVG) is an unavoidable consequence of the error-prone replication of RNA viruses. In recent years, a particular class of DVGs, those containing long deletions or genome rearrangements, has gain interest due to their potential therapeutic and biotechnological applications. Identifying such DVGs in high-throughput sequencing (HTS) data has become an interesting computational problem. Several algorithms have been proposed to accomplish this goal, though all incur false positives, a problem of practical interest if such DVGs have to be synthetized and tested in the laboratory. We present a metasearch tool, DVGfinder, that wraps the two most commonly used DVG search algorithms in a single workflow for the identification of the DVGs in HTS data. DVGfinder processes the results of ViReMa-a and DI-tector and uses a gradient boosting classifier machine learning algorithm to reduce the number of false-positive events. The program also generates output files in user-friendly HTML format, which can help users to explore the DVGs identified in the sample. We evaluated the performance of DVGfinder compared to the two search algorithms used separately and found that it slightly improves sensitivities for low-coverage synthetic HTS data and DI-tector precision for high-coverage samples. The metasearch program also showed higher sensitivity on a real sample for which a set of copy-backs were previously validated.

Recombination in Positive-Strand RNA Viruses

RNA recombination is a major driver of genetic shifts tightly linked to the evolution of RNA viruses. Genomic recombination contributes substantially to the emergence of new viral lineages, expansion in host tropism, adaptations to new environments, and virulence and pathogenesis. Here, we review some of the recent progress that has advanced our understanding of recombination in positive-strand RNA viruses, including recombination triggers and the mechanisms behind them. The study of RNA recombination aids in predicting the probability and outcome of viral recombination events, and in the design of viruses with reduced recombination frequency as candidates for the development of live attenuated vaccines. Surveillance of viral recombination should remain a priority in the detection of emergent viral strains, a goal that can only be accomplished by expanding our understanding of how these events are triggered and regulated.

Next-generation sequencing: A new avenue to understand viral RNA-protein interactions

The genomes of RNA viruses present an astonishing source of both sequence and structural diversity. From intracellular viral RNA-host interfaces to interactions between the RNA genome and structural proteins in virus particles themselves, almost the entire viral lifecycle is accompanied by a myriad of RNA-protein interactions that are required to fulfill their replicative potential. It is therefore important to characterize such rich and dynamic collections of viral RNA-protein interactions to understand virus evolution and their adaptation to their hosts and environment. Recent advances in next-generation sequencing technologies have allowed the characterization of viral RNA-protein interactions, including both transient and conserved interactions, where molecular and structural approaches have fallen short. In this review, we will provide a methodological overview of the high-throughput techniques used to study viral RNA-protein interactions, their biochemical mechanisms, and how they evolved from classical methods as well as one another. We will discuss how different techniques have fueled virus research to characterize how viral RNA and proteins interact, both locally and on a global scale. Finally, we will present examples on how these techniques influence the studies of clinically important pathogens such as HIV-1 and SARS-CoV-2.

A Comparison of Seegene Technologies Novaplex SARS-CoV-2 Variants I, II, and IV Assays with Spike Gene Sequencing for Detection of Known Severe Acute Respiratory Syndrome Coronavirus 2 Variants

Tracking new and emerging severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants has become increasingly important for public health responses, primarily because of variant-dependent transmission, disease severity, and treatment decisions. This evaluation compared Seegene Technologies Novaplex SARS-CoV-2 Variants I, II, and IV (I,II&IV) assays to detect known SARS-CoV-2 variants using traditional spike gene Sanger sequencing results as the gold standard reference. Both RNA extraction and extraction-free protocols were assessed. A total of 156 samples were included in this study. There was 100% (109/109) overall agreement (95% CI, 96.7%-100%) between the spike gene sequencing and the I,II&IV results using extracted RNA for the variants included in the Novaplex assay menus. The RNA extraction-free method was 91.7% (143/156) as sensitive (95% CI, 86.2%-95.5%) as the traditional RNA extraction method. Using the extraction-free method on samples with higher cycle threshold values (>30) resulted in some mutations not being detected, presumably due to lower nucleic acid concentrations in the original samples. In conclusion, the I,II&IV assays provide an accurate, rapid, and less labor-intensive method for detecting SARS-CoV-2 and identifying known variants of interest and concern. The RNA extraction-free method for samples with cycle threshold of <30 could be cost-effective for surveillance purposes. However, spike gene sequencing retains the advantage of detecting more and new variants.

Mouse-adapted SARS-CoV-2 protects animals from lethal SARS-CoV challenge

The emergence of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has resulted in a pandemic causing significant damage to public health and the economy. Efforts to understand the mechanisms of Coronavirus Disease 2019 (COVID-19) have been hampered by the lack of robust mouse models. To overcome this barrier, we used a reverse genetic system to generate a mouse-adapted strain of SARS-CoV-2. Incorporating key mutations found in SARS-CoV-2 variants, this model recapitulates critical elements of human infection including viral replication in the lung, immune cell infiltration, and significant in vivo disease. Importantly, mouse adaptation of SARS-CoV-2 does not impair replication in human airway cells and maintains antigenicity similar to human SARS-CoV-2 strains. Coupled with the incorporation of mutations found in variants of concern, CMA3p20 offers several advantages over other mouse-adapted SARS-CoV-2 strains. Using this model, we demonstrate that SARS-CoV-2-infected mice are protected from lethal challenge with the original Severe Acute Respiratory Syndrome Coronavirus (SARS-CoV), suggesting immunity from heterologous Coronavirus (CoV) strains. Together, the results highlight the use of this mouse model for further study of SARS-CoV-2 infection and disease.

Nucleocapsid mutations in SARS-CoV-2 augment replication and pathogenesis

While SARS-CoV-2 continues to adapt for human infection and transmission, genetic variation outside of the spike gene remains largely unexplored. This study investigates a highly variable region at residues 203–205 in the SARS-CoV-2 nucleocapsid protein. Recreating a mutation found in the alpha and omicron variants in an early pandemic (WA-1) background, we find that the R203K+G204R mutation is sufficient to enhance replication, fitness, and pathogenesis of SARS-CoV-2. The R203K+G204R mutant corresponds with increased viral RNA and protein both in vitro and in vivo. Importantly, the R203K+G204R mutation increases nucleocapsid phosphorylation and confers resistance to inhibition of the GSK-3 kinase, providing a molecular basis for increased virus replication. Notably, analogous alanine substitutions at positions 203+204 also increase SARS-CoV-2 replication and augment phosphorylation, suggesting that infection is enhanced through ablation of the ancestral ‘RG’ motif. Overall, these results demonstrate that variant mutations outside spike are key components in SARS-CoV-2’s continued adaptation to human infection.

SARS-CoV-2 Subgenomic RNAs: Characterization, Utility, and Perspectives

SARS-CoV-2, the etiologic agent at the root of the ongoing COVID-19 pandemic, harbors a large RNA genome from which a tiered ensemble of subgenomic RNAs (sgRNAs) is generated. Comprehensive definition and investigation of these RNA products are important for understanding SARS-CoV-2 pathogenesis. This review summarizes the recent progress on SARS-CoV-2 sgRNA identification, characterization, and application as a viral replication marker. The significance of these findings and potential future research areas of interest are discussed.