Differential genome-wide profiling of tandem 3' UTRs among human breast cancer and normal cells by high-throughput sequencing

Dynamic landscape of tandem 3' UTRs during zebrafish development

Tandem 3' untranslated regions (UTRs), produced by alternative polyadenylation (APA) in the terminal exon of a gene, could have critical roles in regulating gene networks. Here we profiled tandem poly(A) events on a genome-wide scale during the embryonic development of zebrafish (Danio rerio) using a recently developed SAPAS method. We showed that 43% of the expressed protein-coding genes have tandem 3' UTRs. The average 3' UTR length follows a V-shaped dynamic pattern during early embryogenesis, in which the 3' UTRs are first shortened at zygotic genome activation, and then quickly lengthened during gastrulation. Over 4000 genes are found to switch tandem APA sites, and the distinct functional roles of these genes are indicated by Gene Ontology analysis. Three families of cis-elements, including miR-430 seed, U-rich element, and canonical poly(A) signal, are enriched in 3' UTR-shortened/lengthened genes in a stage-specific manner, suggesting temporal regulation coordinated by APA and trans-acting factors. Our results highlight the regulatory role of tandem 3' UTR control in early embryogenesis and suggest that APA may represent a new epigenetic paradigm of physiological regulations.

Estrogen-induced upregulation and 3'-UTR shortening of CDC6

3′-Untranslated region (UTR) shortening of mRNAs via alternative polyadenylation (APA) has important ramifications for gene expression. By using proximal APA sites and switching to shorter 3′-UTRs, proliferating cells avoid miRNA-mediated repression. Such APA and 3′-UTR shortening events may explain the basis of some of the proto-oncogene activation cases observed in cancer cells. In this study, we investigated whether 17 β-estradiol (E2), a potent proliferation signal, induces APA and 3′-UTR shortening to activate proto-oncogenes in estrogen receptor positive (ER+) breast cancers. Our initial probe based screen of independent expression arrays suggested upregulation and 3′-UTR shortening of an essential regulator of DNA replication, CDC6 (cell division cycle 6), upon E2 treatment. We further confirmed the E2- and ER-dependent upregulation and 3′UTR shortening of CDC6, which lead to increased CDC6 protein levels and higher BrdU incorporation. Consequently, miRNA binding predictions and dual luciferase assays suggested that 3′-UTR shortening of CDC6 was a mechanism to avoid 3′-UTR-dependent negative regulations. Hence, we demonstrated CDC6 APA induction by the proliferative effect of E2 in ER+ cells and provided new insights into the complex regulation of APA. E2-induced APA is likely to be an important but previously overlooked mechanism of E2-responsive gene expression.

Landscape of transcription in human cells

Eukaryotic cells make many types of primary and processed RNAs that are found either in specific subcellular compartments or throughout the cells. A complete catalogue of these RNAs is not yet available and their characteristic subcellular localizations are also poorly understood. Because RNA represents the direct output of the genetic information encoded by genomes and a significant proportion of a cell's regulatory capabilities are focused on its synthesis, processing, transport, modification and translation, the generation of such a catalogue is crucial for understanding genome function. Here we report evidence that three-quarters of the human genome is capable of being transcribed, as well as observations about the range and levels of expression, localization, processing fates, regulatory regions and modifications of almost all currently annotated and thousands of previously unannotated RNAs. These observations, taken together, prompt a redefinition of the concept of a gene.

An in-depth map of polyadenylation sites in cancer

We present a comprehensive map of over 1 million polyadenylation sites and quantify their usage in major cancers and tumor cell lines using direct RNA sequencing. We built the Expression and Polyadenylation Database to enable the visualization of the polyadenylation maps in various cancers and to facilitate the discovery of novel genes and gene isoforms that are potentially important to tumorigenesis. Analyses of polyadenylation sites indicate that a large fraction (∼30%) of mRNAs contain alternative polyadenylation sites in their 3' untranslated regions, independent of the cell type. The shortest 3' untranslated region isoforms are preferentially upregulated in cancer tissues, genome-wide. Candidate targets of alternative polyadenylation-mediated upregulation of short isoforms include POLR2K, and signaling cascades of cell-cell and cell-extracellular matrix contact, particularly involving regulators of Rho GTPases. Polyadenylation maps also helped to improve 3' untranslated region annotations and identify candidate regulatory marks such as sequence motifs, H3K36Me3 and Pabpc1 that are isoform dependent and occur in a position-specific manner. In summary, these results highlight the need to go beyond monitoring only the cumulative transcript levels for a gene, to separately analysing the expression of its RNA isoforms.

