Nucleotide-resolution profiling of RNA recombination in the encapsidated genome of a eukaryotic RNA virus by next-generation sequencing

Hecatomb: an integrated software platform for viral metagenomics

BACKGROUND

Modern sequencing technologies offer extraordinary opportunities for virus discovery and virome analysis. Annotation of viral sequences from metagenomic data requires a complex series of steps to ensure accurate annotation of individual reads and assembled contigs. In addition, varying study designs will require project-specific statistical analyses.

FINDINGS

Here we introduce Hecatomb, a bioinformatic platform coordinating commonly used tasks required for virome analysis. Hecatomb means "a great sacrifice." In this setting, Hecatomb is "sacrificing" false-positive viral annotations using extensive quality control and tiered-database searches. Hecatomb processes metagenomic data obtained from both short- and long-read sequencing technologies, providing annotations to individual sequences and assembled contigs. Results are provided in commonly used data formats useful for downstream analysis. Here we demonstrate the functionality of Hecatomb through the reanalysis of a primate enteric and a novel coral reef virome.

CONCLUSION

Hecatomb provides an integrated platform to manage many commonly used steps for virome characterization, including rigorous quality control, host removal, and both read- and contig-based analysis. Each step is managed using the Snakemake workflow manager with dependency management using Conda. Hecatomb outputs several tables properly formatted for immediate use within popular data analysis and visualization tools, enabling effective data interpretation for a variety of study designs. Hecatomb is hosted on GitHub (github.com/shandley/hecatomb) and is available for installation from Bioconda and PyPI.

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.

Covariation of viral recombination with single nucleotide variants during virus evolution revealed by CoVaMa

Adaptation of viruses to their environments occurs through the acquisition of both novel single-nucleotide variants (SNV) and recombination events including insertions, deletions, and duplications. The co-occurrence of SNVs in individual viral genomes during their evolution has been well-described. However, unlike covariation of SNVs, studying the correlation between recombination events with each other or with SNVs has been hampered by their inherent genetic complexity and a lack of bioinformatic tools. Here, we expanded our previously reported CoVaMa pipeline (v0.1) to measure linkage disequilibrium between recombination events and SNVs within both short-read and long-read sequencing datasets. We demonstrate this approach using long-read nanopore sequencing data acquired from Flock House virus (FHV) serially passaged in vitro. We found SNVs that were either correlated or anti-correlated with large genomic deletions generated by nonhomologous recombination that give rise to Defective-RNAs. We also analyzed NGS data from longitudinal HIV samples derived from a patient undergoing antiretroviral therapy who proceeded to virological failure. We found correlations between insertions in the p6Gag and mutations in Gag cleavage sites. This report confirms previous findings and provides insights on novel associations between SNVs and specific recombination events within the viral genome and their role in viral evolution.

Defective Interfering Particles of Negative-Strand RNA Viruses

Viral defective interfering particles (DIPs) were intensely studied several decades ago but research waned leaving open many critical questions. New technologies and other advances led to a resurgence in DIP studies for negative-strand RNA viruses. While DIPs have long been recognized, their exact contribution to the outcome of acute or persistent viral infections has remained elusive. Recent studies have identified defective viral genomes (DVGs) in human infections, including respiratory syncytial virus and influenza, and growing evidence indicates that DVGs influence disease severity and may contribute to viral persistence. Further, several studies have advanced our understanding of key viral and host factors that regulate DIP formation and activity. Here we review these discoveries and highlight key questions moving forward.

ClickSeq: Replacing Fragmentation and Enzymatic Ligation with Click-Chemistry to Prevent Sequence Chimeras

We recently reported a fragmentation-free method for the synthesis of Next-Generation Sequencing libraries called "ClickSeq" that uses biorthogonal click-chemistry in place of enzymes for the ligation of sequencing adaptors. We found that this approach dramatically reduces artifactual chimera formation, allowing the study of rare recombination events that include viral replication intermediates and defective-interfering viral RNAs. ClickSeq illustrates how robust, bio-orthogonal chemistry can be harnessed in vitro to capture and dissect complex biological processes. Here, we describe an updated protocol for the synthesis of "ClickSeq" libraries.

