Emergence of within-host SARS-CoV-2 recombinant genome after coinfection by Gamma and Delta variants

Since the first reports of patients coinfected by two genetically-distinct lineages of SARS-CoV-2, the scientific community raised concerns about the recombination of intra-host viral RNA sequences as a possible mechanism underlying the emergence of novel variants. Indeed, this phenomenon occurs at a relatively high frequency among betacoronaviruses. Nevertheless, the few existing studies about recombination between genetically-distinct lineages of SARS-CoV-2 are restricted to detect the inter-host dissemination of genomes post-recombination events. However, the high genomic similarity between the current co-circulating lineages challenges the identification of these events. Here, we report the first case of intra-host SARS-CoV-2 recombination during a coinfection by the variants of concern (VOC) AY.33 (Delta) and P.1 (Gamma) supported by sequencing reads harboring a mosaic of lineage-defining mutations. By using next-generation sequencing reads intersecting regions that simultaneously overlap lineage-defining mutations from Gamma and Delta, we were able to identify a total of six recombinant regions across the SARS-CoV-2 genome within a sample. Four of them mapped in the spike gene and two in the nucleocapsid gene. We detected mosaic reads harboring a combination of lineage-defining mutations from each VOC. To our knowledge, this is the first report of intra-host RNA-RNA recombination between two lineages of SARS-CoV-2, which can represent a threat to public health management during the COVID-19 pandemic due to the possibility of the emergence of viruses with recombinant phenotypes.

Archaeal tRNA-Splicing Endonuclease as an Effector for RNA Recombination and Novel Trans-Splicing Pathways in Eukaryotes

We have characterized a homodimeric tRNA endonuclease from the euryarchaeota Ferroplasma acidarmanus (FERAC), a facultative anaerobe which can grow at temperatures ranging from 35 to 42 °C. This enzyme, contrary to the eukaryal tRNA endonucleases and the homotetrameric Methanocaldococcus jannaschii (METJA) homologs, is able to cleave minimal BHB (bulge–helix–bulge) substrates at 30 °C. The expression of this enzyme in Schizosaccharomyces pombe (SCHPO) enables the use of its properties as effectors by inserting BHB motif introns into hairpin loops normally seen in mRNA transcripts. In addition, the FERAC endonuclease can create proteins with new functionalities through the recombination of protein domains.

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.

Characterization of a Cytopathogenic Reporter CSFV

Cytopathogenic (cp) pestiviruses frequently emerge in cattle that are persistently infected with the bovine viral diarrhea virus (BVDV) as a consequence of RNA recombination and mutation. They induce apoptosis in infected tissue cultures, are highly attenuated in the immunocompetent host, and unable to establish persistent infections after diaplacental infections. Cp strains of BVDV have been used as naturally attenuated live vaccines and for species-specific plaque reduction tests for the indirect serological detection of BVDV. Here, we present a genetically engineered cp strain of the classical swine fever virus (CSFV). Cytopathogenicity of the strain was induced by the insertion of ubiquitin embedded in a large NS3 to NS4B duplication. The CSFV RNA genome was stabilized by the inactivation of the NS2 autoprotease, hindering the deletion of the insertion and the reversion to a wild-type genome. Additional insertion of a mCherry gene at the 5′-end of the E2 gene allowed fluorescence-verified plaque reduction assays for CSFV, thus providing a novel, cost-efficient diagnostic tool. This genetically stabilized cp CSFV strain could be further used as a basis for potential new modified live vaccines. Taken together, we applied reverse genetics to rationally fixate a typical cp NS3 duplication in a CSFV genome.

Participation of Multifunctional RNA in Replication, Recombination and Regulation of Endogenous Plant Pararetroviruses (EPRVs)

Pararetroviruses, taxon Caulimoviridae, are typical of retroelements with reverse transcriptase and share a common origin with retroviruses and LTR retrotransposons, presumably dating back 1.6 billion years and illustrating the transition from an RNA to a DNA world. After transcription of the viral genome in the host nucleus, viral DNA synthesis occurs in the cytoplasm on the generated terminally redundant RNA including inter- and intra-molecule recombination steps rather than relying on nuclear DNA replication. RNA recombination events between an ancestral genomic retroelement with exogenous RNA viruses were seminal in pararetrovirus evolution resulting in horizontal transmission and episomal replication. Instead of active integration, pararetroviruses use the host DNA repair machinery to prevail in genomes of angiosperms, gymnosperms and ferns. Pararetrovirus integration – leading to Endogenous ParaRetroViruses, EPRVs – by illegitimate recombination can happen if their sequences instead of homologous host genomic sequences on the sister chromatid (during mitosis) or homologous chromosome (during meiosis) are used as template. Multiple layers of RNA interference exist regulating episomal and chromosomal forms of the pararetrovirus. Pararetroviruses have evolved suppressors against this plant defense in the arms race during co-evolution which can result in deregulation of plant genes. Small RNAs serve as signaling molecules for Transcriptional and Post-Transcriptional Gene Silencing (TGS, PTGS) pathways. Different populations of small RNAs comprising 21–24 nt and 18–30 nt in length have been reported for Citrus, Fritillaria, Musa, Petunia, Solanum and Beta. Recombination and RNA interference are driving forces for evolution and regulation of EPRVs.