Genetics and regulatory impact of alternative polyadenylation in human B-lymphoblastoid cells

Gene expression varies widely between individuals of a population, and regulatory change can underlie phenotypes of evolutionary and biomedical relevance. A key question in the field is how DNA sequence variants impact gene expression, with most mechanistic studies to date focused on the effects of genetic change on regulatory regions upstream of protein-coding sequence. By contrast, the role of RNA 3'-end processing in regulatory variation remains largely unknown, owing in part to the challenge of identifying functional elements in 3' untranslated regions. In this work, we conducted a genomic survey of transcript ends in lymphoblastoid cells from genetically distinct human individuals. Our analysis mapped the cis-regulatory architecture of 3' gene ends, finding that transcript end positions did not fall randomly in untranslated regions, but rather preferentially flanked the locations of 3' regulatory elements, including miRNA sites. The usage of these transcript length forms and motifs varied across human individuals, and polymorphisms in polyadenylation signals and other 3' motifs were significant predictors of expression levels of the genes in which they lay. Independent single-gene experiments confirmed the effects of polyadenylation variants on steady-state expression of their respective genes, and validated the regulatory function of 3' cis-regulatory sequence elements that mediated expression of these distinct RNA length forms. Focusing on the immune regulator IRF5, we established the effect of natural variation in RNA 3'-end processing on regulatory response to antigen stimulation. Our results underscore the importance of two mechanisms at play in the genetics of 3'-end variation: the usage of distinct 3'-end processing signals and the effects of 3' sequence elements that determine transcript fate. Our findings suggest that the strategy of integrating observed 3'-end positions with inferred 3' regulatory motifs will prove to be a critical tool in continued efforts to interpret human genome variation.

Alternative cleavage and polyadenylation during colorectal cancer development

PURPOSE

Alternative cleavage and polyadenylation (APA) of mRNAs is a phenomenon that alters 3'-untranslated region length leading to altered posttranscriptional regulation of gene expression. Changing APA patterns have been shown to result in misregulation of genes involved in carcinogenesis; therefore, we hypothesized that altered APA contributes to progression of colorectal cancer, and that measurement of APA may lead to discovery of novel biomarkers.

EXPERIMENTAL DESIGN

We used next-generation sequencing to directly measure global patterns of APA changes during colorectal carcinoma progression in 15 human patient samples. Results were validated in a larger cohort of 50 patients, including 5 normal/carcinoma pairs from individuals.

RESULTS

