Use Of Alternative Polyadenylation Sites Research Articles

An integrative model for alternative polyadenylation, IntMAP, delineates mTOR-modulated endoplasmic reticulum stress response.

3′-untranslated regions (UTRs) can vary through the use of alternative polyadenylation sites during pre-mRNA processing. Multiple publically available pipelines combining high profiling technologies and bioinformatics tools have been developed to catalog changes in 3′-UTR lengths. In our recent RNA-seq experiments using cells with hyper-activated mammalian target of rapamycin (mTOR), we found that cellular mTOR activation leads to transcriptome-wide alternative polyadenylation (APA), resulting in the activation of multiple cellular pathways. Here, we developed a novel bioinformatics algorithm, IntMAP, which integrates RNA-Seq and PolyA Site (PAS)-Seq data for a comprehensive characterization of APA events. By applying IntMAP to the datasets from cells with hyper-activated mTOR, we identified novel APA events that could otherwise not be identified by either profiling method alone. Several transcription factors including Cebpg (CCAAT/enhancer binding protein gamma) were among the newly discovered APA transcripts, indicating that diverse transcriptional networks may be regulated by mTOR-coordinated APA. The prevention of APA in Cebpg using the CRISPR/cas9-mediated genome editing tool showed that mTOR-driven 3′-UTR shortening in Cebpg is critical in protecting cells from endoplasmic reticulum (ER) stress. Taken together, we present IntMAP as a new bioinformatics algorithm for APA analysis by which we expand our understanding of the physiological role of mTOR-coordinated APA events to ER stress response. IntMAP toolbox is available at http://compbio.cs.umn.edu/IntMAP/.

Alternative polyadenylation drives genome-to-phenome information detours in the AMPK\u03b11 and AMPK\u03b12 knockout mice

Currently available mouse knockout (KO) lines remain largely uncharacterized for genome-to-phenome (G2P) information flows. Here we test our hypothesis that altered myogenesis seen in AMPKα1- and AMPKα2-KO mice is caused by use of alternative polyadenylation sites (APSs). AMPKα1 and AMPKα2 are two α subunits of adenosine monophosphate-activated protein kinase (AMPK), which serves as a cellular sensor in regulation of many biological events. A total of 56,483 APSs were derived from gastrocnemius muscles. The differentially expressed APSs (DE-APSs) that were down-regulated tended to be distal. The DE-APSs that were related to reduced and increased muscle mass were down-regulated in AMPKα1-KO mice, but up-regulated in AMPKα2-KO mice, respectively. Five genes: Car3 (carbonic anhydrase 3), Mylk4 (myosin light chain kinase family, member 4), Neb (nebulin), Obscn (obscurin) and Pfkm (phosphofructokinase, muscle) utilized different APSs with potentially antagonistic effects on muscle function. Overall, gene knockout triggers genome plasticity via use of APSs, completing the G2P processes. However, gene-based analysis failed to reach such a resolution. Therefore, we propose that alternative transcripts are minimal functional units in genomes and the traditional central dogma concept should be now examined under a systems biology approach.

Circadian processes in the RNA life cycle.

The circadian clock drives daily rhythms of multiple physiological processes, allowing organisms to anticipate and adjust to periodic changes in environmental conditions. These physiological rhythms are associated with robust oscillations in the expression of at least 30% of expressed genes. While the ability for the endogenous timekeeping system to generate a 24-hr cycle is a cell-autonomous mechanism based on negative autoregulatory feedback loops of transcription and translation involving core-clock genes and their protein products, it is now increasingly evident that additional mechanisms also govern the circadian oscillations of clock-controlled genes. Such mechanisms can take place post-transcriptionally during the course of the RNA life cycle. It has been shown that many steps during RNA processing are regulated in a circadian manner, thus contributing to circadian gene expression. These steps include mRNA capping, alternative splicing, changes in splicing efficiency, and changes in RNA stability controlled by the tail length of polyadenylation or the use of alternative polyadenylation sites. RNA transport can also follow a circadian pattern, with a circadian nuclear retention driven by rhythmic expression within the nucleus of particular bodies (the paraspeckles) and circadian export to the cytoplasm driven by rhythmic proteins acting like cargo. Finally, RNA degradation may also follow a circadian pattern through the rhythmic involvement of miRNAs. In this review, we summarize the current knowledge of the post-transcriptional circadian mechanisms known to play a prominent role in shaping circadian gene expression in mammals. This article is categorized under: RNA Processing > Splicing Regulation/Alternative Splicing RNA Processing > RNA Editing and Modification RNA Export and Localization > Nuclear Export/Import.

