SARS-CoV-2 is a betacoronavirus responsible for the COVID-19 pandemic. Although the SARS-CoV-2 genome was reported recently, its transcriptomic architecture is unknown. Utilizing two complementary sequencing techniques, we present a high-resolution map of the SARS-CoV-2 transcriptome and epitranscriptome. DNA nanoball sequencing shows that the transcriptome is highly complex owing to numerous discontinuous transcription events. In addition to the canonical genomic and 9 subgenomic RNAs, SARS-CoV-2 produces transcripts encoding unknown ORFs with fusion, deletion, and/or frameshift. Using nanopore direct RNA sequencing, we further find at least 41 RNA modification sites on viral transcripts, with the most frequent motif, AAGAA. Modified RNAs have shorter poly(A) tails than unmodified RNAs, suggesting a link between the modification and the 3 0 tail. Functional investigation of the unknown transcripts and RNA modifications discovered in this study will open new directions to our understanding of the life cycle and pathogenicity of SARS-CoV-2. ll

Global investigation of the 3' extremity of mRNA (3'-terminome), despite its importance in gene regulation, has not been feasible due to technical challenges associated with homopolymeric sequences and relative paucity of mRNA. We here develop a method, TAIL-seq, to sequence the very end of mRNA molecules. TAIL-seq allows us to measure poly(A) tail length at the genomic scale. Median poly(A) length is 50-100 nt in HeLa and NIH 3T3 cells. Poly(A) length correlates with mRNA half-life, but not with translational efficiency. Surprisingly, we discover widespread uridylation and guanylation at the downstream of poly(A) tail. The U tails are generally attached to short poly(A) tails (<25 nt), while the G tails are found mainly on longer poly(A) tails (>40 nt), implicating their generic roles in mRNA stability control. TAIL-seq is a potent tool to dissect dynamic control of mRNA turnover and translational control, and to discover unforeseen features of RNA cleavage and tailing.

A hallmark of RNA silencing is a class of approximately 22-nucleotide RNAs that are processed from double-stranded RNA precursors by Dicer. Accurate processing by Dicer is crucial for the functionality of microRNAs (miRNAs). The current model posits that Dicer selects cleavage sites by measuring a set distance from the 3′ overhang of the double-stranded RNA terminus. Here we report that human Dicer anchors not only the 3′ end but also the 5′ end, with the cleavage site determined mainly by the distance (~22 nucleotides) from the 5′ end (5′ counting rule). This cleavage requires a 5′-terminal phosphate group. Further, we identify a novel basic motif (5′ pocket) in human Dicer that recognizes the 5′-phosphorylated end. The 5′ counting rule and the 5′ anchoring residues are conserved in Drosophila Dicer-1, but not in Giardia Dicer. Mutations in the 5′ pocket reduce processing efficiency and alter cleavage sites in vitro. Consistently, miRNA biogenesis is perturbed in vivo when Dicer-null embryonic stem cells are replenished with the 5′-pocket mutant. Thus, 5′-end recognition by Dicer is important for precise and effective biogenesis of miRNAs. Insights from this study should also afford practical benefits to the design of small hairpin RNAs.

RNase III Drosha initiates microRNA (miRNA) maturation by cleaving a primary miRNA transcript and releasing a pre-miRNA with a 2 nt 3' overhang. Dicer recognizes the 2 nt 3' overhang structure to selectively process pre-miRNAs. Here, we find that, unlike prototypic pre-miRNAs (group I), group II pre-miRNAs acquire a shorter (1 nt) 3' overhang from Drosha processing and therefore require a 3'-end mono-uridylation for Dicer processing. The majority of let-7 and miR-105 belong to group II. We identify TUT7/ZCCHC6, TUT4/ZCCHC11, and TUT2/PAPD4/GLD2 as the terminal uridylyl transferases responsible for pre-miRNA mono-uridylation. The TUTs act specifically on dsRNAs with a 1 nt 3' overhang, thereby creating a 2 nt 3' overhang. Depletion of TUTs reduces let-7 levels and disrupts let-7 function. Although the let-7 suppressor, Lin28, induces inhibitory oligo-uridylation in embryonic stem cells, mono-uridylation occurs in somatic cells lacking Lin28 to promote let-7 biogenesis. Our study reveals functional duality of uridylation and introduces TUT7/4/2 as components of the miRNA biogenesis pathway.