Computational analysis of 3'-ends of ESTs shows four classes of alternative polyadenylation in human, mouse, and rat

Deciphering Estrus Expression in Gilts: The Role of Alternative Polyadenylation and LincRNAs in Reproductive Transcriptomics

The fertility rate and litter size of female pigs are critically affected by the expression of estrus. The objective of this study was to elucidate the regulatory mechanisms of estrus expression by analyzing the differential expression of genes and long intergenic non-coding RNAs (lincRNA), as well as the utilization of alternative polyadenylation (APA) sites, in the vulva and vagina during the estrus and diestrus stages of Large White and indigenous Chinese Mi gilts. Our study revealed that the number of differentially expressed genes (DEG) in the vulva was less than that in the vagina, and the DEGs in the vulva were enriched in pathways such as "neural" pathways and steroid hormone responses, including the "Calcium signaling pathway" and "Oxytocin signaling pathway". The DEGs in the vagina were enriched in the "Metabolic pathways" and "VEGF signaling pathway". Furthermore, 27 and 21 differentially expressed lincRNAs (DEL), whose target genes were enriched in the "Endocrine resistance" pathway, were identified in the vulva and vagina, respectively. Additionally, we observed that 63 and 618 transcripts of the 3'-untranslated region (3'-UTR) were lengthened during estrus in the vulva and vagina, respectively. Interestingly, the genes undergoing APA events in the vulva exhibited species-specific enrichment in neural or steroid-related pathways, whereas those in the vagina were enriched in apoptosis or autophagy-related pathways. Further bioinformatic analysis of these lengthened 3'-UTRs revealed the presence of multiple miRNAs binding sites and cytoplasmic polyadenylation element (CPE) regulatory aspects. In particular, we identified more than 10 CPEs in the validated lengthened 3'-UTRs of the NFIX, PCNX4, CEP162 and ABHD2 genes using RT-qPCR. These findings demonstrated the involvement of APA and lincRNAs in the regulation of estrus expression in female pigs, providing new insights into the molecular mechanisms underlying estrus expression in pigs.

Single-cell and spatial mapping Identify cell types and signaling Networks in the human ureter

Tissue engineering offers a promising treatment strategy for ureteral strictures, but its success requires an in-depth understanding of the architecture, cellular heterogeneity, and signaling pathways underlying tissue regeneration. Here, we define and spatially map cell populations within the human ureter using single-cell RNA sequencing, spatial gene expression, and immunofluorescence approaches. We focus on the stromal and urothelial cell populations to enumerate the distinct cell types composing the human ureter and infer potential cell-cell communication networks underpinning the bi-directional crosstalk between these compartments. Furthermore, we analyze and experimentally validate the importance of the sonic hedgehog (SHH) signaling pathway in adult progenitor cell maintenance. The SHH-expressing basal cells support organoid generation in vitro and accurately predict the differentiation trajectory from basal progenitor cells to terminally differentiated umbrella cells. Our results highlight the essential processes involved in adult ureter tissue homeostasis and provide a blueprint for guiding ureter tissue engineering.

Alternative splicing is a Sorghum bicolor defense response to fungal infection

This study provides new insights that alternative splicing participates with transcriptional control in defense responses to Colletotrichum sublineola in sorghum In eukaryotic organisms, alternative splicing (AS) is an important post-transcriptional mechanism to generate multiple transcript isoforms from a single gene. Protein variants translated from splicing isoforms may have altered molecular characteristics in signal transduction and metabolic activities. However, which transcript isoforms will be translated into proteins and the biological functions of the resulting proteoforms are yet to be identified. Sorghum is one of the five major cereal crops, but its production is severely affected by fungal diseases. For example, sorghum anthracnose caused by Colletotrichum sublineola greatly reduces grain yield and biomass production. In this study, next-generation sequencing technology was used to analyze C. sublineola-inoculated sorghum seedlings compared with mock-inoculated control. It was identified that AS regulation may be as important as traditional transcriptional control during defense responses to fungal infection. Moreover, several genes involved in flavonoid and phenylpropanoid biosynthetic pathways were found to undergo multiple AS modifications. Further analysis demonstrated that non-conventional targets of both 5'- and 3'-splice sites were alternatively used in response to C. sublineola infection. Splicing factors were also affected at both transcriptional and post-transcriptional levels. As the first transcriptome report on C. sublineola infected sorghum, our work also suggested that AS plays crucial functions in defense responses to fungal invasion.

Alternative splicing and translation play important roles in hypoxic germination in rice

Post-transcriptional mechanisms (PTMs), including alternative splicing (AS) and alternative translation initiation (ATI), may explain the diversity of proteins involved in plant development and stress responses. Transcriptional regulation is important during the hypoxic germination of rice seeds, but the potential roles of PTMs in this process have not been characterized. We used a combination of proteomics and RNA sequencing to discover how AS and ATI contribute to plant responses to hypoxia. In total, 10 253 intron-containing genes were identified. Of these, ~1741 differentially expressed AS (DAS) events from 811 genes were identified in hypoxia-treated seeds compared with controls. Over 95% of these were not present in the list of differentially expressed genes. In particular, regulatory pathways such as the spliceosome, ribosome, endoplasmic reticulum protein processing and export, proteasome, phagosome, oxidative phosphorylation, and mRNA surveillance showed substantial AS changes under hypoxia, suggesting that AS responses are largely independent of transcriptional regulation. Considerable AS changes were identified, including the preferential usage of some non-conventional splice sites and enrichment of splicing factors in the DAS data sets. Taken together, these results not only demonstrate that AS and ATI function during hypoxic germination but they have also allowed the identification of numerous novel proteins/peptides produced via ATI.

Analyses of alternative polyadenylation: from old school biochemistry to high-throughput technologies

Alternations in usage of polyadenylation sites during transcription termination yield transcript isoforms from a gene. Recent findings of transcriptome-wide alternative polyadenylation (APA) as a molecular response to changes in biology position APA not only as a molecular event of early transcriptional termination but also as a cellular regulatory step affecting various biological pathways. With the development of high-throughput profiling technologies at a single nucleotide level and their applications targeted to the 3'-end of mRNAs, dynamics in the landscape of mRNA 3'-end is measureable at a global scale. In this review, methods and technologies that have been adopted to study APA events are discussed. In addition, various bioinformatics algorithms for APA isoform analysis using publicly available RNA-seq datasets are introduced. [BMB Reports 2017; 50(4): 201-207].

Omni-PolyA: a method and tool for accurate recognition of Poly(A) signals in human genomic DNA

BACKGROUND

Polyadenylation is a critical stage of RNA processing during the formation of mature mRNA, and is present in most of the known eukaryote protein-coding transcripts and many long non-coding RNAs. The correct identification of poly(A) signals (PAS) not only helps to elucidate the 3'-end genomic boundaries of a transcribed DNA region and gene regulatory mechanisms but also gives insight into the multiple transcript isoforms resulting from alternative PAS. Although progress has been made in the in-silico prediction of genomic signals, the recognition of PAS in DNA genomic sequences remains a challenge.

RESULTS

In this study, we analyzed human genomic DNA sequences for the 12 most common PAS variants. Our analysis has identified a set of features that helps in the recognition of true PAS, which may be involved in the regulation of the polyadenylation process. The proposed features, in combination with a recognition model, resulted in a novel method and tool, Omni-PolyA. Omni-PolyA combines several machine learning techniques such as different classifiers in a tree-like decision structure and genetic algorithms for deriving a robust classification model. We performed a comparison between results obtained by state-of-the-art methods, deep neural networks, and Omni-PolyA. Results show that Omni-PolyA significantly reduced the average classification error rate by 35.37% in the prediction of the 12 considered PAS variants relative to the state-of-the-art results.

