Identification of RNA recognition elements in the Saccharomyces cerevisiae transcriptome

A non-canonical Puf3p-binding sequence regulates CAT5/COQ7 mRNA under both fermentable and respiratory conditions in budding yeast

The Saccharomyces cerevisiae uses a highly glycolytic metabolism, if glucose is available, through appropriately suppressing mitochondrial functions except for some of them such as Fe/S cluster biogenesis. Puf3p, a Pumillio family protein, plays a pivotal role in modulating mitochondrial activity, especially during fermentation, by destabilizing its target mRNAs and/or by repressing their translation. Puf3p preferentially binds to 8-nt conserved binding sequences in the 3'-UTR of nuclear-encoded mitochondrial (nc-mitochondrial) mRNAs, leading to broad effects on gene expression under fermentable conditions. To further explore how Puf3p post-transcriptionally regulates nc-mitochondrial mRNAs in response to cell growth conditions, we initially focused on nc-mitochondrial mRNAs known to be enriched in monosomes in a glucose-rich environment. We unexpectedly found that one of the monosome-enriched mRNAs, CAT5/COQ7 mRNA, directly interacts with Puf3p through its non-canonical Puf3p binding sequence, which is generally less considered as a Puf3p binding site. Western blot analysis showed that Puf3p represses translation of Cat5p, regardless of culture in fermentable or respiratory medium. In vitro binding assay confirmed Puf3p's direct interaction with CAT5 mRNA via this non-canonical Puf3p-binding site. Although cat5 mutants of the non-canonical Puf3p-binding site grow normally, Cat5p expression is altered, indicating that CAT5 mRNA is a bona fide Puf3p target with additional regulatory factors acting through this sequence. Unlike other yeast PUF proteins, Puf3p uniquely regulates Cat5p by destabilizing mRNA and repressing translation, shedding new light on an unknown part of the Puf3p regulatory network. Given that pathological variants of human COQ7 lead to CoQ10 deficiency and yeast cat5Δ can be complemented by hCOQ7, our findings may also offer some insights into clinical aspects of COQ7-related disorders.

Dynamical control enables the formation of demixed biomolecular condensates

Cellular matter can be organized into compositionally distinct biomolecular condensates. For example, in Ashbya gossypii, the RNA-binding protein Whi3 forms distinct condensates with different RNA molecules. Using criteria derived from a physical framework for explaining how compositionally distinct condensates can form spontaneously via thermodynamic considerations, we find that condensates in vitro form mainly via heterotypic interactions in binary mixtures of Whi3 and RNA. However, within these condensates, RNA molecules become dynamically arrested. As a result, in ternary systems, simultaneous additions of Whi3 and pairs of distinct RNA molecules lead to well-mixed condensates, whereas delayed addition of an RNA component results in compositional distinctness. Therefore, compositional identities of condensates can be achieved via dynamical control, being driven, at least partially, by the dynamical arrest of RNA molecules. Finally, we show that synchronizing the production of different RNAs leads to more well-mixed, as opposed to compositionally distinct condensates in vivo.

Oxidative stress induces coordinated remodeling of RNA-enzyme interactions

RNA-binding proteins (RBPs) are key post-transcriptional regulators that play a substantial role during stress adaptation. Recent proteome-wide surveys have uncovered a large number of new and “unconventional” RBPs such as metabolic enzymes, yet little is known about the reconfiguration of the RNA-binding proteome (RBPome) and RNA-enzyme interactions in response to cellular stress. Here, we applied RNA-interactome capture to monitor the dynamics of the mRBPome upon mild oxidative stress in the yeast Saccharomyces cerevisiae. Among the 257 proteins that significantly changed RNA associations, we observed the coordinated remodeling of RNA-binding enzymes — particularly of the central carbon metabolism — that complemented known metabolic responses. Furthermore, we recognized the propensity for paralogous specific alterations of enzyme-RNA interactions. Our results suggest coordinated cross talk between RNA-enzyme interactions and intermediary metabolism to maintain the physiological and molecular balance upon oxidative stress, perhaps through specialization of paralogous during evolution.

•

Oxidative stress induces the rearrangement of 257 proteins on polyadenylated RNAs

•

Coordinated response of RNA-enzyme interactions and metabolism

•

Yeast RNA-binding enzymes are paralog specific

•

Integration of three different mass spectrometry analysis tools

Stress-Induced Translation Inhibition through Rapid Displacement of Scanning Initiation Factors

Cellular responses to environmental stress are frequently mediated by RNA-binding proteins (RBPs). Here, we examined global RBP dynamics in Saccharomyces cerevisiae in response to glucose starvation and heat shock. Each stress induced rapid remodeling of the RNA-protein interactome without corresponding changes in RBP abundance. Consistent with general translation shutdown, ribosomal proteins contacting the mRNA showed decreased RNA association. Among translation components, RNA association was most reduced for initiation factors involved in 40S scanning (eukaryotic initiation factor 4A [eIF4A], eIF4B, and Ded1), indicating a common mechanism of translational repression. In unstressed cells, eIF4A, eIF4B, and Ded1 primarily targeted the 5′ ends of mRNAs. Following glucose withdrawal, 5′ binding was abolished within 30 s, explaining the rapid translation shutdown, but mRNAs remained stable. Heat shock induced progressive loss of 5′ RNA binding by initiation factors over ∼16 min and provoked mRNA degradation, particularly for translation-related factors, mediated by Xrn1. Taken together, these results reveal mechanisms underlying translational control of gene expression during stress.

•

A quantitative proteomic approach reveals global stress-induced changes in RNA binding

•

Translation shutdown is driven by rapid loss of mRNA binding by key initiation factors

•

Heat shock induces general mRNA degradation facilitated by Xrn1

Widespread Transcriptional Readthrough Caused by Nab2 Depletion Leads to Chimeric Transcripts with Retained Introns

Nascent RNA sequencing has revealed that pre-mRNA splicing can occur shortly after introns emerge from RNA polymerase II (RNA Pol II). Differences in co-transcriptional splicing profiles suggest regulation by cis- and/or trans-acting factors. Here, we use single-molecule intron tracking (SMIT) to identify a cohort of regulators by machine learning in budding yeast. Of these, Nab2 displays reduced co-transcriptional splicing when depleted. Unexpectedly, these splicing defects are attributable to aberrant "intrusive" transcriptional readthrough from upstream genes, as revealed by long-read sequencing. Transcripts that originate from the intron-containing gene's own transcription start site (TSS) are efficiently spliced, indicating no direct role of Nab2 in splicing per se. This work highlights the coupling between transcription, splicing, and 3' end formation in the context of gene organization along chromosomes. We conclude that Nab2 is required for proper 3' end processing, which ensures gene-specific control of co-transcriptional RNA processing.

Gcn2 eIF2α kinase mediates combinatorial translational regulation through nucleotide motifs and uORFs in target mRNAs

The protein kinase Gcn2 is a central transducer of nutritional stress signaling important for stress adaptation by normal cells and the survival of cancer cells. In response to nutrient deprivation, Gcn2 phosphorylates eIF2α, thereby repressing general translation while enhancing translation of specific mRNAs with upstream ORFs (uORFs) situated in their 5'-leader regions. Here we performed genome-wide measurements of mRNA translation during histidine starvation in fission yeast Schizosaccharomyces pombe. Polysome analyses were combined with microarray measurements to identify gene transcripts whose translation was up-regulated in response to the stress in a Gcn2-dependent manner. We determined that translation is reprogrammed to enhance RNA metabolism and chromatin regulation and repress ribosome synthesis. Interestingly, translation of intron-containing mRNAs was up-regulated. The products of the regulated genes include additional eIF2α kinase Hri2 amplifying the stress signaling and Gcn5 histone acetyl transferase and transcription factors, together altering genome-wide transcription. Unique dipeptide-coding uORFs and nucleotide motifs, such as '5'-UGA(C/G)GG-3', are found in 5' leader regions of regulated genes and shown to be responsible for translational control.

