In order to obtain a systems-level understanding of a complex biological system, detailed proteome information is essential. Despite great progress in proteomics technologies, thorough interrogation of the proteome from quantity-limited biological samples is hampered by inefficiencies during processing. To address these challenges, here we introduce a novel protocol using paramagnetic beads, termed Single-Pot Solid-Phase-enhanced Sample Preparation (SP3). SP3 provides a rapid and unbiased means of proteomic sample preparation in a single tube that facilitates ultrasensitive analysis by outperforming existing protocols in terms of efficiency, scalability, speed, throughput, and flexibility. To illustrate these benefits, characterization of 1,000 HeLa cells and single Drosophila embryos is used to establish that SP3 provides an enhanced platform for profiling proteomes derived from sub-microgram amounts of material. These data present a first view of developmental stage-specific proteome dynamics in Drosophila at a single-embryo resolution, permitting characterization of inter-individual expression variation. Together, the findings of this work position SP3 as a superior protocol that facilitates exciting new directions in multiple areas of proteomics ranging from developmental biology to clinical applications.
Recent studies have suggested that a sub-complex of RNA polymerase II composed of Rpb4 and Rpb7 couples the nuclear and cytoplasmic stages of gene expression by associating with newly made mRNAs in the nucleus, and contributing to their translation and degradation in the cytoplasm. Here we show by yeast two hybrid and co-immunoprecipitation experiments, followed by ribosome fractionation and fluorescent microscopy, that a subunit of the Ccr4-Not complex, Not5, is essential in the nucleus for the cytoplasmic functions of Rpb4. Not5 interacts with Rpb4; it is required for the presence of Rpb4 in polysomes, for interaction of Rpb4 with the translation initiation factor eIF3 and for association of Rpb4 with mRNAs. We find that Rpb7 presence in the cytoplasm and polysomes is much less significant than that of Rpb4, and that it does not depend upon Not5. Hence Not5-dependence unlinks the cytoplasmic functions of Rpb4 and Rpb7. We additionally determine with RNA immunoprecipitation and native gel analysis that Not5 is needed in the cytoplasm for the co-translational assembly of RNA polymerase II. This stems from the importance of Not5 for the association of the R2TP Hsp90 co-chaperone with polysomes translating RPB1 mRNA to protect newly synthesized Rpb1 from aggregation. Hence taken together our results show that Not5 interconnects translation and transcription.
Understanding the relationship between DNA sequence variation and phenotypic variation in complex or quantitative traits is one of the major challenges in modern biology. We are witnessing a deluge of DNA sequence information and association studies of genetic polymorphisms with phenotypes of interest in families and populations. In addition, it has become clear that large portions of eukaryotic genomes beyond protein-coding genes are transcribed, generating numerous noncoding RNA (ncRNA) molecules whose functions remain mostly unknown.DNA oligonucleotide microarrays constitute a powerful technology for studying the expression of genes in different organisms. The Saccharomyces cerevisiae tiling array presents a significant advance over previous array-based platforms. It has a high density of overlapping probes that start on average every 8 bp along each strand of the genome, enabling precise definition of transcript structure. Furthermore, the array includes probes specific for the polymorphic positions of another, distantly related yeast strain, allowing accurate measurement of allele-specific expression in a hybrid of the two strains. This technology thus allows high-resolution, quantitative, strand- and allele-specific measurements of transcription from a full eukaryotic genome. In this chapter, we describe the methods for extracting RNA, synthesizing first-strand cDNA, fragmenting, and labeling of samples for hybridization to the tiling array. Combining genome-wide information on variation in DNA sequence with variation in transcript structure and levels promises to increase our understanding of the genotype-to-phenotype relationship.
In eukaryotic cells, the production of cellular energy requires close interplay between nuclear and mitochondrial genomes. The mitochondrial genome is essential in that it encodes several genes involved in oxidative phosphorylation. Each cell contains several mitochondrial genome copies and mitochondrial DNA recombination is a widespread process occurring in plants, fungi, protists, and invertebrates. Saccharomyces cerevisiae has proved to be an excellent model to dissect mitochondrial biology. Several studies have focused on DNA recombination in this organelle, yet mostly relied on reporter genes or artificial systems. However, no complete mitochondrial recombination map has been released for any eukaryote so far. In the present work, we sequenced pools of diploids originating from a cross between two different S. cerevisiae strains to detect recombination events. This strategy allowed us to generate the first genome-wide map of recombination for yeast mitochondrial DNA. We demonstrated that recombination events are enriched in specific hotspots preferentially localized in non-protein-coding regions. Additionally, comparison of the recombination profiles of two different crosses showed that the genetic background affects hotspot localization and recombination rates. Finally, to gain insights into the mechanisms involved in mitochondrial recombination, we assessed the impact of individual depletion of four genes previously associated with this process. Deletion of NTG1 and MGT1 did not substantially influence the recombination landscape, alluding to the potential presence of additional regulatory factors. Our findings also revealed the loss of large mitochondrial DNA regions in the absence of MHR1, suggesting a pivotal role for Mhr1 in mitochondrial genome maintenance during mating. This study provides a comprehensive overview of mitochondrial DNA recombination in yeast and thus paves the way for future mechanistic studies of mitochondrial recombination and genome maintenance.
Hundreds of transcript isoforms with varying boundaries and alternative regulatory signals are transcribed from the genome, even in a genetically homogeneous population of cells. To study this transcriptional heterogeneity, we developed transcript isoform sequencing (TIF-seq), a method that allows the genome-wide profiling of full-length transcript isoforms defined by their exact 5' and 3' boundaries. TIF-seq entails the generation of full-length cDNA libraries, followed by their circularization and the sequencing of the junction fragments spanning the 5' and 3' transcript ends. By determining the respective co-occurrence of start and end sites of individual transcript molecules, TIF-seq can distinguish variations that conventional approaches for mapping single ends cannot, such as short abortive transcripts, bicistronic messages and overlapping transcripts that differ in lengths. The TIF-seq protocol we describe here can be applied to any eukaryotic organism (e.g., yeast, human), and it requires 6-10 d for generating TIF-seq libraries, 10 d for sequencing and 2-3 d for analysis.
The genetic basis of heritable traits has been studied for decades. Although recent mapping efforts have elucidated genetic determinants of transcript levels, mapping of protein abundance has lagged. Here, we analyze levels of 4084 GFP-tagged yeast proteins in the progeny of a cross between a laboratory and a wild strain using flow cytometry and high-content microscopy. The genotype of trans variants contributed little to protein level variation between individual cells but explained >50% of the variance in the population's average protein abundance for half of the GFP fusions tested. To map trans-acting factors responsible, we performed flow sorting and bulk segregant analysis of 25 proteins, finding a median of five protein quantitative trait loci (pQTLs) per GFP fusion. Further, we find that cis-acting variants predominate; the genotype of a gene and its surrounding region had a large effect on protein level six times more frequently than the rest of the genome combined. We present evidence for both shared and independent genetic control of transcript and protein abundance: More than half of the expression QTLs (eQTLs) contribute to changes in protein levels of regulated genes, but several pQTLs do not affect their cognate transcript levels. Allele replacements of genes known to underlie trans eQTL hotspots confirmed the correlation of effects on mRNA and protein levels. This study represents the first genome-scale measurement of genetic contribution to protein levels in single cells and populations, identifies more than a hundred trans pQTLs, and validates the propagation of effects associated with transcript variation to protein abundance.
Genomic analysis has revealed the existence of a large number of long noncoding RNAs (lncRNAs) with different functions in a variety of organisms, including yeast. Cells display dramatic changes of gene expression upon environmental changes. Upon osmostress, hundreds of stress-responsive genes are induced by the stress-activated protein kinase (SAPK) p38/Hog1. Using whole-genome tiling arrays, we found that Hog1 induces a set of lncRNAs upon stress. One of the genes expressing a Hog1-dependent lncRNA in antisense orientation is CDC28, the cyclin-dependent kinase 1 (CDK1) that controls the cell cycle in yeast. Cdc28 lncRNA mediates the establishment of gene looping and the relocalization of Hog1 and RSC from the 3' UTR to the +1 nucleosome to induce CDC28 expression. The increase in the levels of Cdc28 results in cells able to reenter the cell cycle more efficiently after stress. This may represent a general mechanism to prime expression of genes needed after stresses are alleviated.
Most genomes, including yeast Saccharomyces cerevisiae, are pervasively transcribed producing numerous non-coding RNAs, many of which are unstable and eliminated by nuclear or cytoplasmic surveillance pathways. We previously showed that accumulation of PHO84 antisense RNA (asRNA), in cells lacking the nuclear exosome component Rrp6, is paralleled by repression of sense transcription in a process dependent on the Hda1 histone deacetylase (HDAC) and the H3K4 histone methyl transferase Set1. Here we investigate this process genome-wide and measure the whole transcriptome of various histone modification mutants in a ?rrp6 strain using tiling arrays. We confirm widespread occurrence of potentially antisense-dependent gene regulation and identify three functionally distinct classes of genes that accumulate asRNAs in the absence of Rrp6. These classes differ in whether the genes are silenced by the asRNA and whether the silencing is HDACs and histone methyl transferase-dependent. Among the distinguishing features of asRNAs with regulatory potential, we identify weak early termination by Nrd1/Nab3/Sen1, extension of the asRNA into the open reading frame promoter and dependence of the silencing capacity on Set1 and the HDACs Hda1 and Rpd3 particularly at promoters undergoing extensive chromatin remodelling. Finally, depending on the efficiency of Nrd1/Nab3/Sen1 early termination, asRNA levels are modulated and their capability of silencing is changed.
For a unicellular, nonmotile organism like Saccharomyces cerevisiae, carbon sources act as nutrients and as signaling molecules; consequently, these sources affect various fitness parameters, including growth. It is therefore advantageous for yeast strains to adapt their growth to carbon source variation. The ability of a given genotype to manifest different phenotypes in varying environments is known as phenotypic plasticity. To identify quantitative trait loci (QTL) that drive plasticity in growth, two growth parameters (growth rate and biomass) were measured for a set of meiotic recombinants of two genetically divergent yeast strains grown in different carbon sources. To identify QTL contributing to plasticity across pairs of environments, gene-environment interaction mapping was performed, which identified several QTL that have a differential effect across environments, some of which act antagonistically across pairs of environments. Multi-QTL analysis identified loci interacting with previously known growth affecting QTL as well as novel two-QTL interactions that affect growth. A QTL that had no significant independent effect was found to alter growth rate and biomass for several carbon sources through two-QTL interactions. Our study demonstrates that environment-specific epistatic interactions contribute to the growth plasticity in yeast. We propose that a targeted scan for epistatic interactions, such as the one described here, can help unravel mechanisms regulating phenotypic plasticity.
The modern evolutionary synthesis assumes that mutations occur at random, independently of the environment in which they confer an advantage. However, there are indications that cells facing challenging conditions can adapt rapidly, utilizing processes beyond selection of pre-existing genetic variation. Here, we show that a strong regulatory challenge can induce mutations in many independent yeast cells, in the absence of general mutagenesis. Whole genome sequencing of cell lineages reveals a repertoire of independent mutations within a single lineage that arose only after the cells were exposed to the challenging environment, while other cells in the same lineage adapted without any mutation in their genomes. Thus, our experiments uncovered multiple alternative routes for heritable adaptation that were all induced in the same lineage during a short time period. Our results demonstrate the existence of adaptation mechanisms beyond random mutation, suggesting a tight connection between physiological and genetic processes.
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.
Dissecting the molecular basis of quantitative traits is a significant challenge, and is essential for understanding complex diseases. Even in model organisms, precisely determining causative genes and their interactions has remained elusive, due in part to difficulty in narrowing intervals to single genes, and in detecting epistasis or linked quantitative trait loci. These difficulties are exacerbated by limitations in experimental design, such as low numbers of analyzed individuals, and polymorphisms between parental genomes. We address these challenges by applying three independent high-throughput approaches for QTL mapping to map the genetic variants underlying eleven phenotypes in two genetically distant Saccharomyces cerevisiae strains, namely: 1) individual analysis of over 700 meiotic segregants, 2) bulk segregant analysis, and 3) reciprocal hemizygosity analysis, a new genome-wide method we developed. We identified differences in the performance of each approach and, by combining them, identified eight polymorphic genes that affect eight different phenotypes: colony shape, flocculation, growth on non-fermentable carbon sources, and resistance to drugs, salt, and heat. Our results demonstrate the power of individual segregant analysis to dissect quantitative trait loci and address the underestimated contribution of interactions between variants. We also reveal confounding factors like mutations and aneuploidy in pooled approaches, providing valuable lessons for future designs of complex trait mapping studies.
Antisense transcription, which was initially considered by many as transcriptional noise, is increasingly being recognized as an important regulator of gene expression. It is widespread among all kingdoms of life and has been shown to influence - either through the act of transcription or through the non-coding RNA that is produced - almost all stages of gene expression, from transcription and translation to RNA degradation. Antisense transcription can function as a fast evolving regulatory switch and a modular scaffold for protein complexes, and it can rewire regulatory networks. The genomic arrangement of antisense RNAs opposite sense genes indicates that they might be part of self-regulatory circuits that allow genes to regulate their own expression.
The majority of disease-associated variants lie outside protein-coding regions, suggesting a link between variation in regulatory regions and disease predisposition. We studied differences in chromatin states using five histone modifications, cohesin, and CTCF in lymphoblastoid lines from 19 individuals of diverse ancestry. We found extensive signal variation in regulatory regions, which often switch between active and repressed states across individuals. Enhancer activity is particularly diverse among individuals, whereas gene expression remains relatively stable. Chromatin variability shows genetic inheritance in trios, correlates with genetic variation and population divergence, and is associated with disruptions of transcription factor binding motifs. Overall, our results provide insights into chromatin variation among humans.
Alternative usage of exons provides genomes with plasticity to produce different transcripts from the same gene, modulating the function, localization, and life cycle of gene products. It affects most human genes. For a limited number of cases, alternative functions and tissue-specific roles are known. However, recent high-throughput sequencing studies have suggested that much alternative isoform usage across tissues is nonconserved, raising the question of the extent of its functional importance. We address this question in a genome-wide manner by analyzing the transcriptomes of five tissues for six primate species, focusing on exons that are 1:1 orthologous in all six species. Our results support a model in which differential usage of exons has two major modes: First, most of the exons show only weak differences, which are dominated by interspecies variability and may reflect neutral drift and noisy splicing. These cases dominate the genome-wide view and explain why conservation appears to be so limited. Second, however, a sizeable minority of exons show strong differences between tissues, which are mostly conserved. We identified a core set of 3,800 exons from 1,643 genes that show conservation of strongly tissue-dependent usage patterns from human at least to macaque. This set is enriched for exons encoding protein-disordered regions and untranslated regions. Our findings support the theory that isoform regulation is an important target of evolution in primates, and our method provides a powerful tool for discovering potentially functional tissue-dependent isoforms.
Unraveling the molecular processes that lead from genotype to phenotype is crucial for the understanding and effective treatment of genetic diseases. Knowledge of the causative genetic defect most often does not enable treatment; therefore, causal intermediates between genotype and phenotype constitute valuable candidates for molecular intervention points that can be therapeutically targeted. Mapping genetic determinants of gene expression levels (also known as expression quantitative trait loci or eQTL studies) is frequently used for this purpose, yet distinguishing causation from correlation remains a significant challenge. Here, we address this challenge using extensive, multi-environment gene expression and fitness profiling of hundreds of genetically diverse yeast strains, in order to identify truly causal intermediate genes that condition fitness in a given environment. Using functional genomics assays, we show that the predictive power of eQTL studies for inferring causal intermediate genes is poor unless performed across multiple environments. Surprisingly, although the effects of genotype on fitness depended strongly on environment, causal intermediates could be most reliably predicted from genetic effects on expression present in all environments. Our results indicate a mechanism explaining this apparent paradox, whereby immediate molecular consequences of genetic variation are shared across environments, and environment-dependent phenotypic effects result from downstream integration of environmental signals. We developed a statistical model to predict causal intermediates that leverages this insight, yielding over 400 transcripts, for the majority of which we experimentally validated their role in conditioning fitness. Our findings have implications for the design and analysis of clinical omics studies aimed at discovering personalized targets for molecular intervention, suggesting that inferring causation in a single cellular context can benefit from molecular profiling in multiple contexts.
Frequently during evolution, new phenotypes evolved due to novelty in gene regulation, such as that caused by genome rewiring. This has been demonstrated by comparing common regulatory sequences among species and by identifying single regulatory mutations that are associated with new phenotypes. However, while a single mutation changes a single element, gene regulation is accomplished by a regulatory network involving multiple interactive elements. Therefore, to better understand regulatory evolution, we have studied how mutations contributed to the adaptation of cells to a regulatory challenge. We created a synthetic genome rewiring in yeast cells, challenged their gene regulation, and studied their adaptation. HIS3, an essential enzyme for histidine biosynthesis, was placed exclusively under a GAL promoter, which is induced by galactose and strongly repressed in glucose. Such rewired cells were faced with significant regulatory challenges in a repressive glucose medium. We identified several independent mutations in elements of the GAL system associated with the rapid adaptation of cells, such as the repressor GAL80 and the binding sites of the activator GAL4. Consistent with the extraordinarily high rate of cell adaptation, new regulation emerged during adaptation via multiple trajectories, including those involving mutations in elements of the GAL system. The new regulation of HIS3 tuned its expression according to histidine requirements with or without these significant mutations, indicating that additional factors participated in this regulation and that the regulatory network could reorganize in multiple ways to accommodate different mutations. This study, therefore, stresses network plasticity as an important property for regulatory adaptation and evolution.
Active human promoters produce promoter-upstream transcripts (PROMPTs). Why these RNAs are coupled to decay, whereas their neighboring promoter-downstream mRNAs are not, is unknown. Here high-throughput sequencing demonstrates that PROMPTs generally initiate in the antisense direction closely upstream of the transcription start sites (TSSs) of their associated genes. PROMPT TSSs share features with mRNA-producing TSSs, including stalled RNA polymerase II (RNAPII) and the production of small TSS-associated RNAs. Notably, motif analyses around PROMPT 3 ends reveal polyadenylation (pA)-like signals. Mutagenesis studies demonstrate that PROMPT pA signals are functional but linked to RNA degradation. Moreover, pA signals are under-represented in promoter-downstream versus promoter-upstream regions, thus allowing for more efficient RNAPII progress in the sense direction from gene promoters. We conclude that asymmetric sequence distribution around human gene promoters serves to provide a directional RNA output from an otherwise bidirectional transcription process.
HeLa is the most widely used model cell line for studying human cellular and molecular biology. To date, no genomic reference for this cell line has been released, and experiments have relied on the human reference genome. Effective design and interpretation of molecular genetic studies performed using HeLa cells require accurate genomic information. Here we present a detailed genomic and transcriptomic characterization of a HeLa cell line. We performed DNA and RNA sequencing of a HeLa Kyoto cell line and analyzed its mutational portfolio and gene expression profile. Segmentation of the genome according to copy number revealed a remarkably high level of aneuploidy and numerous large structural variants at unprecedented resolution. Some of the extensive genomic rearrangements are indicative of catastrophic chromosome shattering, known as chromothripsis. Our analysis of the HeLa gene expression profile revealed that several pathways, including cell cycle and DNA repair, exhibit significantly different expression patterns from those in normal human tissues. Our results provide the first detailed account of genomic variants in the HeLa genome, yielding insight into their impact on gene expression and cellular function as well as their origins. This study underscores the importance of accounting for the strikingly aberrant characteristics of HeLa cells when designing and interpreting experiments, and has implications for the use of HeLa as a model of human biology.
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.
Living systems may have evolved probabilistic bet hedging strategies that generate cell-to-cell phenotypic diversity in anticipation of environmental catastrophes, as opposed to adaptation via a deterministic response to environmental changes. Evolution of bet hedging assumes that genotypes segregating in natural populations modulate the level of intraclonal diversity, which so far has largely remained hypothetical. Using a fluorescent P(met17)-GFP reporter, we mapped four genetic loci conferring to a wild yeast strain an elevated cell-to-cell variability in the expression of MET17, a gene regulated by the methionine pathway. A frameshift mutation in the Erc1p transmembrane transporter, probably resulting from a release of laboratory strains from negative selection, reduced P(met17)-GFP expression variability. At a second locus, cis-regulatory polymorphisms increased mean expression of the Mup1p methionine permease, causing increased expression variability in trans. These results demonstrate that an expression quantitative trait locus (eQTL) can simultaneously have a deterministic effect in cis and a probabilistic effect in trans. Our observations indicate that the evolution of transmembrane transporter genes can tune intraclonal variation and may therefore be implicated in both reactive and anticipatory strategies of adaptation.
Owing to their preeminent biological functions, the repertoire of expressed RNA-binding proteins (RBPs) and their activity states are highly informative about cellular systems. We have developed a novel and unbiased technique, called interactome capture, for identifying the active RBPs of cultured cells. By making use of in vivo UV cross-linking of RBPs to polyadenylated RNAs, covalently bound proteins are captured with oligo(dT) magnetic beads. After stringent washes, the mRNA interactome is determined by quantitative mass spectrometry (MS). The protocol takes 3 working days for analysis of single proteins by western blotting, and about 2 weeks for the determination of complete cellular mRNA interactomes by MS. The most important advantage of interactome capture over other in vitro and in silico approaches is that only RBPs bound to RNA in a physiological environment are identified. When applied to HeLa cells, interactome capture revealed hundreds of novel RBPs. Interactome capture can also be broadly used to compare different biological states, including metabolic stress, cell cycle, differentiation, development or the response to drugs.
The throughput of next-generation sequencing machines has increased dramatically over the last few years; yet the cost and time for library preparation have not changed proportionally, thus representing the main bottleneck for sequencing large numbers of samples. Here we present an economical, high-throughput library preparation method for the Illumina platform, comprising a 96-well based method for DNA isolation for yeast cells, a low-cost DNA shearing alternative, and adapter ligation using heat inactivation of enzymes instead of bead cleanups.
Anopheles gambiae mosquitoes are the major vectors of human malaria parasites. However, mosquitoes are not passive hosts for parasites, actively limiting their development in vivo. Our current understanding of the mosquito antiparasitic response is mostly based on the phenotypic analysis of gene knockdowns obtained by RNA interference (RNAi), through the injection or transfection of long dsRNAs in adult mosquitoes or cultured cells, respectively. Recently, RNAi has been extended to silence specifically one allele of a given gene in a heterozygous context, thus allowing to compare the contribution of different alleles to a phenotype in the same genetic background.
The use of alternative poly(A) sites is common and affects the post-transcriptional fate of mRNA, including its stability, subcellular localization and translation. Here, we present a method to identify poly(A) sites in a genome-wide and strand-specific manner. This method, termed 3T-fill, initially fills in the poly(A) stretch with unlabeled dTTPs, allowing sequencing to start directly after the poly(A) tail into the 3-untranslated regions (UTR). Our comparative analysis demonstrates that it outperforms existing protocols in quality and throughput and accurately quantifies RNA levels as only one read is produced from each transcript. We use this method to characterize the diversity of polyadenylation in Saccharomyces cerevisiae, showing that alternative RNA molecules are present even in a genetically identical cell population. Finally, we observe that overlap of convergent 3-UTRs is frequent but sharply limited by coding regions, suggesting factors that restrict compression of the yeast genome.
Transcription termination in Saccharomyces cerevisiae can be performed by at least two distinct pathways and is influenced by the phosphorylation status of the carboxy-terminal domain (CTD) of RNA polymerase II (Pol II). Late termination of mRNAs is performed by the CPF/CF complex, the recruitment of which is dependent on CTD-Ser2 phosphorylation (Ser2P). Early termination of shorter cryptic unstable transcripts (CUTs) and small nucleolar/nuclear RNAs (sno/snRNAs) is performed by the Nrd1-Nab3-Sen1 (NNS) complex that binds phosphorylated CTD-Ser5 (Ser5P) via the CTD-interacting domain (CID) of Nrd1p. In this study, mutants of the different termination pathways were compared by genome-wide expression analysis. Surprisingly, the expression changes observed upon loss of the CTD-Ser2 kinase Ctk1p are more similar to those derived from alterations in the Ser5P-dependent NNS pathway, than from loss of CTD-Ser2P binding factors. Tiling array analysis of ctk1? cells reveals readthrough at snoRNAs, at many cryptic unstable transcripts (CUTs) and stable uncharacterized transcripts (SUTs), but only at some mRNAs. Despite the suggested predominant role in termination of mRNAs, we observed that a CTK1 deletion or a Pol II CTD mutant lacking all Ser2 positions does not result in a global mRNA termination defect. Rather, termination defects in these strains are widely observed at NNS-dependent genes. These results indicate that Ctk1p and Ser2 CTD phosphorylation have a wide impact in termination of small non-coding RNAs but only affect a subset of mRNA coding genes.
The MYST HAT Sas2 is part of the SAS-I complex that acetylates histone H4 lysine 16 (H4 K16Ac) and blocks the propagation of heterochromatin at the telomeres of Saccharomyces cerevisiae. In this study, we investigated Sas2-mediated H4 K16Ac on a genome-wide scale. Interestingly, H4 K16Ac loss in sas2? cells outside of the telomeric regions showed a distinctive pattern in that there was a pronounced decrease of H4 K16Ac within the majority of open reading frames (ORFs), but little change in intergenic regions. Furthermore, regions of low histone H3 exchange and low H3 K56 acetylation showed the most pronounced loss of H4 K16Ac in sas2?, indicating that Sas2 deposited this modification on chromatin independently of histone exchange. In agreement with the effect of Sas2 within ORFs, sas2? caused resistance to 6-azauracil, indicating a positive effect on transcription elongation in the absence of H4 K16Ac. In summary, our data suggest that Sas2-dependent H4 K16Ac is deposited into chromatin independently of transcription and histone exchange, and that it has an inhibitory effect on the ability of PolII to travel through the body of the gene.
In the last decade, it became clear that transcription goes far beyond that of protein-coding genes. Most RNA molecules are transcribed from intergenic regions or introns and exhibit much variability in size, expression level, secondary structure, and evolutionary conservation. While for several types of non-coding RNAs some cellular functions have been reported, like for micro-RNAs and small nucleolar RNAs, for most others no indications of function or regulation have so far been found. Therefore, the RNA population inside a cell is diverse and cryptic and, thus, demands powerful methods to study its composition, abundance, and structure. DNA oligonucleotide microarrays have proven to be of great utility to study transcription of genes in various organisms. Recently, due to advancement in microarray technology, tiling microarrays that extend transcription measurement to genomic regions beyond protein-coding genes were designed for several species. The Saccharomyces cerevisiae yeast tiling array contains overlapping probes across the full genomic sequence, with consecutive probes starting every 8 bp on average on each strand, enabling strand-specific measurement of transcription from a full eukaryotic genome. Here, we describe the methods used to extract yeast RNA, convert it into first-strand cDNA, fragment, and label it for hybridization to the tiling array. This protocol will enable researchers not only to study which genes are expressed and to what levels, but also to identify non-coding RNAs and to study the structure of transcripts including their untranslated regions, alternative start, stop, and processing sites. This information will allow understanding their roles inside cells.
Due to the lack of relevant animal models, development of effective treatments for human mitochondrial diseases has been limited. Here we establish a rapid, yeast-based assay to screen for drugs active against human inherited mitochondrial diseases affecting ATP synthase, in particular NARP (neuropathy, ataxia, and retinitis pigmentosa) syndrome. This method is based on the conservation of mitochondrial function from yeast to human, on the unique ability of yeast to survive without production of ATP by oxidative phosphorylation, and on the amenability of the yeast mitochondrial genome to site-directed mutagenesis. Our method identifies chlorhexidine by screening a chemical library and oleate through a candidate approach. We show that these molecules rescue a number of phenotypes resulting from mutations affecting ATP synthase in yeast. These compounds are also active on human cybrid cells derived from NARP patients. These results validate our method as an effective high-throughput screening approach to identify drugs active in the treatment of human ATP synthase disorders and suggest that this type of method could be applied to other mitochondrial diseases.
Ribonuclease P (RNase P) is an essential endoribonuclease that catalyzes the cleavage of the 5 leader of pre-tRNAs. In addition, a growing number of non-tRNA substrates have been identified in various organisms. RNase P varies in composition, as bacterial RNase P contains a catalytic RNA core and one protein subunit, while eukaryotic nuclear RNase P retains the catalytic RNA but has at least nine protein subunits. The additional eukaryotic protein subunits most likely provide additional functionality to RNase P, with one possibility being additional RNA recognition capabilities. To investigate the possible range of additional RNase P substrates in vivo, a strand-specific, high-density microarray was used to analyze what RNA accumulates with a mutation in the catalytic RNA subunit of nuclear RNase P in Saccharomyces cerevisiae. A wide variety of noncoding RNAs were shown to accumulate, suggesting that nuclear RNase P participates in the turnover of normally unstable nuclear RNAs. In some cases, the accumulated noncoding RNAs were shown to be antisense to transcripts that commensurately decreased in abundance. Pre-mRNAs containing introns also accumulated broadly, consistent with either compromised splicing or failure to efficiently turn over pre-mRNAs that do not enter the splicing pathway. Taken together with the high complexity of the nuclear RNase P holoenzyme and its relatively nonspecific capacity to bind and cleave mixed sequence RNAs, these data suggest that nuclear RNase P facilitates turnover of nuclear RNAs in addition to its role in pre-tRNA biogenesis.
Several studies have shown that promoters of protein-coding genes are origins of pervasive non-coding RNA transcription and can initiate transcription in both directions. However, only recently have researchers begun to elucidate the functional implications of this bidirectionality and non-coding RNA production. Increasing evidence indicates that non-coding transcription at promoters influences the expression of protein-coding genes, revealing a new layer of transcriptional regulation. This regulation acts at multiple levels, from modifying local chromatin to enabling regional signal spreading and more distal regulation. Moreover, the bidirectional activity of a promoter is regulated at multiple points during transcription, giving rise to diverse types of transcripts.
Genome-wide transcription profiling has revealed extensive expression of non-coding RNAs antisense to genes, yet their functions, if any, remain to be understood. In this study, we perform a systematic analysis of sense-antisense expression in response to genetic and environmental changes in yeast. We find that antisense expression is associated with genes of larger expression variability. This is characterized by more switching off at low levels of expression for genes with antisense compared to genes without, yet similar expression at maximal induction. By disrupting antisense transcription, we demonstrate that antisense expression confers an on-off switch on gene regulation for the SUR7 gene. Consistent with this, genes that must respond in a switch-like manner, such as stress-response and environment-specific genes, are enriched for antisense expression. In addition, our data provide evidence that antisense expression initiated from bidirectional promoters enables the spreading of regulatory signals from one locus to neighbouring genes. These results indicate a general regulatory effect of antisense expression on sense genes and emphasize the importance of antisense-initiating regions downstream of genes in models of gene regulation.
When mismatches in heteroduplex DNA formed during meiotic recombination are left unrepaired, post-meiotic segregation of the two mismatched alleles occurs during the ensuing round of mitosis. This gives rise to somatic mosaicism in multicellular organisms and leads to unexpected allelic combinations among progeny. Despite its implications for inheritance, post-meiotic segregation has been studied at only a few loci.
Sen1 of S. cerevisiae is a known component of the NRD complex implicated in transcription termination of nonpolyadenylated as well as some polyadenylated RNA polymerase II transcripts. We now show that Sen1 helicase possesses a wider function by restricting the occurrence of RNA:DNA hybrids that may naturally form during transcription, when nascent RNA hybridizes to DNA prior to its packaging into RNA protein complexes. These hybrids displace the nontranscribed strand and create R loop structures. Loss of Sen1 results in transient R loop accumulation and so elicits transcription-associated recombination. SEN1 genetically interacts with DNA repair genes, suggesting that R loop resolution requires proteins involved in homologous recombination. Based on these findings, we propose that R loop formation is a frequent event during transcription and a key function of Sen1 is to prevent their accumulation and associated genome instability.
Budding yeast noncoding RNAs (ncRNAs) are pervasively transcribed during mitosis, and some regulate mitotic protein-coding genes. However, little is known about ncRNA expression during meiotic development. Using high-resolution profiling we identified an extensive meiotic ncRNA expression program interlaced with the protein-coding transcriptome via sense/antisense transcript pairs, bidirectional promoters, and ncRNAs that overlap the regulatory regions of genes. Meiotic unannotated transcripts (MUTs) are mitotic targets of the conserved exosome component Rrp6, which itself is degraded after the onset of meiosis when MUTs and other ncRNAs accumulate in successive waves. Diploid cells lacking Rrp6 fail to initiate premeiotic DNA replication normally and cannot undergo efficient meiotic development. The present study demonstrates a unique function for budding yeast Rrp6 in degrading distinct classes of meiotically induced ncRNAs during vegetative growth and the onset of meiosis and thus points to a critical role of differential ncRNA expression in the execution of a conserved developmental program.
Accurate estimates of mutation rates provide critical information to analyze genome evolution and organism fitness. We used whole-genome DNA sequencing, pulse-field gel electrophoresis, and comparative genome hybridization to determine mutation rates in diploid vegetative and meiotic mutation accumulation lines of Saccharomyces cerevisiae. The vegetative lines underwent only mitotic divisions while the meiotic lines underwent a meiotic cycle every ?20 vegetative divisions. Similar base substitution rates were estimated for both lines. Given our experimental design, these measures indicated that the meiotic mutation rate is within the range of being equal to zero to being 55-fold higher than the vegetative rate. Mutations detected in vegetative lines were all heterozygous while those in meiotic lines were homozygous. A quantitative analysis of intra-tetrad mating events in the meiotic lines showed that inter-spore mating is primarily responsible for rapidly fixing mutations to homozygosity as well as for removing mutations. We did not observe 1-2 nt insertion/deletion (in-del) mutations in any of the sequenced lines and only one structural variant in a non-telomeric location was found. However, a large number of structural variations in subtelomeric sequences were seen in both vegetative and meiotic lines that did not affect viability. Our results indicate that the diploid yeast nuclear genome is remarkably stable during the vegetative and meiotic cell cycles and support the hypothesis that peripheral regions of chromosomes are more dynamic than gene-rich central sections where structural rearrangements could be deleterious. This work also provides an improved estimate for the mutational load carried by diploid organisms.
Fungal pathogens can be lethal, especially among immunocompromised populations, such as patients with AIDS and recipients of tissue transplantation or chemotherapy. Prolonged usage of antifungal reagents can lead to drug resistance and treatment failure. Understanding mechanisms that underlie drug resistance by pathogenic microorganisms is thus vital for dealing with this emerging issue. In this study, we show that dramatic sequence changes in PDR5, an ABC (ATP-binding cassette) efflux transporter protein gene in an opportunistic fungal pathogen, caused the organism to become hypersensitive to azole, a widely used antifungal drug. Surprisingly, the same mutations conferred growth advantages to the organism on polyenes, which are also commonly used antimycotics. Our results indicate that Pdr5p might be important for ergosterol homeostasis. The observed remarkable sequence divergence in the PDR5 gene in yeast strain YJM789 may represent an interesting case of adaptive loss of gene function with significant clinical implications.
Epigenomes commonly refer to the sequence of presence/absence of specific epigenetic marks along eukaryotic chromatin. Complete histone-borne epigenomes have now been described at single-nucleosome resolution from various organisms, tissues, developmental stages, or diseases, yet their intra-species natural variation has never been investigated. We describe here that the epigenomic sequence of histone H3 acetylation at Lysine 14 (H3K14ac) differs greatly between two unrelated strains of the yeast Saccharomyces cerevisiae. Using single-nucleosome chromatin immunoprecipitation and mapping, we interrogated 58,694 nucleosomes and found that 5,442 of them differed in their level of H3K14 acetylation, at a false discovery rate (FDR) of 0.0001. These Single Nucleosome Epi-Polymorphisms (SNEPs) were enriched at regulatory sites and conserved non-coding DNA sequences. Surprisingly, higher acetylation in one strain did not imply higher expression of the relevant gene. However, SNEPs were enriched in genes of high transcriptional variability and one SNEP was associated with the strength of gene activation upon stimulation. Our observations suggest a high level of inter-individual epigenomic variation in natural populations, with essential questions on the origin of this diversity and its relevance to gene x environment interactions.
Extensive transcription of non-coding RNAs has been detected in eukaryotic genomes and is thought to constitute an additional layer in the regulation of gene expression. Despite this role, their transcription through the cell cycle has not been studied; genome-wide approaches have only focused on protein-coding genes. To explore the complex transcriptome architecture underlying the budding yeast cell cycle, we used 8 bp tiling arrays to generate a 5 minute-resolution, strand-specific expression atlas of the whole genome.
Variation in transcriptional regulation is thought to be a major cause of phenotypic diversity. Although widespread differences in gene expression among individuals of a species have been observed, studies to examine the variability of transcription factor binding on a global scale have not been performed, and thus the extent and underlying genetic basis of transcription factor binding diversity is unknown. By mapping differences in transcription factor binding among individuals, here we present the genetic basis of such variation on a genome-wide scale. Whole-genome Ste12-binding profiles were determined using chromatin immunoprecipitation coupled with DNA sequencing in pheromone-treated cells of 43 segregants of a cross between two highly diverged yeast strains and their parental lines. We identified extensive Ste12-binding variation among individuals, and mapped underlying cis- and trans-acting loci responsible for such variation. We showed that most transcription factor binding variation is cis-linked, and that many variations are associated with polymorphisms residing in the binding motifs of Ste12 as well as those of several proposed Ste12 cofactors. We also identified two trans-factors, AMN1 and FLO8, that modulate Ste12 binding to promoters of more than ten genes under alpha-factor treatment. Neither of these two genes was previously known to regulate Ste12, and we suggest that they may be mediators of gene activity and phenotypic diversity. Ste12 binding strongly correlates with gene expression for more than 200 genes, indicating that binding variation is functional. Many of the variable-bound genes are involved in cell wall organization and biogenesis. Overall, these studies identified genetic regulators of molecular diversity among individuals and provide new insights into mechanisms of gene regulation.
The ability of Anopheles gambiae mosquitoes to transmit Plasmodium parasites is highly variable between individuals. However, the genetic basis of this variability has remained unknown. We combined genome-wide mapping and reciprocal allele-specific RNA interference (rasRNAi) to identify the genomic locus that confers resistance to malaria parasites and demonstrated that polymorphisms in a single gene encoding the antiparasitic thioester-containing protein 1 (TEP1) explain a substantial part of the variability in parasite killing. The link between TEP1 alleles and resistance to malaria may offer new tools for controlling malaria transmission. The successful application of rasRNAi in Anopheles suggests that it could also be applied to other organisms where RNAi is feasible to dissect complex phenotypes to the level of individual quantitative trait alleles.
Homology-dependent gene silencing, a phenomenon described as cosuppression in plants, depends on siRNAs. We provide evidence that in Saccharomyces cerevisiae, which is missing the RNAi machinery, protein coding gene cosuppression exists. Indeed, introduction of an additional copy of PHO84 on a plasmid or within the genome results in the cosilencing of both the transgene and the endogenous gene. This repression is transcriptional and position-independent and requires trans-acting antisense RNAs. Antisense RNAs induce transcriptional gene silencing both in cis and in trans, and the two pathways differ by the implication of the Hda1/2/3 complex. We also show that trans-silencing is influenced by the Set1 histone methyltransferase, which promotes antisense RNA production. Finally we show that although antisense-mediated cis-silencing occurs in other genes, trans-silencing so far depends on features specific to PHO84. All together our data highlight the importance of noncoding RNAs in mediating RNAi-independent transcriptional gene silencing.
Microarrays provide an accurate and cost-effective method for genotyping large numbers of individuals at high resolution. The resulting data permit the identification of loci at which genetic variation is associated with quantitative traits, or fine mapping of meiotic recombination, which is a key determinant of genetic diversity among individuals. Several issues inherent to short oligonucleotide arrays -- cross-hybridization, or variability in probe response to target -- have the potential to produce genotyping errors. There is a need for improved statistical methods for array-based genotyping.
Recent reports have shown that most of the genome is transcribed and that transcription frequently occurs concurrently on both DNA strands. In diploid genomes, the expression level of each allele conditions the degree to which sequence polymorphisms affect the phenotype. It is thus essential to quantify expression in an allele- and strand-specific manner. Using a custom-designed tiling array and a new computational approach, we piloted measuring allele- and strand-specific expression in yeast. Confident quantitative estimates of allele-specific expression were obtained for about half of the coding and non-coding transcripts of a heterozygous yeast strain, of which 371 transcripts (13%) showed significant allelic differential expression greater than 1.5-fold. The data revealed complex allelic differential expression on opposite strands. Furthermore, combining allele-specific expression with linkage mapping enabled identifying allelic variants that act in cis and in trans to regulate allelic expression in the heterozygous strain. Our results provide the first high-resolution analysis of differential expression on all four strands of an eukaryotic genome.
Pervasive and hidden transcription is widespread in eukaryotes, but its global level, the mechanisms from which it originates and its functional significance are unclear. Cryptic unstable transcripts (CUTs) were recently described as a principal class of RNA polymerase II transcripts in Saccharomyces cerevisiae. These transcripts are targeted for degradation immediately after synthesis by the action of the Nrd1-exosome-TRAMP complexes. Although CUT degradation mechanisms have been analysed in detail, the genome-wide distribution at the nucleotide resolution and the prevalence of CUTs are unknown. Here we report the first high-resolution genomic map of CUTs in yeast, revealing a class of potentially functional CUTs and the intrinsic bidirectional nature of eukaryotic promoters. An RNA fraction highly enriched in CUTs was analysed by a 3 Long-SAGE (serial analysis of gene expression) approach adapted to deep sequencing. The resulting detailed genomic map of CUTs revealed that they derive from extremely widespread and very well defined transcription units and do not result from unspecific transcriptional noise. Moreover, the transcription of CUTs predominantly arises within nucleosome-free regions, most of which correspond to promoter regions of bona fide genes. Some of the CUTs start upstream from messenger RNAs and overlap their 5 end. Our study of glycolysis genes, as well as recent results from the literature, indicate that such concurrent transcription is potentially associated with regulatory mechanisms. Our data reveal numerous new CUTs with such a potential regulatory role. However, most of the identified CUTs corresponded to transcripts divergent from the promoter regions of genes, indicating that they represent by-products of divergent transcription occurring at many and possibly most promoters. Eukaryotic promoter regions are thus intrinsically bidirectional, a fundamental property that escaped previous analyses because in most cases divergent transcription generates short-lived unstable transcripts present at very low steady-state levels.
Genome-wide pervasive transcription has been reported in many eukaryotic organisms, revealing a highly interleaved transcriptome organization that involves hundreds of previously unknown non-coding RNAs. These recently identified transcripts either exist stably in cells (stable unannotated transcripts, SUTs) or are rapidly degraded by the RNA surveillance pathway (cryptic unstable transcripts, CUTs). One characteristic of pervasive transcription is the extensive overlap of SUTs and CUTs with previously annotated features, which prompts questions regarding how these transcripts are generated, and whether they exert function. Single-gene studies have shown that transcription of SUTs and CUTs can be functional, through mechanisms involving the generated RNAs or their generation itself. So far, a complete transcriptome architecture including SUTs and CUTs has not been described in any organism. Knowledge about the position and genome-wide arrangement of these transcripts will be instrumental in understanding their function. Here we provide a comprehensive analysis of these transcripts in the context of multiple conditions, a mutant of the exosome machinery and different strain backgrounds of Saccharomyces cerevisiae. We show that both SUTs and CUTs display distinct patterns of distribution at specific locations. Most of the newly identified transcripts initiate from nucleosome-free regions (NFRs) associated with the promoters of other transcripts (mostly protein-coding genes), or from NFRs at the 3 ends of protein-coding genes. Likewise, about half of all coding transcripts initiate from NFRs associated with promoters of other transcripts. These data change our view of how a genome is transcribed, indicating that bidirectionality is an inherent feature of promoters. Such an arrangement of divergent and overlapping transcripts may provide a mechanism for local spreading of regulatory signals-that is, coupling the transcriptional regulation of neighbouring genes by means of transcriptional interference or histone modification.
Antisense noncoding transcripts, genes-within-genes, and convergent gene pairs are prevalent among eukaryotes. The existence of such transcription units raises the question of what happens when RNA polymerase II (RNAPII) molecules collide head-to-head. Here we use a combination of biochemical and genetic approaches in yeast to show that polymerases transcribing opposite DNA strands cannot bypass each other. RNAPII stops but does not dissociate upon head-to-head collision in vitro, suggesting that opposing polymerases represent insurmountable obstacles for each other. Head-to-head collision in vivo also results in RNAPII stopping, and removal of collided RNAPII from the DNA template can be achieved via ubiquitylation-directed proteolysis. Indeed, in cells lacking efficient RNAPII polyubiquitylation, the half-life of collided polymerases increases, so that they can be detected between convergent genes. These results provide insight into fundamental mechanisms of gene traffic control and point to an unexplored effect of antisense transcription on gene regulation via polymerase collision.
Natural populations are known to differ not only in DNA but also in their chromatin-associated epigenetic marks. When such inter-individual epigenomic differences (or "epi-polymorphisms") are observed, their stability is usually not known: they may or may not be reprogrammed over time or upon environmental changes. In addition, their origin may be purely epigenetic, or they may result from regulatory variation encoded in the DNA. Studying epi-polymorphisms requires, therefore, an assessment of their nature and stability. Here we estimate the stability of yeast epi-polymorphisms of chromatin acetylation, and we provide a genome-by-epigenome map of their genetic control. A transient epi-drug treatment was able to reprogram acetylation variation at more than one thousand nucleosomes, whereas a similar amount of variation persisted, distinguishing "labile" from "persistent" epi-polymorphisms. Hundreds of genetic loci underlied acetylation variation at 2,418 nucleosomes either locally (in cis) or distantly (in trans), and this genetic control overlapped only partially with the genetic control of gene expression. Trans-acting regulators were not necessarily associated with genes coding for chromatin modifying enzymes. Strikingly, "labile" and "persistent" epi-polymorphisms were associated with poor and strong genetic control, respectively, showing that genetic modifiers contribute to persistence. These results estimate the amount of natural epigenomic variation that can be lost after transient environmental exposures, and they reveal the complex genetic architecture of the DNA-encoded determinism of chromatin epi-polymorphisms. Our observations provide a basis for the development of population epigenetics.
Eukaryotic genomes are extensively transcribed, forming both messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). ncRNAs made by RNA polymerase II often initiate from bidirectional promoters (nucleosome-depleted chromatin) that synthesize mRNA and ncRNA in opposite directions. We demonstrate that, by adopting a gene-loop conformation, actively transcribed mRNA encoding genes restrict divergent transcription of ncRNAs. Because gene-loop formation depends on a protein factor (Ssu72) that coassociates with both the promoter and the terminator, the inactivation of Ssu72 leads to increased synthesis of promoter-associated divergent ncRNAs, referred to as Ssu72-restricted transcripts (SRTs). Similarly, inactivation of individual gene loops by gene mutation enhances SRT synthesis. We demonstrate that gene-loop conformation enforces transcriptional directionality on otherwise bidirectional promoters.
The exosome is a complex involved in the maturation of rRNA and sn-snoRNA, in the degradation of short-lived noncoding RNAs, and in the quality control of RNAs produced in mutants. It contains two catalytic subunits, Rrp6p and Dis3p, whose specific functions are not fully understood. We analyzed the transcriptome of combinations of Rrp6p and Dis3p catalytic mutants by high-resolution tiling arrays. We show that Dis3p and Rrp6p have both overlapping and specific roles in degrading distinct classes of substrates. We found that transcripts derived from more than half of intron-containing genes are degraded before splicing. Surprisingly, we also show that the exosome degrades large amounts of tRNA precursors despite the absence of processing defects. These results underscore the notion that large amounts of RNAs produced in wild-type cells are discarded before entering functional pathways and suggest that kinetic competition with degradation proofreads the efficiency and accuracy of processing.
The Set3 histone deacetylase complex (Set3C) binds histone H3 dimethylated at lysine 4 (H3K4me2) to mediate deacetylation of histones in 5-transcribed regions. To discern how Set3C affects gene expression, genome-wide transcription was analyzed in yeast undergoing a series of carbon source shifts. Deleting SET3 primarily caused changes during transition periods, as genes were induced or repressed. Surprisingly, a majority of Set3-affected genes are overlapped by noncoding RNA (ncRNA) transcription. Many Set3-repressed genes have H3K4me2 instead of me3 over promoter regions, due to either reduced H3K4me3 or ncRNA transcription from distal or antisense promoters. Set3C also represses internal cryptic promoters, but in different regions of genes than the Set2/Rpd3S pathway. Finally, Set3C stimulates some genes by repressing an overlapping antagonistic antisense transcript. These results show that overlapping noncoding transcription can fine-tune gene expression, not via the ncRNA but by depositing H3K4me2 to recruit the Set3C deacetylase.
Alternative polyadenylation site usage gives rise to variation in 3 ends of transcripts in diverse organisms ranging from yeast to human. Accurate mapping of polyadenylation sites of transcripts is of major biological importance, since the length of the 3UTR can have a strong influence on transcript stability, localization, and translation. However, reads generated using total mRNA sequencing mostly lack the very 3 end of transcripts. Here, we present a method that allows simultaneous analysis of alternative 3 ends and transcriptome dynamics at high throughput. By using transcripts produced in vitro, the high precision of end mapping during the protocol can be controlled. This method is illustrated here for budding yeast. However, this method can be applied to any natural or artificially polyadenylated RNA.
Only a few genes remain in the mitochondrial genome retained by every eukaryotic organism that carry out essential functions and are implicated in severe diseases. Experimentally relocating these few genes to the nucleus therefore has both therapeutic and evolutionary implications. Numerous unproductive attempts have been made to do so, with a total of only 5 successes across all organisms. We have taken a novel approach to relocating mitochondrial genes that utilizes naturally nuclear versions from other organisms. We demonstrate this approach on subunit 9/c of ATP synthase, successfully relocating this gene for the first time in any organism by expressing the ATP9 genes from Podospora anserina in Saccharomyces cerevisiae. This study substantiates the role of protein structure in mitochondrial gene transfer: expression of chimeric constructs reveals that the P. anserina proteins can be correctly imported into mitochondria due to reduced hydrophobicity of the first transmembrane segment. Nuclear expression of ATP9, while permitting almost fully functional oxidative phosphorylation, perturbs many cellular properties, including cellular morphology, and activates the heat shock response. Altogether, our study establishes a novel strategy for allotopic expression of mitochondrial genes, demonstrates the complex adaptations required to relocate ATP9, and indicates a reason that this gene was only transferred to the nucleus during the evolution of multicellular organisms.
RNA sequencing is becoming a standard for expression profiling experiments and many tools have been developed in the past few years to analyze RNA-Seq data. Numerous Bioconductor packages are available for next-generation sequencing data loading in R, e.g. ShortRead and Rsamtools as well as to perform differential gene expression analyses, e.g. DESeq and edgeR. However, the processing tasks lying in between these require the precise interplay of many Bioconductor packages, e.g. Biostrings, IRanges or external solutions are to be sought.
Poly(A) (pA) tail binding proteins (PABPs) control mRNA polyadenylation, stability, and translation. In a purified system, S. cerevisiae PABPs, Pab1p and Nab2p, are individually sufficient to provide normal pA tail length. However, it is unknown how this occurs in more complex environments. Here we find that the nuclear exosome subunit Rrp6p counteracts the in vitro and in vivo extension of mature pA tails by the noncanonical pA polymerase Trf4p. Moreover, PABP loading onto nascent pA tails is controlled by Rrp6p; while Pab1p is the major PABP, Nab2p only associates in the absence of Rrp6p. This is because Rrp6p can interact with Nab2p and displace it from pA tails, potentially leading to RNA turnover, as evidenced for certain pre-mRNAs. We suggest that a nuclear mRNP surveillance step involves targeting of Rrp6p by Nab2p-bound pA-tailed RNPs and that pre-mRNA abundance is regulated at this level.
RNA-binding proteins (RBPs) determine RNA fate from synthesis to decay. Employing two complementary protocols for covalent UV crosslinking of RBPs to RNA, we describe a systematic, unbiased, and comprehensive approach, termed "interactome capture," to define the mRNA interactome of proliferating human HeLa cells. We identify 860 proteins that qualify as RBPs by biochemical and statistical criteria, adding more than 300 RBPs to those previously known and shedding light on RBPs in disease, RNA-binding enzymes of intermediary metabolism, RNA-binding kinases, and RNA-binding architectures. Unexpectedly, we find that many proteins of the HeLa mRNA interactome are highly intrinsically disordered and enriched in short repetitive amino acid motifs. Interactome capture is broadly applicable to study mRNA interactome composition and dynamics in varied biological settings.
Related JoVE Video
Journal of Visualized Experiments
What is Visualize?
JoVE Visualize is a tool created to match the last 5 years of PubMed publications to methods in JoVE's video library.
How does it work?
We use abstracts found on PubMed and match them to JoVE videos to create a list of 10 to 30 related methods videos.
Video X seems to be unrelated to Abstract Y...
In developing our video relationships, we compare around 5 million PubMed articles to our library of over 4,500 methods videos. In some cases the language used in the PubMed abstracts makes matching that content to a JoVE video difficult. In other cases, there happens not to be any content in our video library that is relevant to the topic of a given abstract. In these cases, our algorithms are trying their best to display videos with relevant content, which can sometimes result in matched videos with only a slight relation.