We discovered numerous genes presenting progressive changes in APA. Genes undergoing untranslated region (3'UTR) shortening were enriched for functional groups such as cell-cycle and nucleic acid-binding and processing factors, and those undergoing 3'UTR lengthening or alternative 3'UTR usage were enriched for categories such as cell-cell adhesion and extracellular matrix. We found indications that APA changes result from differential processing of transcripts because of increased expression of cleavage and polyadenylation factors. Quantitative PCR analysis in a larger series of human patient samples, including matched pairs, confirmed APA changes in DMKN, PDXK, and PPIE genes.

CONCLUSIONS

Our results suggest that genes undergoing altered APA during human cancer progression may be useful novel biomarkers and potentially targeted for disease prevention and treatment. We propose that the strategy presented here may be broadly useful in discovery of novel biomarkers for other types of cancer and human disease.

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.

E2F mediates enhanced alternative polyadenylation in proliferation

Background

The majority of mammalian genes contain multiple poly(A) sites in their 3' UTRs. Alternative cleavage and polyadenylation are emerging as an important layer of gene regulation as they generate transcript isoforms that differ in their 3' UTRs, thereby modulating genes' response to 3' UTR-mediated regulation. Enhanced cleavage at 3' UTR proximal poly(A) sites resulting in global 3' UTR shortening was recently linked to proliferation and cancer. However, mechanisms that regulate this enhanced alternative polyadenylation are unknown.

Results

Here, we explored, on a transcriptome-wide scale, alternative polyadenylation events associated with cellular proliferation and neoplastic transformation. We applied a deep-sequencing technique for identification and quantification of poly(A) sites to two human cellular models, each examined under proliferative, arrested and transformed states. In both cell systems we observed global 3' UTR shortening associated with proliferation, a link that was markedly stronger than the association with transformation. Furthermore, we found that proliferation is also associated with enhanced cleavage at intronic poly(A) sites. Last, we found that the expression level of the set of genes that encode for 3'-end processing proteins is globally elevated in proliferation, and that E2F transcription factors contribute to this regulation.

Conclusions

Our results comprehensively identify alternative polyadenylation events associated with cellular proliferation and transformation, and demonstrate that the enhanced alternative polyadenylation in proliferative conditions results not only in global 3' UTR shortening but also in enhanced premature cleavage in introns. Our results also indicate that E2F-mediated co-transcriptional regulation of 3'-end processing genes is one of the mechanisms that links enhanced alternative polyadenylation to proliferation.

Analysis of C. elegans intestinal gene expression and polyadenylation by fluorescence-activated nuclei sorting and 3'-end-seq

Despite the many advantages of Caenorhabditis elegans, biochemical approaches to study tissue-specific gene expression in post-embryonic stages are challenging. Here, we report a novel experimental approach for efficient determination of tissue-specific transcriptomes involving the rapid release and purification of nuclei from major tissues of post-embryonic animals by fluorescence-activated nuclei sorting (FANS), followed by deep sequencing of linearly amplified 3'-end regions of transcripts (3'-end-seq). We employed these approaches to compile the transcriptome of the developed C. elegans intestine and used this to analyse tissue-specific cleavage and polyadenylation. In agreement with intestinal-specific gene expression, highly expressed genes have enriched GATA-elements in their promoter regions and their functional properties are associated with processes that are characteristic for the intestine. We systematically mapped pre-mRNA cleavage and polyadenylation sites, or polyA sites, including more than 3000 sites that have previously not been identified. The detailed analysis of the 3'-ends of the nuclear mRNA revealed widespread alternative polyA site use (APA) in intestinally expressed genes. Importantly, we found that intestinal polyA sites that undergo APA tend to have U-rich and/or A-rich upstream auxiliary elements that may contribute to the regulation of 3'-end formation in the intestine.

Perspectives of integrative cancer genomics in next generation sequencing era

The explosive development of genomics technologies including microarrays and next generation sequencing (NGS) has provided comprehensive maps of cancer genomes, including the expression of mRNAs and microRNAs, DNA copy numbers, sequence variations, and epigenetic changes. These genome-wide profiles of the genetic aberrations could reveal the candidates for diagnostic and/or prognostic biomarkers as well as mechanistic insights into tumor development and progression. Recent efforts to establish the huge cancer genome compendium and integrative omics analyses, so-called "integromics", have extended our understanding on the cancer genome, showing its daunting complexity and heterogeneity. However, the challenges of the structured integration, sharing, and interpretation of the big omics data still remain to be resolved. Here, we review several issues raised in cancer omics data analysis, including NGS, focusing particularly on the study design and analysis strategies. This might be helpful to understand the current trends and strategies of the rapidly evolving cancer genomics research.

Inverted Alu dsRNA structures do not affect localization but can alter translation efficiency of human mRNAs independent of RNA editing

With over one million copies, Alu elements are the most abundant repetitive elements in the human genome. When transcribed, interaction between two Alus that are in opposite orientation gives rise to double-stranded RNA (dsRNA). Although the presence of dsRNA in the cell was previously thought to only occur during viral infection, it is now known that cells express many endogenous small dsRNAs, such as short interfering RNA (siRNAs) and microRNA (miRNAs), which regulate gene expression. It is possible that long dsRNA structures formed from Alu elements influence gene expression. Here, we report that human mRNAs containing inverted Alu elements are present in the mammalian cytoplasm. The presence of these long intramolecular dsRNA structures within 3′-UTRs decreases translational efficiency, and although the structures undergo extensive editing in vivo, the effects on translation are independent of the presence of inosine. As inverted Alus are predicted to reside in >5% of human protein-coding genes, these intramolecular dsRNA structures are important regulators of gene expression.

Genome-wide determination of a broad ESRP-regulated posttranscriptional network by high-throughput sequencing