CONCLUSIONS

The results of our study demonstrate that Omni-PolyA is currently the most accurate model for the prediction of PAS in human and can serve as a useful complement to other PAS recognition methods. Omni-PolyA is publicly available as an online tool accessible at www.cbrc.kaust.edu.sa/omnipolya/ .

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.

Processing and transcriptome expansion at the mRNA 3' end in health and disease: finding the right end

The human transcriptome is highly dynamic, with each cell type, tissue, and organ system expressing an ensemble of transcript isoforms that give rise to considerable diversity. Apart from alternative splicing affecting the "body" of the transcripts, extensive transcriptome diversification occurs at the 3' end. Transcripts differing at the 3' end can have profound physiological effects by encoding proteins with distinct functions or regulatory properties or by affecting the mRNA fate via the inclusion or exclusion of regulatory elements (such as miRNA or protein binding sites). Importantly, the dynamic regulation at the 3' end is associated with various (patho)physiological processes, including the immune regulation but also tumorigenesis. Here, we recapitulate the mechanisms of constitutive mRNA 3' end processing and review the current understanding of the dynamically regulated diversity at the transcriptome 3' end. We illustrate the medical importance by presenting examples that are associated with perturbations of this process and indicate resulting implications for molecular diagnostics as well as potentially arising novel therapeutic strategies.

CFIm25 regulates glutaminase alternative terminal exon definition to modulate miR-23 function

Alternative polyadenylation (APA) and alternative splicing (AS) provide mRNAs with the means to avoid microRNA repression through selective shortening or differential usage of 3'UTRs. The two glutaminase (GLS) mRNA isoforms, termed KGA and GAC, contain distinct 3'UTRs with the KGA isoform subject to repression by miR-23. We show that depletion of the APA regulator CFIm25 causes a strong shift to the usage of a proximal poly(A) site within the KGA 3'UTR and also alters splicing to favor exclusion of the GAC 3'UTR. Surprisingly, we observe that while miR-23 is capable of down-regulating the shortened KGA 3'UTR, it has only minor impact on the full-length KGA 3'UTR, demonstrating that additional potent negative regulation of GLS expression exists beyond this single microRNA targeting site. Finally, we show that the apoptosis induced upon down-regulation of the GAC isoform can be alleviated through concurrent reduction in CFIm25 expression, revealing the sensitivity of glutaminase expression to the levels of RNA processing factors. These results exemplify the complex interplay between RNA processing and microRNA repression in controlling glutamine metabolism in cancer cells.

New frontiers in translational control of the cancer genome

The past several years have seen dramatic leaps in our understanding of how gene expression is rewired at the translation level during tumorigenesis to support the transformed phenotype. This work has been driven by an explosion in technological advances and is revealing previously unimagined regulatory mechanisms that dictate functional expression of the cancer genome. In this Review we discuss emerging trends and exciting new discoveries that reveal how this translational circuitry contributes to specific aspects of tumorigenesis and cancer cell function, with a particular focus on recent insights into the role of translational control in the adaptive response to oncogenic stress conditions.

Dynamic Regulation of Tandem 3' Untranslated Regions in Zebrafish Spleen Cells during Immune Response

Alternative polyadenylation (APA) has been found to be involved in tumorigenesis, development, and cell differentiation, as well as in the activation of several subsets of immune cells in vitro. Whether APA takes place in immune responses in vivo is largely unknown. We profiled the variation in tandem 3' untranslated regions (UTRs) in pathogen-challenged zebrafish and identified hundreds of APA genes with ∼ 10% being immune response genes. The detected immune response APA genes were enriched in TLR signaling, apoptosis, and JAK-STAT signaling pathways. A greater number of microRNA target sites and AU-rich elements were found in the extended 3' UTRs than in the common 3' UTRs of these APA genes. Further analysis suggested that microRNA and AU-rich element-mediated posttranscriptional regulation plays an important role in modulating the expression of APA genes. These results indicate that APA is extensively involved in immune responses in vivo, and it may be a potential new paradigm for immune regulation.

Testis-specific products of the Drosophila melanogaster sbr gene, encoding nuclear export factor 1, are necessary for male fertility

The evolutionarily conserved nuclear export factor 1 (NXF1) provides mRNA export from the nucleus to the cytoplasm. We described several testis-specific transcripts of the Drosophila melanogaster nxf1 gene designated “sbr” in this species via different PCR approaches and CAGE-seq analysis. Characteristically, most of them have truncated 3′UTRs compared with those in other organs. In addition to regular transcripts, there are shorter transcripts that begin in intron 3 of the sbr gene. These short, 5′-truncated testis-specific transcripts vary in terms of transcription start site and their ability to exclude or retain the last 237 nucleotides of intron 3 in their 5′UTR. Using an anti-SBR antibody against the C-terminal portion of this protein, we detected the major SBR protein (74 kDa) in all analyzed organs of the fly as well as a new smaller protein (60 kDa) found only in the testes. This protein corresponds to the detected sbr transcripts that start in intron 3, based on its molecular mass. We investigated the sbr12 allele of the sbr gene, which is lethal in homozygous females and causes dominant sterility in heterozygous males. Sequencing of the sbr12 gene allele revealed a 30-bp deletion in exon 9 without a frame shift.Western blot analysiswith an SBR-specific antibody revealed two bands of the expected size in the testes of heterozygous males. Thus, a mutant protein along with the normal protein presents in the testes of lethal allele-bearing flies and the described shorter testis-specific variant of SBR may account for male sterility.

Alternative mRNA transcription, processing, and translation: insights from RNA sequencing

The human transcriptome comprises >80,000 protein-coding transcripts and the estimated number of proteins synthesized from these transcripts is in the range of 250,000 to 1 million. These transcripts and proteins are encoded by less than 20,000 genes, suggesting extensive regulation at the transcriptional, post-transcriptional, and translational level. Here we review how RNA sequencing (RNA-seq) technologies have increased our understanding of the mechanisms that give rise to alternative transcripts and their alternative translation. We highlight four different regulatory processes: alternative transcription initiation, alternative splicing, alternative polyadenylation, and alternative translation initiation. We discuss their transcriptome-wide distribution, their impact on protein expression, their biological relevance, and the possible molecular mechanisms leading to their alternative regulation. We conclude with a discussion of the coordination and the interdependence of these four regulatory layers.

3'LIFE: a functional assay to detect miRNA targets in high-throughput

MicroRNAs (miRNAs) are short non-coding RNAs that regulate gene output at the post-transcriptional level by targeting degenerate elements primarily in 3'untranslated regions (3'UTRs) of mRNAs. Individual miRNAs can regulate networks of hundreds of genes, yet for the majority of miRNAs few, if any, targets are known. Misexpression of miRNAs is also a major contributor to cancer progression, thus there is a critical need to validate miRNA targets in high-throughput to understand miRNAs' contribution to tumorigenesis. Here we introduce a novel high-throughput assay to detect miRNA targets in 3'UTRs, called Luminescent Identification of Functional Elements in 3'UTRs (3'LIFE). We demonstrate the feasibility of 3'LIFE using a data set of 275 human 3'UTRs and two cancer-relevant miRNAs, let-7c and miR-10b, and compare our results to alternative methods to detect miRNA targets throughout the genome. We identify a large number of novel gene targets for these miRNAs, with only 32% of hits being bioinformatically predicted and 27% directed by non-canonical interactions. Functional analysis of target genes reveals consistent roles for each miRNA as either a tumor suppressor (let-7c) or oncogenic miRNA (miR-10b), and preferentially target multiple genes within regulatory networks, suggesting 3'LIFE is a rapid and sensitive method to detect miRNA targets in high-throughput.

