Extensive alternative polyadenylation during zebrafish development

Benchmarking of methods that identify alternative polyadenylation events in single-/multiple-polyadenylation site genes

Alternative polyadenylation (APA) is a widespread post-transcriptional mechanism that diversifies gene expression by generating messenger RNA isoforms with varying 3' untranslated regions. Accurate identification and quantification of transcriptome-wide polyadenylation site (PAS) usage are essential for understanding APA-mediated gene regulation and its biological implications. In this review, we first review the landscape of computational tools developed to identify APA events from RNA sequencing (RNA-seq) data. We then benchmarked five PAS prediction tools and seven APA detection algorithms using five RNA-seq datasets derived from clear cell renal cell carcinoma (ccRCC) and adjacent normal tissues. By evaluating tool performance across genes with either single or multiple PASs, we revealed substantial variation in accuracy, sensitivity, and consistency among the tools. Based on this comparative analysis, we offer practical guidelines for tool selection and propose considerations for improving APA detection accuracy. Additionally, our analysis identified CCNL2 as a candidate gene exhibiting significant APA regulation in ccRCC, highlighting its potential as a disease-associated biomarker.

The Dynamic Landscape of 3'-UTR Alternative Polyadenylation Across Mouse Fetal Development and Anatomy

Alternative cleavage and polyadenylation (APA) in the 3'-untranslated region (3'-UTR) of mRNA produces transcripts with varied 3'-UTR and plays a key role in development and organogenesis. This work characterizes 3'-UTR APA using 85 high-quality RNA-seq datasets encompassing 12 tissue types and eight developmental stages of mouse fetuses. Results show that 46.4% of expressed genes undergo APA in a tissue-specific manner. Changes in polyadenylation site (pAS) usage often operate beyond transcriptional control, revealing APA as an additional layer of gene regulation. Sequence analysis demonstrates that pAS selection, governed by polyadenylation signal strength and adenine preferences, is evolutionarily conserved between mice and humans. Intriguingly, brain tissues display complex 3'-UTR APA dynamics during development, potentially regulated by RNA-binding proteins such as Rbm38, potentially impacting 3' UTR extension by restricting distal pAS usage. These APA events are associated with a depletion of conserved miRNA binding sites and an enrichment of transposable elements within alternative 3' UTRs. To facilitate further research, this work develops APApedia (http://xozhanglab.com/apapedia/), a comprehensive database cataloging identified 3'-UTR APA events, which serves as a valuable resource for the community to study APA in development and tissue-specific regulation. This comprehensive resource aids in deciphering the functional implications of APA in mouse fetal development.

Pumilio differentially binds to mRNA 3' UTR isoforms to regulate localization of synaptic proteins

In neuronal cells, the regulation of RNA is crucial for the spatiotemporal control of gene expression, but how the correct localization, levels, and function of synaptic proteins are achieved is not well understood. In this study, we globally investigate the role of alternative 3' UTRs in regulating RNA localization in the synaptic regions of the Drosophila brain. We identify direct mRNA targets of the translational repressor Pumilio, finding that mRNAs bound by Pumilio encode proteins enriched in synaptosomes. Pumilio differentially binds to RNA isoforms of the same gene, favoring long, neuronal 3' UTRs. These longer 3' UTRs tend to remain in the neuronal soma, whereas shorter UTR isoforms localize to the synapse. In cultured pumilio mutant neurons, axon outgrowth defects are accompanied by mRNA isoform mislocalization, and proteins encoded by these Pumilio target mRNAs display excessive abundance at synaptic boutons. Our study identifies an important mechanism for the spatiotemporal regulation of protein function in neurons.

The maternal-to-zygotic transition: reprogramming of the cytoplasm and nucleus

A fertilized egg is initially transcriptionally silent and relies on maternally provided factors to initiate development. For embryonic development to proceed, the oocyte-inherited cytoplasm and the nuclear chromatin need to be reprogrammed to create a permissive environment for zygotic genome activation (ZGA). During this maternal-to-zygotic transition (MZT), which is conserved in metazoans, transient totipotency is induced and zygotic transcription is initiated to form the blueprint for future development. Recent technological advances have enhanced our understanding of MZT regulation, revealing common themes across species and leading to new fundamental insights about transcription, mRNA decay and translation.

Using Zebrafish Models to Study Epitranscriptomic Regulation of CNS Functions

Epitranscriptomic regulation of cell functions involves multiple post-transcriptional chemical modifications of coding and non-coding RNA that are increasingly recognized in studying human brain disorders. Although rodent models are presently widely used in neuroepitranscriptomic research, the zebrafish (Danio rerio) has emerged as a useful and promising alternative model species. Mounting evidence supports the importance of RNA modifications in zebrafish CNS function, providing additional insights into epitranscriptomic mechanisms underlying a wide range of brain disorders. Here, we discuss recent data on the role of RNA modifications in CNS regulation, with a particular focus on zebrafish models, as well as evaluate current problems, challenges, and future directions of research in this field of molecular neurochemistry.

InPAS: An R/Bioconductor Package for Identifying Novel Polyadenylation Sites and Alternative Polyadenylation from Bulk RNA-seq Data

BACKGROUND

Alternative cleavage and polyadenylation (APA) is a crucial post-transcriptional gene regulation mechanism that regulates gene expression in eukaryotes by increasing the diversity and complexity of both the transcriptome and proteome. Despite the development of more than a dozen experimental methods over the last decade to identify and quantify APA events, widespread adoption of these methods has been limited by technical, financial, and time constraints. Consequently, APA remains poorly understood in most eukaryotes. However, RNA sequencing (RNA-seq) technology has revolutionized transcriptome profiling and recent studies have shown that RNA-seq data can be leveraged to identify and quantify APA events.

RESULTS

To fully capitalize on the exponentially growing RNA-seq data, we developed InPAS (Identification of Novel alternative PolyAdenylation Sites), an R/Bioconductor package for accurate identification of novel and known cleavage and polyadenylation sites (CPSs), as well as quantification of APA from RNA-seq data of various experimental designs. Compared to other APA analysis tools, InPAS offers several important advantages, including the ability to detect both novel proximal and distal CPSs, to fine tune positions of CPSs using a naïve Bayes classifier based on flanking sequence features, and to identify APA events from RNA-seq data of complex experimental designs using linear models. We benchmarked the performance of InPAS and other leading tools using simulated and experimental RNA-seq data with matched 3'-end RNA-seq data. Our results reveal that InPAS frequently outperforms existing tools in terms of precision, sensitivity, and specificity. Furthermore, we demonstrate its scalability and versatility by applying it to large, diverse RNA-seq datasets.

CONCLUSIONS

InPAS is an efficient and robust tool for identifying and quantifying APA events using readily accessible conventional RNA-seq data. Its versatility opens doors to explore APA regulation across diverse eukaryotic systems with various experimental designs. We believe that InPAS will drive APA research forward, deepening our understanding of its role in regulating gene expression, and potentially leading to the discovery of biomarkers or therapeutics for diseases.

Alternative Splicing in the Heart: The Therapeutic Potential of Regulating the Regulators