Global 3' UTR shortening has a limited effect on protein abundance in proliferating T cells.

Alternative polyadenylation is a cellular mechanism that generates mRNA isoforms differing in their 3' untranslated regions (3' UTRs). Changes in polyadenylation site usage have been described upon induction of proliferation in resting cells, but the underlying mechanism and functional significance of this phenomenon remain largely unknown. To understand the functional consequences of shortened 3' UTR isoforms in a physiological setting, we used 3' end sequencing and quantitative mass spectrometry to determine polyadenylation site usage, mRNA and protein levels in murine and human naive and activated T cells. Although 3' UTR shortening in proliferating cells is conserved between human and mouse, orthologous genes do not exhibit similar expression of alternative 3' UTR isoforms. We generally find that 3' UTR shortening is not accompanied by a corresponding change in mRNA and protein levels. This suggests that although 3' UTR shortening may lead to changes in the RNA-binding protein interactome, it has limited effects on protein output.

Predominance of spliceosomal complex formation over polyadenylation site selection in TDP-43 autoregulation

TDP-43 is a nuclear protein involved in many aspects of RNA metabolism. To ensure cellular viability, its expression levels within cells must be tightly regulated. We have previously demonstrated that TDP-43 autoregulation occurs through the activation of a normally silent intron in its 3′-UTR sequence that results in the use of alternative polyadenylation sites. In this work, we analyse which is the dominant event in autoregulation: the recognition of the splice sites of 3′-UTR intron 7 or the intrinsic quality of the alternative polyadenylation sites. A panel of minigene constructs was tested for autoregulation functionality, protein production and subcellular messenger RNA localization. Our data clearly indicate that constitutive spliceosome complex formation across intron 7 does not lead to high protein production but, on the contrary, to lower TDP-43 messenger RNA and protein levels. This is due to altered nucleocytoplasmic distribution of the RNA that is mostly retained in the nucleus and degraded. This study provides a novel in-depth characterization of how RNA binding proteins can autoregulate their own levels within cells, an essential regulatory process in maintaining cellular viability.

The use of alternative polyadenylation sites renders integrin β1 (Itgb1) mRNA isoforms with differential stability during mammary gland development

Integrins are heterodimeric cell-surface adhesion receptors that play a critical role in tissue development. Characterization of the full-length mRNA encoding the β1 subunit (Itgb1) revealed an alternative functional cleavage and polyadenylation site that yields a new Itgb1 mRNA isoform 578 bp shorter than that previously reported. Using a variety of experimental and bioinformatic approaches, we found that the two Itgb1 isoforms are expressed at different levels in a variety of mouse tissues, including the mammary gland, where they are differentially regulated at successive developmental stages. The longer mRNA species is prevelant during lactation, whereas the shorter is induced after weaning. In 3D cultures, where expression of integrin β1 protein is required for normal formation of acini, experimental blockade of the longer isoform induced enhanced expression of the shorter species which allowed normal morphological mammary differentiation. The short isoform lacks AU-rich motifs and miRNA target sequences that are potentially implicated in the regulation of mRNA stability and translation efficiency. We further determined that the AU-binding protein HuR appears to selectively stabilize the longer isoform in the mammary gland. In summary, the results of the present study identify a new regulatory instance involved in the fine-tuning of Itgb1 expression during mammary gland development and function.

Genome-wide mapping of polyadenylation sites in fission yeast reveals widespread alternative polyadenylation

