UTR Length Research Articles

Differential Expression of Sox11 and Bdnf mRNA Isoforms in the Injured and Regenerating Nervous Systems.

In both the central nervous system (CNS) and the peripheral nervous system (PNS), axonal injury induces changes in neuronal gene expression. In the PNS, a relatively well-characterized alteration in transcriptional activation is known to promote axonal regeneration. This transcriptional cascade includes the neurotrophin Bdnf and the transcription factor Sox11. Although both molecules act to facilitate successful axon regeneration in the PNS, this process does not occur in the CNS. The present study examines the differential expression of Sox11 and Bdnf mRNA isoforms in the PNS and CNS using three experimental paradigms at different time points: (i) the acutely injured CNS (retina after optic nerve crush) and PNS (dorsal root ganglion after sciatic nerve crush), (ii) a CNS regeneration model (retina after optic nerve crush and induced regeneration); and (iii) the retina during a chronic form of central neurodegeneration (the DBA/2J glaucoma model). We find an initial increase of Sox11 in both PNS and CNS after injury; however, the expression of Bdnf isoforms is higher in the PNS relative to the CNS. Sustained upregulation of Sox11 is seen in the injured retina following regeneration treatment, while the expression of two Bdnf mRNA isoforms is suppressed. Furthermore, two isoforms of Sox11 with different 3′UTR lengths are present in the retina, and the long isoform is specifically upregulated in later stages of glaucoma. These results provide insight into the molecular cascades active during axonal injury and regeneration in mammalian neurons.

Alterations in mRNA 3′ UTR Isoform Abundance Accompany Gene Expression Changes in Human Huntington’s Disease Brains

The huntingtin gene has two mRNA isoforms that differ in their 3' UTR length. The relationship of these isoforms with Huntington's disease is not established. We provide evidence that the abundance of huntingtin 3' UTR isoforms differs between patient and control neural stem cells, fibroblasts, motor cortex, and cerebellum. Huntingtin 3' UTR isoforms, including a mid-3' UTR isoform, have different localizations, half-lives, polyA tail lengths, microRNA sites, and RNA-binding protein sites. Isoform shifts in Huntington's disease motor cortex are not limited to huntingtin; 11% of alternatively polyadenylated genes change the abundance of their 3' UTR isoforms. Altered expression of RNA-binding proteins may be associated with aberrant isoform abundance; knockdown of the RNA-binding protein CNOT6 in control fibroblasts leads to huntingtin isoform differences similar to those in disease fibroblasts. These findings demonstrate that mRNA 3'UTR isoform changes are a feature of molecular pathology in the Huntington's disease brain.

Regulation by 3'-Untranslated Regions.

3'-untranslated regions (3'-UTRs) are the noncoding parts of mRNAs. Compared to yeast, in humans, median 3'-UTR length has expanded approximately tenfold alongside an increased generation of alternative 3'-UTR isoforms. In contrast, the number of coding genes, as well as coding region length, has remained similar. This suggests an important role for 3'-UTRs in the biology of higher organisms. 3'-UTRs are best known to regulate diverse fates of mRNAs, including degradation, translation, and localization, but they can also function like long noncoding or small RNAs, as has been shown for whole 3'-UTRs as well as for cleaved fragments. Furthermore, 3'-UTRs determine the fate of proteins through the regulation of protein-protein interactions. They facilitate cotranslational protein complex formation, which establishes a role for 3'-UTRs as evolved eukaryotic operons. Whereas bacterial operons promote the interaction of two subunits, 3'-UTRs enable the formation of protein complexes with diverse compositions. All of these 3'-UTR functions are accomplished by effector proteins that are recruited by RNA-binding proteins that bind to 3'-UTR cis-elements. In summary, 3'-UTRs seem to be major players in gene regulation that enable local functions, compartmentalization, and cooperativity, which makes them important tools for the regulation of phenotypic diversity of higher organisms.