Alternative splicing allows a single gene to produce a variety of protein isoforms. Changes in splicing isoform usage characterize virtually every stage of the differentiation process and define the physiological differences between cardiomyocytes with different function, at different stages of development, and pathological function. Recent identification of cardiac splicing factors provided insights into the mechanisms underlying alternative splicing and revealed how these splicing factors impact functional properties of the heart. Alterations of the splicing of sarcomeric genes, cell signaling proteins, and ion channels have been associated with the development of pathological conditions such as cardiomyopathy and arrhythmia. RBM20, RBM24, PTBP1, RBFOX, and QKI play key roles in cardiac development and pathology. A better understanding of their regulation will yield insights into healthy cardiac development and inform the development of molecular therapeutics.

Full-Length Transcriptome Profile of Apis cerana Revealed by Nanopore Sequencing

The Asian honey bee (Apis cerana) plays a crucial role in providing abundant bee products and in maintaining ecological balance. Despite the availability of the genomic sequence of the Asian honey bee, its transcriptomic information remains largely incomplete. To address this issue, here we constructed three pooled RNA samples from the queen, drone, and worker bees of A. cerana and performed full-length RNA sequencing using Nanopore single-molecule sequencing technology. Ultimately, we obtained 160,811 full-length transcript sequences from 19,859 genes, with 141,189 being novel transcripts, of which 130,367 were functionally annotated. We detected 520, 324, and 1823 specifically expressed transcripts in the queen, worker, and drone bees, respectively. Furthermore, we identified 38,799 alternative splicing (AS) events from 5710 genes, 44,243 alternative polyadenylation (APA) sites from 1649 gene loci, 88,187 simple sequence repeats (SSRs), and 17,387 long noncoding RNAs (lncRNAs). Leveraging these transcripts as references, we identified 6672, 7795, and 6804 differentially expressed transcripts (DETs) in comparisons of queen ovaries vs drone testes, worker ovaries vs drone testes, and worker ovaries vs queen ovaries, respectively. Our research results provide a comprehensive set of reference transcript datasets for Apis cerana, offering important sequence information for further exploration of its gene functions.

Dynamic Alternative Polyadenylation during Litopenaeus Vannamei Metamorphosis Development

As an important mechanism in the post-transcriptional regulation of eukaryotic gene expression, alternative polyadenylation (APA) plays a key role in biological processes such as cell proliferation and differentiation. However, the role and dynamic pattern of APA during Litopenaeus vannamei metamorphosis are poorly understood. Here, RNA-seq data covering from the embryo to the maturation (16 time points) of L. vannamei were utilized. We identified 247 differentially expressed APA events between early and adult stages, and through fuzzy mean clustering analysis, we discovered five dynamic APA patterns. Among them, the gradual elongation of the 3'UTR is the major APA pattern that changes over time, and its genes are enriched in the pathways of protein and energy metabolism. Finally, we constructed mRNA-miRNA and PPI networks and detected several central miRNAs that may regulate L. vannamei development. Our results revealed the complex APA mechanisms in L. vannamei metamorphosis, shedding new light on post-transcriptional regulation of crustacean metamorphosis.

The molecular mechanisms underpinning maternal mRNA dormancy

A large number of mRNAs of maternal origin are produced during oogenesis and deposited in the oocyte. Since transcription stops at the onset of meiosis during oogenesis and does not resume until later in embryogenesis, maternal mRNAs are the only templates for protein synthesis during this period. To ensure that a protein is made in the right place at the right time, the translation of maternal mRNAs must be activated at a specific stage of development. Here we summarize our current understanding of the sophisticated mechanisms that contribute to the temporal repression of maternal mRNAs, termed maternal mRNA dormancy. We discuss mechanisms at the level of the RNA itself, such as the regulation of polyadenine tail length and RNA modifications, as well as at the level of RNA-binding proteins, which often block the assembly of translation initiation complexes at the 5' end of an mRNA or recruit mRNAs to specific subcellular compartments. We also review microRNAs and other mechanisms that contribute to repressing translation, such as ribosome dormancy. Importantly, the mechanisms responsible for mRNA dormancy during the oocyte-to-embryo transition are also relevant to cellular quiescence in other biological contexts.

Zebrafish tsc1 and cxcl12a increase susceptibility to mycobacterial infection

Regulation of host miRNA expression is a contested node that controls the host immune response to mycobacterial infection. The host must counter subversive efforts of pathogenic mycobacteria to launch a protective immune response. Here, we examine the role of miR-126 in the zebrafish-Mycobacterium marinum infection model and identify a protective role for infection-induced miR-126 through multiple effector pathways. We identified a putative link between miR-126 and the tsc1a and cxcl12a/ccl2/ccr2 signalling axes resulting in the suppression of non-tnfa expressing macrophage accumulation at early M. marinum granulomas. Mechanistically, we found a detrimental effect of tsc1a expression that renders zebrafish embryos susceptible to higher bacterial burden and increased cell death via mTOR inhibition. We found that macrophage recruitment driven by the cxcl12a/ccl2/ccr2 signalling axis was at the expense of the recruitment of classically activated tnfa-expressing macrophages and increased cell death around granulomas. Together, our results delineate putative pathways by which infection-induced miR-126 may shape an effective immune response to M. marinum infection in zebrafish embryos.

miR-430 regulates zygotic mRNA during zebrafish embryogenesis

BACKGROUND

Early embryonic developmental programs are guided by the coordinated interplay between maternally inherited and zygotically manufactured RNAs and proteins. Although these processes happen concomitantly and affecting gene function during this period is bound to affect both pools of mRNAs, it has been challenging to study their expression dynamics separately.

RESULTS

By employing SLAM-seq, a nascent mRNA labeling transcriptomic approach, in a developmental time series we observe that over half of the early zebrafish embryo transcriptome consists of maternal-zygotic genes, emphasizing their pivotal role in early embryogenesis. We provide an hourly resolution of de novo transcriptional activation events and follow nascent mRNA trajectories, finding that most de novo transcriptional events are stable throughout this period. Additionally, by blocking microRNA-430 function, a key post transcriptional regulator during zebrafish embryogenesis, we directly show that it destabilizes hundreds of de novo transcribed mRNAs from pure zygotic as well as maternal-zygotic genes. This unveils a novel miR-430 function during embryogenesis, fine-tuning zygotic gene expression.

CONCLUSION

These insights into zebrafish early embryo transcriptome dynamics emphasize the significance of post-transcriptional regulators in zygotic genome activation. The findings pave the way for future investigations into the coordinated interplay between transcriptional and post-transcriptional landscapes required for the establishment of animal cell identities and functions.

Gene expression investigation of four key regulators of polyadenylation and alternative adenylation in the periphery of late-onset Alzheimer's disease patients

BACKGROUND

Alzheimer's disease (AD) is a genetic and sporadic neurodegenerative disease considered by an archetypal cognitive impairment and a decrease in less common cognitive impairment. Notably, the discovery of goals in this paradigm is still a challenge, and understanding basic mechanisms is an important step toward improving disease management. Polyadenylation (PA) and alternative polyadenylation (APA) are two of the most critical RNA processing stages in 3'UTRs that influence various AD-related genes.