Genome-wide determination of poly(A) sites in Medicago truncatula: evolutionary conservation of alternative poly(A) site choice

Background

Alternative polyadenylation (APA) plays an important role in the post-transcriptional regulation of gene expression. Little is known about how APA sites may evolve in homologous genes in different plant species. To this end, comparative studies of APA sites in different organisms are needed. In this study, a collection of poly(A) sites in Medicago truncatula, a model system for legume plants, has been generated and compared with APA sites in Arabidopsis thaliana.

Results

The poly(A) tags from a deep-sequencing protocol were mapped to the annotated M. truncatula genome, and the identified poly(A) sites used to update the annotations of 14,203 genes. The results show that 64% of M. truncatula genes possess more than one poly(A) site, comparable to the percentages reported for Arabidopsis and rice. In addition, the poly(A) signals associated with M. truncatula genes were similar to those seen in Arabidopsis and other plants. The 3′-UTR lengths are correlated in pairs of orthologous genes between M. truncatula and Arabidopsis. Very little conservation of intronic poly(A) sites was found between Arabidopsis and M. truncatula, which suggests that such sites are likely to be species-specific in plants. In contrast, there is a greater conservation of CDS-localized poly(A) sites in these two species. A sizeable number of M. truncatula antisense poly(A) sites were found. A high percentage of the associated target genes possess Arabidopsis orthologs that are also associated with antisense sites. This is suggestive of important roles for antisense regulation of these target genes.

Conclusions

Our results reveal some distinct patterns of sense and antisense poly(A) sites in Arabidopsis and M. truncatula. In so doing, this study lends insight into general evolutionary trends of alternative polyadenylation in plants.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-615) contains supplementary material, which is available to authorized users.

Microarray meta-analysis of RNA-binding protein functions in alternative polyadenylation

Alternative polyadenylation (APA) is a post-transcriptional mechanism to generate diverse mRNA transcripts with different 3′UTRs from the same gene. In this study, we systematically searched for the APA events with differential expression in public mouse microarray data. Hundreds of genes with over-represented differential APA events and the corresponding experiments were identified. We further revealed that global APA differential expression occurred prevalently in tissues such as brain comparing to peripheral tissues, and biological processes such as development, differentiation and immune responses. Interestingly, we also observed widespread differential APA events in RNA-binding protein (RBP) genes such as Rbm3, Eif4e2 and Elavl1. Given the fact that RBPs are considered as the main regulators of differential APA expression, we constructed a co-expression network between APAs and RBPs using the microarray data. Further incorporation of CLIP-seq data of selected RBPs showed that Nova2 represses and Mbnl1 promotes the polyadenylation of closest poly(A) sites respectively. Altogether, our study is the first microarray meta-analysis in a mammal on the regulation of APA by RBPs that integrated massive mRNA expression data under a wide-range of biological conditions. Finally, we present our results as a comprehensive resource in an online website for the research community.

RNA-binding proteins in regulation of alternative cleavage and polyadenylation

Almost all eukaryotic pre-mRNAs are processed at the 3' end by the cleavage and polyadenylation (C/P) reaction, which preludes termination of transcription and gives rise to the poly(A) tail of mature mRNA. Genomic studies in recent years have indicated that most eukaryotic mRNA genes have multiple cleavage and polyadenylation sites (pAs), leading to alternative cleavage and polyadenylation (APA) products. APA isoforms generally differ in their 3' untranslated regions (3' UTRs), but can also have different coding sequences (CDSs). APA expands the repertoire of transcripts expressed from the genome, and is highly regulated under various physiological and pathological conditions. Growing lines of evidence have shown that RNA-binding proteins (RBPs) play important roles in regulation of APA. Some RBPs are part of the machinery for C/P; others influence pA choice through binding to adjacent regions. In this chapter, we review cis elements and trans factors involved in C/P, the significance of APA, and increasingly elucidated roles of RBPs in APA regulation. We also discuss analysis of APA using transcriptome-wide techniques as well as molecular biology approaches.

Mir-20a regulates in vitro mineralization and BMP signaling pathway by targeting BMP-2 transcript in fish

MicroRNAs (miRNAs) are important regulators of vertebrate development but their role during skeletogenesis remains unknown. In this regard, we investigated the mineralogenic activity of miR-20a, a miRNA associated with osteogenesis, in fish bone-derived cells. Expression of miR-20a was up-regulated during differentiation and its overexpression inhibited mineralization, suggesting a role in fish tissue calcification. In this regard, a conserved miR-20a binding site was identified in bone morphogenetic protein 2 (BMP-2) 3'UTR and its functionality was evidenced through luciferase assays, and further confirmed by western-blot and qPCR. Type II BMP receptor (BMPR2) is also targeted by miR-20a in mammalian systems and evidence was collected for the presence of a binding site in fish sequences. We propose that miR-20a is a regulator of BMP pathway through specific action on BMP-2 and possibly BMPR2. Overexpression of miR-20a was also shown to up-regulate matrix Gla protein (MGP) transcript, a physiological inhibitor of calcification previously found to form a complex with BMP-2. We propose that MGP may play a role in the anti-mineralogenic effect promoted by miR-20a by decreasing availability of BMP-2. This study gives new insights into miRNA-mediated regulation of BMP-2, and sheds light into the potential role of miR-20a as a regulator of skeletogenesis.

Genome-wide analysis of poly(A) site selection in Schizosaccharomyces pombe

Polyadenylation of pre-mRNAs, a critical step in eukaryotic gene expression, is mediated by cis elements collectively called the polyadenylation signal. Genome-wide analysis of such polyadenylation signals was missing in fission yeast, even though it is an important model organism. We demonstrate that the canonical AATAAA motif is the most frequent and functional polyadenylation signal in Schizosaccharomyces pombe. Using analysis of RNA-Seq data sets from cells grown under various physiological conditions, we identify 3' UTRs for nearly 90% of the yeast genes. Heterogeneity of cleavage sites is common, as is alternative polyadenylation within and between conditions. We validated the computationally identified sequence elements likely to promote polyadenylation by functional assays, including qRT-PCR and 3'RACE analysis. The biological importance of the AATAAA motif is underlined by functional analysis of the genes containing it. Furthermore, it has been shown that convergent genes require trans elements, like cohesin for efficient transcription termination. Here we show that convergent genes lacking cohesin (on chromosome 2) are generally associated with longer overlapping mRNA transcripts. Our bioinformatic and experimental genome-wide results are summarized and can be accessed and customized in a user-friendly database Pomb(A).

Nucleosome distribution near the 3' ends of genes in the human genome

By systematic analysis of high-throughput sequencing datasets from the human genome, we found that protein-coding genes have a specific chromatin structure near transcription termination sites relative to non-coding genes, one related to polyadenylation. Nucleosome was depleted near the site of cleavage/polyadenylation (polyA site) regardless of its relative position in the gene. DNA sequence plays an improtant role in nucleosome distribution, and conservative sequence elements and the protein binding to them are major determinants in causing nucleosome depletion near polyA sites. Furthermore, nucleosome occupancy was regulated by gene transcription and RNA polymerase II (RNAPII) occupancy. Our results reveal influences on nucleosome occupancy near polyadenylation sites and constitute evidence indicating that nucleosome distribution regulates 3' end processing of protein-coding genes.

DNA damage induces targeted, genome-wide variation of poly(A) sites in budding yeast

Systemic response to DNA damage and other stresses is a complex process that includes changes in the regulation and activity of nearly all stages of gene expression. One gene regulatory mechanism used by eukaryotes is selection among alternative transcript isoforms that differ in polyadenylation [poly(A)] sites, resulting in changes either to the coding sequence or to portions of the 3' UTR that govern translation, stability, and localization. To determine the extent to which this means of regulation is used in response to DNA damage, we conducted a global analysis of poly(A) site usage in Saccharomyces cerevisiae after exposure to the UV mimetic, 4-nitroquinoline 1-oxide (4NQO). Two thousand thirty-one genes were found to have significant variation in poly(A) site distributions following 4NQO treatment, with a strong bias toward loss of short transcripts, including many with poly(A) sites located within the protein coding sequence (CDS). We further explored one possible mechanism that could contribute to the widespread differences in mRNA isoforms. The change in poly(A) site profile was associated with an inhibition of cleavage and polyadenylation in cell extract and a decrease in the levels of several key subunits in the mRNA 3'-end processing complex. Sequence analysis identified differences in the cis-acting elements that flank putatively suppressed and enhanced poly(A) sites, suggesting a mechanism that could discriminate between variable and constitutive poly(A) sites. Our analysis indicates that variation in mRNA length is an important part of the regulatory response to DNA damage.

Characterization of human NLZ1/ZNF703 identifies conserved domains essential for proper subcellular localization and transcriptional repression

NET family members have recently emerged as important players in the development of multiple structures, from the trachea of fly larvae to the vertebrate eye and human breast cancers. However, their mechanisms of action are still poorly understood, and we lack a detailed characterization of their functional domains, as well as gene expression patterns-particularly in adult mammals. Here, we present a characterization of human NLZ1/ZNF703 (NocA-like zinc finger 1/Zinc finger 703), one of the two human NET family member genes. We show that the gene is ubiquitously expressed in adult human and mouse tissues, that three mRNA species with the same coding sequence are generated by alternative polyadenylation, and that the encoded protein contains six evolutionarily conserved domains, three of which are specific to NET proteins. Finally, we present functional evidence that these domains are necessary for proper subcellular distribution of and transcription repression by the NLZ1 protein, but not for its interaction with Groucho family co-repressors.

Alterations in polyadenylation and its implications for endocrine disease

INTRODUCTION

Polyadenylation is the process in which the pre-mRNA is cleaved at the poly(A) site and a poly(A) tail is added - a process necessary for normal mRNA formation. Genes with multiple poly(A) sites can undergo alternative polyadenylation (APA), producing distinct mRNA isoforms with different 3' untranslated regions (3' UTRs) and in some cases different coding regions. Two thirds of all human genes undergo APA. The efficiency of the polyadenylation process regulates gene expression and APA plays an important part in post-transcriptional regulation, as the 3' UTR contains various cis-elements associated with post-transcriptional regulation, such as target sites for micro-RNAs and RNA-binding proteins. Implications of alterations in polyadenylation for endocrine disease: Alterations in polyadenylation have been found to be causative of neonatal diabetes and IPEX (immune dysfunction, polyendocrinopathy, enteropathy, X-linked) and to be associated with type I and II diabetes, pre-eclampsia, fragile X-associated premature ovarian insufficiency, ectopic Cushing syndrome, and many cancer diseases, including several types of endocrine tumor diseases.