5' untranslated regions: the next regulatory sequence in yeast synthetic biology

When developing industrial biotechnology processes, Saccharomyces cerevisiae (baker's yeast or brewer's yeast) is a popular choice as a microbial host. Many tools have been developed in the fields of synthetic biology and metabolic engineering to introduce heterologous pathways and tune their expression in yeast. Such tools mainly focus on controlling transcription, whereas post-transcriptional regulation is often overlooked. Herein we discuss regulatory elements found in the 5' untranslated region (UTR) and their influence on protein synthesis. We provide not only an overall picture, but also a set of design rules on how to engineer a 5' UTR. The reader is also referred to currently available models that allow gene expression to be tuned predictably using different 5' UTRs.

Linking RNA Sequence, Structure, and Function on Massively Parallel High-Throughput Sequencers

High-throughput sequencing methods have revolutionized our ability to catalog the diversity of RNAs and RNA-protein interactions that can exist in our cells. However, the relationship between RNA sequence, structure, and function is enormously complex, demonstrating the need for methods that can provide quantitative thermodynamic and kinetic measurements of macromolecular interaction with RNA, at a scale commensurate with the sequence diversity of RNA. Here, we discuss a class of methods that extend the core functionality of DNA sequencers to enable high-throughput measurements of RNA folding and RNA-protein interactions. Topics discussed include a description of the method and multiple applications to RNA-binding proteins, riboswitch design and engineering, and RNA tertiary structure energetics.

A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

A Quantitative and Predictive Model for RNA Binding by Human Pumilio Proteins

High-throughput methodologies have enabled routine generation of RNA target sets and sequence motifs for RNA-binding proteins (RBPs). Nevertheless, quantitative approaches are needed to capture the landscape of RNA-RBP interactions responsible for cellular regulation. We have used the RNA-MaP platform to directly measure equilibrium binding for thousands of designed RNAs and to construct a predictive model for RNA recognition by the human Pumilio proteins PUM1 and PUM2. Despite prior findings of linear sequence motifs, our measurements revealed widespread residue flipping and instances of positional coupling. Application of our thermodynamic model to published in vivo crosslinking data reveals quantitative agreement between predicted affinities and in vivo occupancies. Our analyses suggest a thermodynamically driven, continuous Pumilio-binding landscape that is negligibly affected by RNA structure or kinetic factors, such as displacement by ribosomes. This work provides a quantitative foundation for dissecting the cellular behavior of RBPs and cellular features that impact their occupancies.

Novel insights into global translational regulation through Pumilio family RNA-binding protein Puf3p revealed by ribosomal profiling

RNA binding proteins (RBPs) can regulate the stability, localization, and translation of their target mRNAs. Among them, Puf3p is a well-known Pumilio family RBP whose biology has been intensively studied. Nevertheless, the impact of Puf3p on the translational regulation of its downstream genes still remains to be investigated at the genome-wide level. In this study, we combined ribosome profiling and RNA-Seq in budding yeast (Saccharomyces cerevisiae) to investigate Puf3p's functions in translational regulation. Comparison of translational efficiency (TE) between wild-type and puf3Δ strains demonstrates extensive translational modulation in the absence of Puf3p (over 27% genes are affected at the genome level). Besides confirming its known role in regulating mitochondrial metabolism, our data demonstrate that Puf3p serves as a key post-transcriptional regulator of downstream RBPs by regulating their translational efficiencies, indicating a network of interactions among RBPs at the post-transcriptional level. Furthermore, Puf3p switches the balance of translational flux between mitochondrial and cytosolic ribosome biogenesis to adapt to changes in cellular metabolism. In summary, our results indicate that TE can be utilized as an informative index to interrogate the mechanism underlying RBP functions, and provide novel insights into Puf3p's mode-of-action.

Structural basis for the dimerization of Nab2 generated by RNA binding provides insight into its contribution to both poly(A) tail length determination and transcript compaction in Saccharomyces cerevisiae

In Saccharomyces cerevisiae generation of export-competent mRNPs terminates the nuclear phase of the gene expression pathway and facilitates transport to the cytoplasm for translation. Nab2 functions in this process to control both mRNP compaction that facilitates movement through nuclear pore complexes and the length of transcript poly(A) tails. Nab2 has a modular structure that includes seven CCCH Zn fingers that bind to A-rich RNAs and fingers 5–7 are critical for these functions. Here, we demonstrate, using both biophysical and structural methods, that binding A₁₁G RNA induces dimerization of Zn fingers 5–7 mediated by the novel spatial arrangement of the fingers promoting each RNA chain binding two protein chains. The dimerization of Nab2 induced by RNA binding provides a basis for understanding its function in both poly(A) tail length regulation and in the compaction of mature transcripts to facilitate nuclear export.

[Effect of nucleolin silencing on differentiation of rat neural stem cells in vitro and the molecular mechanism]

OBJECTIVE

To investigate the effect of nucleolin silencing on the differentiation of rat neural stem cells (NSCs) and the role of Wnt signaling pathway in mediating such effect.

METHODS

Adenovirus vectors expressing small interfering RNA (siRNA) against nucleolin were constructed, verified, and packaged in HEK293A cells. The adenovirus was then transfected into NSCs isolated from neonatal SD rats and the differentiation of the NSCs was examined by detecting the expressions of neuron specific encloase (NSE) and glial fibrillary acidic protein (GFAP) using immunocytochemistry. The expressions of nucleolin, nestin, Wnt3, and β-catenin in the cells were determined with Western blotting.

RESULTS

Restriction endonuclease and sequencing analysis verified successful construction of the adenoviral vector expressing nucleolin siRNA (nucleolin-siRNA2). Infection of rat NSCs with nucleolin-siRNA2 significantly lowered nucleolin protein expression as compared with that in negative and blank control groups (P<0.05). The percentages of NSE-positive cells and GFAP-positive cells were significantly higher in NSCs infected with nucleolin-siRNA (P<0.01); the infection also resulted in obviously lowered expression of nestin protein and increased expressions of Wnt3 protein and β-catenin nucleoprotein in the cells.

CONCLUSIONS

Nucleolin silencing by adenovirus-mediated RNA interference induces the differentiation of NSCs into neurons and astrocytes, which is related with the activation of Wnt signaling pathway.

Analysis of mRNA deadenylation by multi-protein complexes

Poly(A) tails are found at the 3' end of almost every eukaryotic mRNA and are important for the stability of mRNAs and their translation into proteins. Thus, removal of the poly(A) tail, a process called deadenylation, is critical for regulation of gene expression. Most deadenylation enzymes are components of large multi-protein complexes. Here, we describe an in vitro deadenylation assay developed to study the exonucleolytic activities of the multi-protein Ccr4-Not and Pan2-Pan3 complexes. We discuss how this assay can be used with short synthetic RNAs, as well as longer RNA substrates generated using in vitro transcription. Importantly, quantitation of the reactions allows detailed analyses of deadenylation in the presence and absence of accessory factors, leading to new insights into targeted mRNA decay.

Comprehensive and quantitative mapping of RNA-protein interactions across a transcribed eukaryotic genome