Regulatory elements in the 3′ untranslated regions (UTRs) of eukaryotic mRNAs influence mRNA localization, translation, and stability. 3′-UTR length is determined by the location at which mRNAs are cleaved and polyadenylated. The use of alternative polyadenylation sites is common, and can be regulated in different situations. I present a new method to identify cleavage and polyadenylation sites (CSs) at the genome-wide level. The approach is strand-specific, avoids RNA enzymatic modification steps that can introduce sequence-specific biases, and uses unique molecular identifiers to ensure that all identified CS originates from individual RNA molecules. I applied this method to create the first comprehensive genome-wide map of polyadenylation sites of the fission yeast Schizosaccharomyces pombe, comprising the analysis of 2,021,000 individual mRNAs that defined 8,883 CSs. CSs were identified for 90% of coding genes and 50% of ncRNAs. Alternative polyadenylation was prevalent in both groups, with 41% and 45% of all detected genes, respectively, displaying more than one CS. The specificity of the cleavage reaction was gene-specific, resulting in highly variable levels of heterogeneity in 3′-UTR lengths. Finally, I show that for both coding and non-coding genes, the most common regulatory motif associated with CSs in fission yeast is the canonical human AAUAAA sequence.

CPEB1 coordinates alternative 3'-UTR formation with translational regulation.

More than half of mammalian genes generate multiple messenger RNA isoforms that differ in their 3' untranslated regions (3' UTRs) and therefore in regulatory sequences, often associated with cell proliferation and cancer; however, the mechanisms coordinating alternative 3'-UTR processing for specific mRNA populations remain poorly defined. Here we report that the cytoplasmic polyadenylation element binding protein 1 (CPEB1), an RNA-binding protein that regulates mRNA translation, also controls alternative 3'-UTR processing. CPEB1 shuttles to the nucleus, where it co-localizes with splicing factors and mediates shortening of hundreds of mRNA 3' UTRs, thereby modulating their translation efficiency in the cytoplasm. CPEB1-mediated 3'-UTR shortening correlates with cell proliferation and tumorigenesis. CPEB1 binding to pre-mRNAs not only directs the use of alternative polyadenylation sites, but also changes alternative splicing by preventing U2AF65 recruitment. Our results reveal a novel function of CPEB1 in mediating alternative 3'-UTR processing, which is coordinated with regulation of mRNA translation, through its dual nuclear and cytoplasmic functions.

A high-fat diet increases l-carnitine synthesis through a differential maturation of the Bbox1 mRNAs

l-carnitine is a key molecule in both mitochondrial and peroxisomal lipid metabolisms. l-carnitine is biosynthesized from gamma-butyrobetaine by a reaction catalyzed by the gamma-butyrobetaine hydroxylase (Bbox1). The aim of this work was to identify molecular mechanisms involved in the regulation of l-carnitine biosynthesis and availability. Using 3′ RACE, we identified four alternatively polyadenylated Bbox1 mRNAs in rat liver. We utilized a combination of in vitro experiments using hybrid constructs containing the Bbox1 3′ UTR and in vivo experiments on rat liver mRNAs to reveal specificities in the different Bbox1 mRNA isoforms, especially in terms of polyadenylation efficiency, mRNA stability and translation efficiency. This complex maturation process of the Bbox1 mRNAs in the liver was studied on rats fed a high-fat diet. High-fat diet selectively increased the level of three Bbox1 mRNA isoforms in rat liver and the alternative use of polyadenylation sites contributed to the global increase in Bbox1 enzymatic activity and l-carnitine levels. Our results show that the maturation of Bbox1 mRNAs is nutritionally regulated in the liver through a selective polyadenylation process to adjust l-carnitine biosynthesis to the energy supply.

Poly(A) binding protein nuclear 1 levels affect alternative polyadenylation