Next-generation sequencing library preparation method for identification of RNA viruses on the Ion Torrent Sequencing Platform

Next generation sequencing (NGS) is a powerful tool for the characterization, discovery, and molecular identification of RNA viruses. There were multiple NGS library preparation methods published for strand-specific RNA-seq, but some methods are not suitable for identifying and characterizing RNA viruses. In this study, we report a NGS library preparation method to identify RNA viruses using the Ion Torrent PGM platform. The NGS sequencing adapters were directly inserted into the sequencing library through reverse transcription and polymerase chain reaction, without fragmentation and ligation of nucleic acids. The results show that this method is simple to perform, able to identify multiple species of RNA viruses in clinical samples.

Identification of Interactions between Sindbis Virus Capsid Protein and Cytoplasmic vRNA as Novel Virulence Determinants

Alphaviruses are arthropod-borne viruses that represent a significant threat to public health at a global level. While the formation of alphaviral nucleocapsid cores, consisting of cargo nucleic acid and the viral capsid protein, is an essential molecular process of infection, the precise interactions between the two partners are ill-defined. A CLIP-seq approach was used to screen for candidate sites of interaction between the viral Capsid protein and genomic RNA of Sindbis virus (SINV), a model alphavirus. The data presented in this report indicates that the SINV capsid protein binds to specific viral RNA sequences in the cytoplasm of infected cells, but its interaction with genomic RNA in mature extracellular viral particles is largely non-specific in terms of nucleotide sequence. Mutational analyses of the cytoplasmic viral RNA-capsid interaction sites revealed a functional role for capsid binding early in infection. Interaction site mutants exhibited decreased viral growth kinetics; however, this defect was not a function of decreased particle production. Rather mutation of the cytoplasmic capsid-RNA interaction sites negatively affected the functional capacity of the incoming viral genomic RNAs leading to decreased infectivity. Furthermore, cytoplasmic capsid interaction site mutants are attenuated in a murine model of neurotropic alphavirus infection. Collectively, the findings of this study indicate that the identified cytoplasmic interactions of the viral capsid protein and genomic RNA, while not essential for particle formation, are necessary for genomic RNA function early during infection. This previously unappreciated role of capsid protein during the alphaviral replication cycle also constitutes a novel virulence determinant.

Alphaviruses can cause significant disease in infected individuals; however, our understanding of the molecular interactions that enable infection and contribute to the development of disease is limited. The work detailed in this manuscript characterizes the interaction of a viral RNA-binding protein, Capsid, with the viral genomic RNA. Importantly, these interactions were found to be at specific sites on the genome but not essential for virus assembly. Mutation of the capsid / RNA interaction sites decreased the replication of the virus and the severity of disease in a mouse model of infection. Taken together, these findings identify a previously undiscovered determinant of disease severity, and provide a potential basis for the development of new vaccines.

Parallel ClickSeq and Nanopore sequencing elucidates the rapid evolution of defective-interfering RNAs in Flock House virus

Defective-Interfering RNAs (DI-RNAs) have long been known to play an important role in virus replication and transmission. DI-RNAs emerge during virus passaging in both cell-culture and their hosts as a result of non-homologous RNA recombination. However, the principles of DI-RNA emergence and their subsequent evolution have remained elusive. Using a combination of long- and short-read Next-Generation Sequencing, we have characterized the formation of DI-RNAs during serial passaging of Flock House virus (FHV) in cell-culture over a period of 30 days in order to elucidate the pathways and potential mechanisms of DI-RNA emergence and evolution. For short-read RNAseq, we employed 'ClickSeq' due to its ability to sensitively and confidently detect RNA recombination events with nucleotide resolution. In parallel, we used the Oxford Nanopore Technologies's (ONT) MinION to resolve full-length defective and wild-type viral genomes. Together, these accurately resolve both rare and common RNA recombination events, determine the correlation between recombination events, and quantifies the relative abundance of different DI-RNAs throughout passaging. We observe the formation of a diverse pool of defective RNAs at each stage of viral passaging. However, many of these 'intermediate' species, while present in early stages of passaging, do not accumulate. After approximately 9 days of passaging we observe the rapid accumulation of DI-RNAs with a correlated reduction in specific infectivity and with the Nanopore data find that DI-RNAs are characterized by multiple RNA recombination events. This suggests that intermediate DI-RNA species are not competitive and that multiple recombination events interact epistatically to confer 'mature' DI-RNAs with their selective advantage allowing for their rapid accumulation. Alternatively, it is possible that mature DI-RNA species are generated in a single event involving multiple RNA rearrangements. These insights have important consequences for our understanding of the mechanisms, determinants and limitations in the emergence and evolution of DI-RNAs.

