Poly(A) inclusive RNA isoform sequencing (PAIso-seq) reveals wide-spread non-adenosine residues within RNA poly(A) tails

Beyond simple tails: poly(A) tail-mediated RNA epigenetic regulation

The poly(A) tail is an essential structural component of mRNA required for the latter's stability and translation. Recent technologies have enabled transcriptome-wide profiling of the length and composition of poly(A) tails, shedding light on their overlooked regulatory capacities. Notably, poly(A) tails contain not only adenine but also uracil, cytosine, and guanine residues. These findings strongly suggest that poly(A) tails could encode a wealth of regulatory information, similar to known reversible RNA chemical modifications. This review aims to succinctly summarize our current knowledge on the composition, dynamics, and regulatory functions of RNA poly(A) tails. Given their capacity to carry rich regulatory information beyond the genetic code, we propose the concept of 'poly(A) tail epigenetic information' as a new layer of RNA epigenetic regulation.

Maternal mRNA deadenylation is defective in in vitro matured mouse and human oocytes

Oocyte in vitro maturation is a technique in assisted reproductive technology. Thousands of genes show abnormally high expression in in vitro maturated metaphase II (MII) oocytes compared to those matured in vivo in bovines, mice, and humans. The mechanisms underlying this phenomenon are poorly understood. Here, we use poly(A) inclusive RNA isoform sequencing (PAIso-seq) for profiling the transcriptome-wide poly(A) tails in both in vivo and in vitro matured mouse and human oocytes. Our results demonstrate that the observed increase in maternal mRNA abundance is caused by impaired deadenylation in in vitro MII oocytes. Moreover, the cytoplasmic polyadenylation of dormant Btg4 and Cnot7 mRNAs, which encode key components of deadenylation machinery, is impaired in in vitro MII oocytes, contributing to reduced translation of these deadenylase machinery components and subsequently impaired global maternal mRNA deadenylation. Our findings highlight impaired maternal mRNA deadenylation as a distinct molecular defect in in vitro MII oocytes.

FUS reads histone H3K36me3 to regulate alternative polyadenylation

Complex organisms generate differential gene expression through the same set of DNA sequences in distinct cells. The communication between chromatin and RNA regulates cellular behavior in tissues. However, little is known about how chromatin, especially histone modifications, regulates RNA polyadenylation. In this study, we found that FUS was recruited to chromatin by H3K36me3 at gene bodies. The H3K36me3 recognition of FUS was mediated by the proline residues in the ZNF domain. After these proline residues were mutated or H3K36me3 was abolished, FUS dissociated from chromatin and bound more to RNA, resulting in an increase in polyadenylation sites far from stop codons genome-wide. A proline mutation corresponding to a mutation in amyotrophic lateral sclerosis contributed to the hyperactivation of mitochondria and hyperdifferentiation in mouse embryonic stem cells. These findings reveal that FUS is an H3K36me3 reader protein that links chromatin-mediated alternative polyadenylation to human disease.

Deadenylation kinetics of mixed poly(A) tails at single-nucleotide resolution

Shortening of messenger RNA poly(A) tails, or deadenylation, is a rate-limiting step in mRNA decay and is highly regulated during gene expression. The incorporation of non-adenosines in poly(A) tails, or 'mixed tailing', has been observed in vertebrates and viruses. Here, to quantitate the effect of mixed tails, we mathematically modeled deadenylation reactions at single-nucleotide resolution using an in vitro deadenylation system reconstituted with the complete human CCR4-NOT complex. Applying this model, we assessed the disrupting impact of single guanosine, uridine or cytosine to be equivalent to approximately 6, 8 or 11 adenosines, respectively. CCR4-NOT stalls at the 0, -1 and -2 positions relative to the non-adenosine residue. CAF1 and CCR4 enzyme subunits commonly prefer adenosine but exhibit distinct sequence selectivities and stalling positions. Our study provides an analytical framework to monitor deadenylation and reveals the molecular basis of tail sequence-dependent regulation of mRNA stability.

Control of poly(A)-tail length and translation in vertebrate oocytes and early embryos