PERSPECTIVES

Recent developments in high-throughput sequencing have made it possible to characterize polyadenylation genome-wide. Antisense elements inhibiting or enhancing specific poly(A) site usage can induce desired alterations in polyadenylation, and thus hold the promise of new therapeutic approaches.

SUMMARY

This review gives a detailed description of alterations in polyadenylation in endocrine disease, an overview of the current literature on polyadenylation and summarizes the clinical implications of the current state of research in this field.

Alternative polyadenylation: new insights from global analyses

Recent studies have revealed widespread mRNA alternative polyadenylation (APA) in eukaryotes and its dynamic spatial and temporal regulation. APA not only generates proteomic and functional diversity, but also plays important roles in regulating gene expression. Global deregulation of APA has been demonstrated in a variety of human diseases. Recent exciting advances in the field have been made possible in a large part by high throughput analyses using newly developed experimental tools. Here I review the recent progress in global studies of APA and the insights that have emerged from these and other studies that use more conventional methods.

Plant polyadenylation factors: conservation and variety in the polyadenylation complex in plants

Background

Polyadenylation, an essential step in eukaryotic gene expression, requires both cis-elements and a plethora of trans-acting polyadenylation factors. The polyadenylation factors are largely conserved across mammals and fungi. The conservation seems also extended to plants based on the analyses of Arabidopsis polyadenylation factors. To extend this observation, we systemically identified the orthologs of yeast and human polyadenylation factors from 10 plant species chosen based on both the availability of their genome sequences and their positions in the evolutionary tree, which render them representatives of different plant lineages.

Results

The evolutionary trajectories revealed several interesting features of plant polyadenylation factors. First, the number of genes encoding plant polyadenylation factors was clearly increased from “lower” to “higher” plants. Second, the gene expansion in higher plants was biased to some polyadenylation factors, particularly those involved in RNA binding. Finally, while there are clear commonalities, the differences in the polyadenylation apparatus were obvious across different species, suggesting an ongoing process of evolutionary change. These features lead to a model in which the plant polyadenylation complex consists of a conserved core, which is rather rigid in terms of evolutionary conservation, and a panoply of peripheral subunits, which are less conserved and associated with the core in various combinations, forming a collection of somewhat distinct complex assemblies.

Conclusions

The multiple forms of plant polyadenylation complex, together with the diversified polyA signals may explain the intensive alternative polyadenylation (APA) and its regulatory role in biological functions of higher plants.

Mouse endogenous retroviruses can trigger premature transcriptional termination at a distance

Endogenous retrotransposons have caused extensive genomic variation within mammalian species, but the functional implications of such mobilization are mostly unknown. We mapped thousands of endogenous retrovirus (ERV) germline integrants in highly divergent, previously unsequenced mouse lineages, facilitating a comparison of gene expression in the presence or absence of local insertions. Polymorphic ERVs occur relatively infrequently in gene introns and are particularly depleted from genes involved in embryogenesis or that are highly expressed in embryonic stem cells. Their genomic distribution implies ongoing negative selection due to deleterious effects on gene expression and function. A polymorphic, intronic ERV at Slc15a2 triggers up to 49-fold increases in premature transcriptional termination and up to 39-fold reductions in full-length transcripts in adult mouse tissues, thereby disrupting protein expression and functional activity. Prematurely truncated transcripts also occur at Polr1a, Spon1, and up to ∼5% of other genes when intronic ERV polymorphisms are present. Analysis of expression quantitative trait loci (eQTLs) in recombinant BxD mouse strains demonstrated very strong genetic associations between the polymorphic ERV in cis and disrupted transcript levels. Premature polyadenylation is triggered at genomic distances up to >12.5 kb upstream of the ERV, both in cis and between alleles. The parent of origin of the ERV is associated with variable expression of nonterminated transcripts and differential DNA methylation at its 5'-long terminal repeat. This study defines an unexpectedly strong functional impact of ERVs in disrupting gene transcription at a distance and demonstrates that ongoing retrotransposition can contribute significantly to natural phenotypic diversity.