Host-specific accumulation and temperature effects on the generation of dimeric viral RNA species derived from the S-RNA of members of the Tospovirus genus

Polygonum ringspot virus (PolRSV) is a recently characterized Tospovirus reported in Italy. Northern blot analyses of PolRSV infections in Nicotiana benthamiana and tomato plants showed that a viral RNA species with nearly twice the length of the Small genomic RNA (S-RNA) accumulated abundantly in the former host, but was not detected in the latter. Additional assays confirmed that biogenesis of this novel RNA species was common to all PolRSV isolates tested and also to an isolate of Tomato spotted wilt virus (TSWV). Given its size, we hypothesized that the novel RNA species was a dimer molecule and we confirmed this hypothesis by RNA sequencing (RNAseq) analysis and reverse transcription (RT)-PCR of putative predicted dimer junction sites in RNA extracts of N. benthamiana challenged with PolRSV isolates Plg6 and Plg13/2. We also confirmed that these molecules are derived from head-to-tail dimers and often contain deletions at their junction sites. We named these novel molecules imperfect dimer RNAs (IMPD-RNAs). PolRSV IMPD-RNAs systemic accumulation in a range of host plants was restricted to N. benthamiana and Nicotiana occidentalis. Notably, IMPD-RNAs accumulation was modulated by temperature and their generation was restricted to late stages of systemic infection (12 days post-inoculation) in N. benthamiana. Differently from all other PolRSV isolates used in this study, Plg13/2 generated more IMPD-RNAs coupled with low amounts of genomic S-RNA and maintained them even at 18 °C, besides having lost the ability to infect tomato plants. This is the first characterization of S-RNA dimers for Tospovirus, and of occurrence of dimers of genomic segments at the whole organism level for Bunyaviridae.

CoVaMa: Co-Variation Mapper for disequilibrium analysis of mutant loci in viral populations using next-generation sequence data

Next-Generation Sequencing (NGS) has transformed our understanding of the dynamics and diversity of virus populations for human pathogens and model systems alike. Due to the sensitivity and depth of coverage in NGS, it is possible to measure the frequency of mutations that may be present even at vanishingly low frequencies within the viral population. Here, we describe a simple bioinformatic pipeline called CoVaMa (Co-Variation Mapper) scripted in Python that detects correlated patterns of mutations in a viral sample. Our algorithm takes NGS alignment data and populates large matrices of contingency tables that correspond to every possible pairwise interaction of nucleotides in the viral genome or amino acids in the chosen open reading frame. These tables are then analysed using classical linkage disequilibrium to detect and report evidence of epistasis. We test our analysis with simulated data and then apply the approach to find epistatically linked loci in Flock House Virus genomic RNA grown under controlled cell culture conditions. We also reanalyze NGS data from a large cohort of HIV infected patients and find correlated amino acid substitution events in the protease gene that have arisen in response to anti-viral therapy. This both confirms previous findings and suggests new pairs of interactions within HIV protease. The script is publically available at http://sourceforge.net/projects/covama.

Artifact-Free Quantification and Sequencing of Rare Recombinant Viruses by Using Drop-Based Microfluidics

Recombination is an important driver in the evolution of viruses and thus is key to understanding viral epidemics and improving strategies to prevent future outbreaks. Characterization of rare recombinant subpopulations remains technically challenging because of artifacts such as artificial recombinants, known as chimeras, and amplification bias. To overcome this, we have developed a high-throughput microfluidic technique with a second verification step in order to amplify and sequence single recombinant viruses with high fidelity in picoliter drops. We obtained the first artifact-free estimate of in vitro recombination rate between murine norovirus strains MNV-1 and WU20 co-infecting a cell (P(rec) = 3.3 × 10(-4) ± 2 × 10(-5) ) for a 1205 nt region. Our approach represents a time- and cost-effective improvement over current methods, and can be adapted for genomic studies requiring artifact- and bias-free selective amplification, such as microbial pathogens, or rare cancer cells.