METHODS

In this study, we assessed Cleavage and polyadenylation specificity factors 1 and 6 (CPSF1 and CPSF6), cleavage stimulation factor 1 (CSTF1), and WD Repeat Domain 33 (WDR33) genes expression in the periphery of 50 AD patients and 50 healthy individuals with age and gender-matched by quantitative real-time PCR.

RESULTS

Comparing AD patients with healthy people using expression analysis revealed a substantial increase in CSTF1 (posterior beta = 0.773, adjusted P-value = 0.042). Significant positive correlations were found between CSTF1 and CPSF1 (r = 0.365, P < 0.001), WDR33 (r = 0.506, P < 0.001), and CPSF6 (r = 0.446, P < 0.001) expression levels.

CONCLUSION

Although further research is required to determine their potential contribution to AD, our findings offer a fresh perspective on molecular regulatory pathways associated with AD pathogenic mechanisms associated with PA and APA.

Alternative polyadenylation landscape of longissimus dorsi muscle with high and low intramuscular fat content in cattle

Intramuscular fat content is one of the most important factors affecting beef quality. However, the role of alternative polyadenylation (APA) in intramuscular fat deposition remains unclear. We compared APA events in muscle samples from high and low intramuscular fat (IMF) cattle, based on RNA-seq data. A total of 363 significant APAs were identified. Notably, the number of shortened 3'UTR events exceeded the number of lengthened 3'UTR events, and genes associated with shortened 3'UTR events were enriched in fatty acid metabolism-related pathways. Most APA events had alternative 3'UTR (aUTR) lengths of 200 to 300 bp. As the 3'UTR lengthened, the aUTR also lengthened (R2 = 0.79). These findings indicate that genes with longer 3'UTRs are more likely to be regulated by APA in the muscle of cattle with high IMF. To determine whether the identified APA events drove alterations in the expression of fat deposition-related genes, we analyzed the relationship between APA events and differentially expressed genes and identified several genes critical for fat deposition (e.g., PFKL and SLC1A5). Since miRNAs usually bind to the 3'UTR region of protein-coding genes and affect gene expression, we constructed an miRNA-APA network to detect several key miRNAs that may regulate fat deposition. We identified 10 important miRNAs that affect changes in IMF content, which may be gained (gained miRNA-binding sites) or lost (lost miRNA-binding sites) owing to 187 differential APA events. Our study characterized the APA profiles of cattle with high and low intramuscular fat content and provided further insights into the relationship between APA, miRNA, and fat deposition.

Regulation of alternative splicing and polyadenylation in neurons

Cell-type-specific gene expression is a fundamental feature of multicellular organisms and is achieved by combinations of regulatory strategies. Although cell-restricted transcription is perhaps the most widely studied mechanism, co-transcriptional and post-transcriptional processes are also central to the spatiotemporal control of gene functions. One general category of expression control involves the generation of multiple transcript isoforms from an individual gene, whose balance and cell specificity are frequently tightly regulated via diverse strategies. The nervous system makes particularly extensive use of cell-specific isoforms, specializing the neural function of genes that are expressed more broadly. Here, we review regulatory strategies and RNA-binding proteins that direct neural-specific isoform processing. These include various classes of alternative splicing and alternative polyadenylation events, both of which broadly diversify the neural transcriptome. Importantly, global alterations of splicing and alternative polyadenylation are characteristic of many neural pathologies, and recent genetic studies demonstrate how misregulation of individual neural isoforms can directly cause mutant phenotypes.

Mature mRNA processing that deletes 3' end sequences directs translational activation and embryonic development

Eggs accumulate thousands of translationally repressed mRNAs that are translated into proteins after fertilization to direct diverse developmental processes. However, molecular mechanisms underlying the translation of stored mRNAs after fertilization remain unclear. Here, we report a previously unknown RNA processing of 3' end sequences of mature mRNAs that activates the translation of stored mRNAs. Specifically, 9 to 72 nucleotides at the 3' ends of zebrafish pou5f3 and mouse Pou5f1 mRNAs were deleted in the early stages of development. Reporter assays illustrated the effective translation of the truncated forms of mRNAs. Moreover, promotion and inhibition of the shortening of 3' ends accelerated and attenuated Pou5f3 accumulation, respectively, resulting in defective development. Identification of proteins binding to unprocessed and/or processed mRNAs revealed that mRNA shortening acts as molecular switches. Comprehensive analysis revealed that >250 mRNAs underwent this processing. Therefore, our results provide a molecular principle that triggers the translational activation and directs development.

Extension of mRNA poly(A) tails and 3'UTRs during neuronal differentiation exhibits variable association with post-transcriptional dynamics

Differentiation of neural progenitor cells into mature neuronal phenotypes relies on extensive temporospatial coordination of mRNA expression to support the development of functional brain circuitry. Cleavage and polyadenylation of mRNA has tremendous regulatory capacity through the alteration of mRNA stability and modulation of microRNA (miRNA) function, however the extent of utilization in neuronal development is currently unclear. Here, we employed poly(A) tail sequencing, mRNA sequencing, ribosome profiling and small RNA sequencing to explore the functional relationship between mRNA abundance, translation, poly(A) tail length, alternative polyadenylation (APA) and miRNA expression in an in vitro model of neuronal differentiation. Differential analysis revealed a strong bias towards poly(A) tail and 3'UTR lengthening during differentiation, both of which were positively correlated with changes in mRNA abundance, but not translation. Globally, changes in miRNA expression were predominantly associated with mRNA abundance and translation, however several miRNA-mRNA pairings with potential to regulate poly(A) tail length were identified. Furthermore, 3'UTR lengthening was observed to significantly increase the inclusion of non-conserved miRNA binding sites, potentially enhancing the regulatory capacity of these molecules in mature neuronal cells. Together, our findings suggest poly(A) tail length and APA function as part of a rich post-transcriptional regulatory matrix during neuronal differentiation.

Haematopoietic stem and progenitor cell heterogeneity is inherited from the embryonic endothelium

Definitive haematopoietic stem and progenitor cells (HSPCs) generate erythroid, lymphoid and myeloid lineages. HSPCs are produced in the embryo via transdifferentiation of haemogenic endothelial cells in the aorta-gonad-mesonephros (AGM). HSPCs in the AGM are heterogeneous in differentiation and proliferative output, but how these intrinsic differences are acquired remains unanswered. Here we discovered that loss of microRNA (miR)-128 in zebrafish leads to an expansion of HSPCs in the AGM with different cell cycle states and a skew towards erythroid and lymphoid progenitors. Manipulating miR-128 in differentiating haemogenic endothelial cells, before their transition to HSPCs, recapitulated the lineage skewing in both zebrafish and human pluripotent stem cells. miR-128 promotes Wnt and Notch signalling in the AGM via post-transcriptional repression of the Wnt inhibitor csnk1a1 and the Notch ligand jag1b. De-repression of cskn1a1 resulted in replicative and erythroid-biased HSPCs, whereas de-repression of jag1b resulted in G2/M and lymphoid-biased HSPCs with long-term consequence on the respective blood lineages. We propose that HSPC heterogeneity arises in the AGM endothelium and is programmed in part by Wnt and Notch signalling.

