Extensive alternative polyadenylation during zebrafish development

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.

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.

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.

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.

The landscape of alternative polyadenylation in single cells of the developing mouse embryo

3′ untranslated regions (3′ UTRs) post-transcriptionally regulate mRNA stability, localization, and translation rate. While 3′-UTR isoforms have been globally quantified in limited cell types using bulk measurements, their differential usage among cell types during mammalian development remains poorly characterized. In this study, we examine a dataset comprising ~2 million nuclei spanning E9.5–E13.5 of mouse embryonic development to quantify transcriptome-wide changes in alternative polyadenylation (APA). We observe a global lengthening of 3′ UTRs across embryonic stages in all cell types, although we detect shorter 3′ UTRs in hematopoietic lineages and longer 3′ UTRs in neuronal cell types within each stage. An analysis of RNA-binding protein (RBP) dynamics identifies ELAV-like family members, which are concomitantly induced in neuronal lineages and developmental stages experiencing 3′-UTR lengthening, as putative regulators of APA. By measuring 3′-UTR isoforms in an expansive single cell dataset, our work provides a transcriptome-wide and organism-wide map of the dynamic landscape of alternative polyadenylation during mammalian organogenesis.

Alternative polyadenylation regulates localization, half-life and translation of mRNA isoforms. Here the authors investigate alternative polyadenylation using single cell RNA sequencing data from mouse embryos and identify 3’-UTR isoforms that are regulated across cell types and developmental time.

Transcriptomic analysis to elucidate the effects of high stocking density on grass carp (Ctenopharyngodon idella)

Background

Grass carp (Ctenopharyngodon idella) is one of the most widely cultivated fishes in China. High stocking density can reportedly affect fish growth and immunity. Herein we performed PacBio long-read single-molecule real-time (SMRT) sequencing and Illumina RNA sequencing to evaluate the effects of high stocking density on grass carp transcriptome.

Results

SMRT sequencing led to the identification of 33,773 genes (14,946 known and 18,827 new genes). From the structure analysis, 8,009 genes were detected with alternative splicing events, 10,219 genes showed alternative polyadenylation sites and 15,521 long noncoding RNAs. Further, 1,235, 962, and 213 differentially expressed genes (DEGs) were identified in the intestine, muscle, and brain tissues, respectively. We performed functional enrichment analyses of DEGs, and they were identified to be significantly enriched in nutrient metabolism and immune function. The expression levels of several genes encoding apolipoproteins and activities of enzymes involved in carbohydrate enzymolysis were found to be upregulated in the high stocking density group, indicating that lipid metabolism and carbohydrate decomposition were accelerated. Besides, four isoforms of grass carp major histocompatibility complex class II antigen alpha and beta chains in the aforementioned three tissue was showed at least a 4-fold decrease.

Conclusions

The results suggesting that fish farmed at high stocking densities face issues associated with the metabolism and immune system. To conclude, our results emphasize the importance of maintaining reasonable density in grass carp aquaculture.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07924-4.

The brominated flame retardant, TBCO, impairs oocyte maturation in zebrafish (Danio rerio)

The brominated flame retardant, 1,2,5,6-tetrabromocyclooctane (TBCO), has been shown to decrease fecundity in Japanese medaka (Oryzias latipes) and there is indirect evidence from analysis of the transcriptome and proteome that this effect might be due to impaired oogenesis. An assay for disruption of oocyte maturation by chemical stressors has not been developed in Japanese medaka. Thus, using zebrafish (Danio rerio) as a model, objectives of the present study were to determine whether exposure to TBCO has effects on maturation of oocytes and to investigate potential mechanisms. Sexually mature female zebrafish were given a diet of 35.3 or 628.8 μg TBCO / g food for 14 days after which, stage IV oocytes were isolated to assess maturation in response to maturation inducing hormone. To explore potential molecular mechanisms, abundances of mRNAs of a suite of genes that regulate oocyte maturation were quantified by use of quantitative real-time PCR, and abundances of microRNAs were determined by use of miRNAseq. Ex vivo maturation of oocytes from fish exposed to TBCO was significantly less than maturation of oocytes from control fish. The percentage of oocytes which matured from control fish and those exposed to low and high TBCO were 89, 71, and 67%, respectively. Among the suite of genes known to regulate oocyte maturation, mRNA abundance of insulin like growth factor-3 was decreased by 1.64- and 3.44-fold in stage IV oocytes from females given the low and high concentrations of TBCO, respectively, compared to the control group. Abundances of microRNAs regulating the expression of proteins that regulate oocyte maturation, including processes related to insulin-like growth factor, were significantly different in stage IV oocytes from fish exposed to TBCO. Overall, results of this study indicated that impaired oocyte maturation might be a mechanism of reduced reproductive performance in TBCO-exposed fish. Results also suggested that effects of TBCO on oocyte maturation might be due to molecular perturbations on insulin-like growth factor signaling and expression of microRNAs.