RNA-binding proteins (RBPs) control the fate of nearly every transcript in a cell. However, no existing approach for studying these posttranscriptional gene regulators combines transcriptome-wide throughput and biophysical precision. Here, we describe an assay that accomplishes this. Using commonly available hardware, we built a customizable, open-source platform that leverages the inherent throughput of Illumina technology for direct biophysical measurements. We used the platform to quantitatively measure the binding affinity of the prototypical RBP Vts1 for every transcript in the Saccharomyces cerevisiae genome. The scale and precision of these measurements revealed many previously unknown features of this well-studied RBP. Our transcribed genome array (TGA) assayed both rare and abundant transcripts with equivalent proficiency, revealing hundreds of low-abundance targets missed by previous approaches. These targets regulated diverse biological processes including nutrient sensing and the DNA damage response, and implicated Vts1 in de novo gene "birth." TGA provided single-nucleotide resolution for each binding site and delineated a highly specific sequence and structure motif for Vts1 binding. Changes in transcript levels in vts1Δ cells established the regulatory function of these binding sites. The impact of Vts1 on transcript abundance was largely independent of where it bound within an mRNA, challenging prevailing assumptions about how this RBP drives RNA degradation. TGA thus enables a quantitative description of the relationship between variant RNA structures, affinity, and in vivo phenotype on a transcriptome-wide scale. We anticipate that TGA will provide similarly comprehensive and quantitative insights into the function of virtually any RBP.

Nucleolin-Mediated RNA Localization Regulates Neuron Growth and Cycling Cell Size

How can cells sense their own size to coordinate biosynthesis and metabolism with their growth needs? We recently proposed a motor-dependent bidirectional transport mechanism for axon length and cell size sensing, but the nature of the motor-transported size signals remained elusive. Here, we show that motor-dependent mRNA localization regulates neuronal growth and cycling cell size. We found that the RNA-binding protein nucleolin is associated with importin β1 mRNA in axons. Perturbation of nucleolin association with kinesins reduces its levels in axons, with a concomitant reduction in axonal importin β1 mRNA and protein levels. Strikingly, subcellular sequestration of nucleolin or importin β1 enhances axonal growth and causes a subcellular shift in protein synthesis. Similar findings were obtained in fibroblasts. Thus, subcellular mRNA localization regulates size and growth in both neurons and cycling cells.

Premature termination of GAT1 transcription explains paradoxical negative correlation between nitrogen-responsive mRNA, but constitutive low-level protein production

The first step in executing the genetic program of a cell is production of mRNA. In yeast, almost every gene is transcribed as multiple distinct isoforms, differing at their 5′ and/or 3′ termini. However, the implications and functional significance of the transcriptome-wide diversity of mRNA termini remains largely unexplored. In this paper, we show that the GAT1 gene, encoding a transcriptional activator of nitrogen-responsive catabolic genes, produces a variety of mRNAs differing in their 5′ and 3′ termini. Alternative transcription initiation leads to the constitutive, low level production of 2 full length proteins differing in their N-termini, whereas premature transcriptional termination generates a short, highly nitrogen catabolite repression- (NCR-) sensitive transcript that, as far as we can determine, is not translated under the growth conditions we used, but rather likely protects the cell from excess Gat1.

Translational Capacity of a Cell Is Determined during Transcription Elongation via the Ccr4-Not Complex

The current understanding of gene expression considers transcription and translation to be independent processes. Challenging this notion, we found that translation efficiency is determined during transcription elongation through the imprinting of mRNAs with Not1, the central scaffold of the Ccr4-Not complex. We determined that another subunit of the complex, Not5, defines Not1 binding to specific mRNAs, particularly those produced from ribosomal protein genes. This imprinting mechanism specifically regulates ribosomal protein gene expression, which in turn determines the translational capacity of cells. We validate our model by SILAC and polysome profiling experiments. As a proof of concept, we demonstrate that enhanced translation compensates for transcriptional elongation stress. Taken together, our data indicate that in addition to defining mRNA stability, components of the Ccr4-Not imprinting complex regulate RNA translatability, thus ensuring global gene expression homeostasis.

Mapping the Transcriptome-Wide Landscape of RBP Binding Sites Using gPAR-CLIP-seq: Experimental Procedures

An estimated 5-10 % of protein-coding genes in eukaryotic genomes encode RNA-binding proteins (RBPs). Through dynamic changes in RNA recognition, RBPs posttranscriptionally regulate the biogenesis, transport, inheritance, storage, and degradation of RNAs. Understanding such widespread RBP-mediated posttranscriptional regulatory mechanisms requires comprehensive discovery of the in vivo binding sites of RBPs. Here, we describe the experimental procedures of the gPAR-CLIP-seq (global photoactivatable-ribonucleoside-enhanced cross-linking and precipitation followed by deep sequencing) approach we recently developed for capturing and sequencing regions of the transcriptome bound by RBPs in budding yeast. Unlike the standard PAR-CLIP method, which identifies the bound RNA substrates for a single RBP, the gPAR-CLIP-seq method was developed to isolate and sequence all mRNA sites bound by the cellular "RBPome." The gPAR-CLIP-seq approach is readily applicable to a variety of organisms and cell lines to profile global RNA-protein interactions underlying posttranscriptional gene regulation. The complete landscape of RBP binding sites provides insights to the function of all RNA cis-regulatory elements in an organism and reveals fundamental mechanisms of posttranscriptional gene regulation.

The cellular growth rate controls overall mRNA turnover, and modulates either transcription or degradation rates of particular gene regulons

We analyzed 80 different genomic experiments, and found a positive correlation between both RNA polymerase II transcription and mRNA degradation with growth rates in yeast. Thus, in spite of the marked variation in mRNA turnover, the total mRNA concentration remained approximately constant. Some genes, however, regulated their mRNA concentration by uncoupling mRNA stability from the transcription rate. Ribosome-related genes modulated their transcription rates to increase mRNA levels under fast growth. In contrast, mitochondria-related and stress-induced genes lowered mRNA levels by reducing mRNA stability or the transcription rate, respectively. We also detected these regulations within the heterogeneity of a wild-type cell population growing in optimal conditions. The transcriptomic analysis of sorted microcolonies confirmed that the growth rate dictates alternative expression programs by modulating transcription and mRNA decay.

The regulation of overall mRNA turnover keeps a constant ratio between mRNA decay and the dilution of [mRNA] caused by cellular growth. This regulation minimizes the indiscriminate transmission of mRNAs from mother to daughter cells, and favors the response capacity of the latter to physiological signals and environmental changes. We also conclude that, by uncoupling mRNA synthesis from decay, cells control the mRNA abundance of those gene regulons that characterize fast and slow growth.

Snf1-Dependent Transcription Confers Glucose-Induced Decay upon the mRNA Product

In the yeast Saccharomyces cerevisiae, the switch from respiratory metabolism to fermentation causes rapid decay of transcripts encoding proteins uniquely required for aerobic metabolism. Snf1, the yeast ortholog of AMP-activated protein kinase, has been implicated in this process because inhibiting Snf1 mimics the addition of glucose. In this study, we show that the SNF1-dependent ADH2 promoter, or just the major transcription factor binding site, is sufficient to confer glucose-induced mRNA decay upon heterologous transcripts. SNF1-independent expression from the ADH2 promoter prevented glucose-induced mRNA decay without altering the start site of transcription. SNF1-dependent transcripts are enriched for the binding motif of the RNA binding protein Vts1, an important mediator of mRNA decay and mRNA repression whose expression is correlated with decreased abundance of SNF1-dependent transcripts during the yeast metabolic cycle. However, deletion of VTS1 did not slow the rate of glucose-induced mRNA decay. ADH2 mRNA rapidly dissociated from polysomes after glucose repletion, and sequences bound by RNA binding proteins were enriched in the transcripts from repressed cells. Inhibiting the protein kinase A pathway did not affect glucose-induced decay of ADH2 mRNA. Our results suggest that Snf1 may influence mRNA stability by altering the recruitment activity of the transcription factor Adr1.