Sites of transcription initiation drive mRNA isoform selection

The generation of distinct messenger RNA isoforms through alternative RNA processing modulates the expression and function of genes, often in a cell-type-specific manner. Here, we assess the regulatory relationships between transcription initiation, alternative splicing, and 3' end site selection. Applying long-read sequencing to accurately represent even the longest transcripts from end to end, we quantify mRNA isoforms in Drosophila tissues, including the transcriptionally complex nervous system. We find that in Drosophila heads, as well as in human cerebral organoids, 3' end site choice is globally influenced by the site of transcription initiation (TSS). "Dominant promoters," characterized by specific epigenetic signatures including p300/CBP binding, impose a transcriptional constraint to define splice and polyadenylation variants. In vivo deletion or overexpression of dominant promoters as well as p300/CBP loss disrupted the 3' end expression landscape. Our study demonstrates the crucial impact of TSS choice on the regulation of transcript diversity and tissue identity.

Intron dynamics reveal principles of gene regulation during the maternal-to-zygotic transition

The maternal-to-zygotic transition (MZT) is a conserved embryonic process in animals where developmental control shifts from the maternal to zygotic genome. A key step in this transition is zygotic transcription, and deciphering the MZT requires classifying newly transcribed genes. However, due to current technological limitations, this starting point remains a challenge for studying many species. Here, we present an alternative approach that characterizes transcriptome changes based solely on RNA-seq data. By combining intron-mapping reads and transcript-level quantification, we characterized transcriptome dynamics during the Drosophila melanogaster MZT. Our approach provides an accessible platform to investigate transcriptome dynamics that can be applied to the MZT in nonmodel organisms. In addition to classifying zygotically transcribed genes, our analysis revealed that over 300 genes express different maternal and zygotic transcript isoforms due to alternative splicing, polyadenylation, and promoter usage. The vast majority of these zygotic isoforms have the potential to be subject to different regulatory control, and over two-thirds encode different proteins. Thus, our analysis reveals an additional layer of regulation during the MZT, where new zygotic transcripts can generate additional proteome diversity.

Alternative 3'-untranslated regions regulate high-salt tolerance of Spartina alterniflora

High-salt stress continues to challenge the growth and survival of many plants. Alternative polyadenylation (APA) produces mRNAs with different 3'-untranslated regions (3' UTRs) to regulate gene expression at the post-transcriptional level. However, the roles of alternative 3' UTRs in response to salt stress remain elusive. Here, we report the function of alternative 3' UTRs in response to high-salt stress in S. alterniflora (Spartina alterniflora), a monocotyledonous halophyte tolerant of high-salt environments. We found that high-salt stress induced global APA dynamics, and ∼42% of APA genes responded to salt stress. High-salt stress led to 3' UTR lengthening of 207 transcripts through increasing the usage of distal poly(A) sites. Transcripts with alternative 3' UTRs were mainly enriched in salt stress-related ion transporters. Alternative 3' UTRs of HIGH-AFFINITY K+ TRANSPORTER 1 (SaHKT1) increased RNA stability and protein synthesis in vivo. Regulatory AU-rich elements were identified in alternative 3' UTRs, boosting the protein level of SaHKT1. RNAi-knock-down experiments revealed that the biogenesis of 3' UTR lengthening in SaHKT1 was controlled by the poly(A) factor CLEAVAGE AND POLYADENYLATION SPECIFICITY FACTOR 30 (SaCPSF30). Over-expression of SaHKT1 with an alternative 3' UTR in rice (Oryza sativa) protoplasts increased mRNA accumulation of salt-tolerance genes in an AU-rich element-dependent manner. These results suggest that mRNA 3' UTR lengthening is a potential mechanism in response to high-salt stress. These results also reveal complex regulatory roles of alternative 3' UTRs coupling APA and regulatory elements at the post-transcriptional level in plants.

A Global Survey of the Full-Length Transcriptome of Apis mellifera by Single-Molecule Long-Read Sequencing

As important pollinators, honey bees play a crucial role in both maintaining the ecological balance and providing products for humans. Although several versions of the western honey bee genome have already been published, its transcriptome information still needs to be refined. In this study, PacBio single-molecule sequencing technology was used to sequence the full-length transcriptome of mixed samples from many developmental time points and tissues of A. mellifera queens, workers and drones. A total of 116,535 transcripts corresponding to 30,045 genes were obtained. Of these, 92,477 transcripts were annotated. Compared to the annotated genes and transcripts on the reference genome, 18,915 gene loci and 96,176 transcripts were newly identified. From these transcripts, 136,554 alternative splicing (AS) events, 23,376 alternative polyadenylation (APA) sites and 21,813 lncRNAs were detected. In addition, based on the full-length transcripts, we identified many differentially expressed transcripts (DETs) between queen, worker and drone. Our results provide a complete set of reference transcripts for A. mellifera that dramatically expand our understanding of the complexity and diversity of the honey bee transcriptome.

Regulation of neuronal RNA signatures by ELAV/Hu proteins

The RNA-binding proteins encoded by the highly conserved elav/Hu gene family, found in all metazoans, regulate the expression of a wide range of genes, at both the co-transcriptional and posttranscriptional level. Nervous-system-specific ELAV/Hu proteins are prominent for their essential role in neuron differentiation, and mutations have been associated with human neurodevelopmental and neurodegenerative diseases. Drosophila ELAV, the founding member of the protein family, mediates the synthesis of neuronal RNA signatures by promoting alternative splicing and alternative polyadenylation of hundreds of genes. The recent identification of ELAV's direct RNA targets revealed the protein's central role in shaping the neuronal transcriptome, and highlighted the importance of neuronal transcript signatures for neuron maintenance and organism survival. Animals have evolved multiple cellular mechanisms to ensure robustness of ELAV/Hu function. In Drosophila, elav autoregulates in a 3'UTR-dependent manner to maintain optimal protein levels. A complete absence of ELAV causes the activation and nuclear localization of the normally cytoplasmic paralogue FNE, in a process termed EXon-Activated functional Rescue (EXAR). Other species, including mammals, seem to utilize different strategies, such as protein redundancy, to maintain ELAV protein function and effectively safeguard the identity of the neuronal transcriptome. This article is categorized under: RNA Processing > 3' End Processing RNA in Disease and Development > RNA in Development RNA Interactions with Proteins and Other Molecules > Protein-RNA Interactions: Functional Implications.

Transcriptional buffering and 3'UTR lengthening are shaped during human neurodevelopment by shifts in mRNA stability and microRNA load