DRNA-seq transcriptional profiling of the FK506 biosynthetic gene cluster in Streptomyces tsukubaensis NRRL18488 and general analysis of the transcriptome

ABSTRACTFK506 (tacrolimus) is a valuable immunosuppressant produced by several Streptomyces strains. In the genome of the wild type producer Streptomyces tsukubaensis NRRL18488, FK506 biosynthesis is encoded by a gene cluster that spans 83.5 (kb). A whole transcriptome differential shotgun sequencing (dRNA-seq) of S. tsukubaensis was performed to analyze transcription at 2 different time points; before and during active FK506 production. In total, 8,914 transcription start sites were identified in either condition, which enabled precise determination of the 5′-UTR length of the corresponding transcripts as well as the identification of 2 consensus sequence motifs in the promoter regions. The transcription start sites of all gene operons within the FK506 cluster were identified, including 3 examples of leaderless RNA transcripts. These data provide detailed insight into the transcription of the FK506 biosynthetic gene cluster to support future regulatory studies, genetic manipulation, and industrial production.

Abstract 3374: APA isoform diversity in triple negative breast cancers

Abstract Alternative polyadenylation (APA) plays a role in gene expression regulation generally by shortening of 3'UTRs and relieving microRNA-mediated repression. Therefore, APA is gaining increased attention as a potential mechanism to activate oncogenes. Owing to high proliferative indices of triple negative breast cancers (TNBCs), we hypothesized APA to cause 3'UTR length changes in this aggressive subgroup of breast cancers. Our probe-based meta-analysis approach identified 3'UTR length alterations where the significant majority was shortening events (70%, 113 of 165) of mostly proliferation-related transcripts in over 500 TNBC patients compared with normal breast tissue. Representative shortening events correlated with increased protein levels and relapse free survival of patients, suggesting functional significance of isoform variability. To begin addressing the underlying mechanisms of 3’UTR shortening, we turned to APA machinery proteins. We detected variable expression of APA machinery proteins in different breast cancer subtypes but CSTF2 (cleavage stimulation factor 2) has the most prominent overexpression in breast cancer cells. Therefore, among potential regulators of 3’UTR shortening, we further investigated the role of CSTF2 in proximal polyA signal selection. Because some of the TNBC patients are EGFR positive, we found EGF treatment to cause increased CSTF2 levels. Higher CSTF2 levels indeed correlated with further shortening of the 3'UTRs. Accordingly, RNAi-induced silencing of CSTF2 decreased the proliferative rate of cancer cells. Therefore, our integrated approach revealed a pattern of 3'UTR length changes in TNBC patients and a potential link between APA and EGF signaling. Further studies are underway to investigate the mechanistic link between EGF signaling and regulation of 3’UTR lengths. (This work is funded by TUBITAK 112S478 and 114Z884) Citation Format: Ayse Elif Erson-Bensan, Begum H. Akman, Merve Oyken, Hizlan H. Agus, Esra Yavuz, Murat Erdem, Tolga Can. APA isoform diversity in triple negative breast cancers [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr 3374. doi:10.1158/1538-7445.AM2017-3374

Abstract LB-324: Methylation-associated silencing of miR-193a-3p promotes ovarian cancer aggressiveness via targeting GRB7