Gbp2 interacts with THO/TREX through a novel type of RRM domain

Metazoan SR and SR-like proteins are important regulatory factors in RNA splicing, export, translation and RNA decay. We determined the NMR structures and nucleic acid interaction modes of Gbp2 and Hrb1, two paralogous budding yeast proteins with similarities to mammalian SR proteins. Gbp2 RRM1 and RRM2 recognise preferentially RNAs containing the core motif GGUG. Sequence selectivity resides in a non-canonical interface in RRM2 that is highly related to the SRSF1 pseudoRRM. The atypical Gbp2/Hrb1 C-terminal RRM domains (RRM3) do not interact with RNA/DNA, likely because of their novel N-terminal extensions that block the canonical RNA binding interface. Instead, we discovered that RRM3 is crucial for interaction with the THO/TREX complex and identified key residues essential for this interaction. Moreover, Gbp2 interacts genetically with Tho2 as the double deletion shows a synthetic phenotype and preventing Gbp2 interaction with the THO/TREX complex partly supresses gene expression defect associated with inactivation of the latter complex. These findings provide structural and functional insights into the contribution of SR-like proteins in the post-transcriptional control of gene expression.

Evolutionary Conservation and Diversification of Puf RNA Binding Proteins and Their mRNA Targets

Reprogramming of a gene’s expression pattern by acquisition and loss of sequences recognized by specific regulatory RNA binding proteins may be a major mechanism in the evolution of biological regulatory programs. We identified that RNA targets of Puf3 orthologs have been conserved over 100–500 million years of evolution in five eukaryotic lineages. Focusing on Puf proteins and their targets across 80 fungi, we constructed a parsimonious model for their evolutionary history. This model entails extensive and coordinated changes in the Puf targets as well as changes in the number of Puf genes and alterations of RNA binding specificity including that: 1) Binding of Puf3 to more than 200 RNAs whose protein products are predominantly involved in the production and organization of mitochondrial complexes predates the origin of budding yeasts and filamentous fungi and was maintained for 500 million years, throughout the evolution of budding yeast. 2) In filamentous fungi, remarkably, more than 150 of the ancestral Puf3 targets were gained by Puf4, with one lineage maintaining both Puf3 and Puf4 as regulators and a sister lineage losing Puf3 as a regulator of these RNAs. The decrease in gene expression of these mRNAs upon deletion of Puf4 in filamentous fungi (N. crassa) in contrast to the increase upon Puf3 deletion in budding yeast (S. cerevisiae) suggests that the output of the RNA regulatory network is different with Puf4 in filamentous fungi than with Puf3 in budding yeast. 3) The coregulated Puf4 target set in filamentous fungi expanded to include mitochondrial genes involved in the tricarboxylic acid (TCA) cycle and other nuclear-encoded RNAs with mitochondrial function not bound by Puf3 in budding yeast, observations that provide additional evidence for substantial rewiring of post-transcriptional regulation. 4) Puf3 also expanded and diversified its targets in filamentous fungi, gaining interactions with the mRNAs encoding the mitochondrial electron transport chain (ETC) complex I as well as hundreds of other mRNAs with nonmitochondrial functions. The many concerted and conserved changes in the RNA targets of Puf proteins strongly support an extensive role of RNA binding proteins in coordinating gene expression, as originally proposed by Keene. Rewiring of Puf-coordinated mRNA targets and transcriptional control of the same genes occurred at different points in evolution, suggesting that there have been distinct adaptations via RNA binding proteins and transcription factors. The changes in Puf targets and in the Puf proteins indicate an integral involvement of RNA binding proteins and their RNA targets in the adaptation, reprogramming, and function of gene expression.

A map of the evolutionary history of Puf proteins and their RNA targets shows that reprogramming of global gene expression programs via adaptive mutations that affect protein-RNA interactions is an important source of biological diversity.

We set out to trace the evolutionary history of an RNA binding protein and how its interactions with targets change over evolution. Identifying this natural history is a step toward understanding the critical differences between organisms and how gene expression programs are rewired during evolution. Using bioinformatics and experimental approaches, we broadly surveyed the evolution of binding targets of a particular family of RNA binding proteins—the Puf proteins, whose protein sequences and target RNA sequences are relatively well-characterized—across 99 eukaryotic species. We found five groups of species in which targets have been conserved for at least 100 million years and then took advantage of genome sequences from a large number of fungal species to deeply investigate the conservation and changes in Puf proteins and their RNA targets. Our analyses identified multiple and extensive reconfigurations during the natural history of fungi and suggest that RNA binding proteins and their RNA targets are profoundly involved in evolutionary reprogramming of gene expression and help define distinct programs unique to each organism. Continuing to uncover the natural history of RNA binding proteins and their interactions will provide a unique window into the gene expression programs of present day species and point to new ways to engineer gene expression programs.

Integration of Posttranscriptional Gene Networks into Metabolic Adaptation and Biofilm Maturation in Candida albicans

The yeast Candida albicans is a human commensal and opportunistic pathogen. Although both commensalism and pathogenesis depend on metabolic adaptation, the regulatory pathways that mediate metabolic processes in C. albicans are incompletely defined. For example, metabolic change is a major feature that distinguishes community growth of C. albicans in biofilms compared to suspension cultures, but how metabolic adaptation is functionally interfaced with the structural and gene regulatory changes that drive biofilm maturation remains to be fully understood. We show here that the RNA binding protein Puf3 regulates a posttranscriptional mRNA network in C. albicans that impacts on mitochondrial biogenesis, and provide the first functional data suggesting evolutionary rewiring of posttranscriptional gene regulation between the model yeast Saccharomyces cerevisiae and C. albicans. A proportion of the Puf3 mRNA network is differentially expressed in biofilms, and by using a mutant in the mRNA deadenylase CCR4 (the enzyme recruited to mRNAs by Puf3 to control transcript stability) we show that posttranscriptional regulation is important for mitochondrial regulation in biofilms. Inactivation of CCR4 or dis-regulation of mitochondrial activity led to altered biofilm structure and over-production of extracellular matrix material. The extracellular matrix is critical for antifungal resistance and immune evasion, and yet of all biofilm maturation pathways extracellular matrix biogenesis is the least understood. We propose a model in which the hypoxic biofilm environment is sensed by regulators such as Ccr4 to orchestrate metabolic adaptation, as well as the regulation of extracellular matrix production by impacting on the expression of matrix-related cell wall genes. Therefore metabolic changes in biofilms might be intimately linked to a key biofilm maturation mechanism that ultimately results in untreatable fungal disease.

Metabolism is a master regulator of cell biology, including gene regulation, developmental switches and cellular life-death decisions, with the mitochondrion playing a central role in eukaryotes. For the yeast Candida albicans mitochondrial functions have been implicated in host-pathogen interactions, but the regulatory mechanism that control mitochondrial biogenesis are poorly described. We identified the RNA binding protein Puf3 as a new mitochondrial regulator in C. albicans, and show that posttranscriptional regulation and mitochondrial function have important roles during community growth in biofilms. Perturbation of mitochondrial activity or inactivation of a key posttranscriptional regulator, CCR4, led to changes in biofilm maturation, shedding light on the interface between metabolic reprogramming and biofilm developmental pathways. We illuminate a new mechanism that regulates extracellular matrix production, an essential biofilm feature that mediates the notorious drug resistance and immune evasion properties of the biofilm growth mode.

