Whereas DNA viruses are known to be abundant, diverse, and commonly key ecosystem players, RNA viruses are insufficiently studied outside disease settings. In this study, we analyzed ≈28 terabases of Global Ocean RNA sequences to expand Earth’s RNA virus catalogs and their taxonomy, investigate their evolutionary origins, and assess their marine biogeography from pole to pole. Using new approaches to optimize discovery and classification, we identified RNA viruses that necessitate substantive revisions of taxonomy (doubling phyla and adding >50% new classes) and evolutionary understanding. “Species”-rank abundance determination revealed that viruses of the new phyla “ Taraviricota ,” a missing link in early RNA virus evolution, and “ Arctiviricota ” are widespread and dominant in the oceans. These efforts provide foundational knowledge critical to integrating RNA viruses into ecological and epidemiological models.

RNA silencing is a potent means of antiviral defense in plants and animals. A hallmark of this defense response is the production of 21- to 24-nucleotide viral small RNAs via mechanisms that remain to be fully understood. Many viruses encode suppressors of RNA silencing, and some viral RNAs function directly as silencing suppressors as counterdefense. The occurrence of viroid-specific small RNAs in infected plants suggests that viroids can trigger RNA silencing in a host, raising the question of how these noncoding and unencapsidated RNAs survive cellular RNA-silencing systems. We address this question by characterizing the production of small RNAs of Potato spindle tuber viroid (srPSTVds) and investigating how PSTVd responds to RNA silencing. Our molecular and biochemical studies provide evidence that srPSTVds were derived mostly from the secondary structure of viroid RNAs. Replication of PSTVd was resistant to RNA silencing, although the srPSTVds were biologically active in guiding RNA-induced silencing complex (RISC)-mediated cleavage, as shown with a sensor system. Further analyses showed that without possessing or triggering silencing suppressor activities, the PSTVd secondary structure played a critical role in resistance to RISC-mediated cleavage. These findings support the hypothesis that some infectious RNAs may have evolved specific secondary structures as an effective means to evade RNA silencing in addition to encoding silencing suppressor activities. Our results should have important implications in further studies on RNA-based mechanisms of host-pathogen interactions and the biological constraints that shape the evolution of infectious RNA structures.

The formation of secondary structures by a random RNA sequence is studied as a model system for the sequence-structure problem omnipresent in biopolymers. Several toy energy models are introduced to allow detailed analytical and numerical studies. First, a two-replica calculation is performed. By mapping the two-replica problem to the denaturation of a single homogeneous RNA molecule in six-dimensional embedding space, we show that sequence disorder is perturbatively irrelevant, i.e., an RNA molecule with weak sequence disorder is in a molten phase where many secondary structures with comparable total energy coexist. A numerical study of various models at high temperature reproduces behaviors characteristic of the molten phase. On the other hand, a scaling argument based on the external statistics of rare regions can be constructed to show that the low-temperature phase is unstable to sequence disorder. We performed a detailed numerical study of the low-temperature phase using the droplet theory as a guide, and characterized the statistics of large-scale, low-energy excitations of the secondary structures from the ground state structure. We find the excitation energy to grow very slowly (i.e., logarithmically) with the length scale of the excitation, suggesting the existence of a marginal glass phase. The transition between the low-temperature glass phase and the high-temperature molten phase is also characterized numerically. It is revealed by a change in the coefficient of the logarithmic excitation energy, from being disorder dominated to being entropy dominated.

CRISPR/Cas9 genome editing has been proposed as a therapeutic treatment for HIV-1 infection. CRISPR/Cas9 induced double strand breaks (DSBs) targeted to the integrated viral genome have been shown to decrease production of progeny virus. Unfortunately HIV-1 evolves rapidly and may readily produce CRISPR/Cas9 resistant strains. Here we used next-generation sequencing to characterize HIV-1 strains that developed resistance to six different CRISPR/Cas9 guide RNAs (gRNAs). Reverse transcriptase (RT) derived base substitution mutations were commonly found at sites encoding unpaired bases of RNA stem-loop structures. In addition to RT mutations, insertion and/or deletion (indel) mutations were common. Indels localized to the CRISPR/Cas9 cleavage site were major contributors to CRISPR gRNA resistance. While most indels at non-coding regions were a single base pair, 3 base pair indels were observed when a coding region of HIV-1 was targeted. The DSB repair event may preserve the HIV-1 reading frame, while destroying CRISPR gRNA homology. HIV-1 may be successfully edited by CRISPR/Cas9, but the virus remains competent for replication and resistant to further CRISPR/Cas9 targeting at that site. These observations strongly suggest that host DSB repair at CRISPR/Cas9 cleavage sites is a novel and important pathway that may contribute to HIV-1 therapeutic resistance.