Regulation of cytoplasmic mRNA decay

Discoveries made over the past 20 years highlight the importance of mRNA decay as a means of modulating gene expression and thereby protein production. Up until recently, studies largely focused on identifying cis-acting sequences that serve as mRNA stability or instability elements, the proteins that bind these elements, how the process of translation influences mRNA decay and the ribonucleases that catalyse decay. Now, current studies have begun to elucidate how the decay process is regulated. This Review examines our current understanding of how mammalian cell mRNA decay is controlled by different signalling pathways and lays out a framework for future research.

Neuron-specific ELAV/Hu proteins suppress HuR mRNA during neuronal differentiation by alternative polyadenylation

The ubiquitously expressed RNA-binding protein HuR increases the stability and translation of mRNAs encoding growth regulatory proteins that promote proliferation in a variety of cell types. However, the three neuron-specific ELAV/Hu proteins, HuB, HuC and HuD, while binding to the same types of mRNAs, are required instead for neuronal differentiation, and it becomes difficult to reconcile these contrary functions when all four Hu proteins are expressed in the same neuron. HuR mRNA exists as three alternatively polyadenylated variants, a 1.5-kb testes-specific mRNA isoform, a ubiquitous 2.4-kb isoform and a 6.0-kb isoform that we now show is induced during neuronal differentiation and appears to be neuron-specific. This 6.0-kb neuron-specific mRNA isoform is inherently less stable and produces less HuR protein than the ubiquitous 2.4-kb mRNA. Furthermore, we show that neuronal HuB, HuC and HuD, as well as HuR itself, can bind at the 2.4-kb mRNA polyadenylation site, and when overexpressed can affect alternative polyadenylation to generate an extended HuR 3'-UTR that is translationally suppressed. We propose that the regulation of HuR protein expression by alternative polyadenylation allows neurons to post-transcriptionally regulate mRNAs-encoding factors required for proliferation versus differentiation to facilitate neuronal differentiation.

Reconstitution of CF IA from overexpressed subunits reveals stoichiometry and provides insights into molecular topology

Yeast cleavage factor I (CF I) is an essential complex of five proteins that binds signal sequences at the 3' end of yeast mRNA. CF I is required for correct positioning of a larger protein complex, CPF, which contains the catalytic subunits executing mRNA cleavage and polyadenylation. CF I is composed of two parts, CF IA and Hrp1. The CF IA has only four subunits, Rna14, Rna15, Pcf11, and Clp1, but the structural organization has not been fully established. Using biochemical and biophysical methods, we demonstrate that CF IA can be reconstituted from bacterially expressed proteins and that it has 2:2:1:1 stoichiometry of its four proteins, respectively. We also describe mutations that disrupt the dimer interface of Rna14 while preserving the other subunit interactions. On the basis of our results and existing interaction data, we present a topological model for heterohexameric CF IA and its association with RNA and Hrp1.

Signals for pre-mRNA cleavage and polyadenylation

Pre-mRNA cleavage and polyadenylation is an essential step for 3' end formation of almost all protein-coding transcripts in eukaryotes. The reaction, involving cleavage of nascent mRNA followed by addition of a polyadenylate or poly(A) tail, is controlled by cis-acting elements in the pre-mRNA surrounding the cleavage site. Experimental and bioinformatic studies in the past three decades have elucidated conserved and divergent elements across eukaryotes, from yeast to human. Here we review histories and current models of these elements in a broad range of species.

Transcriptional activity regulates alternative cleavage and polyadenylation

Transcriptomic and epigenomic data, as well as reporter and nuclear run-on assays collectively show that transcriptional activity regulates the relative abundance of alternative polyadenylation isoforms, indicating general coupling of 3′ end processing to transcription.

Using RNA-seq and exon array data for a large number of human and mouse tissues and cells, we identified a general correlation between relative expression of alternative polyadenylation (APA) isoforms and gene expression level: short 3′UTR isoforms are relatively more abundant when genes are highly expressed whereas long 3′UTR isoforms are relatively more abundant when genes are lowly expressed.

Using reporter assays with different promoters, we found that induction of transcription leads to more usage of promoter-proximal polyA sites, suggesting modulation of 3′ end processing efficiency by transcriptional activity. Global analysis and reporter-based assays further revealed that regulation of polyA site choice by transcription takes place when genes are regulated under different cell conditions.

Using global and reporter-based nuclear run-on assays, we found that highly expressed genes tend to have more RNA polymerase II pausing at promoter-proximal polyA sites, as compared with lowly expressed genes, supporting the notion that the efficiency of 3′ end processing is coupled to transcriptional activity.

Highly expressed genes have a lower nucleosome level but higher H3K4me3 and H3K36me3 levels at promoter-proximal polyA sites relative to distal ones, as compared with lowly expressed genes, indicating that transcriptional activity impacts 3′ end processing and regulation of APA leaves epigenetic signatures.

Genes containing multiple pre-mRNA cleavage and polyadenylation sites, or polyA sites, express mRNA isoforms with variable 3′ untranslated regions (UTRs). By systematic analysis of human and mouse transcriptomes, we found that short 3′UTR isoforms are relatively more abundant when genes are highly expressed whereas long 3′UTR isoforms are relatively more abundant when genes are lowly expressed. Reporter assays indicated that polyA site choice can be modulated by transcriptional activity through the gene promoter. Using global and reporter-based nuclear run-on assays, we found that RNA polymerase II is more likely to pause at the polyA site of highly expressed genes than that of lowly expressed ones. Moreover, highly expressed genes tend to have a lower level of nucleosome but higher H3K4me3 and H3K36me3 levels at promoter-proximal polyA sites relative to distal ones. Taken together, our results indicate that polyA site usage is generally coupled to transcriptional activity, leading to regulation of alternative polyadenylation by transcription.

Misfolded human tRNA isodecoder binds and neutralizes a 3' UTR-embedded Alu element

Several classes of small noncoding RNAs are key players in cellular metabolism including mRNA decoding, RNA processing, and mRNA stability. Here we show that a tRNA(Asp) isodecoder, corresponding to a human tRNA-derived sequence, binds to an embedded Alu RNA element contained in the 3' UTR of the human aspartyl-tRNA synthetase mRNA. This interaction between two well-known classes of RNA molecules, tRNA and Alu RNA, is driven by an unexpected structural motif and induces a global rearrangement of the 3' UTR. Besides, this 3' UTR contains two functional polyadenylation signals. We propose a model where the tRNA/Alu interaction would modulate the accessibility of the two alternative polyadenylation sites and regulate the stability of the mRNA. This unique regulation mechanism would link gene expression to RNA polymerase III transcription and may have implications in a primate-specific signal pathway.

Differential usage of transcriptional start sites and polyadenylation sites in FMR1 premutation alleles

5′- and 3′-untranslated regions (UTRs) are important regulators of gene expression and play key roles in disease progression and susceptibility. The 5′-UTR of the fragile X mental retardation 1 (FMR1) gene contains a CGG repeat element that is expanded (>200 CGG repeats; full mutation) and methylated in fragile X syndrome (FXS), the most common form of inherited intellectual disability (ID) and known cause of autism. Significant phenotypic involvement has also emerged in some individuals with the premutation (55–200 CGG repeats), including fragile X-associated premature ovarian insufficiency (FXPOI) in females, and the neurodegenerative disorder, fragile X-associated tremor/ataxia syndrome (FXTAS), in older adult carriers. Here, we show that FMR1 mRNA in human and mouse brain is expressed as a combination of multiple isoforms that use alternative transcriptional start sites and different polyadenylation sites. Furthermore, we have identified a novel human transcription start site used in brain but not in lymphoblastoid cells, and have detected FMR1 isoforms generated through the use of both canonical and non-canonical polyadenylation signals. Importantly, in both human and mouse, a specific regulation of the UTRs is observed in brain of FMR1 premutation alleles, suggesting that the transcript variants may play a role in premutation-related pathologies.