A conserved mechanism of TOR-dependent RCK-mediated mRNA degradation regulates autophagy

Autophagy is an essential eukaryotic pathway requiring tight regulation to maintain homeostasis and preclude disease. Using yeast and mammalian cells, we report a conserved mechanism of autophagy regulation by RNA helicase RCK family members in association with the decapping enzyme Dcp2. Under nutrient-replete conditions, Dcp2 undergoes TOR-dependent phosphorylation and associates with RCK members to form a complex with autophagy-related (ATG) mRNA transcripts, leading to decapping, degradation and autophagy suppression. Simultaneous with the induction of ATG mRNA synthesis, starvation reverses the process, facilitating ATG mRNA accumulation and autophagy induction. This conserved post-transcriptional mechanism modulates fungal virulence and the mammalian inflammasome, the latter providing mechanistic insight into autoimmunity reported in a patient with a PIK3CD/p110δ gain-of-function mutation. We propose a dynamic model wherein RCK family members, in conjunction with Dcp2, function in controlling ATG mRNA stability to govern autophagy, which in turn modulates vital cellular processes affecting inflammation and microbial pathogenesis.

The human RBPome: from genes and proteins to human disease

RNA binding proteins (RBPs) play a central role in mediating post transcriptional regulation of genes. However less is understood about them and their regulatory mechanisms. In this study, we construct a catalogue of 1344 experimentally confirmed RBPs. The domain architecture of RBPs enabled us to classify them into three groups - Classical (29%), Non-classical (19%) and unclassified (52%). A higher percentage of proteins with unclassified domains reveals the presence of various uncharacterised motifs that can potentially bind RNA. RBPs were found to be highly disordered compared to Non-RBPs (p<2.2e-16, Fisher's exact test), suggestive of a dynamic regulatory role of RBPs in cellular signalling and homeostasis. Evolutionary analysis in 62 different species showed that RBPs are highly conserved compared to Non-RBPs (p<2.2e-16, Wilcox-test), reflecting the conservation of various biological processes like mRNA splicing and ribosome biogenesis. The expression patterns of RBPs from human proteome map revealed that ~40% of them are ubiquitously expressed and ~60% are tissue-specific. RBPs were also seen to be highly associated with several neurological disorders, cancer and inflammatory diseases. Anatomical contexts like B cells, T-cells, foetal liver and foetal brain were found to be strongly enriched for RBPs, implying a prominent role of RBPs in immune responses and different developmental stages. The catalogue and meta-analysis presented here should form a foundation for furthering our understanding of RBPs and the cellular networks they control, in years to come. This article is part of a Special Issue entitled: Proteomics in India.

Termination of Transcription of Short Noncoding RNAs by RNA Polymerase II

The RNA polymerase II transcription cycle is often divided into three major stages: initiation, elongation, and termination. Research over the last decade has blurred these divisions and emphasized the tightly regulated transitions that occur as RNA polymerase II synthesizes a transcript from start to finish. Transcription termination, the process that marks the end of transcription elongation, is regulated by proteins that interact with the polymerase, nascent transcript, and/or chromatin template. The failure to terminate transcription can cause accumulation of aberrant transcripts and interfere with transcription at downstream genes. Here, we review the mechanism, regulation, and physiological impact of a termination pathway that targets small noncoding transcripts produced by RNA polymerase II. We emphasize the Nrd1-Nab3-Sen1 pathway in yeast, in which the process has been extensively studied. The importance of understanding small RNA termination pathways is underscored by the need to control noncoding transcription in eukaryotic genomes.

Co-transcriptional mRNP formation is coordinated within a molecular mRNP packaging station in S. cerevisiae

In eukaryotes, the messenger RNA (mRNA), the blueprint of a protein‐coding gene, is processed and packaged into a messenger ribonucleoprotein particle (mRNP) by mRNA‐binding proteins in the nucleus. The steps of mRNP formation – transcription, processing, packaging, and the orchestrated release of the export‐competent mRNP from the site of transcription for nuclear mRNA export – are tightly coupled to ensure a highly efficient and regulated process. The importance of highly accurate nuclear mRNP formation is illustrated by the fact that mutations in components of this pathway lead to cellular inviability or to severe diseases in metazoans. We hypothesize that efficient mRNP formation is realized by a molecular mRNP packaging station, which is built by several recruitment platforms and coordinates the individual steps of mRNP formation.

Global mRNA selection mechanisms for translation initiation

Background

The selection and regulation of individual mRNAs for translation initiation from a competing pool of mRNA are poorly understood processes. The closed loop complex, comprising eIF4E, eIF4G and PABP, and its regulation by 4E-BPs are perceived to be key players. Using RIP-seq, we aimed to evaluate the role in gene regulation of the closed loop complex and 4E-BP regulation across the entire yeast transcriptome.

Results

We find that there are distinct populations of mRNAs with coherent properties: one mRNA pool contains many ribosomal protein mRNAs and is enriched specifically with all of the closed loop translation initiation components. This class likely represents mRNAs that rely heavily on the closed loop complex for protein synthesis. Other heavily translated mRNAs are apparently under-represented with most closed loop components except Pab1p. Combined with data showing a close correlation between Pab1p interaction and levels of translation, these data suggest that Pab1p is important for the translation of these mRNAs in a closed loop independent manner. We also identify a translational regulatory mechanism for the 4E-BPs; these appear to self-regulate by inhibiting translation initiation of their own mRNAs.

Conclusions

Overall, we show that mRNA selection for translation initiation is not as uniformly regimented as previously anticipated. Components of the closed loop complex are highly relevant for many mRNAs, but some heavily translated mRNAs interact poorly with this machinery. Therefore, alternative, possibly Pab1p-dependent mechanisms likely exist to load ribosomes effectively onto mRNAs. Finally, these studies identify and characterize a complex self-regulatory circuit for the yeast 4E-BPs.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-014-0559-z) contains supplementary material, which is available to authorized users.

High-throughput characterization of protein-RNA interactions

RNA-binding proteins (RBPs) are important regulators of eukaryotic gene expression. Genomes typically encode dozens to hundreds of proteins containing RNA-binding domains, which collectively recognize diverse RNA sequences and structures. Recent advances in high-throughput methods for assaying the targets of RBPs in vitro and in vivo allow large-scale derivation of RNA-binding motifs as well as determination of RNA–protein interactions in living cells. In parallel, many computational methods have been developed to analyze and interpret these data. The interplay between RNA secondary structure and RBP binding has also been a growing theme. Integrating RNA–protein interaction data with observations of post-transcriptional regulation will enhance our understanding of the roles of these important proteins.

Systematic analysis of the role of RNA-binding proteins in the regulation of RNA stability

mRNA half-lives are transcript-specific and vary over a range of more than 100-fold in eukaryotic cells. mRNA stabilities can be regulated by sequence-specific RNA-binding proteins (RBPs), which bind to regulatory sequence elements and modulate the interaction of the mRNA with the cellular RNA degradation machinery. However, it is unclear if this kind of regulation is sufficient to explain the large range of mRNA stabilities. To address this question, we examined the transcriptome of 74 Schizosaccharomyces pombe strains carrying deletions in non-essential genes encoding predicted RBPs (86% of all such genes). We identified 25 strains that displayed changes in the levels of between 4 and 104 mRNAs. The putative targets of these RBPs formed biologically coherent groups, defining regulons involved in cell separation, ribosome biogenesis, meiotic progression, stress responses and mitochondrial function. Moreover, mRNAs in these groups were enriched in specific sequence motifs in their coding sequences and untranslated regions, suggesting that they are coregulated at the posttranscriptional level. We performed genome-wide RNA stability measurements for several RBP mutants, and confirmed that the altered mRNA levels were caused by changes in their stabilities. Although RBPs regulate the decay rates of multiple regulons, only 16% of all S. pombe mRNAs were affected in any of the 74 deletion strains. This suggests that other players or mechanisms are required to generate the observed range of RNA half-lives of a eukaryotic transcriptome.