ClickSeq: Fragmentation-Free Next-Generation Sequencing via Click Ligation of Adaptors to Stochastically Terminated 3'-Azido cDNAs

We present a simple method called "ClickSeq" for NGS (next-generation sequencing) library synthesis that uses click chemistry rather than enzymatic reactions for the ligation of Illumina sequencing adaptors. In ClickSeq, randomly primed reverse transcription reactions are supplemented with azido-2',3'-dideoxynucleotides that randomly terminate DNA synthesis and release 3'-azido-blocked cDNA fragments in a process akin to dideoxy-Sanger sequencing. Purified fragments are "click ligated" via copper-catalyzed alkyne-azide cycloaddition to DNA oligos modified with a 5'-alkyne group. This generates ssDNA molecules containing an unnatural triazole-linked DNA backbone that is sufficiently biocompatible for PCR amplification to generate a cDNA library for RNAseq. Here, we analyze viral RNAs and mRNA to demonstrate that ClickSeq produces unbiased NGS libraries with low error rates comparable to standard methods. Importantly, ClickSeq is robust against common artifacts of NGS such as chimera formation and artifactual recombination with fewer than 3 aberrant events detected per million reads.

Deep sequencing of virus-derived small interfering RNAs and RNA from viral particles shows highly similar mutational landscapes of a plant virus population

RNA viruses exist within a host as a population of mutant sequences, often referred to as quasispecies. Within a host, sequences of RNA viruses constitute several distinct but interconnected pools, such as RNA packed in viral particles, double-stranded RNA, and virus-derived small interfering RNAs. We aimed to test if the same representation of within-host viral population structure could be obtained by sequencing different viral sequence pools. Using ultradeep Illumina sequencing, the diversity of two coexisting Potato virus Y sequence pools present within a plant was investigated: RNA isolated from viral particles and virus-derived small interfering RNAs (the derivatives of a plant RNA silencing mechanism). The mutational landscape of the within-host virus population was highly similar between both pools, with no notable hotspots across the viral genome. Notably, all of the single-nucleotide polymorphisms with a frequency of higher than 1.6% were found in both pools. Some unique single-nucleotide polymorphisms (SNPs) with very low frequencies were found in each of the pools, with more of them occurring in the small RNA (sRNA) pool, possibly arising through genetic drift in localized virus populations within a plant and the errors introduced during the amplification of silencing signal. Sequencing of the viral particle pool enhanced the efficiency of consensus viral genome sequence reconstruction. Nonhomologous recombinations were commonly detected in the viral particle pool, with a hot spot in the 3' untranslated and coat protein regions of the genome. We stress that they present an important but often overlooked aspect of virus population diversity.

IMPORTANCE

This study is the most comprehensive whole-genome characterization of a within-plant virus population to date and the first study comparing diversity of different pools of viral sequences within a host. We show that both virus-derived small RNAs and RNA from viral particles could be used for diversity assessment of within-plant virus population, since they show a highly congruent portrayal of the virus mutational landscape within a plant. The study is an important baseline for future studies of virus population dynamics, for example, during the adaptation to a new host. The comparison of the two virus sequence enrichment techniques, sequencing of virus-derived small interfering RNAs and RNA from purified viral particles, shows the strength of the latter for the detection of recombinant viral genomes and reconstruction of complete consensus viral genome sequence.

Discovery of functional genomic motifs in viruses with ViReMa-a Virus Recombination Mapper-for analysis of next-generation sequencing data

