General amino acid control in fission yeast is regulated by a nonconserved transcription factor, with functions analogous to Gcn4/Atf4.

Eukaryotes respond to amino acid starvation by enhancing the translation of mRNAs encoding b-ZIP family transcription factors (GCN4 in Saccharomyces cerevisiae and ATF4 in mammals), which launch transcriptional programs to counter this stress. This pathway involves phosphorylation of the eIF2 translation factor by Gcn2-protein kinases and is regulated by upstream ORFs (uORFs) in the GCN4/ATF4 5' leaders. Here, we present evidence that the transcription factors that mediate this response are not evolutionarily conserved. Although cells of the fission yeast Schizosaccharomyces pombe respond transcriptionally to amino acid starvation, they lack clear Gcn4 and Atf4 orthologs. We used ribosome profiling to identify mediators of this response in S. pombe, looking for transcription factors that behave like GCN4 We discovered a transcription factor (Fil1) translationally induced by amino acid starvation in a 5' leader and Gcn2-dependent manner. Like Gcn4, Fil1 is required for the transcriptional response to amino acid starvation, and Gcn4 and Fil1 regulate similar genes. Despite their similarities in regulation, function, and targets, Fil1 and Gcn4 belong to different transcription factor families (GATA and b-ZIP, respectively). Thus, the same functions are performed by nonorthologous proteins under similar regulation. These results highlight the plasticity of transcriptional networks, which maintain conserved principles with nonconserved regulators.

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.

Cotranslational protein-RNA associations predict protein-protein interactions.

BACKGROUND: Most cellular proteins function as part of stable protein complexes. We recently showed that around 38% of proteins associate with mRNAs that encode interacting proteins, reflecting the cotranslational formation of the complex between the bait protein and the nascent peptides encoded by the interacting mRNAs. Here we hypothesise that these cotranslational protein-mRNA associations can be used to predict protein-protein interactions. RESULTS: We found that the fission yeast Exo2 protein, which encodes an exonuclease of the XRN1 family, coimmunoprecipitates with the eti1 mRNA, which codes for a protein of unknown function and uninformative sequence. Based on this protein-mRNA association, we predicted that the Exo2 and Eti1 protein are part of the same complex, and confirmed this hypothesis by coimmunoprecipitation and colocalization of the proteins. Similarly, we show that the cotranslational interaction between the Sty1 MAP kinase and the cip2 mRNA, which encodes an RNA-binding protein, predicts a complex between Sty1 and Cip2. CONCLUSIONS: Our results demonstrate that cotranslational protein-mRNA associations can be used to identify new components of protein complexes.

Widespread cotranslational formation of protein complexes.

Most cellular processes are conducted by multi-protein complexes. However, little is known about how these complexes are assembled. In particular, it is not known if they are formed while one or more members of the complexes are being translated (cotranslational assembly). We took a genomic approach to address this question, by systematically identifying mRNAs associated with specific proteins. In a sample of 31 proteins from Schizosaccharomyces pombe that did not contain RNA-binding domains, we found that ∼38% copurify with mRNAs that encode interacting proteins. For example, the cyclin-dependent kinase Cdc2p associates with the rum1 and cdc18 mRNAs, which encode, respectively, an inhibitor of Cdc2p kinase activity and an essential regulator of DNA replication. Both proteins interact with Cdc2p and are key cell cycle regulators. We obtained analogous results with proteins with different structures and cellular functions (kinesins, protein kinases, transcription factors, proteasome components, etc.). We showed that copurification of a bait protein and of specific mRNAs was dependent on the presence of the proteins encoded by the interacting mRNAs and on polysomal integrity. These results indicate that these observed associations reflect the cotranslational interaction between the bait and the nascent proteins encoded by the interacting mRNAs. Therefore, we show that the cotranslational formation of protein-protein interactions is a widespread phenomenon.

Sharing and archiving nucleic acid structure mapping data.

Nucleic acids are particularly amenable to structural characterization using chemical and enzymatic probes. Each individual structure mapping experiment reveals specific information about the structure and/or dynamics of the nucleic acid. Currently, there is no simple approach for making these data publically available in a standardized format. We therefore developed a standard for reporting the results of single nucleotide resolution nucleic acid structure mapping experiments, or SNRNASMs. We propose a schema for sharing nucleic acid chemical probing data that uses generic public servers for storing, retrieving, and searching the data. We have also developed a consistent nomenclature (ontology) within the Ontology of Biomedical Investigations (OBI), which provides unique identifiers (termed persistent URLs, or PURLs) for classifying the data. Links to standardized data sets shared using our proposed format along with a tutorial and links to templates can be found at http://snrnasm.bio.unc.edu.