Coupling mRNA synthesis and decay

What has been will be again, what has been done will be done again; there is nothing new under the sun. -Ecclesiastes 1:9 (New International Version) Posttranscriptional regulation of gene expression has an important role in defining the phenotypic characteristics of an organism. Well-defined steps in mRNA metabolism that occur in the nucleus-capping, splicing, and polyadenylation-are mechanistically linked to the process of transcription. Recent evidence suggests another link between RNA polymerase II (Pol II) and a posttranscriptional process that occurs in the cytoplasm-mRNA decay. This conclusion appears to represent a conundrum. How could mRNA synthesis in the nucleus and mRNA decay in the cytoplasm be mechanistically linked? After a brief overview of mRNA processing, we will review the recent evidence for transcription-coupled mRNA decay and the possible involvement of Snf1, the Saccharomyces cerevisiae ortholog of AMP-activated protein kinase, in this process.

Contributions of transcription and mRNA decay to gene expression dynamics of fission yeast in response to oxidative stress

The cooperation of transcriptional and post-transcriptional levels of control to shape gene regulation is only partially understood. Here we show that a combination of two simple and non-invasive genomic techniques, coupled with kinetic mathematical modeling, afford insight into the intricate dynamics of RNA regulation in response to oxidative stress in the fission yeast Schizosaccharomyces pombe. This study reveals a dominant role of transcriptional regulation in response to stress, but also points to the first minutes after stress induction as a critical time when the coordinated control of mRNA turnover can support the control of transcription for rapid gene regulation. In addition, we uncover specialized gene expression strategies associated with distinct functional gene groups, such as simultaneous transcriptional repression and mRNA destabilization for genes encoding ribosomal proteins, delayed mRNA destabilization with varying contribution of transcription for ribosome biogenesis genes, dominant roles of mRNA stabilization for genes functioning in protein degradation, and adjustment of both transcription and mRNA turnover during the adaptation to stress. We also show that genes regulated independently of the bZIP transcription factor Atf1p are predominantly controlled by mRNA turnover, and identify putative cis-regulatory sequences that are associated with different gene expression strategies during the stress response. This study highlights the intricate and multi-faceted interplay between transcription and RNA turnover during the dynamic regulatory response to stress.

Harnessing natural sequence variation to dissect posttranscriptional regulatory networks in yeast

Understanding how genomic variation influences phenotypic variation through the molecular networks of the cell is one of the central challenges of biology. Transcriptional regulation has received much attention, but equally important is the posttranscriptional regulation of mRNA stability. Here we applied a systems genetics approach to dissect posttranscriptional regulatory networks in the budding yeast Saccharomyces cerevisiae. Quantitative sequence-to-affinity models were built from high-throughput in vivo RNA binding protein (RBP) binding data for 15 yeast RBPs. Integration of these models with genome-wide mRNA expression data allowed us to estimate protein-level RBP regulatory activity for individual segregants from a genetic cross between two yeast strains. Treating these activities as a quantitative trait, we mapped trans-acting loci (activity quantitative trait loci, or aQTLs) that act via posttranscriptional regulation of transcript stability. We predicted and experimentally confirmed that a coding polymorphism at the IRA2 locus modulates Puf4p activity. Our results also indicate that Puf3p activity is modulated by distinct loci, depending on whether it acts via the 5′ or the 3′ untranslated region of its target mRNAs. Together, our results validate a general strategy for dissecting the connectivity between posttranscriptional regulators and their upstream signaling pathways.

Alternative polyadenylation diversifies post-transcriptional regulation by selective RNA-protein interactions

Recent research has uncovered extensive variability in the boundaries of transcript isoforms, yet the functional consequences of this variation remain largely unexplored. Here, we systematically discriminate between the molecular phenotypes of overlapping coding and non‐coding transcriptional events from each genic locus using a novel genome‐wide, nucleotide‐resolution technique to quantify the half‐lives of 3′ transcript isoforms in yeast. Our results reveal widespread differences in stability among isoforms for hundreds of genes in a single condition, and that variation of even a single nucleotide in the 3′ untranslated region (UTR) can affect transcript stability. While previous instances of negative associations between 3′ UTR length and transcript stability have been reported, here, we find that shorter isoforms are not necessarily more stable. We demonstrate the role of RNA‐protein interactions in conditioning isoform‐specific stability, showing that PUF3 binds and destabilizes specific polyadenylation isoforms. Our findings indicate that although the functional elements of a gene are encoded in DNA sequence, the selective incorporation of these elements into RNA through transcript boundary variation allows a single gene to have diverse functional consequences.

Global regulation of mRNA translation and stability in the early Drosophila embryo by the Smaug RNA-binding protein

Background

Smaug is an RNA-binding protein that induces the degradation and represses the translation of mRNAs in the early Drosophila embryo. Smaug has two identified direct target mRNAs that it differentially regulates: nanos and Hsp83. Smaug represses the translation of nanos mRNA but has only a modest effect on its stability, whereas it destabilizes Hsp83 mRNA but has no detectable effect on Hsp83 translation. Smaug is required to destabilize more than one thousand mRNAs in the early embryo, but whether these transcripts represent direct targets of Smaug is unclear and the extent of Smaug-mediated translational repression is unknown.

Results

To gain a panoramic view of Smaug function in the early embryo, we identified mRNAs that are bound to Smaug using RNA co-immunoprecipitation followed by hybridization to DNA microarrays. We also identified mRNAs that are translationally repressed by Smaug using polysome gradients and microarrays. Comparison of the bound mRNAs to those that are translationally repressed by Smaug and those that require Smaug for their degradation suggests that a large fraction of Smaug’s target mRNAs are both translationally repressed and degraded by Smaug. Smaug directly regulates components of the TRiC/CCT chaperonin, the proteasome regulatory particle and lipid droplets, as well as many metabolic enzymes, including several glycolytic enzymes.

Conclusions

Smaug plays a direct and global role in regulating the translation and stability of a large fraction of the mRNAs in the early Drosophila embryo, and has unanticipated functions in control of protein folding and degradation, lipid droplet function and metabolism.

Computational analysis of conserved RNA secondary structure in transcriptomes and genomes

Transcriptomics experiments and computational predictions both enable systematic discovery of new functional RNAs. However, many putative noncoding transcripts arise instead from artifacts and biological noise, and current computational prediction methods have high false positive rates. I discuss prospects for improving computational methods for analyzing and identifying functional RNAs, with a focus on detecting signatures of conserved RNA secondary structure. An interesting new front is the application of chemical and enzymatic experiments that probe RNA structure on a transcriptome-wide scale. I review several proposed approaches for incorporating structure probing data into the computational prediction of RNA secondary structure. Using probabilistic inference formalisms, I show how all these approaches can be unified in a well-principled framework, which in turn allows RNA probing data to be easily integrated into a wide range of analyses that depend on RNA secondary structure inference. Such analyses include homology search and genome-wide detection of new structural RNAs.

Effects of the yeast RNA-binding protein Whi3 on the half-life and abundance of CLN3 mRNA and other targets