The contribution of mRNA half-life is commonly overlooked when examining changes in mRNA abundance during development. mRNA levels of some genes are regulated by transcription rate only, but others may be regulated by mRNA half-life only shifts. Furthermore, transcriptional buffering is predicted when changes in transcription rates have compensating shifts in mRNA half-life resulting in no change to steady-state levels. Likewise, transcriptional boosting should result when changes in transcription rate are accompanied by amplifying half-life shifts. During neurodevelopment there is widespread 3'UTR lengthening that could be shaped by differential shifts in the stability of existing short or long 3'UTR transcript isoforms. We measured transcription rate and mRNA half-life changes during induced human Pluripotent Stem Cell (iPSC)-derived neuronal development using RATE-seq. During transitions to progenitor and neuron stages, transcriptional buffering occurred in up to 50%, and transcriptional boosting in up to 15%, of genes with changed transcription rates. The remaining changes occurred by transcription rate only or mRNA half-life only shifts. Average mRNA half-life decreased two-fold in neurons relative to iPSCs. Short gene isoforms were more destabilized in neurons and thereby increased the average 3'UTR length. Small RNA sequencing captured an increase in microRNA copy number per cell during neurodevelopment. We propose that mRNA destabilization and 3'UTR lengthening are driven in part by an increase in microRNA load in neurons. Our findings identify mRNA stability mechanisms in human neurodevelopment that regulate gene and isoform level abundance and provide a precedent for similar post-transcriptional regulatory events as other tissues develop.

SLAMseq resolves the kinetics of maternal and zygotic gene expression during early zebrafish embryogenesis

The maternal-to-zygotic transition (MZT) is a key developmental process in metazoan embryos that involves the activation of zygotic transcription (ZGA) and degradation of maternal transcripts. We employed metabolic mRNA sequencing (SLAMseq) to deconvolute the compound embryonic transcriptome in zebrafish. While mitochondrial zygotic transcripts prevail prior to MZT, we uncover the spurious transcription of hundreds of short and intron-poor genes as early as the 2-cell stage. Upon ZGA, most zygotic transcripts originate from thousands of maternal-zygotic (MZ) genes that are transcribed at rates comparable to those of hundreds of purely zygotic genes and replenish maternal mRNAs at distinct timescales. Rapid replacement of MZ transcripts involves transcript decay features unrelated to major maternal degradation pathways and promotes de novo synthesis of the core gene expression machinery by increasing poly(A)-tail length and translation efficiency. SLAMseq hence provides insights into the timescales, molecular features, and regulation of MZT during zebrafish embryogenesis.

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.

Regulation and function of alternative polyadenylation in development and differentiation

Alternative processing of nascent mRNAs is widespread in eukaryotic organisms and greatly impacts the output of gene expression. Specifically, alternative cleavage and polyadenylation (APA) is a co-transcriptional molecular process that switches the polyadenylation site (PAS) at which a nascent mRNA is cleaved, resulting in mRNA isoforms with different 3'UTR length and content. APA can potentially affect mRNA translation efficiency, localization, stability, and mRNA seeded protein-protein interactions. APA naturally occurs during development and cellular differentiation, with around 70% of human genes displaying APA in particular tissues and cell types. For example, neurons tend to express mRNAs with long 3'UTRs due to preferential processing at PASs more distal than other PASs used in other cell types. In addition, changes in APA mark a variety of pathological states, including many types of cancer, in which mRNAs are preferentially cleaved at more proximal PASs, causing expression of mRNA isoforms with short 3'UTRs. Although APA has been widely reported, both the function of APA in development and the mechanisms that regulate the choice of 3'end cut sites in normal and pathogenic conditions are still poorly understood. In this review, we summarize current understanding of how APA is regulated during development and cellular differentiation and how the resulting change in 3'UTR content affects multiple aspects of gene expression. With APA being a widespread phenomenon, the advent of cutting-edge scientific techniques and the pressing need for in-vivo studies, there has never been a better time to delve into the intricate mechanisms of alternative cleavage and polyadenylation.

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.

Expression Profile of Sphingosine Kinase 1 Isoforms in Human Cancer Tissues and Cells: Importance and Clinical Relevance of the Neglected 1b-Isoform

Background

Overexpression of sphingosine kinase 1 (SphK1) is casually associated with many types of cancer, and inhibitors of SphK1 sensitize tumors to chemotherapy. SphK1 is expressed as two major isoforms, SphK1a and SphK1b. To date, no information has been reported on the SphK1 isoform expression profile and its clinical relevance.

Objective

The objective is to examine the expression profile of the SphK1a and SPhK1b isoforms in human cancer and noncancer tissues and cell lines and explore their clinical relevance.

Methods

We used PCR to qualitatively examine the expression profile of these two isoforms in breast, liver, and prostate cancer tissues plus paired adjacent tissues and in 11 cancer and normal cell lines (breast, cervical, bone, prostate, colon, brain, mesothelioma tumor and benign, and human kidney cells).

Results

We found that SphK1a was ubiquitously expressed in all cancer cells and tissues tested; in contrast, SphK1b was only expressed in selective cell types in breast, prostate, and lung cancer.

Conclusions

Our data suggest that SphK1a is important for generic SphK1/S1P functions, and SphK1b mediates specialized and/or unique pathways in a specific type of tissue and could be a biomarker for cancer. This discovery is important for future SphK1-related cancer research and may have clinical implications in drug development associated with SphK1-directed cancer treatment.

Gene product diversity: adaptive or not?

One gene does not equal one RNA or protein. The genomic revolution has revealed numerous different RNA and protein molecules that can be produced from one gene, such as circular RNAs generated by back-splicing, proteins with residues mismatching the genomic encoding because of RNA editing, and proteins extended in the C terminus via stop codon readthrough in translation. Are these diverse products results of exquisite gene regulations or imprecise biological processes? While there are cases where the gene product diversity appears beneficial, genome-scale patterns suggest that much of this diversity arises from nonadaptive, molecular errors. This finding has important implications for studying the functions of diverse gene products and for understanding the fundamental properties and evolution of cellular life.

A critical developmental window for ELAV/Hu-dependent mRNA signatures at the onset of neuronal differentiation

Cell-type-specific gene regulatory programs are essential for cell differentiation and function. In animal neurons, the highly conserved ELAV/Hu family of proteins promotes alternative splicing and polyadenylation of mRNA precursors to create unique neuronal transcript isoforms. Here, we assess transcriptome profiles and neurogenesis success in Drosophila models engineered to express differing levels of ELAV activity in the course of development. We show that the ELAV-mediated establishment of a subset of neuronal mRNA isoforms at the onset of neuron differentiation constitutes a developmental bottleneck that cannot be overcome later by the nuclear activation of the paralog found in neurons (FNE). Loss of ELAV function outside of that critical time window results in neurological defects. We find that FNE, when activated early enough, can restore ELAV-dependent neuronal mRNA isoforms and fully rescue development. Our findings demonstrate the essential role of robust cellular strategies to maintain ELAV activity and intact neuronal signatures in neurogenesis and neuronal function.

A Survey on Methods for Predicting Polyadenylation Sites from DNA Sequences, Bulk RNA-seq, and Single-cell RNA-seq

