Analyzing viral epitranscriptomes using nanopore direct RNA sequencing.

RNA modifications are a common occurrence across all domains of life. Several chemical modifications, including N6-methyladenosine, have also been found in viral transcripts and viral RNA genomes. Some of the modifications increase the viral replication efficiency while also helping the virus to evade the host immune system. Nonetheless, there are numerous examples in which the host's RNA modification enzymes function as antiviral factors. Although established methods like MeRIP-seq and miCLIP can provide a transcriptome- wide overview of how viral RNA is modified, it is difficult to distinguish between the complex overlapping viral transcript isoforms using the short read-based techniques. Nanopore direct RNA sequencing (DRS) provides both long reads and direct signal readings, which may carry information about the modifications. Here, we describe a refined protocol for analyzing the RNA modifications in viral transcriptomes using nanopore technology.

Functional and molecular dissection of HCMV long non-coding RNAs.

Small, compact genomes confer a selective advantage to viruses, yet human cytomegalovirus (HCMV) expresses the long non-coding RNAs (lncRNAs); RNA1.2, RNA2.7, RNA4.9, and RNA5.0. Little is known about the function of these lncRNAs in the virus life cycle. Here, we dissected the functional and molecular landscape of HCMV lncRNAs. We found that HCMV lncRNAs occupy ~ 30% and 50-60% of total and poly(A)+viral transcriptome, respectively, throughout virus life cycle. RNA1.2, RNA2.7, and RNA4.9, the three abundantly expressed lncRNAs, appear to be essential in all infection states. Among these three lncRNAs, depletion of RNA2.7 and RNA4.9 results in the greatest defect in maintaining latent reservoir and promoting lytic replication, respectively. Moreover, we delineated the global post-transcriptional nature of HCMV lncRNAs by nanopore direct RNA sequencing and interactome analysis. We revealed that the lncRNAs are modified with N6-methyladenosine (m6A) and interact with m6A readers in all infection states. In-depth analysis demonstrated that m6A machineries stabilize HCMV lncRNAs, which could account for the overwhelming abundance of viral lncRNAs. Our study lays the groundwork for understanding the viral lncRNA-mediated regulation of host-virus interaction throughout the HCMV life cycle.

Partitioning RNAs by length improves transcriptome reconstruction from short-read RNA-seq data.

The accuracy of methods for assembling transcripts from short-read RNA sequencing data is limited by the lack of long-range information. Here we introduce Ladder-seq, an approach that separates transcripts according to their lengths before sequencing and uses the additional information to improve the quantification and assembly of transcripts. Using simulated data, we show that a kallisto algorithm extended to process Ladder-seq data quantifies transcripts of complex genes with substantially higher accuracy than conventional kallisto. For reference-based assembly, a tailored scheme based on the StringTie2 algorithm reconstructs a single transcript with 30.8% higher precision than its conventional counterpart and is more than 30% more sensitive for complex genes. For de novo assembly, a similar scheme based on the Trinity algorithm correctly assembles 78% more transcripts than conventional Trinity while improving precision by 78%. In experimental data, Ladder-seq reveals 40% more genes harboring isoform switches compared to conventional RNA sequencing and unveils widespread changes in isoform usage upon m6A depletion by Mettl14 knockout.