Abstract The human growth factor receptor-bound protein-7 (GRB7) is a pivotal mediator involved in receptor tyrosine kinase signaling and governing diverse cellular processes. We and others have reported that GRB7 is frequently upregulated and associated with the progression of human cancers. Clinicopathological analysis has shown that the overexpressed GRB7 is correlated with high grade and metastatic ovarian cancers. However, the molecular mechanisms leading to the upregulation of GRB7 remain largely unknown. In this study, we propose that the epigenetic modification of GRB7 at the post-transcriptional level may be a crucial factor leading to GRB7 upregulation in ovarian cancers. We found that the protein level of GRB7 was much higher than its mRNA levels in a subset of ovarian cancer cell lines, indicating aberrant post-transcriptional alterations in GRB7. Using an in silico study, we identified miR-193a-3p and its isoform, miR-193b-3p, by specifically targeting the conserved site (GGCCAGT) at position 332-338 in the 3’ UTR (length: 387) of the human GRB7 gene. QPCR and western blot analyses revealed that only the expression levels of miR-193a-3p but not miR-193b-3p were inversely correlated in ovarian cancer cell lines accompanied with high levels of GRB7 protein levels. Further biochemical analyses using luciferase reporter assays in combination with mutational analyses identified miR-193a-3p specifically targeting GRB7 3’UTR, confirming miR-193a-3p, instead of its isoform miR-193b-3p, has a more important role in GRB7 post-transcriptional regulation in ovarian cancer cells. Constitutive expression or knockdown of miR-193a-3p could significantly alter GRB7 expression, as well as oncogenic capacities of ovarian cancer cells such as cell growth, cell migration and cell invasion. Importantly, we observed that there was a significant stepwise decrease of miR-193a-3p expression from low to high grade ovarian cancers, and was inversely correlated with GRB7 expressions. Intriguingly, treatment with the DNA methyl transferase inhibitor (5’-Aza-dc) could restore the expression of miR-193a-3p, suggesting the epigenetic inhibition of miR-193a-3p is clinically relevant. Importantly, such epigenetic alteration may be the clue leading to overexpression of GRB7 in ovarian cancer cells. Taken together, our findings suggest that epigenetically silencing of miR-193a-3p upregulates GRB7 and contribute to ovarian cancer aggressiveness in tumor progression. Citation Format: Kangmei Chen, Hextan YS Ngan, David W. Chan. Methylation-associated silencing of miR-193a-3p promotes ovarian cancer aggressiveness via targeting GRB7 [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2017; 2017 Apr 1-5; Washington, DC. Philadelphia (PA): AACR; Cancer Res 2017;77(13 Suppl):Abstract nr LB-324. doi:10.1158/1538-7445.AM2017-LB-324

MicroRNAs control mRNA fate by compartmentalization based on 3\u2032 UTR length in male germ cells

BackgroundPost-transcriptional regulation of gene expression can be achieved through the control of mRNA stability, cytoplasmic compartmentalization, 3′ UTR length and translational efficacy. Spermiogenesis, a process through which haploid male germ cells differentiate into spermatozoa, represents an ideal model for studying post-transcriptional regulation in vivo because it involves a large number of transcripts that are physically sequestered in ribonucleoprotein particles (RNPs) and thus subjected to delayed translation. To explore how small RNAs regulate mRNA fate, we conducted RNA-Seq analyses to determine not only the levels of both mRNAs and small noncoding RNAs, but also their cytoplasmic compartmentalization during spermiogenesis.ResultAmong all small noncoding RNAs studied, miRNAs displayed the most dynamic changes in both abundance and subcytoplasmic localization. mRNAs with shorter 3′ UTRs became increasingly enriched in RNPs from pachytene spermatocytes to round spermatids, and the enrichment of shorter 3′ UTR mRNAs in RNPs coincided with newly synthesized miRNAs that target these mRNAs at sites closer to the stop codon. In contrast, the translocation of longer 3′ UTR mRNAs from RNPs to polysomes correlated with the production of new miRNAs that target these mRNAs at sites distal to the stop codon.ConclusionsmiRNAs appear to control cytoplasmic compartmentalization of mRNAs based on 3′ UTR length. Our data suggest that transcripts with longer 3′ UTRs tend to contain distal miRNA binding sites and are thus targeted to polysomes for translation followed by degradation. In contrast, those with shorter 3′ UTRs only possess proximal miRNA binding sites, which, therefore, are targeted into RNPs for enrichment and delayed translation.

Multi-platform analysis reveals a complex transcriptome architecture of a circovirus

In this study, we used Pacific Biosciences RS II long-read and Illumina HiScanSQ short-read sequencing technologies for the characterization of porcine circovirus type 1 (PCV-1) transcripts. Our aim was to identify novel RNA molecules and transcript isoforms, as well as to determine the exact 5′- and 3′-end sequences of previously described transcripts with single base-pair accuracy. We discovered a novel 3′-UTR length isoform of the Cap transcript, and a non-spliced Cap transcript variant. Additionally, our analysis has revealed a 3′-UTR isoform of Rep and two 5′-UTR isoforms of Rep’ transcripts, and a novel splice variant of the longer Rep’ transcript. We also explored two novel long transcripts, one with a previously identified splice site, and a formerly undetected mRNA of ORF3. Altogether, our methods have identified nine novel RNA molecules, doubling the size of PCV-1 transcriptome that had been known before. Additionally, our investigations revealed an intricate pattern of transcript overlapping, which might produce transcriptional interference between the transcriptional machineries of adjacent genes, and thereby may potentially play a role in the regulation of gene expression in circoviruses.

Regulation of AGR2 expression via 3’UTR shortening

One recently discussed general mechanism affecting gene expression is 3’-untranslated region (3’UTR) length. Events such as shortening, translocation or loss of 3’UTRs are observed within oncogenes and are proposed to associate with increased expression. Thus, increased efforts are being made to understand constitutive and differential transcript 3’end formation. Investigation of AGR2 mRNA revealed a direct impact of its 3’UTR length on AGR2 expression. In silico analyses identified several regulatory sequences within the distal part of AGR2 mRNA that may regulate 3’UTR length and associated protein levels. Short 3’UTRs were observed in a panel of AGR2-positive cancer cell lines and in human breast cancer specimens, in which more extensive 3’UTR shortening correlated with increased AGR2 protein levels. AGR2 is an important member of PI3K/AKT signalling pathway, which along with the proposed involvement of mTOR in the regulation of alternative polyadenylation, prompted us to study the role of mTOR in relation to AGR2 mRNA 3’UTR shortening. A direct impact of mTOR signalling on AGR2 3’UTR shortening associated with increased protein synthesis was found, which led to the identification of a novel molecular mechanism involved in upregulation of AGR2 levels in mTOR-activated cells via modulating the 3’UTR length of AGR2 mRNA.

Molecular Regulation of Alternative Polyadenylation (APA) within the Drosophila Nervous System

Alternative polyadenylation (APA) is a widespread gene regulatory mechanism that generates mRNAs with different 3′-ends, allowing them to interact with different sets of RNA regulators such as microRNAs and RNA-binding proteins. Recent studies have shown that during development, neural tissues produce mRNAs with particularly long 3′UTRs, suggesting that such extensions might be important for neural development and function. Despite this, the mechanisms underlying neural APA are not well understood. Here, we investigate this problem within the Drosophila nervous system, focusing on the roles played by general cleavage and polyadenylation factors (CPA factors). In particular, we examine the model that modulations in CPA factor concentration may affect APA during development. For this, we first analyse the expression of the Drosophila orthologues of all mammalian CPA factors and note that their expression decreases during embryogenesis. In contrast to this global developmental decrease in CPA factor expression, we see that cleavage factor I (CFI) expression is actually elevated in the late embryonic central nervous system, suggesting that CFI might play a special role in neural tissues. To test this, we use the UAS/Gal4 system to deplete CFI proteins from neural tissue and observe that in this condition, multiple genes switch their APA patterns, demonstrating a role of CFI in APA control during Drosophila neural development. Furthermore, analysis of genes with 3′UTR extensions of different length leads us to suggest a novel relation between 3′UTR length and sensitivity to CPA factor expression. Our work thus contributes to the understanding of the mechanisms of APA control within the developing central nervous system.

The role of alternative polyadenylation in the antiviral innate immune response

Alternative polyadenylation (APA) is an important regulatory mechanism of gene functions in many biological processes. However, the extent of 3′ UTR variation and the function of APA during the innate antiviral immune response are unclear. Here, we show genome-wide poly(A) sites switch and average 3′ UTR length shortens gradually in response to vesicular stomatitis virus (VSV) infection in macrophages. Genes with APA and mRNA abundance change are enriched in immune-related categories such as the Toll-like receptor, RIG-I-like receptor, JAK-STAT and apoptosis-related signalling pathways. The expression of 3′ processing factors is down-regulated upon VSV infection. When the core 3′ processing factors are knocked down, viral replication is affected. Thus, our study reports the annotation of genes with APA in antiviral immunity and highlights the roles of 3′ processing factors on 3′ UTR variation upon viral infection.

Effects of Kasugamycin on the Translatome of Escherichia coli

It is long known that Kasugamycin inhibits translation of canonical transcripts containing a 5’-UTR with a Shine Dalgarno (SD) motif, but not that of leaderless transcripts. To gain a global overview of the influence of Kasugamycin on translation efficiencies, the changes of the translatome of Escherichia coli induced by a 10 minutes Kasugamycin treatment were quantified. The effect of Kasugamycin differed widely, 102 transcripts were at least twofold more sensitive to Kasugamycin than average, and 137 transcripts were at least twofold more resistant, and there was a more than 100-fold difference between the most resistant and the most sensitive transcript. The 5’-ends of 19 transcripts were determined from treated and untreated cultures, but Kasugamycin resistance did neither correlate with the presence or absence of a SD motif, nor with differences in 5’-UTR lengths or GC content. RNA Structure Logos were generated for the 102 Kasugamycin-sensitive and for the 137 resistant transcripts. For both groups a short Shine Dalgarno (SD) motif was retrieved, but no specific motifs associated with resistance or sensitivity could be found. Notably, this was also true for the region -3 to -1 upstream of the start codon and the presence of an extended SD motif, which had been proposed to result in Kasugamycin resistance. Comparison of the translatome results with the database RegulonDB showed that the transcript with the highest resistance was leaderless, but no further leaderless transcripts were among the resistant transcripts. Unexpectedly, it was found that translational coupling might be a novel feature that is associated with Kasugamycin resistance. Taken together, Kasugamycin has a profound effect on translational efficiencies of E. coli transcripts, but the mechanism of action is different than previously described.

A Novel Strategy to Identify mRNA 3′ UTR and mRNA Level Difference in Crohn’s Disease

Abstract Background and Aims: Crohn’s disease (CD) is a chronic and refractory intestinal inflammatory disease. The 3′ untranslated region (3′ UTR) of mRNA transcript can regulate its own translational process post-transciptionally, and its role in CD remains unclear. The aim of this study was to identify the critical genes influencing CD phenotype via the regulation of mRNA 3′ UTR. Methods: Isoform Structural Change Model (IsoSCM) was used to evaluate the length of 3′ UTR of the whole transciptome from a RNA-seq based online CD database. Correlations between 3′ UTR length status and mRNA level were analyzed by Spearman coefficients. Correlated genes list was overlapped with published critical CD related gene list. Univariate and multivariate analysis were performed to identify the genes associated with CD phenotype. Results: Compared with normal control, 3′ UTR of 34192 genes were shortened and 57389 were lengthened among 432784 genes in CD patients. The 3′ UTR changes of 255 genes were found significantly correlated with their mRNA level. Furthermore, 8 overlapped genes were shared by above 255 correlated genes and known CD related gene list (ABCB1, FAF1, TUT2, IL27, JAK2, LTB, NAGLU and PTPN2). The mRNA level of ABCB1 and JAK2, and shortened 3′ UTR length of JAK2 transcript were found to be correlated with deep ulcer in CD patients by univariate and multivariate analysis. Conclusions: The 3′ UTR changes of CD related genes were confirmed to be related to the phenotype of CD. Our findings may provide further insight into the epigenetic mechanisms and therapeutic strategies for CD.

Cellular Localization and Regulation of Expression of the PLET1 Gene in Porcine Placenta.

The placenta expressed transcript 1 (PLET1) gene, which is expressed in placentas of pigs and mice, has been found to have a potential role in trophoblast cell fate decision in mice. Results of this study showed that the porcine PLET1 mRNA and protein were expressed exclusively in trophoblast cells on Days 15, 26, 50, and 95 of gestation (gestation length in the pig is 114 days), indicating that the PLET1 could be a useful marker for porcine trophoblast cells. Additionally, PLET1 protein was found to be redistributed from cytoplasm to the apical side of trophoblast cells as gestation progresses, which suggests a role of PLET1 in the establishment of a stable trophoblast and endometrial epithelial layers. In addition, two transcripts that differ in the 3′ UTR length but encode identical protein were identified to be generated by the alternative cleavage and polyadenylation (APA), and the expression of PLET1-L transcript was significantly upregulated in porcine placentas as gestation progresses. Furthermore, we demonstrated the interaction between the miR-365-3p and PLET1 gene using luciferase assay system. Our findings imply an important role of PLET1 in the placental development in pigs.

Genome-wide map of RNA degradation kinetics patterns in dendritic cells after LPS stimulation facilitates identification of primary sequence and secondary structure motifs in mRNAs.

BackgroundImmune cells have to change their gene expression patterns dynamically in response to external stimuli such as lipopolysaccharide (LPS). The gene expression is regulated at multiple steps in eukaryotic cells, in which control of RNA levels at both the transcriptional level and the post-transcriptional level plays important role. Impairment of the control leads to aberrant immune responses such as excessive or impaired production of cytokines. However, genome-wide studies focusing on the post-transcriptional control were relatively rare until recently. Moreover, several RNA cis elements and RNA-binding proteins have been found to be involved in the process, but our general understanding remains poor, partly because identification of regulatory RNA motifs is very challenging in spite of its importance. We took advantage of genome-wide measurement of RNA degradation in combination with estimation of degradation kinetics by qualitative approach, and performed de novo prediction of RNA sequence and structure motifs.MethodsTo classify genes by their RNA degradation kinetics, we first measured RNA degradation time course in mouse dendritic cells after LPS stimulation and the time courses were clustered to estimate degradation kinetics and to find patterns in the kinetics. Then genes were clustered by their similarity in degradation kinetics patterns. The 3′ UTR sequences of a cluster was subjected to de novo sequence or structure motif prediction.ResultsThe quick degradation kinetics was found to be strongly associated with lower gene expression level, immediate regulation (both induction and repression) of gene expression level, and longer 3′ UTR length. De novo sequence motif prediction found AU-rich element-like and TTP-binding sequence-like motifs which are enriched in quickly degrading genes. De novo structure motif prediction found a known functional motif, namely stem-loop structure containing sequence bound by RNA-binding protein Roquin and Regnase-1, as well as unknown motifs.ConclusionsThe current study indicated that degradation kinetics patterns lead to classification different from that by gene expression and the differential classification facilitates identification of functional motifs. Identification of novel motif candidates implied post-transcriptional controls different from that by known pairs of RNA-binding protein and RNA motif.Electronic supplementary materialThe online version of this article (doi:10.1186/s12864-016-3325-7) contains supplementary material, which is available to authorized users.

KDM5 lysine demethylases are involved in maintenance of 3'UTR length.

The complexity by which cells regulate gene and protein expression is multifaceted and intricate. Regulation of 3' untranslated region (UTR) processing of mRNA has been shown to play a critical role in development and disease. However, the process by which cells select alternative mRNA forms is not well understood. We discovered that the Saccharomyces cerevisiae lysine demethylase, Jhd2 (also known as KDM5), recruits 3'UTR processing machinery and promotes alteration of 3'UTR length for some genes in a demethylase-dependent manner. Interaction of Jhd2 with both chromatin and RNA suggests that Jhd2 affects selection of polyadenylation sites through a transcription-coupled mechanism. Furthermore, its mammalian homolog KDM5B (also known as JARID1B or PLU1), but not KDM5A (also known as JARID1A or RBP2), promotes shortening of CCND1 transcript in breast cancer cells. Consistent with these results, KDM5B expression correlates with shortened CCND1 in human breast tumor tissues. In contrast, both KDM5A and KDM5B are involved in the lengthening of DICER1. Our findings suggest both a novel role for this family of demethylases and a novel targetable mechanism for 3'UTR processing.

Roar: detecting alternative polyadenylation with standard mRNA sequencing libraries.

BackgroundPost-transcriptional regulation is a complex mechanism that plays a central role in defining multiple cellular identities starting from a common genome. Modifications in the length of 3’UTRs have been found to play an important role in this context, since alternative 3’ UTRs could lead to differences for example in regulation by microRNAs and cellular localization of the transcripts thus altering their fate.ResultsWe propose a strategy to identify the genes undergoing regulation of 3’ UTR length using RNA sequencing data obtained from standard libraries, thus widely applicable to data originally obtained to perform classical differential expression analyses. We decided to exploit previously annotated APA sites from public databases, in contrast with other approaches recently proposed in which the location of the APA site is inferred from the data together with the relative abundance of the isoforms.We demonstrate the reliability of our method by comparing it to the results of other microarray based or specific RNA-seq libraries methods and show that using APA sites databases results in higher sensitivity compared to de novo site prediction approach.ConclusionsWe implemented the algorithm in a Bioconductor package to facilitate its broad usage in the scientific community. The ability of this approach to detect shortening from libraries with a number of reads comparable to that needed for differential expression analyses makes it useful for investigating if alternative polyadenylation is relevant in a certain biological process without requiring specific experimental assays.Electronic supplementary materialThe online version of this article (doi:10.1186/s12859-016-1254-8) contains supplementary material, which is available to authorized users.

Identifying Aspects of the Post-Transcriptional Program Governing the Proteome of the Green Alga Micromonas pusilla.

Micromonas is a unicellular motile alga within the Prasinophyceae, a green algal group that is related to land plants. This picoeukaryote (<2 μm diameter) is widespread in the marine environment but is not well understood at the cellular level. Here, we examine shifts in mRNA and protein expression over the course of the day-night cycle using triplicated mid-exponential, nutrient replete cultures of Micromonas pusilla CCMP1545. Samples were collected at key transition points during the diel cycle for evaluation using high-throughput LC-MS proteomics. In conjunction, matched mRNA samples from the same time points were sequenced using pair-ended directional Illumina RNA-Seq to investigate the dynamics and relationship between the mRNA and protein expression programs of M. pusilla. Similar to a prior study of the marine cyanobacterium Prochlorococcus, we found significant divergence in the mRNA and proteomics expression dynamics in response to the light:dark cycle. Additionally, expressional responses of genes and the proteins they encoded could also be variable within the same metabolic pathway, such as we observed in the oxygenic photosynthesis pathway. A regression framework was used to predict protein levels from both mRNA expression and gene-specific sequence-based features. Several features in the genome sequence were found to influence protein abundance including codon usage as well as 3’ UTR length and structure. Collectively, our studies provide insights into the regulation of the proteome over a diel cycle as well as the relationships between transcriptional and translational programs in the widespread marine green alga Micromonas.

A comprehensive analysis of 3' end sequencing data sets reveals novel polyadenylation signals and the repressive role of heterogeneous ribonucleoprotein C on cleavage and polyadenylation.

Alternative polyadenylation (APA) is a general mechanism of transcript diversification in mammals, which has been recently linked to proliferative states and cancer. Different 3′ untranslated region (3′ UTR) isoforms interact with different RNA-binding proteins (RBPs), which modify the stability, translation, and subcellular localization of the corresponding transcripts. Although the heterogeneity of pre-mRNA 3′ end processing has been established with high-throughput approaches, the mechanisms that underlie systematic changes in 3′ UTR lengths remain to be characterized. Through a uniform analysis of a large number of 3′ end sequencing data sets, we have uncovered 18 signals, six of which are novel, whose positioning with respect to pre-mRNA cleavage sites indicates a role in pre-mRNA 3′ end processing in both mouse and human. With 3′ end sequencing we have demonstrated that the heterogeneous ribonucleoprotein C (HNRNPC), which binds the poly(U) motif whose frequency also peaks in the vicinity of polyadenylation (poly(A)) sites, has a genome-wide effect on poly(A) site usage. HNRNPC-regulated 3′ UTRs are enriched in ELAV-like RBP 1 (ELAVL1) binding sites and include those of the CD47 gene, which participate in the recently discovered mechanism of 3′ UTR–dependent protein localization (UDPL). Our study thus establishes an up-to-date, high-confidence catalog of 3′ end processing sites and poly(A) signals, and it uncovers an important role of HNRNPC in regulating 3′ end processing. It further suggests that U-rich elements mediate interactions with multiple RBPs that regulate different stages in a transcript's life cycle.

727 - Expression of PTCH1b tumor suppressor gene is controlled by different 5’-untranslated region cis-regulatory elements

Background The PTCH1 tumor suppressor gene encodes for a 12-pass transmembrane receptor with a negative regulatory role in Hedgehog signaling pathway. The PTCH1 germline mutations cause Gorlin syndrome, disorder characterized by developmental abnormalities and tumor susceptibility. Most of the malformations are caused by PTCH1 haploinsufficiency, indicative of fine-tuning needed to properly regulate the activity of pathway which is involved in pathogenesis of various tumors. Since in patients with different tumors and healthy controls we identified 5 to 8 CGG repeats located 4 bases upstream of a translation initiation site, in this study we wanted to examine how 5’ untranslated region (5’UTR) regulates the expression of PTCH1 transcript 1b. Material and methods Various in silico tools and databases were used to reveal all potential cis-regulatory elements in the PTCH1b 5’UTR. Since PTCH1b transcript has two different-sized 5’UTRs, we built pGL3-P-based plasmids by inserting upstream of firefly luciferase gene either 188- or 300bp-long 5’UTR, each harboring 5 to 8 CGG repeats. As the last 76bp of 5’UTR were predicted as an internal ribosome entry site (IRES), we constructed bicistronic pRuF vectors by cloning each PTCH1b 5’UTR between Renilla and firefly luciferase gene. Reporter gene assays and qPCR were performed in transfected MCF-7, HCT 116 and HEK 293T cells. Results Dual luciferase assays showed that shorter 5’UTR significantly increased reporter activity, with a subtle reduction with increased number of repeats. Longer 5’UTR led to much reduced reporter activity, without a difference among repeats. Luciferase mRNA quantification showed that both 5’UTR lengths significantly increased transcription. Site-directed mutagenesis proved hypothesis that 2 potential upstream open reading frames, contained in first 112bp of longer 5’UTR, might account for this severe reduction in reporter activity. Both 5’UTR lengths significantly increased firefly luciferase activity of pRuF vectors (proved by PCR this is not due to alternative splicing), with no difference among repeats. Firefly luciferase activity was significantly reduced when predicted IRES motif was removed from pRuF plasmid. Firefly/Renilla luciferase mRNA ratios were the same as for empty vector, indicating that observed higher firefly luciferase activity with equal mRNA levels should be due to a post-transcriptional level of regulation, i.e., cap-independent translation of firefly luciferase gene. Conclusions All our results point to the exceptionally complex and so far unexplored role of 5’UTR in the regulation of PTCH1b expression, whilst the existence of an IRES motif would enable Ptch1 protein to be synthesized under conditions when the general level of protein synthesis is reduced, such as in hypoxia, which is known activator of Hedgehog signaling pathway.