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.

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

High-throughput genomics of SARS-CoV-2 is essential to characterize virus evolution and to identify adaptations that affect pathogenicity or transmission. While single-nucleotide variations (SNVs) are commonly considered as driving virus adaption, RNA recombination events that delete or insert nucleic acid sequences are also critical. Whole genome targeting sequencing of SARS-CoV-2 is typically achieved using pairs of primers to generate cDNA amplicons suitable for next-generation sequencing (NGS). However, paired-primer approaches impose constraints on where primers can be designed, how many amplicons are synthesized and requires multiple PCR reactions with non-overlapping primer pools. This imparts sensitivity to underlying SNVs and fails to resolve RNA recombination junctions that are not flanked by primer pairs. To address these limitations, we have designed an approach called 'Tiled-ClickSeq', which uses hundreds of tiled-primers spaced evenly along the virus genome in a single reverse-transcription reaction. The other end of the cDNA amplicon is generated by azido-nucleotides that stochastically terminate cDNA synthesis, removing the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended to the cDNA fragment using click-chemistry and a PCR reaction generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5'UTR, at high depth and specificity to the virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and clinical samples to simultaneously characterize minority variants, sub-genomic mRNAs (sgmRNAs), structural variants (SVs) and D-RNAs. Tiled-ClickSeq therefore provides a convenient and robust platform for SARS-CoV-2 genomics that captures the full range of RNA species in a single, simple assay.

Long Noncoding RNAs in Plant Viroids and Viruses: A Review.

Infectious long-noncoding (lnc) RNAs related to plants can be of both viral and non-viral origin. Viroids are infectious plant lncRNAs that are not related to viruses and carry the circular, single-stranded, non-coding RNAs that replicate with host enzymatic activities via a rolling circle mechanism. Viroids interact with host processes in complex ways, emerging as one of the most productive tools for studying the functions of lncRNAs. Defective (D) RNAs, another category of lnc RNAs, are found in a variety of plant RNA viruses, most of which are noncoding. These are derived from and are replicated by the helper virus. D RNA-virus interactions evolve into mutually beneficial combinations, enhancing virus fitness via competitive advantages of moderated symptoms. Yet the satellite RNAs are single-stranded and include either large linear protein-coding ss RNAs, small linear ss RNAs, or small circular ss RNAs (virusoids). The satellite RNAs lack sequence homology to the helper virus, but unlike viroids need a helper virus to replicate and encapsidate. They can attenuate symptoms via RNA silencing and enhancement of host defense, but some can be lethal as RNA silencing suppressor antagonists. Moreover, selected viruses produce lncRNAs by incomplete degradation of genomic RNAs. They do not replicate but may impact viral infection, gene regulation, and cellular functions. Finally, the host plant lncRNAs can also contribute during plant-virus interactions, inducing plant defense and the regulation of gene expression, often in conjunction with micro and/or circRNAs.

Defective RNA particles derived from Tomato black ring virus genome interfere with the replication of parental virus.

Tomato black ring virus (TBRV) is the only member of the Nepovirus genus that is known to form defective RNA particles (D RNAs) during replication. Here, de novo generation of D RNAs was observed during prolonged passages of TBRV isolates originated from Solanum lycopersicum and Lactuca sativa in Chenopodium quinoa plants. D RNAs of about 500 nt derived by a single deletion in the RNA1 molecule and contained a portion of the 5' untranslated region and viral replicase, and almost the entire 3' non-coding region. Short regions of sequence complementarity were found at the 5' and 3' junction borders, which can facilitate formation of the D RNAs. Moreover, in this study we analyzed the effects of D RNAs on TBRV replication and symptoms development of infected plants. C. quinoa, S. lycopersicum, Nicotiana tabacum, and L. sativa were infected with the original TBRV isolates (TBRV-D RNA) and those containing additional D RNA particles (TBRV + D RNA). The viral accumulation in particular hosts was measured up to 28 days post inoculation by RT-qPCR. Statistical analyses revealed that D RNAs interfere with TBRV replication and thus should be referred to as defective interfering particles. The magnitude of the interference effect depends on the interplay between TBRV isolate and host species.