Alternative polyadenylation events differ dramatically between Tongcheng and Large White pigs in response to PRRSV infection

Alternative polyadenylation (APA) is a widespread post-transcriptional regulation mechanism that increases the biological complexity of transcriptome and proteome. However, it is unclear whether APA regulation plays a role in genetic resistance to porcine reproductive and respiratory syndrome virus (PRRSV). Here, we reported genome-wide APA regulation of porcine alveolar macrophages in PRRSV-resistant Tongcheng (TC) pigs and PRRSV-susceptible Large White (LW) pigs upon PRRSV infection. Using 3' mRNA sequencing strategy, we detected 75 981 high-quality APA sites in porcine alveolar macrophages of TC and LW pigs. Furthermore, 1202 and 1089 differentially expressed APA sites, as well as 79 and 117 untranslated region-APA switching genes were identified in TC pigs and LW pigs upon PRRSV infection respectively. The APA events in TC pigs and LW pigs were involved in different biological pathways, while APA events in TC pigs are directly associated with the immune response to PRRSV infection. In addition, we identified genetic variations affecting polyadenylation signal between TC pigs and LW pigs. These findings would provide helpful information on APA regulation for further understanding of genetic resistance to PRRSV.

LABRAT reveals association of alternative polyadenylation with transcript localization, RNA binding protein expression, transcription speed, and cancer survival

BACKGROUND

The sequence content of the 3' UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging.

RESULTS

To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated.

CONCLUSIONS

We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.

A CRISPR RNA-binding protein screen reveals regulators of RUNX1 isoform generation

The proper balance of hematopoietic stem cell (HSC) self-renewal and differentiation is critical for normal hematopoiesis and is disrupted in hematologic malignancy. Among regulators of HSC fate, transcription factors have a well-defined central role, and mutations promote malignant transformation. More recently, studies have illuminated the importance of posttranscriptional regulation by RNA-binding proteins (RBPs) in hematopoiesis and leukemia development. However, the RBPs involved and the breadth of regulation are only beginning to be elucidated. Furthermore, the intersection between posttranscriptional regulation and hematopoietic transcription factor function is poorly understood. Here, we studied the posttranscriptional regulation of RUNX1, a key hematopoietic transcription factor. Alternative polyadenylation (APA) of RUNX1 produces functionally antagonistic protein isoforms (RUNX1a vs RUNX1b/c) that mediate HSC self-renewal vs differentiation, an RNA-processing event that is dysregulated in malignancy. Consequently, RBPs that regulate this event directly contribute to healthy and aberrant hematopoiesis. We modeled RUNX1 APA using a split GFP minigene reporter and confirmed the sensitivity of our model to detect changes in RNA processing. We used this reporter in a clustered regularly interspaced short palindromic repeats (CRISPR) screen consisting of single guide RNAs exclusively targeting RBPs and uncovered HNRNPA1 and KHDRBS1 as antagonistic regulators of RUNX1a isoform generation. Overall, our study provides mechanistic insight into the posttranscriptional regulation of a key hematopoietic transcription factor and identifies RBPs that may have widespread and important functions in hematopoiesis.

Generation of 3'UTR knockout cell lines by CRISPR/Cas9-mediated genome editing

In addition to the protein code, messenger RNAs (mRNAs) also contain untranslated regions (UTRs). 3'UTRs span the region between the translational stop codon and the poly(A) tail. Sequence elements located in 3'UTRs are essential for pre-mRNA processing. 3'UTRs also contain elements that can regulate protein abundance, localization, and function. At least half of all human genes use alternative cleavage and polyadenylation (APA) to further diversify the regulatory potential of protein functions. Traditional gene editing approaches are designed to disrupt the production of functional proteins. Here, we describe a method that allows investigators to manipulate 3'UTR sequences of endogenous genes for both single- 3'UTR and multi-3'UTR genes. As 3'UTRs can regulate individual functions of proteins, techniques to manipulate 3'UTRs at endogenous gene loci will help to disentangle multi-functionality of proteins. Furthermore, the ability to directly examine the impact of gene regulatory elements in 3'UTRs will provide further insights into their functional significance.