The choice for a polyadenylation site determines the length of the 3′-untranslated region (3′-UTRs) of an mRNA. Inclusion or exclusion of regulatory sequences in the 3′-UTR may ultimately affect gene expression levels. Poly(A) binding protein nuclear 1 (PABPN1) is involved in polyadenylation of pre-mRNAs. An alanine repeat expansion in PABPN1 (exp-PABPN1) causes oculopharyngeal muscular dystrophy (OPMD). We hypothesized that previously observed disturbed gene expression patterns in OPMD muscles may have been the result of an effect of PABPN1 on alternative polyadenylation, influencing mRNA stability, localization and translation. A single molecule polyadenylation site sequencing method was developed to explore polyadenylation site usage on a genome-wide level in mice overexpressing exp-PABPN1. We identified 2012 transcripts with altered polyadenylation site usage. In the far majority, more proximal alternative polyadenylation sites were used, resulting in shorter 3′-UTRs. 3′-UTR shortening was generally associated with increased expression. Similar changes in polyadenylation site usage were observed after knockdown or overexpression of expanded but not wild-type PABPN1 in cultured myogenic cells. Our data indicate that PABPN1 is important for polyadenylation site selection and that reduced availability of functional PABPN1 in OPMD muscles results in use of alternative polyadenylation sites, leading to large-scale deregulation of gene expression.

The Long, the Short, and the Micro: A PolyA Tale of Pax3 in Satellite Cells

The use of alternative polyadenylation sites is emerging as an important regulator of gene expression. In this issue of Cell Stem Cell, Boutet et al. (2012) report that alternative 3'UTRs of the Pax3 transcript restrict its expression to axial satellite cells through miR-mediated targeting of one of the isoforms.

A PASport to Cellular Proliferation

The use of alternative polyadenylation sites produces mRNA isoforms with different 3' untranslated regions. A recent report in Science (Sandberg et al., 2008) suggests that alternative polyadenylation is connected to microRNA-mediated regulation of gene expression as part of a global program for cellular proliferation.

Recent integrations of mammalianHmg retropseudogenes

We propose that select retropseudogenes of the high mobility group nonhistone chromosomal protein genes have recently integrated into mammalian genomes on the basis of the high sequence identity of the copies to the cDNA sequences derived from the original genes. These include the Hmg1 gene family in mice and the Hmgn2 family in humans. We investigated orthologous loci of several strains and species of Mus for presence or absence of apparently young Hmg1 retropseudogenes. Three of four analysed elements were specific to Mus musculus, two of which were not fixed, indicative of recent evolutionary origins. Additionally, we datamined a presumptive subfamily (Hmgz) of mouse Hmg1, but only identified one true element in the GenBank database, which is not consistent with a separate subfamily status. Two of four analysed Hmgn2 retropseudogenes were specific for the human genome, whereas a third was identified in human, chimpanzee and gorilla genomes, and a fourth additionally found in orangutan but absent in African green monkey. Flanking target-site duplications were consistent with LINE integration sites supporting LINE machinery for their mechanism of amplification. The human Hmgn2 retropseudogenes were full length, whereas the mouse Hmg1 elements were either full length or 3'-truncated at specific positions, most plausibly the result of use of alternative polyadenylation sites. The nature of their recent amplification success in relation to other retropseudogenes is unclear, although availability of a large number of transcripts during gametogenesis may be a reason. It is apparent that retropseudogenes continue to shape mammalian genomes, and may provide insight into the process of retrotransposition, as well as offer potential use as phylogenetic markers.

Brain-specific promoter and polyadenylation sites of the -adducin pre-mRNA generate an unusually long 3'-UTR

Adducins are a family of membrane skeleton proteins composed of α-, β- and γ-subunits that promote actin and spectrin association in erythrocytes. The α- and γ-subunits are expressed ubiquitously, while the β-subunit is found in brain and erythropoietic tissues. The brain β-adducin protein is similar in size to that of spleen, but the mRNA transcript is a brain-specific one that has not been yet characterized, having an estimated length of 8–9 kb instead of the 3–4 kb of spleen mRNA. Here, we show the molecular basis for these differences by determining the structure of the brain-specific β-adducin transcript in rats, mice and humans. We identified a brain-specific promoter in rodents that, apparently, was not conserved in humans. In addition, we present evidence that the brain-mRNAs are formed by a common mechanism consisting in the tissue-specific use of alternative polyadenylation sites generating unusually long 3′-untranslated region of up to 6.6 kb. This hypothesis is supported by the presence of highly-conserved regions flanking the brain-specific polyadenylation site that suggest the involvement of these sequences in the translational regulation, stability and/or subcellular localization of the β-adducin transcript in the brain.