Nonhierarchical ribonucleoprotein assembly suggests a strain-propagation model for protein-facilitated RNA folding.

Proteins play diverse and critical roles in cellular ribonucleoproteins (RNPs) including promoting formation of and stabilizing active RNA conformations. Yet, the conformational changes required to convert large RNAs into active RNPs have proven difficult to characterize fully. Here we use high-resolution approaches to monitor both local nucleotide flexibility and solvent accessibility for nearly all nucleotides in the bI3 group I intron RNP in four assembly states: the free RNA, maturase-bound RNA, Mrs1-bound RNA, and the complete six-component holocomplex. The free RNA is misfolded relative to the secondary structure required for splicing. The maturase and Mrs1 proteins each stabilized long-range tertiary interactions, but neither protein alone induced folding into the functional secondary structure. In contrast, simultaneous binding by both proteins results in large secondary structure rearrangements in the RNA and yielded the catalytically active group I intron structure. Secondary and tertiary folding of the RNA component of the bI3 RNP are thus not independent: RNA folding is strongly nonhierarchical. These results emphasize that protein-mediated stabilization of RNA tertiary interactions functions to pull the secondary structure into an energetically disfavored, but functional, conformation and emphasize a new role for facilitator proteins in RNP assembly.

Proteomic analysis of meiosis and characterization of novel short open reading frames in the fission yeast Schizosaccharomyces pombe.

Meiosis is the process by which haploid gametes are produced from diploid precursor cells. We used stable isotope labeling by amino acids in cell culture (SILAC) to characterize the meiotic proteome in the fission yeast Schizosaccharomyces pombe. We compared relative levels of proteins extracted from cells harvested around meiosis I with those of meiosis II, and proteins from premeiotic S phase with the interval between meiotic divisions, when S phase is absent. Our proteome datasets revealed peptides corresponding to short open reading frames (sORFs) that have been previously identified by ribosome profiling as new translated regions. We verified expression of selected sORFs by Western blotting and analyzed the phenotype of deletion mutants. Our data provide a resource for studying meiosis that may help understand differences between meiosis I and meiosis II and how S phase is suppressed between the two meiotic divisions.

SHAPE analysis of long-range interactions reveals extensive and thermodynamically preferred misfolding in a fragile group I intron RNA.

Most functional RNAs require proteins to facilitate formation of their active structures. In the case of the yeast bI3 group I intron, splicing requires binding by two proteins, the intron-encoded bI3 maturase and the nuclear encoded Mrs1. Here, we use selective 2'-hydroxyl acylation analyzed by primer extension (SHAPE) chemistry coupled with analysis of point mutants to map long-range interactions in this RNA. This analysis reveals two critical features of the free RNA state. First, the catalytic intron is separated from the flanking exons via a stable anchoring helix. This anchoring helix creates an autonomous structural domain for the intron and functions to prevent misfolding with the flanking exons. Second, the thermodynamically most stable structure for the free RNA is not consistent with the catalytically active conformation as phylogenetically conserved elements form stable, non-native structures. These results highlight a fragile bI3 RNA for which binding by protein cofactors functions to promote extensive secondary structure rearrangements that are an obligatory prerequisite for forming the catalytically active tertiary structure.

Effects of cycloheximide on the interpretation of ribosome profiling experiments in Schizosaccharomyces pombe.

Stress conditions lead to global and gene-specific changes in RNA translation. Ribosome profiling experiments have identified genome-wide alterations in the distribution of ribosomes along mRNAs. However, it is contentious whether these changes reflect real responses, or whether they are artefacts caused by the use of inhibitors of translation (notably cycloheximide). To address this issue we performed ribosome profiling with the fission yeast Schizosaccharomyces pombe under conditions of exponential growth (unstressed) and nitrogen starvation (nutritional stress), and both in the presence and absence of cycloheximide. We examined several aspects of the translational response, including density of ribosomal footprints on coding sequences, 5' leader ribosomal densities, distribution of ribosomes along coding sequences, and ribosome codon occupancies. Cycloheximide had minor effects on overall ribosome density, which affected mostly mRNAs encoding ribosomal proteins. Nitrogen starvation caused an accumulation of ribosomes on 5' leaders in both cycloheximide-treated and untreated cells. By contrast, stress-induced ribosome accumulation on the 5' side of coding sequences was cycloheximide-dependent. Finally, codon occupancy showed strong positive correlations in cycloheximide-treated and untreated cells. Our results demonstrate that cycloheximide does influence some of the results of ribosome profiling experiments, although it is not clear if this effect is always artefactual.

