Pervasive isoform-specific translational regulation via alternative transcription start sites in mammals

Selective translational usage of TSS and core promoters revealed by translatome sequencing

Background

In mammals, fine-tuned regulation of gene expression leads to transcription initiation from diverse transcription start sites (TSSs) and multiple core promoters. Although polysome association is a critical step in translation, whether polysome selectively uses TSSs and core promoters and how this could impact translation remains elusive.

Results

In this study, we used CAGE followed by deep sequencing to globally profile the transcript 5′ isoforms in the translatome and transcriptome of human HEK293 cells at single-nucleotide resolution. By comparing the two profiles, we identified the 5′ isoforms preferentially used in translatome and revealed a widespread selective usage of TSSs (32.0%) and core promoters (48.7%) by polysome. We discovered the transcription initiation patterns and the sequence characteristics that were highly correlated with polysome selection. We further identified 5804 genes significantly enriched or depleted in translatome and showed that polysome selection was an important contributing factor to the abundance of related gene products. Moreover, after comparison with public transcriptome CAGE data from 180 human tissues and primary cells, we raised a question on whether it is a widely adopted mechanism to regulate translation efficiency by changing the transcription initiation sites on the transcription level in cells of different conditions.

Conclusions

Using HEK293 cells as a model, we delineated an indirect selection toward TSSs and core promoters by the translation machinery. Our findings lend additional evidence for a much closer coordination between transcription and translation, warranting future translatome studies in more cell types and conditions to develop a more intricate regulatory model for gene expression.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5650-0) contains supplementary material, which is available to authorized users.

Evidence for an Integrated Gene Repression Mechanism Based on mRNA Isoform Toggling in Human Cells

We recently described an unconventional mode of gene regulation in budding yeast by which transcriptional and translational interference collaborate to down-regulate protein expression. Developmentally timed transcriptional interference inhibited production of a well translated mRNA isoform and resulted in the production of an mRNA isoform containing inhibitory upstream open reading frames (uORFs) that prevented translation of the main ORF. Transcriptional interference and uORF-based translational repression are established mechanisms outside of yeast, but whether this type of integrated regulation was conserved was unknown. Here we find that, indeed, a similar type of regulation occurs at the locus for the human oncogene MDM2 We observe evidence of transcriptional interference between the two MDM2 promoters, which produce a poorly translated distal promoter-derived uORF-containing mRNA isoform and a well-translated proximal promoter-derived transcript. Down-regulation of distal promoter activity markedly up-regulates proximal promoter-driven expression and results in local reduction of histone H3K36 trimethylation. Moreover, we observe that this transcript toggling between the two MDM2 isoforms naturally occurs during human embryonic stem cell differentiation programs.

Pervasive, Coordinated Protein-Level Changes Driven by Transcript Isoform Switching during Meiosis

To better understand the gene regulatory mechanisms that program developmental processes, we carried out simultaneous genome-wide measurements of mRNA, translation, and protein through meiotic differentiation in budding yeast. Surprisingly, we observed that the levels of several hundred mRNAs are anti-correlated with their corresponding protein products. We show that rather than arising from canonical forms of gene regulatory control, the regulation of at least 380 such cases, or over 8% of all measured genes, involves temporally regulated switching between production of a canonical, translatable transcript and a 5' extended isoform that is not efficiently translated into protein. By this pervasive mechanism for the modulation of protein levels through a natural developmental program, a single transcription factor can coordinately activate and repress protein synthesis for distinct sets of genes. The distinction is not based on whether or not an mRNA is induced but rather on the type of transcript produced.

A novel thyroid hormone receptor isoform, TRβ2-46, promotes SKP2 expression and retinoblastoma cell proliferation