Alternative mRNA Polyadenylation Can Potentially Affect Detection of Gene Expression by Affymetrix GeneChip?? Arrays

DNA microarrays have been widely used to examine gene expression. The Affymetrix GeneChip is one of the most commonly used platforms, employing DNA probes of 25 nucleotides designed to hybridise to different regions of target mRNA. The targeted region is often biased toward the 3' end of mRNA, which can lead to biases in detection. A large number of mammalian genes can undergo alternative polyadenylation under different cellular conditions. Multiple polyadenylation sites can lead to variable transcripts with different hybridisation properties. Here, we surveyed probes on human, mouse and rat GeneChip arrays and found that the detection of a significant proportion of mRNAs can potentially be affected by alternative polyadenylation. This could lead to inaccurate interpretation of GeneChip data when the changes of expression values actually result from alternative use of polyadenylation sites.

Complexity and species variation of the kidney-type glutaminase gene

Increased expression of rat kidney-type glutaminase (KGA) during metabolic acidosis results from selective mRNA stabilization. This process is mediated by an 8-base AU-sequence that functions as a pH-response element (pHRE). LLC-PK1-FBPase+ cells, a pH-responsive porcine kidney cell line, express four distinct GA mRNAs. RNase H mapping indicated that three of the GA mRNAs are generated by use of alternative polyadenylation sites and are homologs of the rat KGA mRNA, while the fourth contains a different COOH-terminal coding and 3'-untranslated sequence. PCR cloning and sequencing established that the latter GA mRNA is the homolog of the human GAC mRNA. A rat GAC cDNA was also cloned from a rat kidney library. The 3'-untranslated regions of the GAC mRNAs, but not the porcine or human KGA mRNAs, contain identifiable pHREs. The human KGA gene spans 82 kb and is composed of 19 exons. The unique sequence from the hGAC cDNA is contained in a single exon. Thus in humans, alternative splicing of the initial transcript could produce two GA mRNAs, only one of which may be increased during acidosis.

The serotonin 5-HT2B receptor from the puffer fish Tetraodon fluviatilis: cDNA cloning, genomic organization and alternatively spliced variants

We cloned the 5-HT2B serotonin receptor from the puffer fish Tetraodon fluviatilis. Two cDNAs differing in length because of the use of alternative polyadenylation sites were isolated. We partly characterized the genomic organization of the 5-HT2B gene and we found two introns conserved in position between the puffer fish and mammals. In addition, four splice variants which would generate truncated forms of the receptor were detected.

HnRNP F influences binding of a 64-kilodalton subunit of cleavage stimulation factor to mRNA precursors in mouse B cells.

Previous studies on the regulation of polyadenylation of the immunoglobulin (Ig) heavy-chain pre-mRNA argued for trans-acting modifiers of the cleavage-polyadenylation reaction operating differentially during B-cell developmental stages. Using four complementary approaches, we demonstrate that a change in the level of hnRNP F is an important determinant in the regulated use of alternative polyadenylation sites between memory and plasma stage B cells. First, by Western analyses of cellular proteins, the ratio of hnRNP F to H or H' was found to be higher in memory B cells than in plasma cells. In memory B cells the activity of CstF-64 binding to pre-mRNA, but not its amount, was reduced. Second, examination of the complexes formed on input pre-mRNA in nuclear extracts revealed large assemblages containing hnRNP H, H', and F but deficient in CstF-64 in memory B-cell extracts but not in plasma cells. Formation of these large complexes is dependent on the region downstream of the AAUAAA in pre-mRNA, suggesting that CstF-64 and the hnRNPs compete for a similar region. Third, using a recombinant protein we showed that hnRNP F could bind to the region downstream of a poly(A) site, block CstF-64 association with RNA, and inhibit the cleavage reaction. Fourth, overexpression of recombinant hnRNP F in plasma cells resulted in a decrease in the endogenous Ig heavy-chain mRNA secretory form-to-membrane ratio. These results demonstrate that mammalian hnRNP F can act as a negative regulator in the pre-mRNA cleavage reaction and that increased expression of F in memory B cells contributes to the suppression of the Ig heavy-chain secretory poly(A) site.