AnABlast: a new in silico strategy for the genome-wide search of novel genes and fossil regions.

Genome annotation, assisted by computer programs, is one of the great advances in modern biology. Nevertheless, the in silico identification of small and complex coding sequences is still challenging. We observed that amino acid sequences inferred from coding-but rarely from non-coding-DNA sequences accumulated alignments in low-stringency BLAST searches, suggesting that this alignments accumulation could be used to highlight coding regions in sequenced DNA. To investigate this possibility, we developed a computer program (AnABlast) that generates profiles of accumulated alignments in query amino acid sequences using a low-stringency BLAST strategy. To validate this approach, all six-frame translations of DNA sequences between every two annotated exons of the fission yeast genome were analysed with AnABlast. AnABlast-generated profiles identified three new copies of known genes, and four new genes supported by experimental evidence. New pseudogenes, ancestral carboxyl- and amino-terminal subtractions, complex gene rearrangements, and ancient fragments of mitDNA and of bacterial origin, were also inferred. Thus, this novel in silico approach provides a powerful tool to uncover new genes, as well as fossil-coding sequences, thus providing insight into the evolutionary history of annotated genomes.

The translational landscape of fission-yeast meiosis and sporulation.

Sexual development in Schizosaccharomyces pombe culminates in meiosis and sporulation. We used ribosome profiling to investigate the translational landscape of this process. We show that the translation efficiency of hundreds of genes is regulated in complex patterns, often correlating with changes in RNA levels. Ribosome-protected fragments show a three-nucleotide periodicity that identifies translated sequences and their reading frame. Using this property, we identified 46 new translated genes and found that 24% of noncoding RNAs are actively translated. We also detected 19 nested antisense genes, in which both DNA strands encode translated mRNAs. Finally, we identified 1,735 translated upstream open reading frames (ORFs) in leader sequences. In S. pombe, in contrast with Saccharomyces cerevisiae, sexual development is not accompanied by large increases in upstream ORF use, thus suggesting that this is an organism-specific adaptation, not a general feature of developmental processes.

The Mrs1 splicing factor binds the bI3 group I intron at each of two tetraloop-receptor motifs.

Most large ribozymes require protein cofactors in order to function efficiently. The yeast mitochondrial bI3 group I intron requires two proteins for efficient splicing, Mrs1 and the bI3 maturase. Mrs1 has evolved from DNA junction resolvases to function as an RNA cofactor for at least two group I introns; however, the RNA binding site and the mechanism by which Mrs1 facilitates splicing were unknown. Here we use high-throughput RNA structure analysis to show that Mrs1 binds a ubiquitous RNA tertiary structure motif, the GNRA tetraloop-receptor interaction, at two sites in the bI3 RNA. Mrs1 also interacts at similar tetraloop-receptor elements, as well as other structures, in the self-folding Azoarcus group I intron and in the RNase P enzyme. Thus, Mrs1 recognizes general features found in the tetraloop-receptor motif. Identification of the two Mrs1 binding sites now makes it possible to create a model of the complete six-component bI3 ribonucleoprotein. All protein cofactors bind at the periphery of the RNA such that every long-range RNA tertiary interaction is stabilized by protein binding, involving either Mrs1 or the bI3 maturase. This work emphasizes the strong evolutionary pressure to bolster RNA tertiary structure with RNA-binding interactions as seen in the ribosome, spliceosome, and other large RNA machines.

High-throughput SHAPE and hydroxyl radical analysis of RNA structure and ribonucleoprotein assembly.

RNA folds to form complex structures vital to many cellular functions. Proteins facilitate RNA folding at both the secondary and tertiary structure levels. An absolute prerequisite for understanding RNA folding and ribonucleoprotein (RNP) assembly reactions is a complete understanding of the RNA structure at each stage of the folding or assembly process. Here we provide a guide for comprehensive and high-throughput analysis of RNA secondary and tertiary structure using SHAPE and hydroxyl radical footprinting. As an example of the strong and sometimes surprising conclusions that can emerge from high-throughput analysis of RNA folding and RNP assembly, we summarize the structure of the bI3 group I intron RNA in four distinct states. Dramatic structural rearrangements occur in both secondary and tertiary structure as the RNA folds from the free state to the active, six-component, RNP complex. As high-throughput and high-resolution approaches are applied broadly to large protein-RNA complexes, other proteins previously viewed as making simple contributions to RNA folding are also likely to be found to exert multifaceted, long-range, cooperative, and nonadditive effects on RNA folding. These protein-induced contributions add another level of control, and potential regulatory function, in RNP complexes.