Retinoblastoma is a childhood retinal tumor that develops from cone photoreceptor precursors in response to inactivating RB1 mutations and loss of functional RB protein. The cone precursor's response to RB loss involves cell type-specific signaling circuitry that helps to drive tumorigenesis. One component of the cone precursor circuitry, the thyroid hormone receptor β2 (TRβ2), enables the aberrant proliferation of diverse RB-deficient cells in part by opposing the down-regulation of S-phase kinase-associated protein 2 (SKP2) by the more widely expressed and tumor-suppressive TRβ1. However, it is unclear how TRβ2 opposes TRβ1 to enable SKP2 expression and cell proliferation. Here, we show that in human retinoblastoma cells TRβ2 mRNA encodes two TRβ2 protein isoforms: a predominantly cytoplasmic 54-kDa protein (TRβ2-54) corresponding to the well-characterized full-length murine Trβ2 and an N-terminally truncated and exclusively cytoplasmic 46-kDa protein (TRβ2-46) that starts at Met-79. Whereas TRβ2 knockdown decreased SKP2 expression and impaired retinoblastoma cell cycle progression, re-expression of TRβ2-46 but not TRβ2-54 stabilized SKP2 and restored proliferation to an extent similar to that of ectopic SKP2 restoration. We conclude that TRβ2-46 is an oncogenic thyroid hormone receptor isoform that promotes SKP2 expression and SKP2-dependent retinoblastoma cell proliferation.

Translatomics: The Global View of Translation

In all kingdoms of life, proteins are synthesized by ribosomes in a process referred to as translation. The amplitude of translational regulation exceeds the sum of transcription, mRNA degradation and protein degradation. Therefore, it is essential to investigate translation in a global scale. Like the other “omics”-methods, translatomics investigates the totality of the components in the translation process, including but not limited to translating mRNAs, ribosomes, tRNAs, regulatory RNAs and nascent polypeptide chains. Technical advances in recent years have brought breakthroughs in the investigation of these components at global scale, both for their composition and dynamics. These methods have been applied in a rapidly increasing number of studies to reveal multifaceted aspects of translation control. The process of translation is not restricted to the conversion of mRNA coding sequences into polypeptide chains, it also controls the composition of the proteome in a delicate and responsive way. Therefore, translatomics has extended its unique and innovative power to many fields including proteomics, cancer research, bacterial stress response, biological rhythmicity and plant biology. Rational design in translation can enhance recombinant protein production for thousands of times. This brief review summarizes the main state-of-the-art methods of translatomics, highlights recent discoveries made in this field and introduces applications of translatomics on basic biological and biomedical research.

Post-transcriptional Processing of mRNA in Neurons: The Vestiges of the RNA World Drive Transcriptome Diversity

Neurons are morphologically complex cells that rely on the compartmentalization of protein expression to develop and maintain their extraordinary cytoarchitecture. This formidable task is achieved, at least in part, by targeting mRNA to subcellular compartments where they are rapidly translated. mRNA transcripts are the conveyor of genetic information from DNA to the translational machinery, however, they are also endowed with additional functions linked to both the coding sequence (open reading frame, or ORF) and the flanking 5′ and 3′ untranslated regions (UTRs), that may harbor coding-independent functions. In this review, we will highlight recent evidences supporting new coding-dependent and -independent functions of mRNA and discuss how nuclear and cytoplasmic post-transcriptional modifications of mRNA contribute to localization and translation in mammalian cells with specific emphasis on neurons. We also describe recently developed techniques that can be employed to study RNA dynamics at subcellular level in eukaryotic cells in developing and regenerating neurons.

SEASTAR: systematic evaluation of alternative transcription start sites in RNA

Alternative first exons diversify the transcriptomes of eukaryotes by producing variants of the 5' Untranslated Regions (5'UTRs) and N-terminal coding sequences. Accurate transcriptome-wide detection of alternative first exons typically requires specialized experimental approaches that are designed to identify the 5' ends of transcripts. We developed a computational pipeline SEASTAR that identifies first exons from RNA-seq data alone then quantifies and compares alternative first exon usage across multiple biological conditions. The exons inferred by SEASTAR coincide with transcription start sites identified directly by CAGE experiments and bear epigenetic hallmarks of active promoters. To determine if differential usage of alternative first exons can yield insights into the mechanism controlling gene expression, we applied SEASTAR to an RNA-seq dataset that tracked the reprogramming of mouse fibroblasts into induced pluripotent stem cells. We observed dynamic temporal changes in the usage of alternative first exons, along with correlated changes in transcription factor expression. Using a combined sequence motif and gene set enrichment analysis we identified N-Myc as a regulator of alternative first exon usage in the pluripotent state. Our results demonstrate that SEASTAR can leverage the available RNA-seq data to gain insights into the control of gene expression and alternative transcript variation in eukaryotic transcriptomes.