Personality traits' trending to influence outcome in bipolar disorder

Exon trapping was performed from a partial cosmid, PAC, and P1 clone contig from human chromosome 21 betweenMX1and 21qter to identify genes that may be involved in the pathogenesis of Down syndrome or several of the genetic diseases that map to chromosome 21q22.3. One 19-bp exon showed identity to three ESTs. The complete sequence of the EST clones, RT-PCR, and cDNA library screening were used to determine the full-length cDNA sequence of 2.2 kb with an open reading frame of 256-amino-acids. The putative 256-amino-acid peptide has homology with a hypotheticalCaehorhabditis elegansprotein of unknown function. Northern blot analysis of this gene, termed C21orf2 (chromosome 21 open reading frame 2), revealed two ubiquitously expressed mRNAs of 2.2 and 1.2 kb produced by use of alternative polyadenylation sites. Hybridization of the EST clones to a cosmid contig in chromosome 21q22.3 mapped C21orf2 just distal toPFKL,a critical mapping region for several genetic diseases. Comparison to publicly available genomic sequence, and additional data, revealed that the gene is split into seven exons over 10.5 kb, further refining the mapping position to only 1.2 kb distal toPFKLwith the direction of transcription toward the centromere. The 5′UTR is contiguous withD21S400,and intron 2 contains a 52-bp VNTR polymorphism. Given its mapping position, C21orf2 is a candidate for involvement in disorders including autoimmune polyglandular disease type I (also called autoimmune polyendocrinopathy candidiasis ectodermal dystrophy or APECED) and the autosomal nonsyndromic deafness loci, DFNB8 and DFNB10. Mutation analysis using sequencing of RT-PCR and genomic DNA-derived PCR products, SSCP, and Southern and Northern blot analyses in APECED patients excluded C21orf2 as the gene for APECED.

Assessing the feasibility of an allelic association between bipolar affective disorder and region 21q22.3

Exon trapping was performed from a partial cosmid, PAC, and P1 clone contig from human chromosome 21 betweenMX1and 21qter to identify genes that may be involved in the pathogenesis of Down syndrome or several of the genetic diseases that map to chromosome 21q22.3. One 19-bp exon showed identity to three ESTs. The complete sequence of the EST clones, RT-PCR, and cDNA library screening were used to determine the full-length cDNA sequence of 2.2 kb with an open reading frame of 256-amino-acids. The putative 256-amino-acid peptide has homology with a hypotheticalCaehorhabditis elegansprotein of unknown function. Northern blot analysis of this gene, termed C21orf2 (chromosome 21 open reading frame 2), revealed two ubiquitously expressed mRNAs of 2.2 and 1.2 kb produced by use of alternative polyadenylation sites. Hybridization of the EST clones to a cosmid contig in chromosome 21q22.3 mapped C21orf2 just distal toPFKL,a critical mapping region for several genetic diseases. Comparison to publicly available genomic sequence, and additional data, revealed that the gene is split into seven exons over 10.5 kb, further refining the mapping position to only 1.2 kb distal toPFKLwith the direction of transcription toward the centromere. The 5′UTR is contiguous withD21S400,and intron 2 contains a 52-bp VNTR polymorphism. Given its mapping position, C21orf2 is a candidate for involvement in disorders including autoimmune polyglandular disease type I (also called autoimmune polyendocrinopathy candidiasis ectodermal dystrophy or APECED) and the autosomal nonsyndromic deafness loci, DFNB8 and DFNB10. Mutation analysis using sequencing of RT-PCR and genomic DNA-derived PCR products, SSCP, and Southern and Northern blot analyses in APECED patients excluded C21orf2 as the gene for APECED.