Tissue-specific alternative splicing is achieved through the coordinated assembly of RNA binding proteins at specific sites to enhance or silence splicing at nearby splice sites. We used high-throughput sequencing (RNA-Seq) to investigate the complete spectrum of alternative splicing events that are regulated by the epithelium-specific splicing regulatory proteins ESRP1 and ESRP2. We also combined this analysis with direct RNA sequencing (DRS) to reveal ESRP-mediated regulation of alternative polyadenylation. To define binding motifs that mediate direct regulation of splicing and polyadenylation by ESRP, SELEX-Seq analysis was performed, coupling traditional SELEX with high-throughput sequencing. Identification and scoring of high-affinity ESRP1 binding motifs within ESRP target genes allowed the generation of RNA maps that define the position-dependent activity of the ESRPs in regulating cassette exons and alternative 3' ends. These extensive analyses provide a comprehensive picture of the functions of the ESRPs in an epithelial posttranscriptional gene expression program.

Global microRNA level regulation of EGFR-driven cell-cycle protein network in breast cancer

A genome-wide microRNA (miRNome) screen coupled with high-throughput monitoring of protein levels reveals complex, modular miRNA regulation of the EGFR-driven cell-cycle network, and identifies new miRNAs that can suppress breast cancer cell proliferation.

We interrogated, for the first time, a mammalian oncogenic signaling network with the miRNome and report the outputs at the protein level.

Whole-genome microRNA (miRNA) effects on a given protein are generally mild, supporting a fine-tuning role for miRNAs, and these effects are dominated by sequence-matching mechanisms.

We developed a novel network-analysis methodology with a bipartite graph model to identify proteins co-regulated by miRNAs. Besides the sequence-based mechanism, our results demonstrated that miRNAs simultaneously regulate several proteins belonging to the same functional module.

We identified three miRNAs, miR-124, miR-147 and miR-193a-3p, as novel tumor suppressors that co-regulate EGFR-driven cell-cycle network proteins, and inhibit cell-cycle progression and proliferation in breast cancer.

Our results demonstrate the potential to steer miRNA research toward the network level, underlining the need for systematic approaches before positioning miRNAs as drugs or drug targets.

The EGFR-driven cell-cycle pathway has been extensively studied due to its pivotal role in breast cancer proliferation and pathogenesis. Although several studies reported regulation of individual pathway components by microRNAs (miRNAs), little is known about how miRNAs coordinate the EGFR protein network on a global miRNA (miRNome) level. Here, we combined a large-scale miRNA screening approach with a high-throughput proteomic readout and network-based data analysis to identify which miRNAs are involved, and to uncover potential regulatory patterns. Our results indicated that the regulation of proteins by miRNAs is dominated by the nucleotide matching mechanism between seed sequences of the miRNAs and 3′-UTR of target genes. Furthermore, the novel network-analysis methodology we developed implied the existence of consistent intrinsic regulatory patterns where miRNAs simultaneously co-regulate several proteins acting in the same functional module. Finally, our approach led us to identify and validate three miRNAs (miR-124, miR-147 and miR-193a-3p) as novel tumor suppressors that co-target EGFR-driven cell-cycle network proteins and inhibit cell-cycle progression and proliferation in breast cancer.

Shortening of 3'UTRs correlates with poor prognosis in breast and lung cancer

A major part of the post-transcriptional regulation of gene expression is affected by trans-acting elements, such as microRNAs, binding the 3' untraslated region (UTR) of their target mRNAs. Proliferating cells partly escape this type of negative regulation by expressing shorter 3' UTRs, depleted of microRNA binding sites, compared to non-proliferating cells. Using large-scale gene expression datasets, we show that a similar phenomenon takes place in breast and lung cancer: tumors expressing shorter 3' UTRs tend to be more aggressive and to result in shorter patient survival. Moreover, we show that a gene expression signature based only on the expression ratio of alternative 3' UTRs is a strong predictor of survival in both tumors. Genes undergoing 3'UTR shortening in aggressive tumors of the two tissues significantly overlap, and several of them are known to be involved in tumor progression. However the pattern of 3' UTR shortening in aggressive tumors in vivo is clearly distinct from analogous patterns involved in proliferation and transformation.

Brain-derived neurotrophic factor expression is repressed during myogenic differentiation by miR-206