Alternative polyadenylation (APA) plays important roles in modulating mRNA stability, translation, and subcellular localization, and contributes extensively to shaping eukaryotic transcriptome complexity and proteome diversity. Identification of poly(A) sites (pAs) on a genome-wide scale is a critical step toward understanding the underlying mechanism of APA-mediated gene regulation. A number of established computational tools have been proposed to predict pAs from diverse genomic data. Here we provided an exhaustive overview of computational approaches for predicting pAs from DNA sequences, bulk RNA sequencing (RNA-seq) data, and single-cell RNA sequencing (scRNA-seq) data. Particularly, we examined several representative tools using bulk RNA-seq and scRNA-seq data from peripheral blood mononuclear cells and put forward operable suggestions on how to assess the reliability of pAs predicted by different tools. We also proposed practical guidelines on choosing appropriate methods applicable to diverse scenarios. Moreover, we discussed in depth the challenges in improving the performance of pA prediction and benchmarking different methods. Additionally, we highlighted outstanding challenges and opportunities using new machine learning and integrative multi-omics techniques, and provided our perspective on how computational methodologies might evolve in the future for non-3' untranslated region, tissue-specific, cross-species, and single-cell pA prediction.

Diverse cell-specific patterns of alternative polyadenylation in Drosophila

Most genes in higher eukaryotes express isoforms with distinct 3' untranslated regions (3' UTRs), generated by alternative polyadenylation (APA). Since 3' UTRs are predominant locations of post-transcriptional regulation, APA can render such programs conditional, and can also alter protein sequences via alternative last exon (ALE) isoforms. We previously used 3'-sequencing from diverse Drosophila samples to define multiple tissue-specific APA landscapes. Here, we exploit comprehensive single nucleus RNA-sequencing data (Fly Cell Atlas) to elucidate cell-type expression of 3' UTRs across >250 adult Drosophila cell types. We reveal the cellular bases of multiple tissue-specific APA/ALE programs, such as 3' UTR lengthening in differentiated neurons and 3' UTR shortening in spermatocytes and spermatids. We trace dynamic 3' UTR patterns across cell lineages, including in the male germline, and discover new APA patterns in the intestinal stem cell lineage. Finally, we correlate expression of RNA binding proteins (RBPs), miRNAs and global levels of cleavage and polyadenylation (CPA) factors in several cell types that exhibit characteristic APA landscapes, yielding candidate regulators of transcriptome complexity. These analyses provide a comprehensive foundation for future investigations of mechanisms and biological impacts of alternative 3' isoforms across the major cell types of this widely-studied model organism.

MiR-202-3p determines embryo viability during mid-blastula transition

Developmental growth is an intricate process involving the coordinated regulation of the expression of various genes, and microRNAs (miRNAs) play crucial roles in diverse processes throughout animal development. The mid-blastula transition (MBT) is a developmental milestone when maternal RNAs are cleared and the zygotic genome programmed asynchronous cell division begins to drive embryogenesis. While mechanisms underlying MBT have been intensively revealed, factors regulating cell proliferation at the transition remain largely unknown. We report here a microRNA, miR-202-3p to be a key factor that determines embryonic fate during MBT in zebrafish. A miR-202-3p antagomir specifically terminated embryo development at the mid-blastula stage. In vivo deletion of the miR-202 locus recapitulated the fatal phenotypes, which were rescued only by miR-202-3p or its precursor. Transcriptome comparison revealed >250 RNAs including both maternal and zygotic origins were dysregulated at MBT in the miR-202^−/− embryos, corresponding with arrays of homeostatic disorders leading to massive apoptosis. A trio of genes: nfkbiaa, perp and mgll, known to be intimately involved with cell proliferation and survival, were identified as direct targets of miR-202-3p. Importantly, over- or under-expression of any of the trio led to developmental delay or termination at the blastula or gastrula stages. Furthermore, nfkbiaa and perp were shown to inter-regulate each other. Thus, miR-202-3p mediates a regulatory network whose components interact closely during MBT to determine embryonic viability and development.

Developmentally regulated alternate 3' end cleavage of nascent transcripts controls dynamic changes in protein expression in an adult stem cell lineage

Alternative polyadenylation (APA) generates transcript isoforms that differ in the position of the 3' cleavage site, resulting in the production of mRNA isoforms with different length 3' UTRs. Although widespread, the role of APA in the biology of cells, tissues, and organisms has been controversial. We identified >500 Drosophila genes that express mRNA isoforms with a long 3' UTR in proliferating spermatogonia but a short 3' UTR in differentiating spermatocytes due to APA. We show that the stage-specific choice of the 3' end cleavage site can be regulated by the arrangement of a canonical polyadenylation signal (PAS) near the distal cleavage site but a variant or no recognizable PAS near the proximal cleavage site. The emergence of transcripts with shorter 3' UTRs in differentiating cells correlated with changes in expression of the encoded proteins, either from off in spermatogonia to on in spermatocytes or vice versa. Polysome gradient fractionation revealed >250 genes where the long 3' UTR versus short 3' UTR mRNA isoforms migrated differently, consistent with dramatic stage-specific changes in translation state. Thus, the developmentally regulated choice of an alternative site at which to make the 3' end cut that terminates nascent transcripts can profoundly affect the suite of proteins expressed as cells advance through sequential steps in a differentiation lineage.

Temporal Dynamic Analysis of Alternative Splicing During Embryonic Development in Zebrafish

Alternative splicing is pervasive in mammalian genomes and involved in embryo development, whereas research on crosstalk of alternative splicing and embryo development was largely restricted to mouse and human and the alternative splicing regulation during embryogenesis in zebrafish remained unclear. We constructed the alternative splicing atlas at 18 time-course stages covering maternal-to-zygotic transition, gastrulation, somitogenesis, pharyngula stages, and post-fertilization in zebrafish. The differential alternative splicing events between different developmental stages were detected. The results indicated that abundance alternative splicing and differential alternative splicing events are dynamically changed and remarkably abundant during the maternal-to-zygotic transition process. Based on gene expression profiles, we found splicing factors are expressed with specificity of developmental stage and largely expressed during the maternal-to-zygotic transition process. The better performance of cluster analysis was achieved based on the inclusion level of alternative splicing. The biological function analysis uncovered the important roles of alternative splicing during embryogenesis. The identification of isoform switches of alternative splicing provided a new insight into mining the regulated mechanism of transcript isoforms, which always is hidden by gene expression. In conclusion, we inferred that alternative splicing activation is synchronized with zygotic genome activation and discovered that alternative splicing is coupled with transcription during embryo development in zebrafish. We also unveiled that the temporal expression dynamics of splicing factors during embryo development, especially co-orthologous splicing factors. Furthermore, we proposed that the inclusion level of alternative splicing events can be employed for cluster analysis as a novel parameter. This work will provide a deeper insight into the regulation of alternative splicing during embryogenesis in zebrafish.

Dysregulated RNA polyadenylation contributes to metabolic impairment in non-alcoholic fatty liver disease