Whi3 is an RNA binding protein known to bind the mRNA of the yeast G1 cyclin gene CLN3. It inhibits CLN3 function, but the mechanism of this inhibition is unclear; in previous studies, Whi3 made no observable difference to CLN3 mRNA levels, translation, or protein abundance. Here, we re-approach this issue using microarrays, RNA-Seq, ribosome profiling, and other methods. By multiple methods, we find that the whi3 mutation causes a small but consistent increase in the abundance of hundreds of mRNAs, including the CLN3 mRNA. The effect on various mRNAs is roughly in proportion to the density of GCAU or UGCAU motifs carried by these mRNAs, which may be a binding site for Whi3. mRNA instability of Whi3 targets may in part depend on a 3′ AU rich element (ARE), AUUUUA. In addition, the whi3 mutation causes a small increase in the translational efficiency of CLN3 mRNA. The increase in CLN3 mRNA half-life and abundance together with the increase in translational efficiency is fully sufficient to explain the small-cell phenotype of whi3 mutants. Under stress conditions, Whi3 becomes a component of P-bodies or stress granules, but Whi3 also acts under non-stress condition, when no P-bodies are visible. We suggest that Whi3 may be a very broadly-acting, but mild, modulator of mRNA stability. In CLN3, Whi3 may bind to the 3′ GCAU motifs to attract the Ccr4-Not complex to promote RNA deadenylation and turnover, and Whi3 may bind to the 5′ GCAU motifs to inhibit translation.

Whi3, an S. cerevisiae RNA-binding protein, is a component of stress granules that regulates levels of its target mRNAs

RNA binding proteins (RBPs) are vital to the regulation of mRNA transcripts, and can alter mRNA localization, degradation, translation, and storage. Whi3 was originally identified in a screen for small cell size mutants, and has since been characterized as an RBP. The identification of Whi3-interacting mRNAs involved in mediating cellular responses to stress suggested that Whi3 might be involved in stress-responsive RNA processing. We show that Whi3 localizes to stress granules in response to glucose deprivation or heat shock. The kinetics and pattern of Whi3 localization in response to a range of temperatures were subtly but distinctly different from those of known components of RNA processing granules. Deletion of Whi3 resulted in an increase in the relative abundance of Whi3 target RNAs, either in the presence or absence of heat shock. Increased levels of the CLN3 mRNA in whi3Δ cells may explain their decreased cell size. Another mRNA target of Whi3 encodes the zinc-responsive transcription factor Zap1, suggesting a role for Whi3 in response to zinc stress. Indeed, we found that whi3Δ cells have enhanced sensitivity to zinc toxicity. Together our results suggest an expanded model for Whi3 function: in addition to its role as a regulator of the cell cycle, Whi3 may have a role in stress-dependent RNA processing and responses to a variety of stress conditions.

Carbon source-dependent alteration of Puf3p activity mediates rapid changes in the stabilities of mRNAs involved in mitochondrial function

The Puf family of RNA-binding proteins regulates gene expression primarily by interacting with the 3′ untranslated region (3′ UTR) of targeted mRNAs and inhibiting translation and/or stimulating decay. Physical association and computational analyses of yeast Puf3p identified >150 potential mRNA targets involved in mitochondrial function. However, only COX17 has been established as a target of Puf3p-mediated deadenylation and decapping. We have identified 10 new targets that are rapidly degraded in a Puf3p-dependent manner. We also observed changes in Puf3p activity in response to environmental conditions. Puf3p promotes rapid degradation of mRNA targets in the fermentable carbon source dextrose. However, Puf3p-mediated decay activity is inhibited in carbon sources that require mitochondrial function for efficient cell growth. In addition, the activity of Puf3p is rapidly altered by changing the carbon source. PUF3 expression is not decreased at the RNA or protein level by different carbon sources and localization is not significantly altered, suggesting that Puf3p activity is regulated posttranslationally. Finally, under conditions when Puf3p is unable to stimulate decay, Puf3p can still bind its target mRNAs. Together, these experiments provide insight into the carbon source-specific control of Puf3p activity and how such alterations allow Puf3p to dynamically regulate mitochondrial function.

Finding the target sites of RNA-binding proteins

RNA–protein interactions differ from DNA–protein interactions because of the central role of RNA secondary structure. Some RNA-binding domains (RBDs) recognize their target sites mainly by their shape and geometry and others are sequence-specific but are sensitive to secondary structure context. A number of small- and large-scale experimental approaches have been developed to measure RNAs associated in vitro and in vivo with RNA-binding proteins (RBPs). Generalizing outside of the experimental conditions tested by these assays requires computational motif finding. Often RBP motif finding is done by adapting DNA motif finding methods; but modeling secondary structure context leads to better recovery of RBP-binding preferences. Genome-wide assessment of mRNA secondary structure has recently become possible, but these data must be combined with computational predictions of secondary structure before they add value in predicting in vivo binding. There are two main approaches to incorporating structural information into motif models: supplementing primary sequence motif models with preferred secondary structure contexts (e.g., MEMERIS and RNAcontext) and directly modeling secondary structure recognized by the RBP using stochastic context-free grammars (e.g., CMfinder and RNApromo). The former better reconstruct known binding preferences for sequence-specific RBPs but are not suitable for modeling RBPs that recognize shape and geometry of RNAs. Future work in RBP motif finding should incorporate interactions between multiple RBDs and multiple RBPs in binding to RNA. WIREs RNA 2014, 5:111–130. doi: 10.1002/wrna.1201

Puf3p induces translational repression of genes linked to oxidative stress

In response to stress, the translation of many mRNAs in yeast can change in a fashion discordant with the general repression of translation. Here, we use machine learning to mine the properties of these mRNAs to determine specific translation control signals. We find a strong association between transcripts acutely translationally repressed under oxidative stress and those associated with the RNA-binding protein Puf3p, a known regulator of cellular mRNAs encoding proteins targeted to mitochondria. Under oxidative stress, a PUF3 deleted strain exhibits more robust growth than wild-type cells and the shift in translation from polysomes to monosomes is attenuated, suggesting puf3Δ cells perceive less stress. In agreement, the ratio of reduced:oxidized glutathione, a major antioxidant and indicator of cellular redox state, is increased in unstressed puf3Δ cells but remains lower under stress. In untreated conditions, Puf3p migrates with polysomes rather than ribosome-free fractions, but this is lost under stress. Finally, reverse transcriptase-polymerase chain reaction (RT-PCR) of Puf3p targets following affinity purification shows Puf3p-mRNA associations are maintained or increased under oxidative stress. Collectively, these results point to Puf3p acting as a translational repressor in a manner exceeding the global translational response, possibly by temporarily limiting synthesis of new mitochondrial proteins as cells adapt to the stress.

A network of interdependent molecular interactions describes a higher order Nrd1-Nab3 complex involved in yeast transcription termination

Nab3 and Nrd1 are yeast heterogeneous nuclear ribonucleoprotein (hnRNP)-like proteins that heterodimerize and bind RNA. Genetic and biochemical evidence reveals that they are integral to the termination of transcription of short non-coding RNAs by RNA polymerase II. Here we define a Nab3 mutation (nab3Δ134) that removes an essential part of the protein's C terminus but nevertheless can rescue, in trans, the phenotype resulting from a mutation in the RNA recognition motif of Nab3. This low complexity region of Nab3 appears intrinsically unstructured and can form a hydrogel in vitro. These data support a model in which multiple Nrd1-Nab3 heterodimers polymerize onto substrate RNA to effect termination, allowing complementation of one mutant Nab3 molecule by another lacking a different function. The self-association property of Nab3 adds to the previously documented interactions between these hnRNP-like proteins, RNA polymerase II, and the nascent transcript, leading to a network of nucleoprotein interactions that define a higher order Nrd1-Nab3 complex. This was underscored from the synthetic phenotypes of yeast strains with pairwise combinations of Nrd1 and Nab3 mutations known to affect their distinct biochemical activities. The mutations included a Nab3 self-association defect, a Nab3-Nrd1 heterodimerization defect, a Nrd1-polymerase II binding defect, and an Nab3-RNA recognition motif mutation. Although no single mutation was lethal, cells with any two mutations were not viable for four such pairings, and a fifth displayed a synthetic growth defect. These data strengthen the idea that a multiplicity of interactions is needed to assemble a higher order Nrd1-Nab3 complex that coats specific nascent RNAs in preparation for termination.