Brain-derived neurotrophic factor (BDNF) is required for efficient skeletal-muscle regeneration and perturbing its expression causes abnormalities in the proliferation and differentiation of skeletal muscle cells. In this study, we investigated the mechanism of BDNF suppression that occurs during myogenic differentiation. BDNF is expressed at the mRNA level as two isoforms that differ in the length of their 3'UTRs as a result of alternative cleavage and polyadenylation. Sequence analysis revealed the presence of three miR-206 target sites in the long BDNF 3'UTR (BDNF-L), whereas only one site was found in the short mRNA BDNF 3'UTR (BDNF-S). miR-206 is known to regulate the differentiation of C2C12 myoblasts and its expression is induced during the transition from myoblasts to myotubes. We thus examined whether miR-206-mediated suppression is responsible for the expression pattern of BDNF during myogenic differentiation. BDNF-L was suppressed to a greater extent than BDNF-S during differentiation of C2C12 myoblasts. Transfection of a miR-206 precursor decreased activity of reporters representative of the BDNF-L 3'UTR, but not BDNF-S 3'UTR, and repressed endogenous BDNF mRNA levels. This suppression was found to be dependent on the presence of multiple miR-206 target sites in the BDNF-L 3'UTR. Conversely, suppression of miR-206 levels resulted in de-repression of BDNF 3'UTR reporter activity and increased endogenous BDNF-L mRNA levels. A receptor for BDNF, p75(NTR) , was also suppressed during differentiation and in response to miR-206, but this appeared to not be entirely mediated via a miR-206 target site its 3'UTR. Based on these observations, BDNF represents a novel target through which miR-206 controls the initiation and maintenance of the differentiated state of muscle cells. These results further suggest that miR-206 might play a role in regulating retrograde signaling of BDNF at the neuromuscular junction.

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.

Ending the message: poly(A) signals then and now

Polyadenylation [poly(A)] signals (PAS) are a defining feature of eukaryotic protein-coding genes. The central sequence motif AAUAAA was identified in the mid-1970s and subsequently shown to require flanking, auxiliary elements for both 3'-end cleavage and polyadenylation of premessenger RNA (pre-mRNA) as well as to promote downstream transcriptional termination. More recent genomic analysis has established the generality of the PAS for eukaryotic mRNA. Evidence for the mechanism of mRNA 3'-end formation is outlined, as is the way this RNA processing reaction communicates with RNA polymerase II to terminate transcription. The widespread phenomenon of alternative poly(A) site usage and how this interrelates with pre-mRNA splicing is then reviewed. This shows that gene expression can be drastically affected by how the message is ended. A central theme of this review is that while genomic analysis provides generality for the importance of PAS selection, detailed mechanistic understanding still requires the direct analysis of specific genes by genetic and biochemical approaches.

HLA-DRB1 may be antagonistically regulated by the coordinately evolved promoter and 3'-UTR under stabilizing selection

HLA-DRB1 is the most polymorphic MHC (major histocompatibility complex) class II gene in human, and plays a crucial role in the development and function of the immune system. Extensive polymorphisms exist in the promoter and 3'-UTR of HLA-DRB1, especially a LTR (Long terminal repeat) element in the promoter, which may be involved in the expression regulation. However, it remains unknown how the polymorphisms in the whole promoter region and 3'-UTR to regulate the gene expression. In this study, we investigated the extensive polymorphisms in the HLA-DRB1 promoter and 3'-UTR, and how these polymorphisms affect the gene expression in both independent and jointly manners. It was observed that most of the haplotypes in the DRB1 promoter and 3'-UTR were clustered into 4 conserved lineages (H1, H2, H3 and H4), and showed high linkage disequilibrium. Compared with H1 and H2 lineage, a LTR element in the promoter of H3 and H4 lineage significantly suppressed the promoter activity, whereas the activity of the linked 3'-UTR increased, leading to no apparent difference in the final expression product between H1/H2 and H3/H4 lineage. Nevertheless, compared with the plasmid with a promoter and 3'-UTR from the same lineage, the recombinant plasmid with a promoter from H2 and a 3'-UTR from H3 showed about double fold increased luciferase activity, Conversely, the recombinant plasmid with a promoter from H3 and a 3'-UTR from H2 resulted in about 2-fold decreased luciferase activity. These results indicate that the promoter and 3'-UTR of HLA-DRB1 may antagonistically regulate the gene expression, which may be subjected to stabilizing selection. These findings may provide a novel insight into the mechanisms of the diseases associated with HLA-DRB1 genes.