Systematic identification of A-to-I RNA editing in zebrafish development and adult organs

A-to-I RNA editing is a common post transcriptional mechanism, mediated by the Adenosine deaminase that acts on RNA (ADAR) enzymes, that increases transcript and protein diversity. The study of RNA editing is limited by the absence of editing maps for most model organisms, hindering the understanding of its impact on various physiological conditions. Here, we mapped the vertebrate developmental landscape of A-to-I RNA editing, and generated the first comprehensive atlas of editing sites in zebrafish. Tens of thousands unique editing events and 149 coding sites were identified with high-accuracy. Some of these edited sites are conserved between zebrafish and humans. Sequence analysis of RNA over seven developmental stages revealed high levels of editing activity in early stages of embryogenesis, when embryos rely on maternal mRNAs and proteins. In contrast to the other organisms studied so far, the highest levels of editing were detected in the zebrafish ovary and testes. This resource can serve as the basis for understanding of the role of editing during zebrafish development and maturity.

Thermostability, Tunability, and Tenacity of RNA as Rubbery Anionic Polymeric Materials in Nanotechnology and Nanomedicine-Specific Cancer Targeting with Undetectable Toxicity

RNA nanotechnology is the bottom-up self-assembly of nanometer-scale architectures, resembling LEGOs, composed mainly of RNA. The ideal building material should be (1) versatile and controllable in shape and stoichiometry, (2) spontaneously self-assemble, and (3) thermodynamically, chemically, and enzymatically stable with a long shelf life. RNA building blocks exhibit each of the above. RNA is a polynucleic acid, making it a polymer, and its negative-charge prevents nonspecific binding to negatively charged cell membranes. The thermostability makes it suitable for logic gates, resistive memory, sensor set-ups, and NEM devices. RNA can be designed and manipulated with a level of simplicity of DNA while displaying versatile structure and enzyme activity of proteins. RNA can fold into single-stranded loops or bulges to serve as mounting dovetails for intermolecular or domain interactions without external linking dowels. RNA nanoparticles display rubber- and amoeba-like properties and are stretchable and shrinkable through multiple repeats, leading to enhanced tumor targeting and fast renal excretion to reduce toxicities. It was predicted in 2014 that RNA would be the third milestone in pharmaceutical drug development. The recent approval of several RNA drugs and COVID-19 mRNA vaccines by FDA suggests that this milestone is being realized. Here, we review the unique properties of RNA nanotechnology, summarize its recent advancements, describe its distinct attributes inside or outside the body and discuss potential applications in nanotechnology, medicine, and material science.

Multi-Dimensional Transcriptome Analysis Reveals Modulation of Cholesterol Metabolism as Highly Integrated Response to Brain Injury

Zebrafish is an attractive model to investigate regeneration of the nervous system. Despite major progress in our understanding of the underlying processes, the transcriptomic changes are largely unknown. We carried out a computational analysis of the transcriptome of the regenerating telencephalon integrating changes in the expression of mRNAs, their splice variants and investigated the putative role of regulatory RNAs in the modulation of these transcriptional changes. Profound changes in the expression of genes and their splice variants engaged in many distinct processes were observed. Differential transcription and splicing are important processes in response to injury of the telencephalon. As exemplified by the coordinated regulation of the cholesterol synthesizing enzymes and transporters, the genome responded to injury of the telencephalon in a multi-tiered manner with distinct and interwoven changes in expression of enzymes, transporters and their regulatory molecules. This coordinated genomic response involved a decrease of the mRNA of the key transcription factor SREBF2, induction of microRNAs (miR-182, miR-155, miR-146, miR-31) targeting cholesterol genes, shifts in abundance of splice variants as well as regulation of long non-coding RNAs. Cholesterol metabolism appears to be switched from synthesis to relocation of cholesterol. Based on our in silico analyses, this switch involves complementary and synergistic inputs by different regulatory principles. Our studies suggest that adaptation of cholesterol metabolism is a key process involved in regeneration of the injured zebrafish brain.