Intragenic transcriptional interference regulates the human immune ligand MICA

Many human genes have tandem promoters driving overlapping transcription, but the value of this distributed promoter configuration is generally unclear. Here we show that MICA, a gene encoding a ligand for the activating immune receptor NKG2D, contains a conserved upstream promoter that expresses a noncoding transcript. Transcription from the upstream promoter represses the downstream standard promoter activity in cis through transcriptional interference. The effect of transcriptional interference depends on the strength of transcription from the upstream promoter and can be described quantitatively by a simple reciprocal repressor function. Transcriptional interference coincides with recruitment at the standard downstream promoter of the FACT histone chaperone complex, which is involved in nucleosomal remodelling during transcription. The mechanism is invoked in the regulation of MICA expression by the physiological inputs interferon-γ and interleukin-4 that act on the upstream promoter. Genome-wide analysis indicates that transcriptional interference between tandem intragenic promoters may constitute a general mechanism with widespread importance in human transcriptional regulation.

Distinct transcriptomic and exomic abnormalities within myelodysplastic syndrome marrow cells

To provide biologic insights into mechanisms underlying myelodysplastic syndromes (MDS) we evaluated the CD34⁺ marrow cells transcriptome using high-throughput RNA sequencing (RNA-Seq). We demonstrated significant differential gene expression profiles (GEPs) between MDS and normal and identified 41 disease classifier genes. Additionally, two main clusters of GEPs distinguished patients based on their major clinical features, particularly between those whose disease remained stable versus patients who transformed into acute myeloid leukemia within 12 months. The genes whose expression was associated with disease outcome were involved in functional pathways and biologic processes highly relevant for MDS. Combined with exomic analysis we identified differential isoform usage of genes in MDS mutational subgroups, with consequent dysregulation of distinct biologic functions. This combination of clinical, transcriptomic and exomic findings provides valuable understanding of mechanisms underlying MDS and its progression to a more aggressive stage and also facilitates prognostic characterization of MDS patients.

The 5' Untranslated Region of the Major Immediate Early mRNA Is Necessary for Efficient Human Cytomegalovirus Replication

The human cytomegalovirus (HCMV) immediate early 1 (IE1) and IE2 proteins are critical regulators of virus replication. Both proteins are needed to efficiently establish lytic infection, and nascent expression of IE1 and IE2 is critical for reactivation from latency. The regulation of IE1 and IE2 protein expression is thus a central event in the outcome of HCMV infection. Transcription of the primary transcript encoding both IE1 and IE2 is well studied, but relatively little is known about the posttranscriptional mechanisms that control IE1 and IE2 protein synthesis. The mRNA 5' untranslated region (5' UTR) plays an important role in regulating mRNA translation. Therefore, to better understand the control of IE1 and IE2 mRNA translation, we examined the role of the shared 5' UTR of the IE1 and IE2 mRNAs (MIE 5' UTR) in regulating translation. In a cell-free system, the MIE 5' UTR repressed translation, as predicted based on its length and sequence composition. However, in transfected cells we found that the MIE 5' UTR increased the expression of a reporter gene and enhanced its association with polysomes, demonstrating that the MIE 5' UTR has a positive role in translation control. We also found that the MIE 5' UTR was necessary for efficient IE1 and IE2 translation during infection. Replacing the MIE 5' UTR with an unstructured sequence of the same length decreased IE1 and IE2 protein expression despite similar levels of IE1 and IE2 mRNA and reduced the association of the IE1 and IE2 mRNAs with polysomes. The wild-type MIE 5'-UTR sequence was also necessary for efficient HCMV replication. Together these data identify the shared 5' UTR of the IE1 and IE2 mRNAs as an important regulator of HCMV lytic replication.IMPORTANCE The HCMV IE1 and IE2 proteins are critical regulators of HCMV replication, both during primary infection and during reactivation from viral latency. Thus, defining factors that regulate IE1 and IE2 expression is important for understanding the molecular events controlling the HCMV replicative cycle. Here we identify a positive role for the MIE 5' UTR in mediating the efficient translation of the IE1 and IE2 mRNAs. This result is an important advance for several reasons. To date, most studies of IE1 and IE2 regulation have focused on defining events that regulate IE1 and IE2 transcription. Our work reveals that in addition to the regulation of transcription, IE1 and IE2 are also regulated at the level of translation. Therefore, this study is important in that it identifies an additional layer of regulation controlling IE1 and IE2 expression and thus HCMV pathogenesis. These translational regulatory events could potentially be targeted by novel antiviral therapeutics that limit IE1 and IE2 mRNA translation and thus inhibit lytic replication or prevent HCMV reactivation.