Characterization and prediction of mRNA polyadenylation sites in human genes

The accurate identification of potential poly(A) sites has contributed to all many studies with regard to alternative polyadenylation. The aim of this study was the development of a machine-learning methodology that will help to discriminate real polyadenylation signals from randomly occurring signals in genomic sequence. Since previous studies have revealed that RNA secondary structure in certain genes has significant impact, the authors tried to computationally pinpoint common structural patterns around the poly(A) sites and to investigate how RNA secondary structure may influence polyadenylation. This involved an initial study on the impact of RNA structure and it was found using motif search tools that hairpin structures might be important. Thus, it was propose that, in addition to the sequence pattern around poly(A) sites, there exists a widespread structural pattern that is also employed during human mRNA polyadenylation. In this study, the authors present a computational model that uses support vector machines to predict human poly(A) sites. The results show that this predictive model has a comparable performance to the current prediction tool. In addition, it was identified common structural patterns associated with polyadenylation using several motif finding programs and this provides new insight into the role of RNA secondary structure plays in polyadenylation.

Probabilistic analysis of probe reliability in differential gene expression studies with short oligonucleotide arrays

Probe defects are a major source of noise in gene expression studies. While existing approaches detect noisy probes based on external information such as genomic alignments, we introduce and validate a targeted probabilistic method for analyzing probe reliability directly from expression data and independently of the noise source. This provides insights into the various sources of probe-level noise and gives tools to guide probe design.

POLYAR, a new computer program for prediction of poly(A) sites in human sequences

BACKGROUND

mRNA polyadenylation is an essential step of pre-mRNA processing in eukaryotes. Accurate prediction of the pre-mRNA 3'-end cleavage/polyadenylation sites is important for defining the gene boundaries and understanding gene expression mechanisms.

RESULTS

28761 human mapped poly(A) sites have been classified into three classes containing different known forms of polyadenylation signal (PAS) or none of them (PAS-strong, PAS-weak and PAS-less, respectively) and a new computer program POLYAR for the prediction of poly(A) sites of each class was developed. In comparison with polya_svm (till date the most accurate computer program for prediction of poly(A) sites) while searching for PAS-strong poly(A) sites in human sequences, POLYAR had a significantly higher prediction sensitivity (80.8% versus 65.7%) and specificity (66.4% versus 51.7%) However, when a similar sort of search was conducted for PAS-weak and PAS-less poly(A) sites, both programs had a very low prediction accuracy, which indicates that our knowledge about factors involved in the determination of the poly(A) sites is not sufficient to identify such polyadenylation regions.

CONCLUSIONS

We present a new classification of polyadenylation sites into three classes and a novel computer program POLYAR for prediction of poly(A) sites/regions of each of the class. In tests, POLYAR shows high accuracy of prediction of the PAS-strong poly(A) sites, though this program's efficiency in searching for PAS-weak and PAS-less poly(A) sites is not very high but is comparable to other available programs. These findings suggest that additional characteristics of such poly(A) sites remain to be elucidated. POLYAR program with a stand-alone version for downloading is available at http://cub.comsats.edu.pk/polyapredict.htm.

Genomic analysis of expressed sequence tags in American black bear Ursus americanus

Background

Species of the bear family (Ursidae) are important organisms for research in molecular evolution, comparative physiology and conservation biology, but relatively little genetic sequence information is available for this group. Here we report the development and analyses of the first large scale Expressed Sequence Tag (EST) resource for the American black bear (Ursus americanus).

Results

Comprehensive analyses of molecular functions, alternative splicing, and tissue-specific expression of 38,757 black bear EST sequences were conducted using the dog genome as a reference. We identified 18 genes, involved in functions such as lipid catabolism, cell cycle, and vesicle-mediated transport, that are showing rapid evolution in the bear lineage Three genes, Phospholamban (PLN), cysteine glycine-rich protein 3 (CSRP3) and Troponin I type 3 (TNNI3), are related to heart contraction, and defects in these genes in humans lead to heart disease. Two genes, biphenyl hydrolase-like (BPHL) and CSRP3, contain positively selected sites in bear. Global analysis of evolution rates of hibernation-related genes in bear showed that they are largely conserved and slowly evolving genes, rather than novel and fast-evolving genes.

Conclusion

We provide a genomic resource for an important mammalian organism and our study sheds new light on the possible functions and evolution of bear genes.

Heterogeneity in mammalian RNA 3' end formation

Precisely directed cleavage and polyadenylation of mRNA is a fundamental part of eukaryotic gene expression. Yet, 3' end heterogeneity has been documented for thousands of mammalian genes, and usage of one cleavage and polyadenylation signal over another has been shown to impact gene expression in many cases. Building upon the rich biochemical and genetic understanding of the 3' end formation, recent genomic studies have begun to suggest that widespread changes in mRNA cleavage and polyadenylation may be a part of large, dynamic gene regulatory programs. In this review, we begin with a modest overview of the studies that defined the mechanisms of mammalian 3' end formation, and then discuss how recent genomic studies intersect with these more traditional approaches, showing that both will be crucial for expanding our understanding of this facet of gene regulation.

An unexpected ending: noncanonical 3' end processing mechanisms

Proper 3' end processing of a nascent transcript is critical for the functionality of the mature RNA. Although it has long been thought that virtually all long RNA polymerase II transcripts terminate in a poly(A) tail that is generated by endonucleolytic cleavage followed by polyadenylation, noncanonical 3' end processing mechanisms have recently been identified at several gene loci. Unexpectedly, enzymes with well-characterized roles in other RNA processing events, such as tRNA biogenesis and pre-mRNA splicing, cleave these nascent transcripts to generate their mature 3' ends despite the presence of nearby polyadenylation signals. In fact, the presence of multiple potential 3' end cleavage sites is the norm at many human genes, and recent work suggests that the choice among sites is regulated during development and in response to cellular cues. It is, therefore, becoming increasing clear that the selection of a proper 3' end cleavage site represents an important step in the regulation of gene expression and that the mature 3' ends of RNA polymerase II transcripts can be generated via multiple mechanisms.

Genetic and molecular analyses of PEG10 reveal new aspects of genomic organization, transcription and translation

The paternally expressed gene PEG10 is a retrotransposon derived gene adapted through mammalian evolution located on human chromosome 7q21. PEG10 codes for at least two proteins, PEG10-RF1 and PEG10-RF1/2, by -1 frameshift translation. Overexpression or reinduced PEG10 expression was seen in malignancies, like hepatocellular carcinoma or B-cell acute and chronic lymphocytic leukemia. PEG10 was also shown to promote adipocyte differentiation. Experimental evidence suggests that the PEG10-RF1 protein is an inhibitor of apoptosis and mediates cell proliferation. Here we present new data on the genomic organization of PEG10 by identifying the major transcription start site, a new splice variant and report the cloning and analysis of 1.9 kb of the PEG10 promoter. Furthermore, we show for the first time that PEG10 translation is initiated at a non-AUG start codon upstream of the previously predicted AUG codon as well as at the AUG codon. The finding that PEG10 translation is initiated at different sides adds a new aspect to the already interesting feature of PEG10's -1 frameshift translation mechanism. It is now important to unravel the cellular functions of the PEG10 protein variants and how they are related to normal or pathological conditions. The generated promoter-reporter constructs can be used for future studies to investigate how PEG10 expression is regulated. In summary, our study provides new data on the genomic organization as well as expression and translation of PEG10, a prerequisite in order to study and understand the role of PEG10 in cancer, embryonic development and normal cell homeostasis.