During oocyte maturation and early embryogenesis, changes in mRNA poly(A)-tail lengths strongly influence translation, but how these tail-length changes are orchestrated has been unclear. Here, we performed tail-length and translational profiling of mRNA reporter libraries (each with millions of 3' UTR sequence variants) in frog oocytes and embryos and in fish embryos. Contrasting to previously proposed cytoplasmic polyadenylation elements (CPEs), we found that a shorter element, UUUUA, together with the polyadenylation signal (PAS), specify cytoplasmic polyadenylation, and we identified contextual features that modulate the activity of both elements. In maturing oocytes, this tail lengthening occurs against a backdrop of global deadenylation and the action of C-rich elements that specify tail-length-independent translational repression. In embryos, cytoplasmic polyadenylation becomes more permissive, and additional elements specify waves of stage-specific deadenylation. Together, these findings largely explain the complex tapestry of tail-length changes observed in early frog and fish development, with strong evidence of conservation in both mice and humans.

Degradation and translation of maternal mRNA for embryogenesis

Maternal mRNAs accumulate during egg growth and must be judiciously degraded or translated to ensure successful development of mammalian embryos. In this review we integrate recent investigations into pathways controlling rapid degradation of maternal mRNAs during the maternal-to-zygotic transition. Degradation is not indiscriminate, and some mRNAs are selectively protected and rapidly translated after fertilization for reprogramming the zygotic genome during early embryogenesis. Oocyte specific cofactors and pathways have been illustrated to control different futures of maternal mRNAs. We discuss mechanisms that control the fate of maternal mRNAs during late oogenesis and after fertilization. Issues to be resolved in current maternal mRNA research are described, and future research directions are proposed.

Self-feedback loop-containing synthetic mRNA switches for controlled microRNA sensing

Synthetic mRNA switches are powerful cell fate manipulation tools that sense cellular input molecules to directly control protein expression at the translational level. The lack of available switch designs that can mimic the natural sophisticated protein regulation is a fundamental issue that limits the application of synthetic mRNA switches. Here we report a new set of synthetic mRNA switches by incorporating self-feedback loop machineries to dynamically control protein expression levels upon sensing cellular microRNAs. We redesigned the coding region of the switch to express output protein along with mRNA regulatory proteins. RNA-binding proteins (RBPs) and RBP-binding RNA motifs (aptamers) guide the regulatory proteins to act on their own mRNAs, enhancing or flattening the effect of microRNA sensing. Importantly, we demonstrated that the switches with the positive feedback feature can enlarge a high-or-low microRNA effect into a nearly all-or-none pattern, substantially boosting the use of synthetic mRNA switches as high-performance microRNA sensors or binary cell regulation tools. We believe these novel mRNA switch designs provide new strategies to construct complex mRNA-based genetic circuits for future molecular sensing and cell engineering.

Poly(A) tale: From A to A; RNA polyadenylation in prokaryotes and eukaryotes

Most eukaryotic mRNAs and different non-coding RNAs undergo a form of 3' end processing known as polyadenylation. Polyadenylation machinery is present in almost all organisms except few species. In bacteria, the machinery has evolved from PNPase, which adds heteropolymeric tails, to a poly(A)-specific polymerase. Differently, a complex machinery for accurate polyadenylation and several non-canonical poly(A) polymerases are developed in eukaryotes. The role of poly(A) tail has also evolved from serving as a degradative signal to a stabilizing modification that also regulates translation. In this review, we discuss poly(A) tail emergence in prokaryotes and its development into a stable, yet dynamic feature at the 3' end of mRNAs in eukaryotes. We also describe how appearance of novel poly(A) polymerases gives cells flexibility to shape poly(A) tail. We explain how poly(A) tail dynamics help regulate cognate RNA metabolism in a context-dependent manner, such as during oocyte maturation. Finally, we describe specific mRNAs in metazoans that bear stem-loops instead of poly(A) tails. We conclude with how recent discoveries about poly(A) tail can be applied to mRNA technology. This article is categorized under: RNA Evolution and Genomics > RNA and Ribonucleoprotein Evolution RNA Processing > 3' End Processing RNA Turnover and Surveillance > Regulation of RNA Stability.

An extended wave of global mRNA deadenylation sets up a switch in translation regulation across the mammalian oocyte-to-embryo transition

Without new transcription, gene expression across the oocyte-to-embryo transition (OET) relies instead on regulation of mRNA poly(A) tails to control translation. However, how tail dynamics shape translation across the OET in mammals remains unclear. We perform long-read RNA sequencing to uncover poly(A) tail lengths across the mouse OET and, incorporating published ribosome profiling data, provide an integrated, transcriptome-wide analysis of poly(A) tails and translation across the entire transition. We uncover an extended wave of global deadenylation during fertilization in which short-tailed, oocyte-deposited mRNAs are translationally activated without polyadenylation through resistance to deadenylation. Subsequently, in the embryo, mRNAs are readenylated and translated in a surge of global polyadenylation. We further identify regulation of poly(A) tail length at the isoform level and stage-specific enrichment of mRNA sequence motifs among regulated transcripts. These data provide insight into the stage-specific mechanisms of poly(A) tail regulation that orchestrate gene expression from oocyte to embryo in mammals.

A genome-wide perspective of the maternal mRNA translation program during oocyte development

Transcriptional and post-transcriptional regulations control gene expression in most cells. However, critical transitions during the development of the female gamete relies exclusively on regulation of mRNA translation in the absence of de novo mRNA synthesis. Specific temporal patterns of maternal mRNA translation are essential for the oocyte progression through meiosis, for generation of a haploid gamete ready for fertilization and for embryo development. In this review, we will discuss how mRNAs are translated during oocyte growth and maturation using mostly a genome-wide perspective. This broad view on how translation is regulated reveals multiple divergent translational control mechanisms required to coordinate protein synthesis with progression through the meiotic cell cycle and with development of a totipotent zygote.

Alternative cleavage and polyadenylation of the Ccnb1 mRNA defines accumulation of cyclin protein during the meiotic cell cycle

Progression through the mitotic and meiotic cell cycle is driven by fluctuations in the levels of cyclins, the regulatory subunits controlling the localization and activity of CDK1 kinases. Cyclin levels are regulated through a precise balance of synthesis and degradation. Here we demonstrate that the synthesis of Cyclin B1 during the oocyte meiotic cell cycle is defined by the selective translation of mRNA variants generated through alternative cleavage and polyadenylation (APA). Using gene editing in mice, we introduced mutations into the proximal and distal polyadenylation elements of the 3' untranslated region (UTR) of the Ccnb1 mRNA. Through in vivo loss-of-function experiments, we demonstrate that the translation of mRNA with a short 3' UTR specifies Cyclin B1 protein levels that set the timing of meiotic re-entry. In contrast, translation directed by a long 3' UTR is necessary to direct Cyclin B1 protein accumulation during the MI/MII transition. These findings establish that the progression through the cell cycle is dependent on the selective translation of multiple mRNA variants generated by APA.

Decoding the Heterogeneity and Specialized Function of Translation Machinery Through Ribosome Profiling in Yeast Mutants of Initiation Factors

The nuanced heterogeneity and specialized functions of translation machinery are increasingly recognized as crucial for precise translational regulation. Here, high-throughput ribosomal profiling (ribo-seq) is used to analyze the specialized roles of eukaryotic initiation factors (eIFs) in the budding yeast. By examining changes in ribosomal distribution across the genome resulting from knockouts of eIF4A, eIF4B, eIF4G1, CAF20, or EAP1, or knockdowns of eIF1, eIF1A, eIF4E, or PAB1, two distinct initiation-factor groups, the "looping" and "scanning" groups are discerned, based on similarities in the ribosomal landscapes their perturbation induced. The study delves into the cis-regulatory sequence features of genes influenced predominantly by each group, revealing that genes more dependent on the looping-group factors generally have shorter transcripts and poly(A) tails. In contrast, genes more dependent on the scanning-group factors often possess upstream open reading frames and exhibit a higher GC content in their 5' untranslated regions. From the ribosomal RNA fragments identified in the ribo-seq data, ribosomal heterogeneity associated with perturbation of specific initiation factors is further identified, suggesting their potential roles in regulating ribosomal components. Collectively, the study illuminates the complexity of translational regulation driven by heterogeneity and specialized functions of translation machinery, presenting potential approaches for targeted gene translation manipulation.

Sequencing of Transcriptome-Wide Poly(A) Tails by PAIso-seq

Poly(A) tails are added to most eukaryotic mRNA and have essential regulatory functions. However, due to its homopolymeric nature, the sequence information in poly(A) tails is challenging to obtain in transcriptome measurement studies. In this chapter, we describe the detailed procedures of poly(A) inclusive full-length RNA isoform-sequencing (PAIso-seq), a method that can measure transcriptome-wide poly(A) tails from as low as nanogram level of total RNA based on the PacBio HiFi sequencing platform. The accurate length and base composition of poly(A) tails can be obtained along with the full-length cDNA.

Improving the RNA velocity approach with single-cell RNA lifecycle (nascent, mature and degrading RNAs) sequencing technologies

We presented an experimental method called FLOUR-seq, which combines BD Rhapsody and nanopore sequencing to detect the RNA lifecycle (including nascent, mature, and degrading RNAs) in cells. Additionally, we updated our HIT-scISOseq V2 to discover a more accurate RNA lifecycle using 10x Chromium and Pacbio sequencing. Most importantly, to explore how single-cell full-length RNA sequencing technologies could help improve the RNA velocity approach, we introduced a new algorithm called 'Region Velocity' to more accurately configure cellular RNA velocity. We applied this algorithm to study spermiogenesis and compared the performance of FLOUR-seq with Pacbio-based HIT-scISOseq V2. Our findings demonstrated that 'Region Velocity' is more suitable for analyzing single-cell full-length RNA data than traditional RNA velocity approaches. These novel methods could be useful for researchers looking to discover full-length RNAs in single cells and comprehensively monitor RNA lifecycle in cells.

RNA modification in mRNA cancer vaccines

RNA modification is manifested as chemically altered nucleotides, widely exists in diverse natural RNAs, and is closely related to RNA structure and function. Currently, mRNA-based vaccines have received great attention and rapid development as novel and mighty fighters against various diseases including cancer. The achievement of RNA vaccines in clinical application is largely attributed to some methodological innovations including the incorporation of modified nucleotides into the synthetic RNA. The selection of optimal RNA modifications aimed at reducing the instability and immunogenicity of RNA molecules is a very critical task to improve the efficacy and safety of mRNA vaccines. This review summarizes the functions of RNA modifications and their application in mRNA vaccines, highlights recent advances of mRNA vaccines in cancer immunotherapy, and provides perspectives for future development of mRNA vaccines in the context of personalized tumor therapy.

Stress responses of plants through transcriptome plasticity by mRNA alternative polyadenylation

The sessile nature of plants confines their responsiveness to changing environmental conditions. Gene expression regulation becomes a paramount mechanism for plants to adjust their physiological and morphological behaviors. Alternative polyadenylation (APA) is known for its capacity to augment transcriptome diversity and plasticity, thereby furnishing an additional set of tools for modulating gene expression. APA has also been demonstrated to exhibit intimate associations with plant stress responses. In this study, we review APA dynamic features and consequences in plants subjected to both biotic and abiotic stresses. These stresses include adverse environmental stresses, and pathogenic attacks, such as cadmium toxicity, high salt, hypoxia, oxidative stress, cold, heat shock, along with bacterial, fungal, and viral infections. We analyzed the overarching research framework employed to elucidate plant APA response and the alignment of polyadenylation site transitions with the modulation of gene expression levels within the ambit of each stress condition. We also proposed a general APA model where transacting factors, including poly(A) factors, epigenetic regulators, RNA m6A modification factors, and phase separation proteins, assume pivotal roles in APA related transcriptome plasticity during stress response in plants.

CLT-seq as a universal homopolymer-sequencing concept reveals poly(A)-tail-tuned ncRNA regulation

Dynamic tuning of the poly(A) tail is a crucial mechanism for controlling translation and stability of eukaryotic mRNA. Achieving a comprehensive understanding of how this regulation occurs requires unbiased abundance quantification of poly(A)-tail transcripts and simple poly(A)-length measurement using high-throughput sequencing platforms. Current methods have limitations due to complicated setups and elaborate library preparation plans. To address this, we introduce central limit theorem (CLT)-managed RNA-seq (CLT-seq), a simple and straightforward homopolymer-sequencing method. In CLT-seq, an anchor-free oligo(dT) primer rapidly binds to and unbinds from anywhere along the poly(A) tail string, leading to position-directed reverse transcription with equal probability. The CLT mechanism enables the synthesized poly(T) lengths, which correspond to the templated segment of the poly(A) tail, to distribute normally. Based on a well-fitted pseudogaussian-derived poly(A)-poly(T) conversion model, the actual poly(A)-tail profile is reconstructed from the acquired poly(T)-length profile through matrix operations. CLT-seq follows a simple procedure without requiring RNA-related pre-treatment, enrichment or selection, and the CLT-shortened poly(T) stretches are more compatible with existing sequencing platforms. This proof-of-concept approach facilitates direct homopolymer base-calling and features unbiased RNA-seq. Therefore, CLT-seq provides unbiased, robust and cost-efficient transcriptome-wide poly(A)-tail profiling. We demonstrate that CLT-seq on the most common Illumina platform delivers reliable poly(A)-tail profiling at a transcriptome-wide scale in human cellular contexts. We find that the poly(A)-tail-tuned ncRNA regulation undergoes a dynamic, complex process similar to mRNA regulation. Overall, CLT-seq offers a simplified, effective and economical approach to investigate poly(A)-tail regulation, with potential implications for understanding gene expression and identifying therapeutic targets.

A unique poly(A) tail profile uncovers the stability and translational activation of TOP transcripts during neuronal differentiation

Cell differentiation is associated with global changes in translational activity. Here, we characterize how mRNA poly(A) tail processing supports this dynamic. We observe that decreased translation during neuronal differentiation of P19 cells correlates with the downregulation of 5'-terminal oligopyrimidine (TOP) transcripts which encode the translational machinery. Despite their downregulation, TOP transcripts remain highly stable and show increased translation as cells differentiate. Changes in TOP mRNA metabolism are reflected by their accumulation with poly(A) tails ∼60-nucleotide (nt) long. The dynamic changes in poly(A) processing can be partially recapitulated by depleting LARP1 or activating the mTOR pathway in undifferentiated cells. Although mTOR-induced accumulation of TOP mRNAs with tails ∼60-nt long does not trigger differentiation, it is associated with reduced proliferation of neuronal progenitors. We propose that while TOP mRNAs are transcriptionally silenced, their post-transcriptional regulation mediated by a specific poly(A) processing ensures an adequate supply of ribosomes to complete differentiation.

CPEB and translational control by cytoplasmic polyadenylation: impact on synaptic plasticity, learning, and memory

The late 1990s were banner years in molecular neuroscience; seminal studies demonstrated that local protein synthesis, at or near synapses, was necessary for synaptic plasticity, the underlying cellular basis of learning and memory [1, 2]. The newly made proteins were proposed to "tag" the stimulated synapse, distinguishing it from naive synapses, thereby forming a cellular memory [3]. Subsequent studies demonstrated that the transport of mRNAs from soma to dendrite was linked with translational unmasking at synapses upon synaptic stimulation. It soon became apparent that one prevalent mechanism governing these events is cytoplasmic polyadenylation, and that among the proteins that control this process, CPEB, plays a central role in synaptic plasticity, and learning and memory. In vertebrates, CPEB is a family of four proteins, all of which regulate translation in the brain, that have partially overlapping functions, but also have unique characteristics and RNA binding properties that make them control different aspects of higher cognitive function. Biochemical analysis of the vertebrate CPEBs demonstrate them to respond to different signaling pathways whose output leads to specific cellular responses. In addition, the different CPEBs, when their functions go awry, result in pathophysiological phenotypes resembling specific human neurological disorders. In this essay, we review key aspects of the vertebrate CPEB proteins and cytoplasmic polyadenylation within the context of brain function.

Single-cell quantification of ribosome occupancy in early mouse development

Translation regulation is critical for early mammalian embryonic development¹. However, previous studies had been restricted to bulk measurements², precluding precise determination of translation regulation including allele-specific analyses. Here, to address this challenge, we developed a novel microfluidic isotachophoresis (ITP) approach, named RIBOsome profiling via ITP (Ribo-ITP), and characterized translation in single oocytes and embryos during early mouse development. We identified differential translation efficiency as a key mechanism regulating genes involved in centrosome organization and N⁶-methyladenosine modification of RNAs. Our high-coverage measurements enabled, to our knowledge, the first analysis of allele-specific ribosome engagement in early development. These led to the discovery of stage-specific differential engagement of zygotic RNAs with ribosomes and reduced translation efficiency of transcripts exhibiting allele-biased expression. By integrating our measurements with proteomics data, we discovered that ribosome occupancy in germinal vesicle-stage oocytes is the predominant determinant of protein abundance in the zygote. The Ribo-ITP approach will enable numerous applications by providing high-coverage and high-resolution ribosome occupancy measurements from ultra-low input samples including single cells.

Germ cell-specific eIF4E1b regulates maternal mRNA translation to ensure zygotic genome activation

Translation of maternal mRNAs is detected before transcription of zygotic genes and is essential for mammalian embryo development. How certain maternal mRNAs are selected for translation instead of degradation and how this burst of translation affects zygotic genome activation remain unknown. Using gene-edited mice, we document that the oocyte-specific eukaryotic translation initiation factor 4E family member 1b (eIF4E1b) is the regulator of maternal mRNA expression that ensures subsequent reprogramming of the zygotic genome. In oocytes, eIF4E1b binds to transcripts encoding translation machinery proteins, chromatin remodelers, and reprogramming factors to promote their translation in zygotes and protect them from degradation. The protein products are thought to establish an open chromatin landscape in one-cell zygotes to enable transcription of genes required for cleavage stage development. Our results define a program for rapid resetting of the zygotic epigenome that is regulated by maternal mRNA expression and provide new insights into the mammalian maternal-to-zygotic transition.

Selective Translation of Maternal mRNA by eIF4E1B Controls Oocyte to Embryo Transition

Maternal messenger ribonucleic acids (mRNAs) are driven by a highly orchestrated scheme of recruitment to polysomes and translational activation. However, selecting and regulating individual mRNAs for the translation from a competitive pool of mRNAs are little-known processes. This research shows that the maternal eukaryotic translation initiation factor 4e1b (Eif4e1b) expresses during the oocyte-to-embryo transition (OET), and maternal deletion of Eif4e1b leads to multiple defects concerning oogenesis and embryonic developmental competence during OET. The linear amplification of complementary deoxyribonucleic acid (cDNA) ends, and sequencing (LACE-seq) is used to identify the distinct subset of mRNA and its CG-rich binding sites within the 5' untranslated region (UTR) targeted by eIF4E1B. The proteomics analyses indicate that eIF4E1B-specific bound genes show stronger downregulation at the protein level, which further verify a group of proteins that plays a crucial role in oocyte maturation and embryonic developmental competence is insufficiently synthesized in Eif4e1b-cKO oocytes during OET. Moreover, the biochemical results in vitro are combined to further confirm the maternal-specific translation activation model assembled by eIF4E1B and 3'UTR-associated mRNA binding proteins. The findings demonstrate the indispensability of eIF4E1B for selective translation activation in mammalian oocytes and provide a potential network regulated by eIF4E1B in OET.

Start of life controlled by poly(A) tail-mediated remodeling

Understanding a remarkable event at the start of life, the oocyte-to-embryo transition (OET), has remained elusive, especially in humans. Using newly developed techniques, Liu et al. showed that human maternal mRNAs undergo global poly(A) tail-mediated remodeling during OET, identified the enzymes involved, and demonstrated the essentiality of remodeling for embryo cleavage.

DIS3L2 ribonuclease degrades terminal-uridylated RNA to ensure oocyte maturation and female fertility

During oocyte development in mice, transcripts accumulate in the growth phase and are subsequently degraded during maturation. At the transition point between growth and maturation, oocytes have an intact nucleus or germinal vesicle (GV), and terminal uridylation labels RNA for degradation in meiosis I. By profiling the transcriptome using single-oocyte long-read PacBio RNA sequencing, we document that a small cohort of mRNAs are polyadenylated after terminal uridylation in GV oocytes [designated uridylated-poly(A) RNA]. Because DIS3L2 ribonuclease is known to degrade uridylated transcripts, we established oocyte-specific Dis3l2 knockout mice (Dis3l2cKO). Upon DIS3L2 depletion, uridylated-poly(A) RNAs remain intact which increases their abundance, and they predominate in the transcriptome of Dis3l2cKO oocytes. The abundance of uridylated-poly(A) RNA in Dis3l2cKO oocytes arises not only from insufficient degradation, but also from the stabilizing effect of subsequent polyadenylation. Uridylated-poly(A) RNAs have shorter poly(A) tails and their translation activity decreases in Dis3l2cKO oocytes. Almost all Dis3l2cKO oocytes arrest at the GV stage, and female mice are infertile. Our study demonstrates multiple fates for RNA after terminal uridylation and highlights the role of DIS3L2 ribonuclease in safeguarding the transcriptome and ensuring female fertility.

Remodeling of maternal mRNA through poly(A) tail orchestrates human oocyte-to-embryo transition

Poly(A)-tail-mediated post-transcriptional regulation of maternal mRNAs is vital in the oocyte-to-embryo transition (OET). Nothing is known about poly(A) tail dynamics during the human OET. Here, we show that poly(A) tail length and internal non-A residues are highly dynamic during the human OET, using poly(A)-inclusive RNA isoform sequencing (PAIso-seq). Unexpectedly, maternal mRNAs undergo global remodeling: after deadenylation or partial degradation into 3'-UTRs, they are re-polyadenylated to produce polyadenylated degradation intermediates, coinciding with massive incorporation of non-A residues, particularly internal long consecutive U residues, into the newly synthesized poly(A) tails. Moreover, TUT4 and TUT7 contribute to the incorporation of these U residues, BTG4-mediated deadenylation produces substrates for maternal mRNA re-polyadenylation, and TENT4A and TENT4B incorporate internal G residues. The maternal mRNA remodeling is further confirmed using PAIso-seq2. Importantly, maternal mRNA remodeling is essential for the first cleavage of human embryos. Together, these findings broaden our understanding of the post-transcriptional regulation of maternal mRNAs during the human OET.

MAEL gene contributes to bovine testicular development through the m5C-mediated splicing

Knowledge of RNA molecules regulating testicular development and spermatogenesis in bulls is essential for elite bull selection and an ideal breeding program. Herein, we performed direct RNA sequencing (DRS) to explore the functional characterization of RNA molecules produced in the testicles of 9 healthy Simmental bulls at three testicular development stages (prepuberty, puberty, and postpuberty). We identified 5,043 differentially expressed genes associated with testicular weight. These genes exhibited more alternative splicing at sexual maturity, particularly alternative 3' (A3) and 5' (A5) splice sites usage and exon skipping (SE). The expression of hub genes in testicular developmental stages was also affected by both m6A and m5C RNA modifications. We found m5C-mediated splicing events significantly (p < 0.05) increased MAEL gene expression at the isoform level, likely promoting spermatogenesis. Our findings highlight the complexity of RNA processing and expression as well as the regulation of transcripts involved in testicular development and spermatogenesis.

Measuring the tail: Methods for poly(A) tail profiling

The 3'-end poly(A) tail is an important and potent feature of most mRNA molecules that affects mRNA fate and translation efficiency. Polyadenylation is a posttranscriptional process that occurs in the nucleus by canonical poly(A) polymerases (PAPs). In some specific instances, the poly(A) tail can also be extended in the cytoplasm by noncanonical poly(A) polymerases (ncPAPs). This epitranscriptomic regulation of mRNA recently became one of the most interesting aspects in the field. Advances in RNA sequencing technologies and software development have allowed the precise measurement of poly(A) tails, identification of new ncPAPs, expansion of the function of known enzymes, discovery and a better understanding of the physiological role of tail heterogeneity, and recognition of a correlation between tail length and RNA translatability. Here, we summarize the development of polyadenylation research methods, including classic low-throughput approaches, Illumina-based genome-wide analysis, and advanced state-of-art techniques that utilize long-read third-generation sequencing with Pacific Biosciences and Oxford Nanopore Technologies platforms. A boost in technical opportunities over recent decades has allowed a better understanding of the regulation of gene expression at the mRNA level. This article is categorized under: RNA Methods > RNA Analyses In Vitro and In Silico.

Nano3P-seq: transcriptome-wide analysis of gene expression and tail dynamics using end-capture nanopore cDNA sequencing

RNA polyadenylation plays a central role in RNA maturation, fate, and stability. In response to developmental cues, polyA tail lengths can vary, affecting the translation efficiency and stability of mRNAs. Here we develop Nanopore 3' end-capture sequencing (Nano3P-seq), a method that relies on nanopore cDNA sequencing to simultaneously quantify RNA abundance, tail composition, and tail length dynamics at per-read resolution. By employing a template-switching-based sequencing protocol, Nano3P-seq can sequence RNA molecule from its 3' end, regardless of its polyadenylation status, without the need for PCR amplification or ligation of RNA adapters. We demonstrate that Nano3P-seq provides quantitative estimates of RNA abundance and tail lengths, and captures a wide diversity of RNA biotypes. We find that, in addition to mRNA and long non-coding RNA, polyA tails can be identified in 16S mitochondrial ribosomal RNA in both mouse and zebrafish models. Moreover, we show that mRNA tail lengths are dynamically regulated during vertebrate embryogenesis at an isoform-specific level, correlating with mRNA decay. Finally, we demonstrate the ability of Nano3P-seq in capturing non-A bases within polyA tails of various lengths, and reveal their distribution during vertebrate embryogenesis. Overall, Nano3P-seq is a simple and robust method for accurately estimating transcript levels, tail lengths, and tail composition heterogeneity in individual reads, with minimal library preparation biases, both in the coding and non-coding transcriptome.

Stage-specific requirement for METTL3-dependent m6A modification during dental pulp stem cell differentiation

BACKGROUND

N⁶-methyladenosine (m⁶A) is the most prevalent epigenetic modification in eukaryotic messenger RNAs and plays a critical role in cell fate transition. However, it remains to be elucidated how m⁶A marks functionally impact the transcriptional cascades that orchestrate stem cell differentiation. The present study focuses on the biological function and mechanism of m⁶A methylation in dental pulp stem cell (DPSC) differentiation.

METHODS

m⁶A RNA immunoprecipitation sequencing was utilized to assess the m⁶A-mRNA landscape during DPSC differentiation. Ectopic transplantation of DPSCs in immunodeficient mice was conducted to verify the in vitro findings. RNA sequencing and m⁶A RNA immunoprecipitation sequencing were combined to identify the candidate targets. RNA immunoprecipitation and RNA/protein stability of Noggin (NOG) were evaluated. The alteration in poly(A) tail was measured by 3'-RACE and poly(A) tail length assays.

RESULTS

We characterized a dynamic m⁶A-mRNA landscape during DPSC mineralization with increasing enrichment in the 3' untranslated region (UTR). Methyltransferase-like 3 (METTL3) was identified as the key m⁶A player, and METTL3 knockdown disrupted functional DPSC differentiation. Moreover, METTL3 overexpression enhanced DPSC mineralization. Increasing m⁶A deposition in the 3' UTR restricted NOG expression, which is required for DPSC mineralization. This stage-specific m⁶A methylation and destabilization of NOG was suppressed by METTL3 knockdown only in differentiated DPSCs. Furthermore, METTL3 promotes the degradation of m⁶A-tagged NOG by shortening the poly(A) tail length in the differentiated stage.

CONCLUSIONS

Our results address an essential role of dynamic m⁶A signaling in the temporal control of DPSC differentiation and provide new insight into epitranscriptomic mechanisms in stem cell-based therapy.

PABPN1 functions as a hub in the assembly of nuclear poly(A) domains that are essential for mouse oocyte development

Growing oocytes store a large amount of maternal mRNA to support the subsequent "maternal-zygotic transition" process. At present, it is not clear how the growing oocytes store and process the newly transcribed mRNA under physiological conditions. In this study, we report non-membrane-bound compartments, nuclear poly(A) domains (NPADs), as the hub for newly transcribed mRNA, in developing mouse oocytes. The RNA binding protein PABPN1 promotes the formation of NPAD through its N-terminal disordered domain and RNA-recognized motif by means of liquid phase separation. Pabpn1-null growing oocytes cannot form NPAD normally in vivo and have defects in stability of oocyte growing-related transcripts and formation of long 3' untranslated region isoform transcripts. Ultimately, Pabpn1^fl/fl;Gdf9-Cre mice are completely sterile with primary ovarian insufficiency. These results demonstrate that NPAD formed by the phase separation properties of PABPN1-mRNA are the hub of the newly transcribed mRNA and essential for the development of oocytes and female reproduction.

An atlas of plant full-length RNA reveals tissue-specific and monocots-dicots conserved regulation of poly(A) tail length

Poly(A) tail is a hallmark of eukaryotic messenger RNA and its length plays an essential role in regulating mRNA metabolism. However, a comprehensive resource for plant poly(A) tail length has yet to be established. Here, we applied a poly(A)-enrichment-free, nanopore-based method to profile full-length RNA with poly(A) tail information in plants. Our atlas contains over 120 million polyadenylated mRNA molecules from seven different tissues of Arabidopsis, as well as the shoot tissue of maize, soybean and rice. In most tissues, the size of plant poly(A) tails shows peaks at approximately 20 and 45 nucleotides, while the poly(A) tails in pollen exhibit a distinct pattern with strong peaks centred at 55 and 80 nucleotides. Moreover, poly(A) tail length is regulated in a gene-specific manner-mRNAs with short half-lives in general have long poly(A) tails, while mRNAs with long half-lives are featured with relatively short poly(A) tails that peak at ~45 nucleotides. Across species, poly(A) tails in the nucleus are almost twice as long as in the cytoplasm. Our comprehensive dataset lays the groundwork for future functional and evolutionary studies on poly(A) tail length regulation in plants.

Drought induces epitranscriptome and proteome changes in stem-differentiating xylem of Populus trichocarpa

Understanding gene expression and regulation requires insights into RNA transcription, processing, modification, and translation. However, the relationship between the epitranscriptome and the proteome under drought stress remains undetermined in poplar (Populus trichocarpa). In this study, we used Nanopore direct RNA sequencing and tandem mass tag-based proteomic analysis to examine epitranscriptomic and proteomic regulation induced by drought treatment in stem-differentiating xylem (SDX). Our results revealed a decreased full-length read ratio under drought treatment and, especially, a decreased association between transcriptome and proteome changes in response to drought. Epitranscriptome analysis of cellulose- and lignin-related genes revealed an increased N6-Methyladenosine (m6A) ratio, which was accompanied by decreased RNA abundance and translation, under drought stress. Interestingly, usage of the distal poly(A) site increased during drought stress. Finally, we found that transcripts of highly expressed genes tend to have shorter poly(A) tail length (PAL), and drought stress increased the percentage of transcripts with long PAL. These findings provide insights into the interplay among m6A, polyadenylation, PAL, and translation under drought stress in P. trichocarpa SDX.

Chemically Modified Poly(A) Analogs Targeting PABP: Structure Activity Relationship and Translation Inhibitory Properties

Poly(A)‐binding protein (PABP) is an essential element of cellular translational machinery. Recent studies have revealed that poly(A) tail modifications can modulate mRNA stability and translational potential, and that oligoadenylate‐derived PABP ligands can act as effective translational inhibitors with potential applications in pain management. Although extensive research has focused on protein‐RNA and protein‐protein interactions involving PABPs, further studies are required to examine the ligand specificity of PABP. In this study, we developed a microscale thermophoresis‐based assay to probe the interactions between PABP and oligoadenylate analogs containing different chemical modifications. Using this method, we evaluated oligoadenylate analogs modified with nucleobase, ribose, and phosphate moieties to identify modification hotspots. In addition, we determined the susceptibility of the modified oligos to CNOT7 to identify those with the potential for increased cellular stability. Consequently, we selected two enzymatically stable oligoadenylate analogs that inhibit translation in rabbit reticulocyte lysates with a higher potency than a previously reported PABP ligand. We believe that the results presented in this study and the implemented methodology can be capitalized upon in the future development of RNA‐based biological tools.

Transcriptome-wide measurement of poly(A) tail length and composition at subnanogram total RNA sensitivity by PAIso-seq

Poly(A) tails are added to the 3' ends of most mRNAs in a non-templated manner and play essential roles in post-transcriptional regulation, including mRNA export, stability and translation. Measuring poly(A) tails is critical for understanding their regulatory roles in almost every aspect of biological and medical studies. Previous methods for analyzing poly(A) tails require large amounts of input RNA (microgram-level total RNA), which limits their application. We recently developed a poly(A) inclusive full-length RNA isoform-sequencing method (PAIso-seq) at single-oocyte-level sensitivity (a single mammalian oocyte contains ~0.5 ng of total RNA) based on PacBio sequencing that enabled accurate measurement of the poly(A) tail length and non-A residues within the body of poly(A) tails along with the full-length cDNA, providing the opportunity to study precious in vivo samples with very limited input material. Here, we describe a detailed protocol for PAIso-seq library preparation from single mouse oocytes or bulk oocyte samples. In addition, we provide a complete bioinformatic pipeline to perform the analysis from the raw data to downstream analysis. The minimum time required is ~14.5 h for PAIso-seq double-stranded cDNA preparation, 2 d for PacBio sequencing in HiFi mode and 8 h for the initial data analysis.

Novel Mobile Phase to Control Charge States and Metal Adducts in the LC/MS for mRNA Characterization Assays

Mass spectrometry is a widely used tool in the characterization of oligonucleotides. This analysis can be challenging due to the large number of possible charge states of oligonucleotides, which can limit the sensitivity of the assay, along with the propensity of oligonucleotides to readily form adducts with free alkali metals. To reduce the adduct formation, oligonucleotides are typically purified with desalting columns prior to analysis. We have developed a mobile phase that gives superior reduction in charge states and adduct formation compared to previously reported methods and, more importantly, obviates the requirement of desalting samples prior to mass spectrometric analysis, significantly decreasing the sample preparation time and amount of RNA required for analysis. We have applied this mobile phase to develop methods to quantify the 5'-capping efficiency and to characterize the polyadenosine (poly(A)) tail of mRNA synthesized in vitro: two critical quality attributes of mRNA therapeutics. Through this, we were able to demonstrate RNA that was co-transcriptionally capped to have capping efficiency equivalent (the percent total molecules that contain a cap) to other reports in the literature using materials that were generated using the same synthesis procedure. Furthermore, by using a mobile phase mixture comprised of hexafluoroisopropanol, triethylammonium acetate, triethylamine, and ethanol, we were able to determine the size distribution of the poly(A) tail in various mRNA samples from DNA templates that ranged from 50 to 150 nt poly(A) and verify that distribution with commercially available RNA standards, successfully demonstrating that this mobile phase composition could be used for characterization assays for both mRNA caps and tails.

PolyAtailor: measuring poly(A) tail length from short-read and long-read sequencing data

The poly(A) tail is a dynamic addition to the eukaryotic mRNA and the change in its length plays an essential role in regulating gene expression through affecting nuclear export, mRNA stability and translation. Only recently high-throughput sequencing strategies began to emerge for transcriptome-wide profiling of poly(A) tail length in diverse developmental stages and organisms. However, there is currently no easy-to-use and universal tool for measuring poly(A) tails in sequencing data from different sequencing protocols. Here we established PolyAtailor, a unified and efficient framework, for identifying and analyzing poly(A) tails from PacBio-based long reads or next generation short reads. PolyAtailor provides two core functions for measuring poly(A) tails, namely Tail_map and Tail_scan, which can be used for profiling tails with or without using a reference genome. Particularly, PolyAtailor can identify all potential tails in a read, providing users with detailed information such as tail position, tail length, tail sequence and tail type. Moreover, PolyAtailor integrates rich functions for poly(A) tail and poly(A) site analyses, such as differential poly(A) length analysis, poly(A) site identification and annotation, and statistics and visualization of base composition in tails. We compared PolyAtailor with three latest methods, FLAMAnalysis, FLEPSeq and PAIsoSeqAnalysis, using data from three sequencing protocols in HeLa samples and Arabidopsis. Results show that PolyAtailor is effective in measuring poly(A) tail length and detecting significance of differential poly(A) length, which achieves much higher sensitivity and accuracy than competing methods. PolyAtailor is available at https://github.com/BMILAB/PolyAtailor.

Uncoupling transcription and translation through miRNA-dependent poly(A) length control in haploid male germ cells

As one of the post-transcriptional regulatory mechanisms, uncoupling of transcription and translation plays an essential role in development and adulthood physiology. However, it remains elusive how thousands of mRNAs get translationally silenced while stability is maintained for hours or even days before translation. In addition to oocytes and neurons, developing spermatids display significant uncoupling of transcription and translation for delayed translation. Therefore, spermiogenesis represents an excellent in vivo model for investigating the mechanism underlying uncoupled transcription and translation. Through full-length poly(A) deep sequencing, we discovered dynamic changes in poly(A) length through deadenylation and re-polyadenylation. Deadenylation appeared to be mediated by microRNAs (miRNAs), and transcripts with shorter poly(A) tails tend to be sequestered into ribonucleoprotein (RNP) granules for translational repression and stabilization. In contrast, re-polyadenylation might allow for translocation of the translationally repressed transcripts from RNP granules to polysomes. Overall, our data suggest that miRNA-dependent poly(A) length control represents a previously unreported mechanism underlying uncoupled translation and transcription in haploid male mouse germ cells.

Ultrasensitive Ribo-seq reveals translational landscapes during mammalian oocyte-to-embryo transition and pre-implantation development

In mammals, translational control plays critical roles during oocyte-to-embryo transition (OET) when transcription ceases. However, the underlying regulatory mechanisms remain challenging to study. Here, using low-input Ribo-seq (Ribo-lite), we investigated translational landscapes during OET using 30-150 mouse oocytes or embryos per stage. Ribo-lite can also accommodate single oocytes. Combining PAIso-seq to interrogate poly(A) tail lengths, we found a global switch of translatome that closely parallels changes of poly(A) tails upon meiotic resumption. Translation activation correlates with polyadenylation and is supported by polyadenylation signal proximal cytoplasmic polyadenylation elements (papCPEs) in 3' untranslated regions. By contrast, translation repression parallels global de-adenylation. The latter includes transcripts containing no CPEs or non-papCPEs, which encode many transcription regulators that are preferentially re-activated before zygotic genome activation. CCR4-NOT, the major de-adenylation complex, and its key adaptor protein BTG4 regulate translation downregulation often independent of RNA decay. BTG4 is not essential for global de-adenylation but is required for selective gene de-adenylation and production of very short-tailed transcripts. In sum, our data reveal intimate interplays among translation, RNA stability and poly(A) tail length regulation underlying mammalian OET.

Detection of rare thalassemia mutations using long-read single-molecule real-time sequencing

Gap- polymerase chain reaction (PCR), reverse dot-blot assay (RDB), real-time PCR based multicolor melting curve analysis (MMCA assay), multiplex ligation-dependent probe amplification (MLPA) and Sanger sequencing are conventional methods to diagnose thalassemia but all of them have limitations. In this study, we applied single-molecule real-time (SMRT) sequencing following multiplex long-range PCR to uncover rare mutations in nine patients and their family members. The patients with different results between Gap-PCR and MMCA assay or with phenotype not matching genotype were included. Using SMRT sequencing, we first identified the carriers with ααα^anti3.7/HKαα, -α^762bpα/αα (chr16:172,648-173,409), αα^fusion/α^QSα (in a trans configuration), two cases with novel gene rearrangements and another case with a novel 341 bp insertion in α-globin gene cluster, respectively. One carrier with --^SEA/ααα^anti4.2, and two carriers with the coexistence of globin variant and an α-globin gene duplication were also found. Most importantly, we could determine two defects in α-globin gene cluster being a cis or trans configuration in a single test. Our results showed that SMRT has great advantages in detection of α-globin gene triplications, rare deletions and determination of a cis or trans configuration. SMRT is a comprehensive and one-step method for thalassemia screening and diagnosis, especially for detection of rare thalassemia mutations.

Cnot8 eliminates naïve regulation networks and is essential for naïve-to-formative pluripotency transition

Mammalian early epiblasts at different phases are characterized by naïve, formative, and primed pluripotency states, involving extensive transcriptome changes. Here, we report that deadenylase Cnot8 of Ccr4-Not complex plays essential roles during the transition from naïve to formative state. Knock out (KO) Cnot8 resulted in early embryonic lethality in mice, but Cnot8 KO embryonic stem cells (ESCs) could be established. Compared with the cells differentiated from normal ESCs, Cnot8 KO cells highly expressed a great many genes during their differentiation into the formative state, including several hundred naïve-like genes enriched in lipid metabolic process and gene expression regulation that may form the naïve regulation networks. Knockdown expression of the selected genes of naïve regulation networks partially rescued the differentiation defects of Cnot8 KO ESCs. Cnot8 depletion led to the deadenylation defects of its targets, increasing their poly(A) tail lengths and half-life, eventually elevating their expression levels. We further found that Cnot8 was involved in the clearance of targets through its deadenylase activity and the binding of Ccr4-Not complex, as well as the interacting with Tob1 and Pabpc1. Our results suggest that Cnot8 eliminates naïve regulation networks through mRNA clearance, and is essential for naïve-to-formative pluripotency transition.

OUP accepted manuscript

Alternative polyadenylation increases transcript diversities at the 3’ end, regulating biological processes including cell differentiation, embryonic development and cancer progression. Here, we present a Bayesian method SCAPE, which enables de novo identification and quantification of polyadenylation (pA) sites at single-cell level by utilizing insert size information. We demonstrated its accuracy and robustness and identified 31 558 sites from 36 mouse organs, 43.8% (13 807) of which were novel. We illustrated that APA isoforms were associated with miRNAs binding and regulated in tissue-, cell type-and tumor-specific manners where no difference was found at gene expression level, providing an extra layer of information for cell clustering. Furthermore, we found genome-wide dynamic changes of APA usage during erythropoiesis and induced pluripotent stem cell (iPSC) differentiation, suggesting APA contributes to the functional flexibility and diversity of single cells. We expect SCAPE to aid the analyses of cellular dynamics and diversities in health and disease.

Landscape of transcription termination in Arabidopsis revealed by single-molecule nascent RNA sequencing

BACKGROUND

The dynamic process of transcription termination produces transient RNA intermediates that are difficult to distinguish from each other via short-read sequencing methods.

RESULTS

Here, we use single-molecule nascent RNA sequencing to characterize the various forms of transient RNAs during termination at genome-wide scale in wildtype Arabidopsis and in atxrn3, fpa, and met1 mutants. Our data reveal a wide range of termination windows among genes, ranging from ~ 50 nt to over 1000 nt. We also observe efficient termination before downstream tRNA genes, suggesting that chromatin structure around the promoter region of tRNA genes may block pol II elongation. 5' Cleaved readthrough transcription in atxrn3 with delayed termination can run into downstream genes to produce normally spliced and polyadenylated mRNAs in the absence of their own transcription initiation. Consistent with previous reports, we also observe long chimeric transcripts with cryptic splicing in fpa mutant; but loss of CG DNA methylation has no obvious impact on termination in the met1 mutant.

CONCLUSIONS

Our method is applicable to establish a comprehensive termination landscape in a broad range of species.

Post-transcriptional regulation in spermatogenesis: all RNA pathways lead to healthy sperm

Multiple RNA pathways are required to produce functional sperm. Here, we review RNA post-transcriptional regulation during spermatogenesis with particular emphasis on the role of 3' end modifications. From early studies in the 1970s, it became clear that spermiogenesis transcripts could be stored for days only to be translated at advanced stages of spermatid differentiation. The transition between the translationally repressed and active states was observed to correlate with the shortening of the transcripts' poly(A) tail, establishing a link between RNA 3' end metabolism and male germ cell differentiation. Since then, numerous RNA metabolic pathways have been implicated not only in the progression through spermatogenesis, but also in the maintenance of genomic integrity. Recent studies have characterized the elusive 3' biogenesis of Piwi-interacting RNAs (piRNAs), identified a critical role for messenger RNA (mRNA) 3' uridylation in meiotic progression, established the mechanisms that destabilize transcripts with long 3' untranslated regions (3'UTRs) in post-mitotic cells, and defined the physiological relevance of RNA exonucleases and deadenylases in male germ cells. In this review, we discuss RNA processing in the male germline in the light of the most recent findings. A brief recollection of different RNA-processing events will aid future studies exploring post-transcriptional regulation in spermatogenesis.

Pervasive 3'-UTR Isoform Switches During Mouse Oocyte Maturation

Oocyte maturation is the foundation for developing healthy individuals of mammals. Upon germinal vesicle breakdown, oocyte meiosis resumes and the synthesis of new transcripts ceases. To quantitatively profile the transcriptomic dynamics after meiotic resumption throughout the oocyte maturation, we generated transcriptome sequencing data with individual mouse oocytes at three main developmental stages: germinal vesicle (GV), metaphase I (MI), and metaphase II (MII). When clustering the sequenced oocytes, results showed that isoform-level expression analysis outperformed gene-level analysis, indicating isoform expression provided extra information that was useful in distinguishing oocyte stages. Comparing transcriptomes of the oocytes at the GV stage and the MII stage, in addition to identification of differentially expressed genes (DEGs), we detected many differentially expressed transcripts (DETs), some of which came from genes that were not identified as DEGs. When breaking down the isoform-level changes into alternative RNA processing events, we found the main source of isoform composition changes was the alternative usage of polyadenylation sites. With detailed analysis focusing on the alternative usage of 3′-UTR isoforms, we identified, out of 3,810 tested genes, 512 (13.7%) exhibiting significant switches of 3′-UTR isoforms during the process of moues oocyte maturation. Altogether, our data and analyses suggest the importance of examining isoform abundance changes during oocyte maturation, and further investigation of the pervasive 3′-UTR isoform switches in the transition may deepen our understanding on the molecular mechanisms underlying mammalian early development.

A comprehensive map of alternative polyadenylation in African American and European American lung cancer patients

Deciphering the post-transcriptional mechanisms (PTM) regulating gene expression is critical to understand the dynamics underlying transcriptomic regulation in cancer. Alternative polyadenylation (APA)-regulation of mRNA 3'UTR length by alternating poly(A) site usage-is a key PTM mechanism whose comprehensive analysis in cancer remains an important open challenge. Here we use a method and analysis pipeline that sequences 3'end-enriched RNA directly to overcome the saturation limitation of traditional 5'-3' based sequencing. We comprehensively map the APA landscape in lung cancer in a cohort of 98 tumor/non-involved tissues derived from European American and African American patients. We identify a global shortening of 3'UTR transcripts in lung cancer, with notable functional implications on the expression of both coding and noncoding genes. We find that APA of non-coding RNA transcripts (long non-coding RNAs and microRNAs) is a recurrent event in lung cancer and discover that the selection of alternative polyA sites is a form of non-coding RNA expression control. Our results indicate that mRNA transcripts from EAs are two times more likely than AAs to undergo APA in lung cancer. Taken together, our findings comprehensively map and identify the important functional role of alternative polyadenylation in determining transcriptomic heterogeneity in lung cancer.

FLEP-seq: simultaneous detection of RNA polymerase II position, splicing status, polyadenylation site and poly(A) tail length at genome-wide scale by single-molecule nascent RNA sequencing

Elongation, splicing and polyadenylation are fundamental steps of transcription, and studying their coordination requires simultaneous monitoring of these dynamic processes on one transcript. We recently developed a full-length nascent RNA sequencing method in the model plant Arabidopsis that simultaneously detects RNA polymerase II position, splicing status, polyadenylation site and poly(A) tail length at genome-wide scale. This method allows calculation of the kinetics of cotranscriptional splicing and detects polyadenylated transcripts with unspliced introns retained at specific positions posttranscriptionally. Here we describe a detailed protocol for this method called FLEP-seq (full-length elongating and polyadenylated RNA sequencing) that is applicable to plants. Library production requires as little as one nanogram of nascent RNA (after rRNA/tRNA removal), and either Nanopore or PacBio platforms can be used for sequencing. We also provide a complete bioinformatic pipeline from raw data processing to downstream analysis. The minimum time required for FLEP-seq, including RNA extraction and library preparation, is 36 h. The subsequent long-read sequencing and initial data analysis ranges between 31 and 40 h, depending on the sequencing platform.

DNA transposons mediate duplications via transposition-independent and -dependent mechanisms in metazoans

Despite long being considered as "junk", transposable elements (TEs) are now accepted as catalysts of evolution. One example is Mutator-like elements (MULEs, one type of terminal inverted repeat DNA TEs, or TIR TEs) capturing sequences as Pack-MULEs in plants. However, their origination mechanism remains perplexing, and whether TIR TEs mediate duplication in animals is almost unexplored. Here we identify 370 Pack-TIRs in 100 animal reference genomes and one Pack-TIR (Ssk-FB4) family in fly populations. We find that single-copy Pack-TIRs are mostly generated via transposition-independent gap filling, and multicopy Pack-TIRs are likely generated by transposition after replication fork switching. We show that a proportion of Pack-TIRs are transcribed and often form chimeras with hosts. We also find that Ssk-FB4s represent a young protein family, as supported by proteomics and signatures of positive selection. Thus, TIR TEs catalyze new gene structures and new genes in animals via both transposition-independent and -dependent mechanisms.

Single-molecule polyadenylated tail sequencing (SM-PAT-Seq) to measure polyA tail lengths transcriptome-wide

Polyadenylation of the 3' end of mRNAs is an important mechanism for regulating their stability and translation. We developed a nucleotide-resolution, transcriptome-wide, single-molecule SM-PAT-Seq method to accurately measure the polyA tail lengths of individual transcripts using long-read sequencing. The method generates cDNA using a double stranded splint adaptor targeting the far 3' end of the polyA tail for first strand synthesis along with random hexamers for second strand synthesis. This straight-forward method yields accurate polyA tail sequence lengths, can identify non-A residues in those tails, and quantitate transcript abundance.