Alternative start and termination sites of transcription drive most transcript isoform differences across human tissues

Most human genes generate multiple transcript isoforms. The differential expression of these isoforms can help specify cell types. Diverse transcript isoforms arise from the use of alternative transcription start sites, polyadenylation sites and splice sites; however, the relative contribution of these processes to isoform diversity in normal human physiology is unclear. To address this question, we investigated cell type-dependent differences in exon usage of over 18 000 protein-coding genes in 23 cell types from 798 samples of the Genotype-Tissue Expression Project. We found that about half of the expressed genes displayed tissue-dependent transcript isoforms. Alternative transcription start and termination sites, rather than alternative splicing, accounted for the majority of tissue-dependent exon usage. We confirmed the widespread tissue-dependent use of alternative transcription start sites in a second, independent dataset, Cap Analysis of Gene Expression data from the FANTOM consortium. Moreover, our results indicate that most tissue-dependent splicing involves untranslated exons and therefore may not increase proteome complexity. Thus, alternative transcription start and termination sites are the principal drivers of transcript isoform diversity across tissues, and may underlie the majority of cell type specific proteomes and functions.

Genome-wide identification and differential analysis of translational initiation

Translation is principally regulated at the initiation stage. The development of the translation initiation (TI) sequencing (TI-seq) technique has enabled the global mapping of TIs and revealed unanticipated complex translational landscapes in metazoans. Despite the wide adoption of TI-seq, there is no computational tool currently available for analyzing TI-seq data. To fill this gap, we develop a comprehensive toolkit named Ribo-TISH, which allows for detecting and quantitatively comparing TIs across conditions from TI-seq data. Ribo-TISH can also predict novel open reading frames (ORFs) from regular ribosome profiling (rRibo-seq) data and outperform several established methods in both computational efficiency and prediction accuracy. Applied to published TI-seq/rRibo-seq data sets, Ribo-TISH uncovers a novel signature of elevated mitochondrial translation during amino-acid deprivation and predicts novel ORFs in 5'UTRs, long noncoding RNAs, and introns. These successful applications demonstrate the power of Ribo-TISH in extracting biological insights from TI-seq/rRibo-seq data.

Light Controls Protein Localization through Phytochrome-Mediated Alternative Promoter Selection

Alternative promoter usage is a proteome-expanding mechanism that allows multiple pre-mRNAs to be transcribed from a single gene. The impact of this mechanism on the proteome and whether it is positively exploited in normal organismal responses remain unclear. We found that the plant photoreceptor phytochrome induces genome-wide changes in alternative promoter selection in Arabidopsis thaliana. Through this mechanism, protein isoforms with different N termini are produced that display light-dependent differences in localization. For instance, shade-grown plants accumulate a cytoplasmic isoform of glycerate kinase (GLYK), an essential photorespiration enzyme that was previously thought to localize exclusively to the chloroplast. Cytoplasmic GLYK constitutes a photorespiratory bypass that alleviates fluctuating light-induced photoinhibition. Therefore, phytochrome controls alternative promoter selection to modulate protein localization in response to changing light conditions. This study suggests that alternative promoter usage represents another ubiquitous layer of gene expression regulation in eukaryotes that contributes to diversification of the proteome.