Molecular mechanisms of eukaryotic pre-mRNA 3' end processing regulation

Messenger RNA (mRNA) 3′ end formation is a nuclear process through which all eukaryotic primary transcripts are endonucleolytically cleaved and most of them acquire a poly(A) tail. This process, which consists in the recognition of defined poly(A) signals of the pre-mRNAs by a large cleavage/polyadenylation machinery, plays a critical role in gene expression. Indeed, the poly(A) tail of a mature mRNA is essential for its functions, including stability, translocation to the cytoplasm and translation. In addition, this process serves as a bridge in the network connecting the different transcription, capping, splicing and export machineries. It also participates in the quantitative and qualitative regulation of gene expression in a variety of biological processes through the selection of single or alternative poly(A) signals in transcription units. A large number of protein factors associates with this machinery to regulate the efficiency and specificity of this process and to mediate its interaction with other nuclear events. Here, we review the eukaryotic 3′ end processing machineries as well as the comprehensive set of regulatory factors and discuss the different molecular mechanisms of 3′ end processing regulation by proposing several overlapping models of regulation.

Reprogramming of 3' untranslated regions of mRNAs by alternative polyadenylation in generation of pluripotent stem cells from different cell types

BACKGROUND

The 3' untranslated regions (3'UTRs) of mRNAs contain cis elements involved in post-transcriptional regulation of gene expression. Over half of all mammalian genes contain multiple polyadenylation sites that lead to different 3'UTRs for a gene. Studies have shown that the alternative polyadenylation (APA) pattern varies across tissues, and is dynamically regulated in proliferating or differentiating cells. Generation of induced pluripotent stem (iPS) cells, in which differentiated cells are reprogrammed to an embryonic stem (ES) cell-like state, has been intensively studied in recent years. However, it is not known how 3'UTRs are regulated during cell reprogramming.

METHODS/MAIN FINDINGS

Using a computational method that robustly examines APA across DNA microarray data sets, we analyzed 3'UTR dynamics in generation of iPS cells from different cell types. We found that 3'UTRs shorten during reprogramming of somatic cells, the extent of which depends on the type of source cell. By contrast, reprogramming of spermatogonial cells involves 3'UTR lengthening. The alternative polyadenylation sites that are highly responsive to change of cell state in generation of iPS cells are also highly regulated during embryonic development in opposite directions. Compared with other sites, they are more conserved, can lead to longer alternative 3'UTRs, and are associated with more cis elements for polyadenylation. Consistently, reprogramming of somatic cells and germ cells involves significant upregulation and downregulation, respectively, of mRNAs encoding polyadenylation factors, and RNA processing is one of the most significantly regulated biological processes during cell reprogramming. Furthermore, genes containing target sites of ES cell-specific microRNAs (miRNAs) in different portions of 3'UTR are distinctively regulated during cell reprogramming, suggesting impact of APA on miRNA targeting.

CONCLUSIONS/SIGNIFICANCE

Taken together, these findings indicate that reprogramming of 3'UTRs by APA, which result from regulation of both general polyadenylation activity and cell type-specific factors and can reset post-transcriptional gene regulatory programs in the cell, is an integral part of iPS cell generation, and the APA pattern can be a good biomarker for cell type and state, useful for sample classification. The results also suggest that perturbation of the mRNA polyadenylation machinery or RNA processing activity may facilitate generation of iPS cells.

Discovery of novel human transcript variants by analysis of intronic single-block EST with polyadenylation site

BACKGROUND

Alternative polyadenylation sites within a gene can lead to alternative transcript variants. Although bioinformatic analysis has been conducted to detect polyadenylation sites using nucleic acid sequences (EST/mRNA) in the public databases, one special type, single-block EST is much less emphasized. This bias leaves a large space to discover novel transcript variants.

RESULTS

In the present study, we identified novel transcript variants in the human genome by detecting intronic polyadenylation sites. Poly(A/T)-tailed ESTs were obtained from single-block ESTs and clustered into 10,844 groups standing for 5,670 genes. Most sites were not found in other alternative splicing databases. To verify that these sites are from expressed transcripts, we analyzed the supporting EST number of each site, blasted representative ESTs against known mRNA sequences, traced terminal sequences from cDNA clones, and compared with the data of Affymetrix tiling array. These analyses confirmed about 84% (9,118/10,844) of the novel alternative transcripts, especially, 33% (3,575/10,844) of the transcripts from 2,704 genes were taken as high-reliability. Additionally, RT-PCR confirmed 38% (10/26) of predicted novel transcript variants.

CONCLUSION

Our results provide evidence for novel transcript variants with intronic poly(A) sites. The expression of these novel variants was confirmed with computational and experimental tools. Our data provide a genome-wide resource for identification of novel human transcript variants with intronic polyadenylation sites, and offer a new view into the mystery of the human transcriptome.

Polypyrimidine tract binding protein prevents activity of an intronic regulatory element that promotes usage of a composite 3'-terminal exon

Alternative splicing of 3'-terminal exons plays a critical role in gene expression by producing mRNA with distinct 3'-untranslated regions that regulate their fate and their expression. The Xenopus alpha-tropomyosin pre-mRNA possesses a composite internal/3'-terminal exon (exon 9A9') that is differentially processed depending on the embryonic tissue. Exon 9A9' is repressed in non-muscle tissue by the polypyrimidine tract binding protein, whereas it is selected as a 3'-terminal or internal exon in myotomal cells and adult striated muscles, respectively. We report here the identification of an intronic regulatory element, designated the upstream terminal exon enhancer (UTE), that is required for the specific usage of exon 9A9' as a 3'-terminal exon in the myotome. We demonstrate that polypyrimidine tract binding protein prevents the activity of UTE in non-muscle cells, whereas a subclass of serine/arginine rich (SR) proteins promotes the selection of exon 9A9' in a UTE-dependent way. Morpholino-targeted blocking of UTE in the embryo strongly reduced the inclusion of exon 9A9' as a 3'-terminal exon in the endogenous mRNA, demonstrating the function of UTE under physiological circumstances. This strategy allowed us to reveal a splicing pathway that generates a mRNA with no in frame stop codon and whose steady-state level is translation-dependent. This result suggests that a non-stop decay mechanism participates in the strict control of the 3'-end processing of the alpha-tropomyosin pre-mRNA.

Mechanistic insights into the link between a polymorphism of the 3'UTR of the SLC7A1 gene and hypertension