Structure and dynamics of SARS-CoV-2 proofreading exoribonuclease ExoN

High-fidelity replication of the large RNA genome of coronaviruses (CoVs) is mediated by a 3'-to-5' exoribonuclease (ExoN) in non-structural protein 14 (nsp14), which excises nucleotides including antiviral drugs mis-incorporated by the low-fidelity viral RNA-dependent RNA polymerase (RdRp) and has also been implicated in viral RNA recombination and resistance to innate immunity. Here we determined a 1.6-Å resolution crystal structure of SARS-CoV-2 ExoN in complex with its essential co-factor, nsp10. The structure shows a highly basic and concave surface flanking the active site, comprising several Lys residues of nsp14 and the N-terminal amino group of nsp10. Modeling suggests that this basic patch binds to the template strand of double-stranded RNA substrates to position the 3' end of the nascent strand in the ExoN active site, which is corroborated by mutational and computational analyses. Molecular dynamics simulations further show remarkable flexibility of multi-domain nsp14 and suggest that nsp10 stabilizes ExoN for substrate RNA-binding to support its exoribonuclease activity. Our high-resolution structure of the SARS-CoV-2 ExoN-nsp10 complex serves as a platform for future development of anti-coronaviral drugs or strategies to attenuate the viral virulence.

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’. Tiled-ClickSeq 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, obviating the need for a paired-primer. A sequencing adaptor containing a Unique Molecular Identifier (UMI) is appended using click-chemistry and a PCR reaction using Illumina adaptors generates a final NGS library. Tiled-ClickSeq provides complete genome coverage, including the 5’UTR, at high depth and specificity to virus on both Illumina and Nanopore NGS platforms. Here, we analyze multiple SARS-CoV-2 isolates and 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.

Reovirus RNA recombination is sequence directed and generates internally deleted defective genome segments during passage

For viruses with segmented genomes, genetic diversity is generated by genetic drift, reassortment, and recombination. Recombination produces RNA populations distinct from full-length gene segments and can influence viral population dynamics, persistence, and host immune responses. Viruses in the Reoviridae family, including rotavirus and mammalian orthoreovirus (reovirus), have been reported to package segments containing rearrangements or internal deletions. Rotaviruses with RNA segments containing rearrangements have been isolated from immunocompromised and immunocompetent children and in vitro following serial passage at relatively high multiplicity. Reoviruses that package small, defective RNA segments have established chronic infections in cells and in mice. However, the mechanism and extent of Reoviridae RNA recombination are undefined. Towards filling this gap in knowledge, we determined the titers and RNA segment profiles for reovirus and rotavirus following serial passage in cultured cells. The viruses exhibited occasional titer reductions characteristic of interference. Reovirus strains frequently accumulated segments that retained 5′ and 3′ terminal sequences and featured large internal deletions, while similarly fragmented segments were rarely detected in rotavirus populations. Using next-generation RNA-sequencing to analyze RNA molecules packaged in purified reovirus particles, we identified distinct recombination sites within individual viral genome segments. Recombination junctions were frequently but not always characterized by short direct sequence repeats upstream and downstream that spanned junction sites. Taken together, these findings suggest that reovirus accumulates defective gene segments featuring internal deletions during passage and undergoes sequence-directed recombination at distinct sites. IMPORTANCE Viruses in the Reoviridae family include important pathogens of humans and other animals and have segmented RNA genomes. Recombination in RNA virus populations can facilitate novel host exploration and increased disease severity. The extent, patterns, and mechanisms of Reoviridae recombination and the functions and effects of recombined RNA products are poorly understood. Here, we provide evidence that mammalian orthoreovirus regularly synthesizes RNA recombination products that retain terminal sequences but contain internal deletions, while rotavirus rarely synthesizes such products. Recombination occurs more frequently at specific sites in the mammalian orthoreovirus genome, and short regions of identical sequence are often detected at junction sites. These findings suggest that mammalian orthoreovirus recombination events are directed in part by RNA sequences. An improved understanding of recombined viral RNA synthesis may enhance our capacity to engineer improved vaccines and virotherapies in the future.