We developed an algorithm named ViReMa (Viral-Recombination-Mapper) to provide a versatile platform for rapid, sensitive and nucleotide-resolution detection of recombination junctions in viral genomes using next-generation sequencing data. Rather than mapping read segments of pre-defined lengths and positions, ViReMa dynamically generates moving read segments. ViReMa initially attempts to align the 5' end of a read to the reference genome(s) with the Bowtie seed-based alignment. A new read segment is then made by either extracting any unaligned nucleotides at the 3' end of the read or by trimming the first nucleotide from the read. This continues iteratively until all portions of the read are either mapped or trimmed. With multiple reference genomes, it is possible to detect virus-to-host or inter-virus recombination. ViReMa is also capable of detecting insertion and substitution events and multiple recombination junctions within a single read. By mapping the distribution of recombination events in the genome of flock house virus, we demonstrate that this information can be used to discover de novo functional motifs located in conserved regions of the viral genome.

Identification and survey of a novel avian coronavirus in ducks

The rapid discovery of novel viruses using next generation sequencing (NGS) technologies including DNA-Seq and RNA-Seq, has greatly expanded our understanding of viral diversity in recent years. The timely identification of novel viruses using NGS technologies is also important for us to control emerging infectious diseases caused by novel viruses. In this study, we identified a novel duck coronavirus (CoV), distinct with chicken infectious bronchitis virus (IBV), using RNA-Seq. The novel duck-specific CoV was a potential novel species within the genus Gammacoronavirus, as indicated by sequences of three regions in the viral 1b gene. We also performed a survey of CoVs in domestic fowls in China using reverse-transcription polymerase chain reaction (RT-PCR), targeting the viral nucleocapsid (N) gene. A total of 102 CoV positives were identified through the survey. Phylogenetic analysis of the viral N sequences suggested that CoVs in domestic fowls have diverged into several region-specific or host-specific clades or subclades in the world, and IBVs can infect ducks, geese and pigeons, although they mainly circulate in chickens. Moreover, this study provided novel data supporting the notion that some host-specific CoVs other than IBVs circulate in ducks, geese and pigeons, and indicated that the novel duck-specific CoV identified through RNA-Seq in this study is genetically closer to some CoVs circulating in wild water fowls. Taken together, this study shed new insight into the diversity, distribution, evolution and control of avian CoVs.

Confessions of an icosahedral virus crystallographer

This is a personal history of my structural studies of icosahedral viruses that evolved from crystallographic studies, to hybrid methods with electron cryo-microscopy and image reconstruction (cryoEM) and then developed further by incorporating a variety of physical methods to augment the high resolution crystallographic studies. It is not meant to be comprehensive, even for my own work, but hopefully provides some perspective on the growth of our understanding of these remarkable biologic assemblies. The goal is to provide a historical perspective for those new to the field and to emphasize the limitations of any one method, even those that provide atomic resolution information about viruses.

Genetic recombination in plant-infecting messenger-sense RNA viruses: overview and research perspectives

RNA recombination is one of the driving forces of genetic variability in (+)-strand RNA viruses. Various types of RNA-RNA crossovers were described including crosses between the same or different viral RNAs or between viral and cellular RNAs. Likewise, a variety of molecular mechanisms are known to support RNA recombination, such as replicative events (based on internal or end-to-end replicase switchings) along with non-replicative joining among RNA fragments of viral and/or cellular origin. Such mechanisms as RNA decay or RNA interference are responsible for RNA fragmentation and trans-esterification reactions which are likely accountable for ligation of RNA fragments. Numerous host factors were found to affect the profiles of viral RNA recombinants and significant differences in recombination frequency were observed among various RNA viruses. Comparative analyses of viral sequences allowed for the development of evolutionary models in order to explain adaptive phenotypic changes and co-evolving sites. Many questions remain to be answered by forthcoming RNA recombination research. (1) How various factors modulate the ability of viral replicase to switch templates, (2) What is the intracellular location of RNA-RNA template switchings, (3) Mechanisms and factors responsible for non-replicative RNA recombination, (4) Mechanisms of integration of RNA viral sequences with cellular genomic DNA, and (5) What is the role of RNA splicing and ribozyme activity. From an evolutionary stand point, it is not known how RNA viruses parasitize new host species via recombination, nor is it obvious what the contribution of RNA recombination is among other RNA modification pathways. We do not understand why the frequency of RNA recombination varies so much among RNA viruses and the status of RNA recombination as a form of sex is not well documented.