Cis-regulatory elements explain most of the mRNA stability variation across genes in yeast

The stability of mRNA is one of the major determinants of gene expression. Although a wealth of sequence elements regulating mRNA stability has been described, their quantitative contributions to half-life are unknown. Here, we built a quantitative model for Saccharomyces cerevisiae based on functional mRNA sequence features that explains 59% of the half-life variation between genes and predicts half-life at a median relative error of 30%. The model revealed a new destabilizing 3' UTR motif, ATATTC, which we functionally validated. Codon usage proves to be the major determinant of mRNA stability. Nonetheless, single-nucleotide variations have the largest effect when occurring on 3' UTR motifs or upstream AUGs. Analyzing mRNA half-life data of 34 knockout strains showed that the effect of codon usage not only requires functional decapping and deadenylation, but also the 5'-to-3' exonuclease Xrn1, the nonsense-mediated decay genes, but not no-go decay. Altogether, this study quantitatively delineates the contributions of mRNA sequence features on stability in yeast, reveals their functional dependencies on degradation pathways, and allows accurate prediction of half-life from mRNA sequence.

Major splice variants and multiple polyadenylation site utilization in mRNAs encoding human translation initiation factors eIF4E1 and eIF4E3 regulate the translational regulators?

Alternative polyadenylation is an important and pervasive mechanism that generates heterogeneous 3'-termini of mRNA and is considered an important regulator of gene expression. We performed bioinformatics analyses of ESTs and the 3'-UTRs of the main transcript splice variants of the translational initiation factor eIF4E1 and its family members, eIF4E2 and eIF4E3. This systematic analysis led to the prediction of new polyadenylation signals. All identified polyadenylation sites were subsequently verified by 3'RACE of transcripts isolated from human lymphoblastic cell lines. This led to the observation that multiple simultaneous polyadenylation site utilization occurs in single cell population. Importantly, we described the use of new polyadenylation site in the eIF4E1 mRNA, which lacked any known polyadenylation signal. The proportion of eIF4E1 transcripts derived from the first two polyadenylation sites in eIF4E1 mRNA achieved 15% in a wide range of cell lines. This result demonstrates the ubiquitous presence of ARE-lacking transcripts, which escape HuR/Auf1-mediated control, the main mechanism of eIF4E1 gene expression regulation. We found many EST clones documenting the significant production of transcript variants 2-4 of eIF4E2 gene that encode proteins with C-termini that were distinct from the mainly studied prototypical isoform A. Similarly, eIF4E3 mRNAs are produced as two main variants with the same very long 3'-UTR with potential for heavy post-transcriptional regulation. We identified sparsely documented transcript variant 1 of eIF4E3 gene in human placenta. eIF4E3 truncated transcript variants were found mainly in brain. We propose to elucidate the minor splice variants of eIF4E2 and eIF4E3 in great detail because they might produce proteins with modified features that fulfill different cellular roles from their major counterparts.

Efficient and Accurate Translation Initiation Directed by TISU Involves RPS3 and RPS10e Binding and Differential Eukaryotic Initiation Factor 1A Regulation