Quantification of mRNA stability of stress-responsive yeast genes following conditional excision of open reading frames

Eukaryotic cells rapidly adjust the levels of mRNAs in response to environmental stress primarily by controlling transcription and mRNA turnover. How different stress conditions influence the fate of stress-responsive mRNAs, however, is relatively poorly understood. This is largely due to the fact that mRNA half-life assays are traditionally based on interventions (e.g., temperature-shifts using temperature-sensitive RNA polymerase II alleles or treatment with general transcription inhibitory drugs), which, rather than blocking, specifically induce transcription of stress-responsive genes. To study the half-lives of the latter suite of mRNAs, we developed and describe here a minimally perturbing alternative method, coined CEO, which is based on discontinuance of transcription following the conditional excision of open reading frames. Using CEO, we confirm that the target of rapamycin complex I (TORC1), a nutrient-activated, central stimulator of eukaryotic cell growth, favors the decay of mRNAs that depend on the stress- and/or nutrient-regulated transcription factors Msn2/4 and Gis1 for their transcription. We further demonstrate that TORC1 controls the stability of these mRNAs via the Rim15-Igo1/2-PP2A^Cdc55 effector branch, which reportedly also controls Gis1 promoter recruitment. These data pinpoint PP2A^Cdc55 as a central node in homo-directional coordination of transcription and post-transcriptional mRNA stabilization of a specific array of nutrient-regulated genes.

Extensive transcriptional heterogeneity revealed by isoform profiling

Transcript function is determined by sequence elements arranged on an individual RNA molecule. Variation in transcripts can affect messenger RNA stability, localization and translation, or produce truncated proteins that differ in localization or function. Given the existence of overlapping, variable transcript isoforms, determining the functional impact of the transcriptome requires identification of full-length transcripts, rather than just the genomic regions that are transcribed. Here, by jointly determining both transcript ends for millions of RNA molecules, we reveal an extensive layer of isoform diversity previously hidden among overlapping RNA molecules. Variation in transcript boundaries seems to be the rule rather than the exception, even within a single population of yeast cells. Over 26 major transcript isoforms per protein-coding gene were expressed in yeast. Hundreds of short coding RNAs and truncated versions of proteins are concomitantly encoded by alternative transcript isoforms, increasing protein diversity. In addition, approximately 70% of genes express alternative isoforms that vary in post-transcriptional regulatory elements, and tandem genes frequently produce overlapping or even bicistronic transcripts. This extensive transcript diversity is generated by a relatively simple eukaryotic genome with limited splicing, and within a genetically homogeneous population of cells. Our findings have implications for genome compaction, evolution and phenotypic diversity between single cells. These data also indicate that isoform diversity as well as RNA abundance should be considered when assessing the functional repertoire of genomes.

Protein targeting to subcellular organelles via MRNA localization

Cells have complex membranous organelles for the compartmentalization and the regulation of most intracellular processes. Organelle biogenesis and maintenance requires newly synthesized proteins, each of which needs to go from the ribosome translating its mRNA to the correct membrane for insertion or transclocation to an a organellar subcompartment. Decades of research have revealed how proteins are targeted to the correct organelle and translocated across one or more organelle membranes ro the compartment where they function. The paradigm examples involve interactions between a peptide sequence in the protein, localization factors, and various membrane embedded translocation machineries. Membrane translocation is either cotranslational or posttranslational depending on the protein and target organelle. Meanwhile research in embryos, neurons and yeast revealed an alternative targeting mechanism in which the mRNA is localized and only then translated to synthesize the protein in the correct location. In these cases, the targeting information is coded by the cis-acting sequences in the mRNA ("Zipcodes") that interact with localization factors and, in many cases, are transported by the molecular motors on the cytoskeletal filaments. Recently, evidence has been found for this "mRNA based" mechanism in organelle protein targeting to endoplasmic reticulum, mitochondria, and the photosynthetic membranes within chloroplasts. Here we review known and potential roles of mRNA localization in protein targeting to and within organelles. This article is part of a Special Issue entitled: Protein Import and Quality Control in Mitochondria and Plastids.

Pervasive and dynamic protein binding sites of the mRNA transcriptome in Saccharomyces cerevisiae

Background

Protein-RNA interactions are integral components of nearly every aspect of biology, including regulation of gene expression, assembly of cellular architectures, and pathogenesis of human diseases. However, studies in the past few decades have only uncovered a small fraction of the vast landscape of the protein-RNA interactome in any organism, and even less is known about the dynamics of protein-RNA interactions under changing developmental and environmental conditions.

Results

Here, we describe the gPAR-CLIP (global photoactivatable-ribonucleoside-enhanced crosslinking and immunopurification) approach for capturing regions of the untranslated, polyadenylated transcriptome bound by RNA-binding proteins (RBPs) in budding yeast. We report over 13,000 RBP crosslinking sites in untranslated regions (UTRs) covering 72% of protein-coding transcripts encoded in the genome, confirming 3' UTRs as major sites for RBP interaction. Comparative genomic analyses reveal that RBP crosslinking sites are highly conserved, and RNA folding predictions indicate that secondary structural elements are constrained by protein binding and may serve as generalizable modes of RNA recognition. Finally, 38% of 3' UTR crosslinking sites show changes in RBP occupancy upon glucose or nitrogen deprivation, with major impacts on metabolic pathways as well as mitochondrial and ribosomal gene expression.

Conclusions

Our study offers an unprecedented view of the pervasiveness and dynamics of protein-RNA interactions in vivo.

Extensive transcript diversity and novel upstream open reading frame regulation in yeast

To understand the diversity of transcripts in yeast (Saccharomyces cerevisiae) we analyzed the transcriptional landscapes for cells grown under 18 different environmental conditions. Each sample was analyzed using RNA-sequencing, and a total of 670,446,084 uniquely mapped reads and 377,263 poly-adenylated end tags were produced. Consistent with previous studies, we find that the majority of yeast genes are expressed under one or more different conditions. By directly comparing the 5′ and 3′ ends of the transcribed regions, we find extensive differences in transcript ends across many conditions, especially those of stationary phase, growth in grape juice, and salt stimulation, suggesting differential choice of transcription start and stop sites is pervasive in yeast. Relative to the exponential growth condition (i.e., YPAD), transcripts differing at the 5′ ends and 3′ ends are predicted to differ in their annotated start codon in 21 genes and their annotated stop codon in 63 genes. Many (431) upstream open reading frames (uORFs) are found in alternate 5′ ends and are significantly enriched in transcripts produced during the salt response. Mutational analysis of five genes with uORFs revealed that two sets of uORFs increase the expression of a reporter construct, indicating a role in activation which had not been reported previously, whereas two other uORFs decreased expression. In addition, RNA binding protein motifs are statistically enriched for alternate ends under many conditions. Overall, these results demonstrate enormous diversity of transcript ends, and that this heterogeneity is regulated under different environmental conditions. Moreover, transcript end diversity has important biological implications for the regulation of gene expression. In addition, our data also serve as a valuable resource for the scientific community.