We previously identified the polymorphism ss52051869 in the 3'UTR of human SLC7A1, and demonstrated that it might participate in the apparent link between altered endothelial function, decreased L-arginine and nitric oxide (NO) metabolism, and a genetic predisposition to essential hypertension. Here, we demonstrate that the major allele contains a consensus sequence for the transcription factor SP1 and binds to SP1, in contrast, the minor allele fails to bind to SP1. Resequencing of the entire SLC7A1 coding sequence failed to find other informative polymorphisms, indicating that ss52051869 plays a key role in the biochemical and clinical association. In conjunction, the short and long variants of the 3'UTR of SLC7A1 contain three and four potential microRNA-122 (miR-122) binding sites, respectively. We found that the minor allele is more frequently associated with SLC7A1 bearing a long 3'UTR, while the major allele is more likely to accompany a short 3'UTR only (P=0.034). As such, reporter genes containing the long 3'UTR from SLC7A1 show much less gene expression than those containing short 3'UTR from SLC7A1, regardless of their allele status (P<0.01), suggesting that an alternative polyadenylation event and/or miRNA-122 binding sites may also play a role in controlling gene expression. It is therefore possible that binding of miR-122 to the 3'UTR may cause the depression of gene expression, contributing to the lesser level of SLC7A1 and the endothelial dysfunction seen in hypertensive subjects. Taken together, these data provide novel insights into the mechanism by which ss52051869 influences SLC7A1 gene expression.

Alternative polyadenylation: a mechanism maximizing transcriptome diversity in higher eukaryotes

Based on comparative genome analyses, the increases in protein-coding gene number could not account for the increases of morphological and behavioral complexity of higher eukaryotes. Transcriptional regulations, alternative splicing and the involvement of non-coding RNA in gene expression regulations have been credited for the drastic increase of transcriptome complexity. However, an emerging theme of another mechanism that contributes to the formation of alternative mRNA 3'-ends is alternative polyadenylation (APA). First, recent studies indicated that APA is a wide spread phenomenon across the transcriptomes of higher eukaryotes and being regulated by developmental and environmental cues. Secondly, our characterization of the Arabidopsis polyadenylation factors suggested that plant polyadenylation has also evolved to regulate the expression of specific genes by means of APA and therefore the specific biological functions. Finally, Phylogenetic analyses of eukaryotic polyadenylation factors from several organisms revealed that the number of polyadenylation factors tends to increase in higher eukaryotes, which provides the potential for their functional differentiation in regulating gene expression through APA. Based on above evidence, we, thus, hypothesize that APA, serving as an additional mechanism, contributes to the complexity of higher eukaryotes.

Identification of gene 3' ends by automated EST cluster analysis

The properties and biology of mRNA transcripts can be affected profoundly by the choice of alternative polyadenylation sites, making definition of the 3' ends of transcripts essential for understanding their regulation. Here we show that 22-52% of sequences in commonly used human and murine "full-length" transcript databases may not currently end at bona fide polyadenylation sites. To identify probable transcript termini over the entire murine and human genomes, we analyzed the EST databases for positional clustering of EST ends. The analysis yielded 58,282 murine- and 86,410 human-candidate polyadenylation sites, of which 75% mapped to 23,091 known murine transcripts and 22,891 known human transcripts. The murine dataset correctly predicted 97% of the 3' ends in a manually curated and experimentally supported benchmark transcript set. Of currently known genes, 15% had no associated prediction and 25% had only a single predicted termination site. The remaining genes had an average of 3-4 alternative polyadenylation sites predicted for each murine or human transcript, respectively. The results are made available in the form of tables and an interactive web site that can be mined for rapid assessment of the validity of 3' ends in existing collections, enumeration of potential alternative 3' polyadenylation sites of known transcripts, direct retrieval of terminal sequences for design of probes, and detection of polyadenylation sites not currently mapped to known genes.

Species-specific cis-regulatory elements in the 3'-untranslated region direct alternative polyadenylation of bone morphogenetic protein 2 mRNA

BMP2 (bone morphogenetic protein 2) is a multifunctional member of the transforming growth factor-beta family of growth factors. Disruption of BMP2 signaling results in developmental defects, cancers, and other diseases. BMP2 mRNAs are alternatively polyadenylated, resulting in mRNAs with distinct 3'-untranslated regions. The longer mRNA contains additional putative binding sites for post-transcriptional regulatory factors, including micro-RNAs. We combined functional assays with computational analyses of emerging genome data to define site- and species-specific polyadenylation determinants. In all mouse and human cell lines tested, shorter mRNAs resulting from using the first polyadenylation signal (PA1) were more abundant than mRNAs from the second signal (PA2). However, the PA1/PA2 usage ratios were 2-3-fold higher in human than in mouse cells. Expression of human BMP2 constructs in mouse cells and mouse constructs in human cells showed that cis-regulatory elements direct species-specific 3' processing of BMP2 transcripts. A 72-nucleotide region downstream of PA2 in the mouse sequence contains two novel cis-acting elements previously hypothesized to regulate polyadenylation in a bioinformatics analysis. Mutations that humanized the mouse-specific elements lowered the affinity for cleavage stimulation factor CstF64 and significantly weakened the PA2 signal relative to the PA1 signal. Thus, we have experimentally defined for the first time cis-regulatory elements that control a species-specific difference in the 3'-end processing of BMP2 and potentially of other genes.

MicroRNA-mediated up-regulation of an alternatively polyadenylated variant of the mouse cytoplasmic {beta}-actin gene

Actin is a major cytoskeletal protein in eukaryotes. Recent studies suggest more diverse functional roles for this protein. Actin mRNA is known to be localized to neuronal synapses and undergoes rapid deadenylation during early developmental stages. However, its 3′-untranslated region (UTR) is not characterized and there are no experimentally determined polyadenylation (polyA) sites in actin mRNA. We have found that the cytoplasmic β-actin (Actb) gene generates two alternative transcripts terminated at tandem polyA sites. We used 3′-RACE, EST end analysis and in situ hybridization to unambiguously establish the existence of two 3′-UTRs of varying length in Actb transcript in mouse neuronal cells. Further analyses showed that these two tandem polyA sites are used in a tissue-specific manner. Although the longer 3′-UTR was expressed at a relatively lower level, it conferred higher translational efficiency to the transcript. The longer transcript harbours a conserved mmu-miR-34a/34b-5p target site. Sequence-specific anti-miRNA molecule, mutations of the miRNA target region in the 3′-UTR resulted in reduced expression. The expression was restored by a mutant miRNA complementary to the mutated target region implying that miR-34 binding to Actb 3′-UTR up-regulates target gene expression. Heterogeneity of the Actb 3′-UTR could shed light on the mechanism of miRNA-mediated regulation of messages in neuronal cells.

Phylogenetic analysis of mRNA polyadenylation sites reveals a role of transposable elements in evolution of the 3'-end of genes

mRNA polyadenylation is an essential step for the maturation of almost all eukaryotic mRNAs, and is tightly coupled with termination of transcription in defining the 3′-end of genes. Large numbers of human and mouse genes harbor alternative polyadenylation sites [poly(A) sites] that lead to mRNA variants containing different 3′-untranslated regions (UTRs) and/or encoding distinct protein sequences. Here, we examined the conservation and divergence of different types of alternative poly(A) sites across human, mouse, rat and chicken. We found that the 3′-most poly(A) sites tend to be more conserved than upstream ones, whereas poly(A) sites located upstream of the 3′-most exon, also termed intronic poly(A) sites, tend to be much less conserved. Genes with longer evolutionary history are more likely to have alternative polyadenylation, suggesting gain of poly(A) sites through evolution. We also found that nonconserved poly(A) sites are associated with transposable elements (TEs) to a much greater extent than conserved ones, albeit less frequently utilized. Different classes of TEs have different characteristics in their association with poly(A) sites via exaptation of TE sequences into polyadenylation elements. Our results establish a conservation pattern for alternative poly(A) sites in several vertebrate species, and indicate that the 3′-end of genes can be dynamically modified by TEs through evolution.