Canonical translation initiation involves ribosomal scanning, but short 5' untranslated region (5'UTR) mRNAs are translated in a scanning-independent manner. The extent and mechanism of scanning-independent translation are not fully understood. Here we report that short 5'UTR mRNAs constitute a substantial fraction of the translatome. Short 5'UTR mRNAs are enriched with TISU (translation initiator of short 5'UTR), a 12-nucleotide element directing efficient scanning-independent translation. Comprehensive mutagenesis revealed that each AUG codon-flanking nucleotide of TISU contributes to translational strength, but only a few are important for accuracy. Using site-specific UV cross-linking of ribosomal complexes assembled on TISU mRNA, we demonstrate specific binding of TISU to ribosomal proteins at the E and A sites. We identified RPS3 as the major TISU binding protein in the 48S complex A site. Upon 80S complex formation, RPS3 interaction is weakened and switched to RPS10e (formerly called RPS10). We further demonstrate that TISU is particularly dependent on eukaryotic initiation factor 1A (eIF1A) which interacts with both RPS3 and RPS10e. Our findings suggest that the cap-recruited ribosome specifically binds the TISU nucleotides at the A and E sites in cooperation with eIF1A to promote scanning arrest.

Identification and prediction of alternative transcription start sites that generate rod photoreceptor-specific transcripts from ubiquitously expressed genes

Transcriptome complexity is substantially increased by the use of multiple transcription start sites for a given gene. By utilizing a rod photoreceptor-specific chromatin signature, and the RefSeq database of established transcription start sites, we have identified essentially all known rod photoreceptor genes as well as a group of novel genes that have a high probability of being expressed in rod photoreceptors. Approximately half of these novel rod genes are transcribed into multiple mRNA and/or protein isoforms through alternative transcriptional start sites (ATSS), only one of which has a rod-specific epigenetic signature and gives rise to a rod transcript. This suggests that, during retina development, some genes use ATSS to regulate cell type and temporal specificity, effectively generating a rod transcript from otherwise ubiquitously expressed genes. Biological confirmation of the relationship between epigenetic signatures and gene expression, as well as comparison of our genome-wide chromatin signature maps with available data sets for retina, namely a ChIP-on-Chip study of Polymerase-II (Pol-II) binding sites, ChIP-Seq studies for NRL- and CRX- binding sites and DHS (University of Washington data, available on UCSC mouse Genome Browser as a part of ENCODE project) fully support our hypothesis and together accurately identify and predict an array of new rod transcripts. The same approach was used to identify a number of TSS that are not currently in RefSeq. Biological conformation of the use of some of these TSS suggests that this method will be valuable for exploring the range of transcriptional complexity in many tissues. Comparison of mouse and human genome-wide data indicates that most of these alternate TSS appear to be present in both species, indicating that our approach can be useful for identification of regulatory regions that might play a role in human retinal disease.

A researcher's guide to the galaxy of IRESs

The idea of internal initiation is frequently exploited to explain the peculiar translation properties or unusual features of some eukaryotic mRNAs. In this review, we summarize the methods and arguments most commonly used to address cases of translation governed by internal ribosome entry sites (IRESs). Frequent mistakes are revealed. We explain why "cap-independent" does not readily mean "IRES-dependent" and why the presence of a long and highly structured 5' untranslated region (5'UTR) or translation under stress conditions cannot be regarded as an argument for appealing to internal initiation. We carefully describe the known pitfalls and limitations of the bicistronic assay and artefacts of some commercially available in vitro translation systems. We explain why plasmid DNA transfection should not be used in IRES studies and which control experiments are unavoidable if someone decides to use it anyway. Finally, we propose a workflow for the validation of IRES activity, including fast and simple experiments based on a single genetic construct with a sequence of interest.

Cap-proximal nucleotides via differential eIF4E binding and alternative promoter usage mediate translational response to energy stress

Transcription start-site (TSS) selection and alternative promoter (AP) usage contribute to gene expression complexity but little is known about their impact on translation. Here we performed TSS mapping of the translatome following energy stress. Assessing the contribution of cap-proximal TSS nucleotides, we found dramatic effect on translation only upon stress. As eIF4E levels were reduced, we determined its binding to capped-RNAs with different initiating nucleotides and found the lowest affinity to 5'cytidine in correlation with the translational stress-response. In addition, the number of differentially translated APs was elevated following stress. These include novel glucose starvation-induced downstream transcripts for the translation regulators eIF4A and Pabp, which are also translationally-induced despite general translational inhibition. The resultant eIF4A protein is N-terminally truncated and acts as eIF4A inhibitor. The induced Pabp isoform has shorter 5'UTR removing an auto-inhibitory element. Our findings uncovered several levels of coordination of transcription and translation responses to energy stress.

