The Role of Alternative Polyadenylation in the Regulation of Subcellular RNA Localization

Alternative polyadenylation (APA) is a widespread and conserved regulatory mechanism that generates diverse 3′ ends on mRNA. APA patterns are often tissue specific and play an important role in cellular processes such as cell proliferation, differentiation, and response to stress. Many APA sites are found in 3′ UTRs, generating mRNA isoforms with different 3′ UTR contents. These alternate 3′ UTR isoforms can change how the transcript is regulated, affecting its stability and translation. Since the subcellular localization of a transcript is often regulated by 3′ UTR sequences, this implies that APA can also change transcript location. However, this connection between APA and RNA localization has only recently been explored. In this review, we discuss the role of APA in mRNA localization across distinct subcellular compartments. We also discuss current challenges and future advancements that will aid our understanding of how APA affects RNA localization and molecular mechanisms that drive these processes.

Mature Neurons sensitivity to oxidative stress is epigenetically programmed by alternative splicing and mRNA stability

Neuronal differentiation is a complex process that entails extensive morphological, transcriptional, metabolic, and functional changes that dictate neuronal lineage commitment. Much less understood is the role that epigenetic and epi-transcriptional reprogramming plays in the process of neuronal differentiation and maturation. To depict the whole landscape of transcriptomics and epigenetic changes during neuronal differentiation and maturation, we differentiated SH-SY5Y cells and performed RNA sequencing on differentiated and undifferentiated cells. 728 differentially expressed genes (DEGs) enriched in synaptic signaling and cell morphogenesis pathways were observed. Moreover, transcriptome-wide mRNA stability profiling revealed that genes with altered stability were exceptionally enriched for redox homeostasis pathways. Mature neurons are known to be highly sensitive to oxidative stress, which is crucial in the pathophysiology of neurodegenerative disease. Our results suggest that this heightened sensitivity is regulated at the mRNA stability level (i.e., epigenetic) rather than at the transcriptional level. Alternative splicing analysis revealed the exon skipping and alternative mRNA isoforms enriched for morphogenesis related pathway. Alternatively, alternative 5 and 3 prime splicing site, intron retention and mutually exclusive exon events exclusively clustered in the translation and translation initiation pathways, suggesting the potential effect of alternative splicing on translation following neuronal maturation. Splice motif analysis revealed enriched motifs for RBPs that regulate various splice types and can be further correlated to distinct phenotypical changes during neuronal differentiation and maturation. Here we present an extensive exploration of the transcriptional and epigenetic changes and their potential association with the process of neuronal differentiation, providing a new insight into understanding the molecular mechanism of neuronal function and behavior.

FMRP regulates mRNAs encoding distinct functions in the cell body and dendrites of CA1 pyramidal neurons

Neurons rely on translation of synaptic mRNAs in order to generate activity-dependent changes in plasticity. Here we develop a strategy combining compartment-specific CLIP and TRAP in conditionally tagged mice to precisely define the ribosome-bound dendritic transcriptome of CA1 pyramidal neurons. We identify CA1 dendritic transcripts with differentially localized mRNA isoforms generated by alternative polyadenylation and alternative splicing, including many which have altered protein-coding capacity. Among dendritic mRNAs, FMRP targets were found to be overrepresented. Cell-type specific FMRP-CLIP and TRAP in microdissected CA1 neuropil revealed 383 dendritic FMRP targets and suggests that FMRP differentially regulates functionally distinct modules in CA1 dendrites and cell bodies. FMRP regulates ~15-20% of mRNAs encoding synaptic functions and 10% of chromatin modulators, in the dendrite and cell body, respectively. In the absence of FMRP, dendritic FMRP targets had increased ribosome association, consistent with a function for FMRP in synaptic translational repression. Conversely, downregulation of FMRP targets involved in chromatin regulation in cell bodies and suggest a role for FMRP in stabilizing mRNAs containing stalled ribosomes in this compartment. Together, the data support a model in which FMRP regulates the translation and expression of synaptic and nuclear proteins within different compartments of a single neuronal cell type.