Pre-mRNA processing is an essential mechanism for the generation of mature mRNA and the regulation of gene expression in eukaryotic cells. While defects in pre-mRNA processing have been implicated in a number of diseases their involvement in metabolic pathologies is still unclear. Here, we show that both alternative splicing and alternative polyadenylation, two major steps in pre-mRNA processing, are significantly altered in non-alcoholic fatty liver disease (NAFLD). Moreover, we find that Serine and Arginine Rich Splicing Factor 10 (SRSF10) binding is enriched adjacent to consensus polyadenylation motifs and its expression is significantly decreased in NAFLD, suggesting a role mediating pre-mRNA dysregulation in this condition. Consistently, inactivation of SRSF10 in mouse and human hepatocytes in vitro, and in mouse liver in vivo, was found to dysregulate polyadenylation of key metabolic genes such as peroxisome proliferator-activated receptor alpha (PPARA) and exacerbate diet-induced metabolic dysfunction. Collectively our work implicates dysregulated pre-mRNA polyadenylation in obesity-induced liver disease and uncovers a novel role for SRSF10 in this process.

An Evolutionarily Conserved AU-Rich Element in the 3' Untranslated Region of a Transcript Misannotated as a Long Noncoding RNA Regulates RNA Stability

One of the primary mechanisms of post-transcriptional gene regulation is the modulation of RNA stability. We recently discovered that LINC00675, a transcript annotated as a long noncoding RNA (lncRNA), is transcriptionally regulated by FOXA1 and encodes a highly conserved small protein that localizes to the endoplasmic reticulum, hence renamed as FORCP (FOXA1-regulated conserved small protein). Here, we show that the endogenous FORCP transcript is rapidly degraded and rendered unstable as a result of 3'UTR-mediated degradation. Surprisingly, although the FORCP transcript is a canonical nonsense-mediated decay (NMD) and microRNA (miRNA) target, we found that it is not degraded by NMD or miRNAs. Targeted deletion of an evolutionarily conserved region in the FORCP 3'UTR using CRISPR/Cas9 significantly increased the stability of the FORCP transcript. Interestingly, this region requires the presence of an immediate downstream 55-nt-long sequence for transcript stability regulation. Functionally, colorectal cancer cells lacking this conserved region expressed from the endogenous FORCP locus displayed decreased proliferation and clonogenicity. These data demonstrate that the FORCP transcript is destabilized via conserved elements within its 3'UTR and emphasize the need to interrogate the function of a given 3'UTR in its native context.

Regulation of the Alternative Neural Transcriptome by ELAV/Hu RNA Binding Proteins

The process of alternative polyadenylation (APA) generates multiple 3' UTR isoforms for a given locus, which can alter regulatory capacity and on occasion change coding potential. APA was initially characterized for a few genes, but in the past decade, has been found to be the rule for metazoan genes. While numerous differences in APA profiles have been catalogued across genetic conditions, perturbations, and diseases, our knowledge of APA mechanisms and biology is far from complete. In this review, we highlight recent findings regarding the role of the conserved ELAV/Hu family of RNA binding proteins (RBPs) in generating the broad landscape of lengthened 3' UTRs that is characteristic of neurons. We relate this to their established roles in alternative splicing, and summarize ongoing directions that will further elucidate the molecular strategies for neural APA, the in vivo functions of ELAV/Hu RBPs, and the phenotypic consequences of these regulatory paradigms in neurons.

Amino Acids and IGF1 Regulation of Fish Muscle Growth Revealed by Transcriptome and microRNAome Integrative Analyses of Pacu (Piaractus mesopotamicus) Myotubes

Amino acids (AA) and IGF1 have been demonstrated to play essential roles in protein synthesis and fish muscle growth. The myoblast cell culture is useful for studying muscle regulation, and omics data have contributed enormously to understanding its molecular biology. However, to our knowledge, no study has performed the large-scale sequencing of fish-cultured muscle cells stimulated with pro-growth signals. In this work, we obtained the transcriptome and microRNAome of pacu (Piaractus mesopotamicus)-cultured myotubes treated with AA or IGF1. We identified 1228 and 534 genes differentially expressed by AA and IGF1. An enrichment analysis showed that AA treatment induced chromosomal changes, mitosis, and muscle differentiation, while IGF1 modulated IGF/PI3K signaling, metabolic alteration, and matrix structure. In addition, potential molecular markers were similarly modulated by both treatments. Muscle-miRNAs (miR-1, -133, -206 and -499) were up-regulated, especially in AA samples, and we identified molecular networks with omics integration. Two pairs of genes and miRNAs demonstrated a high-level relationship, and involvement in myogenesis and muscle growth: marcksb and miR-29b in AA, and mmp14b and miR-338-5p in IGF1. Our work helps to elucidate fish muscle physiology and metabolism, highlights potential molecular markers, and creates a perspective for improvements in aquaculture and in in vitro meat production.

Alternative polyadenylation dysregulation contributes to the differentiation block of acute myeloid leukemia

Posttranscriptional regulation has emerged as a driver for leukemia development and an avenue for therapeutic targeting. Among posttranscriptional processes, alternative polyadenylation (APA) is globally dysregulated across cancer types. However, limited studies have focused on the prevalence and role of APA in myeloid leukemia. Furthermore, it is poorly understood how altered poly(A) site usage of individual genes contributes to malignancy or whether targeting global APA patterns might alter oncogenic potential. In this study, we examined global APA dysregulation in patients with acute myeloid leukemia (AML) by performing 3' region extraction and deep sequencing (3'READS) on a subset of AML patient samples along with healthy hematopoietic stem and progenitor cells (HSPCs) and by analyzing publicly available data from a broad AML patient cohort. We show that patient cells exhibit global 3' untranslated region (UTR) shortening and coding sequence lengthening due to differences in poly(A) site (PAS) usage. Among APA regulators, expression of FIP1L1, one of the core cleavage and polyadenylation factors, correlated with the degree of APA dysregulation in our 3'READS data set. Targeting global APA by FIP1L1 knockdown reversed the global trends seen in patients. Importantly, FIP1L1 knockdown induced differentiation of t(8;21) cells by promoting 3'UTR lengthening and downregulation of the fusion oncoprotein AML1-ETO. In non-t(8;21) cells, FIP1L1 knockdown also promoted differentiation by attenuating mechanistic target of rapamycin complex 1 (mTORC1) signaling and reducing MYC protein levels. Our study provides mechanistic insights into the role of APA in AML pathogenesis and indicates that targeting global APA patterns can overcome the differentiation block in patients with AML.

CRISPR-Mediated Knockout of Long 3' UTR mRNA Isoforms in mESC-Derived Neurons

Alternative cleavage and polyadenylation (APA) is pervasive, occurring for more than 70% of human and mouse genes. Distal poly(A) site selection to generate longer 3′ UTR mRNA isoforms is prevalent in the nervous system, affecting thousands of genes. Here, we establish mouse embryonic stem cell (mESC)-derived neurons (mES-neurons) as a suitable system to study long 3′ UTR isoforms. RNA-seq analysis revealed that mES-neurons show widespread 3′ UTR lengthening that closely resembles APA patterns found in mouse cortex. mESCs are highly amenable to genetic manipulation. We present a method to eliminate long 3′ UTR isoform expression using CRISPR/Cas9 editing. This approach can lead to clones with the desired deletion within several weeks. We demonstrate this strategy on the Mprip gene as a proof-of-principle. To confirm loss of long 3′ UTR expression and the absence of cryptic poly(A) site usage stemming from the CRISPR deletion, we present a simple and cost-efficient targeted long-read RNA-sequencing strategy using the Oxford Nanopore Technologies platform. Using this method, we confirmed specific loss of the Mprip long 3′ UTR isoform. CRISPR gene editing of mESCs thus serves as a highly relevant platform for studying the molecular and cellular functions of long 3′ UTR mRNA isoforms.