DOI:http://dx.doi.org/10.7554/eLife.21907.001

The production of new proteins is a complex process that occurs in two steps known as transcription and translation. During transcription, the cell copies a section of DNA to make molecules of messenger ribonucleic acid (or mRNA for short) in the nucleus of the cell. The mRNA then leaves the nucleus and enters another cell compartment called the cytoplasm, where it serves as a template to make proteins during translation. A mRNA molecule contains a sequence of building blocks known as nucleotides. There are four different types of nucleotides in mRNA and the order they appear in the sequence determines how the protein is built.

Both transcription and translation consume a lot of energy so they are highly regulated and sensitive to environmental changes. However, since transcription and translation happen in different cell compartments, it is not known if and how they are coordinated under stress.

Tamarkin-Ben-Harush et al. studied transcription and translation in mouse cells that were starved of glucose. The experiments show that the identity of the very first nucleotide in the mRNA – which is dictated during transcription – has a dramatic influence on the translation of the mRNA, especially when the cells are starved of glucose. This first nucleotide affects the ability of a protein called eIF4E, which recruits the machinery needed for translation, to bind to the mRNA. The experiments also show that there is a dramatic increase in the number of distinct mRNAs that are transcribed from the same section of DNA but translated in a different way during glucose starvation.

The findings of Tamarkin-Ben-Harush et al. show that transcription and translation are highly coordinated when cells are starved of glucose, allowing the cells to cope with the stress. The next step is to further analyze the data to find out more about how transcription and translation are linked.

DOI:http://dx.doi.org/10.7554/eLife.21907.002

Global analysis of regulatory divergence in the evolution of mouse alternative polyadenylation

Alternative polyadenylation (APA), which is regulated by both cis‐elements and trans‐factors, plays an important role in post‐transcriptional regulation of eukaryotic gene expression. However, comparing to the extensively studied transcription and alternative splicing, the extent of APA divergence during evolution and the relative cis‐ and trans‐contribution remain largely unexplored. To directly address these questions for the first time in mammals, by using deep sequencing‐based methods, we measured APA divergence between C57BL/6J and SPRET/EiJ mouse strains as well as allele‐specific APA pattern in their F1 hybrids. Among the 24,721 polyadenylation sites (pAs) from 7,271 genes expressing multiple pAs, we identified 3,747 pAs showing significant divergence between the two strains. After integrating the allele‐specific data from F1 hybrids, we demonstrated that these events could be predominately attributed to cis‐regulatory effects. Further systematic sequence analysis of the regions in proximity to cis‐divergent pAs revealed that the local RNA secondary structure and a poly(U) tract in the upstream region could negatively modulate the pAs usage.

Pervasive isoform-specific translational regulation via alternative transcription start sites in mammals

Transcription initiated at alternative sites can produce mRNA isoforms with different 5'UTRs, which are potentially subjected to differential translational regulation. However, the prevalence of such isoform-specific translational control across mammalian genomes is currently unknown. By combining polysome profiling with high-throughput mRNA 5' end sequencing, we directly measured the translational status of mRNA isoforms with distinct start sites. Among 9,951 genes expressed in mouse fibroblasts, we identified 4,153 showed significant initiation at multiple sites, of which 745 genes exhibited significant isoform-divergent translation. Systematic analyses of the isoform-specific translation revealed that isoforms with longer 5'UTRs tended to translate less efficiently. Further investigation of cis-elements within 5'UTRs not only provided novel insights into the regulation by known sequence features, but also led to the discovery of novel regulatory sequence motifs. Quantitative models integrating all these features explained over half of the variance in the observed isoform-divergent translation. Overall, our study demonstrated the extensive translational regulation by usage of alternative transcription start sites and offered comprehensive understanding of translational regulation by diverse sequence features embedded in 5'UTRs.