Long-read isoform sequencing reveals tissue-specific isoform expression between active and hibernating brown bears (Ursus arctos)

Understanding hibernation in brown bears (Ursus arctos) can provide insight into some human diseases. During hibernation, brown bears experience periods of insulin resistance, physical inactivity, extreme bradycardia, obesity, and the absence of urine production. These states closely mimic aspects of human diseases such as type 2 diabetes, muscle atrophy, as well as renal and heart failure. The reversibility of these states from hibernation to active season enables the identification of mediators with possible therapeutic value for humans. Recent studies have identified genes and pathways that are differentially expressed between active and hibernation seasons. However, little is known about the role of differential expression of gene isoforms on hibernation physiology. To identify both distinct and novel mRNA isoforms, full-length RNA-sequencing (Iso-Seq) was performed on adipose, skeletal muscle, and liver from three individuals sampled during both active and hibernation seasons. The existing reference annotation was improved by combining it with the Iso-Seq data. Short-read RNA-sequencing data from six individuals was mapped to the new reference annotation to quantify differential isoform usage between tissues and seasons. We identified differentially expressed isoforms in all three tissues, to varying degrees. Adipose had a high level of differential isoform usage with isoform switching, regardless of whether the genes were differentially expressed. Our analyses revealed that differential isoform usage, even in the absence of differential gene expression, is an important mechanism for modulating genes during hibernation. These findings demonstrate the value of isoform expression studies and will serve as the basis for deeper exploration into hibernation biology.

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.

Single-cell characterization of CRISPR-modified transcript isoforms with nanopore sequencing

We developed a single-cell approach to detect CRISPR-modified mRNA transcript structures. This method assesses how genetic variants at splicing sites and splicing factors contribute to alternative mRNA isoforms. We determine how alternative splicing is regulated by editing target exon-intron segments or splicing factors by CRISPR-Cas9 and their consequences on transcriptome profile. Our method combines long-read sequencing to characterize the transcript structure and short-read sequencing to match the single-cell gene expression profiles and gRNA sequence and therefore provides targeted genomic edits and transcript isoform structure detection at single-cell resolution.

Quantification of alternative 3′UTR isoforms from single cell RNA-seq data with scUTRquant