Computational Detection of MicroRNA Targets

MicroRNAs (miRNAs) are small noncoding RNAs that are recognized as posttranscriptional regulators of gene expression. These molecules have been shown to play important roles in several cellular processes. MiRNAs act on their target by guiding the RISC complex and binding to the mRNA molecule. Thus, it is recognized that the function of a miRNA is determined by the function of its target (s). By using high-throughput methodologies, novel miRNAs are being identified, but their functions remain uncharted. Target validation is crucial to properly understand the specific role of a miRNA in a cellular pathway. However, molecular techniques for experimental validation of miRNA-target interaction are expensive, time-consuming, laborious, and can be not accurate in inferring true interactions. Thus, accurate miRNA target predictions are helpful to understand the functions of miRNAs. There are several algorithms proposed for target prediction and databases containing miRNA-target information. However, these available computational tools for prediction still generate a large number of false positives and fail to detect a considerable number of true targets, which indicates the necessity of highly confident approaches to identify bona fide miRNA-target interactions. This chapter focuses on tools and strategies used for miRNA target prediction, by providing practical insights and outlooks.

Sex-biased dynamics of three-spined stickleback (Gasterosteus aculeatus) gene expression patterns

The study of the differences between sexes presents an excellent model to unravel how phenotypic variation is achieved from a similar genetic background. Sticklebacks are of particular interest since evidence of a heteromorphic chromosome pair has not always been detected. The present study investigated sex-biased mRNA and small non-coding RNA (sncRNA) expression patterns in the brain, adipose tissues, and gonads of the three-spined stickleback. The sncRNA analysis indicated that regulatory functions occurred mainly in the gonads. Alleged miRNA-mRNA interactions were established and a mapping bias of differential expressed transcripts towards chromosome 19 was observed. Key players previously shown to control sex determination and differentiation in other fish species but also genes like gapdh were among the transcripts identified. This is the first report in the three-spined stickleback demonstrating tissue-specific expression comprising both mRNA and sncRNA between sexes, emphasizing the importance of mRNA-miRNA interactions as well as new presumed genes not yet identified to have gender-specific roles.

Imaging Organization of RNA Processing within the Nucleus

Within the nucleus, messenger RNA is generated and processed in a highly organized and regulated manner. Messenger RNA processing begins during transcription initiation and continues until the RNA is translated and degraded. Processes such as 5' capping, alternative splicing, and 3' end processing have been studied extensively with biochemical methods and more recently with single-molecule imaging approaches. In this review, we highlight how imaging has helped understand the highly dynamic process of RNA processing. We conclude with open questions and new technological developments that may further our understanding of RNA processing.

CstF64-Induced Shortening of the BID 3'UTR Promotes Esophageal Squamous Cell Carcinoma Progression by Disrupting ceRNA Cross-talk with ZFP36L2

The majority of human genes have multiple polyadenylation sites, which are differentially used through the process of alternative polyadenylation (APA). Dysregulation of APA contributes to numerous diseases, including cancer. However, specific genes subject to APA that impact oncogenesis have not been well characterized, and many cancer APA landscapes remain underexplored. Here, we used dynamic analyses of APA from RNA-seq (DaPars) to define both the 3'UTR APA profile in esophageal squamous cell carcinoma (ESCC) and to identify 3'UTR shortening events that may drive tumor progression. In four distinct squamous cell carcinoma datasets, BID 3'UTRs were recurrently shortened and BID mRNA levels were significantly upregulated. Moreover, system correlation analysis revealed that CstF64 is a candidate upstream regulator of BID 3'UTR length. Mechanistically, a shortened BID 3'UTR promoted proliferation of ESCC cells by disrupting competing endogenous RNA (ceRNA) cross-talk, resulting in downregulation of the tumor suppressor gene ZFP36L2. These in vitro and in vivo results were supported by human patient data whereby 3'UTR shortening of BID and low expression of ZFP36L2 are prognostic factors of survival in ESCC. Collectively, these findings demonstrate that a key ceRNA network is disrupted through APA and promotes ESCC tumor progression.Significance: High-throughput analysis of alternative polyadenylation in esophageal squamous cell carcinoma identifies recurrent shortening of the BID 3'UTR as a driver of disease progression.

IL-6 Signaling Protects Zebrafish Larvae during Staphylococcus epidermidis Infection in a Bath Immersion Model

The host immune responses to Staphylococcus epidermidis, a frequent cause of nosocomial infections, are not well understood. We have established a bath immersion model of this infection in zebrafish (Danio rerio) larvae. Macrophages play a primary role in the host immune response and are involved in clearance of infection in the larvae. S. epidermidis infection results in upregulation of tlr-2 There is marked inflammation characterized by heightened NF-κB signaling and elevation of several proinflammatory cytokines. There is rapid upregulation of il-1b and tnf-a transcripts, whereas an increase in il-6 levels is relatively more delayed. The IL-6 signaling pathway is further amplified by elevation of IL-6 signal transducer (il-6st) levels, which negatively correlates with miRNA dre-miR-142a-5p. Enhanced IL-6 signaling is protective to the host in this model as inhibition of the signaling pathway resulted in increased mortality upon S. epidermidis infection. Our study describes the host immune responses to S. epidermidis infection, establishes the importance of IL-6 signaling, and identifies a potential role of miR-142-5p-il-6st interaction in this infection model.

Alternative polyadenylation: An untapped source for prostate cancer biomarkers and therapeutic targets?

OBJECTIVE

To review alternative polyadenylation (APA) as a mechanism of gene regulation and consider potential roles for APA in prostate cancer (PCa) biology and treatment.

METHODS

An extensive review of mRNA polyadenylation, APA, and PCa literature was performed. This review article introduces APA and its association with human disease, outlines the mechanisms and components of APA, reviews APA in cancer biology, and considers whether APA may contribute to PCa progression and/or produce novel biomarkers and therapeutic targets for PCa.

RESULTS

Eukaryotic mRNA 3'-end cleavage and polyadenylation play a critical role in gene expression. Most human genes encode more than one polyadenylation signal, and produce more than one transcript isoform, through APA. Polyadenylation can occur throughout the gene body to generate transcripts with differing 3'-termini and coding sequence. Differences in 3'-untranslated regions length can modify post-transcriptional gene regulation by microRNAs and RNA binding proteins, and alter mRNA stability, translation efficiency, and subcellular localization. Distinctive APA patterns are associated with human diseases, tissue origins, and changes in cellular proliferation rate and differentiation state. APA events may therefore generate unique mRNA biomarkers or therapeutic targets in certain cancer types or phenotypic states.

CONCLUSIONS

The full extent of cancer-associated and tissue-specific APA events have yet to be defined, and the mechanisms and functional consequences of APA in cancer remain incompletely understood. There is evidence that APA is active in PCa, and that it may be an untapped resource for PCa biomarkers or therapeutic targets.