Although half of human genes use alternative polyadenylation (APA) to generate mRNA isoforms that encode the same protein but differ in their 3′UTRs, most single cell RNA-sequencing (scRNA-seq) pipelines only measure gene expression. Here, we describe an open-access pipeline, called scUTRquant (https://github.com/Mayrlab/scUTRquant), that measures gene and 3′UTR isoform expression from scRNA-seq data obtained from known cell types in any species. scUTRquant-derived gene and 3′UTR transcript counts were validated against standard methods which demonstrated their accuracy. 3′UTR isoform quantification was substantially more reproducible than previous methods. scUTRquant provides an atlas of high-confidence 3′ end cleavage sites at single-nucleotide resolution to allow APA comparison across mouse datasets. Analysis of 120 mouse cell types revealed that during differentiation genes either change their expression or they change their 3′UTR isoform usage. Therefore, we identified thousands of genes with 3′UTR isoform changes that have previously not been implicated in specific biological processes.

The Expression of Oct3/4A mRNA and Not Its Isoforms is Upregulated by the HPV16 E7 Oncoprotein

Oct3/4 a transcription factor is involved in maintaining the characteristics of cancer stem cells. Oct3/4 can be expressed differentially with respect to the progression of CC. In addition, Oct3/4 can give rise to three isoforms by alternative splicing of the mRNA Oct3/4A, Oct3/4B and Oct3/4B1. The aim of this study was to evaluate the mRNA expression from Oct3/4A, Oct3/4B and Oct3/4B1 in low-grade squamous intraepithelial lesion (LSIL), high-grade squamous intraepithelial lesion (HSIL), cervical cancer (CC) samples, and measure the effect of the HPV16 E7 oncoprotein on the mRNA expression from Oct3/4 isoforms in the C-33 A cell line. The expression levels of Oct3/4A, Oct3/4B and Oct3/4B1 mRNA were analyzed by reverse transcription quantitative polymerase chain reaction (RT-qPCR) in patients with LSILs, HSILs and CC. Additionally, C-33 A cells that expressed the HPV16 E7 oncoprotein were established to evaluate the effect of E7 on the expression of Oct3/4 mRNA isoforms. Oct3/4A (p=0.02), Oct3/4B (p=0. 001) and Oct3/4B1 (p<0. 0001) expression is significantly higher in patients with LSIL, HSIL and CC than in woman with non-IL. In the C-33 A cell line, the expression of Oct3/4A mRNA in the presence of the E7 oncoprotein increased compared to that in nontransfected C-33 A cells. Oct3/4B and Oct3/4B1 mRNA were expressed at similar levels among the different groups. These data indicate that only the mRNA of Oct3/4A is upregulated by the HPV16 E7 oncoprotein.

16 Identification of Structural Variants Associated with Mastitis in Holstein Dairy Cows Using Whole Genome Sequencing and RNA-sequencing

For the dairy industry, animal welfare and economic viability are key factors in determining the long-term sustainability of the industry. One challenge in lactating dairy cows is mastitis, as there are significant costs associated with each mastitis case due to treatment, milk loss, and potential cow culling. Differentially expressed (DE) genes, DE mRNA isoforms and DE long non-coding RNA (lncRNA) candidates were previously identified by our group between healthy and mastitic samples from six Holstein dairy cows. In total, 7 candidate genes, 5 mRNA isoforms and 4 lncRNA were targeted as they play an important role due to their association with mastitis or the immune system. The aim of the current research is to identify new structural variants in the transcriptome of these previously identified potential candidates using RNA-Sequencing. These structural variants could occur in both coding and intergenic regions of the genome and potentially impact the amino acid that is transcribed, creating a new functional or non-functional protein. In addition, whole genome sequencing (WGS) analysis was performed using hair samples collected from ten Holstein dairy cows (first lactation). Five of these animals had no previous reports of mastitis and five animals had at least one report of mastitis in her lifespan. The WGS results will be used to identify additional structural variants, especially within introns, of the previously identified candidates that could cause a genetic variation in individuals and an association to their susceptibility to mastitis disease.

MicroSalmon: A Comprehensive, Searchable Resource of Predicted MicroRNA Targets and 3′UTR Cis-Regulatory Elements in the Full-Length Sequenced Atlantic Salmon Transcriptome

Complete 3′UTRs unambiguously assigned to specific mRNA isoforms from the Atlantic salmon full-length (FL) transcriptome were collected into a 3′UTRome. miRNA response elements (MREs) and other cis-regulatory motifs were subsequently predicted and assigned to 3′UTRs of all FL-transcripts. The MicroSalmon GitHub repository provides all results. RNAHybrid and sRNAtoolbox tools predicted the MREs. UTRscan and the Teiresias algorithm predicted other 3′UTR cis-acting motifs, both known vertebrate motifs and putative novel motifs. MicroSalmon provides search programs to retrieve all FL-transcripts targeted by a miRNA (median number 1487), all miRNAs targeting an FL-transcript (median number 27), and other cis-acting motifs. As thousands of FL-transcripts may be targets of each miRNA, additional experimental strategies are necessary to reduce the likely true and relevant targets to a number that may be functionally validated. Low-complexity motifs known to affect mRNA decay in vertebrates were over-represented. Many of these were enriched in the terminal end, while purine- or pyrimidine-rich motifs with unknown functions were enriched immediately downstream of the stop codon. Furthermore, several novel complex motifs were over-represented, indicating conservation and putative function. In conclusion, MicroSalmon is an extensive and useful, searchable resource for study of Atlantic salmon transcript regulation by miRNAs and cis-acting 3′UTR motifs.