B.S. 1994-1997, Molecular Biophysics and Biochemistry, Yale University.
Ph.D. 1997-2001, Genetics, Stanford University. Advisor: Ronald W. Davis.
Postdoctoral research at Stanford Genome Technology Center 2002 and visiting group leader 2003-2013.
Group leader at European Molecular Biology Laboratory (EMBL), 2003-present.
Founding Joint Head (chair) of Genome Biology Unit and Senior Scientist, EMBL 2009-2013.
Co-founder, Sophia Genetics 2011.
ERC Advanced Investigator (twice), 2012-present.
Associate Head of Genome Biology Unit and Senior Scientist, EMBL 2013-2016.
Senior Scientist and Director of Life Science Alliance, EMBL 2016-present.
Professor of Genetics, Stanford University, 2013-present.
Co-Director of the Stanford Genome Technology Center, 2013-present.
Co-Director, Stanford Genome Technology Center (2013 - Present)
Senior Scientist and Director of Life Science Alliance, a Stanford-EMBL collaboration, European Molecular Biology Laboratory (EMBL) (2014 - Present)
Director of Graduate Admissions, Department of Genetics, Stanford University (2013 - Present)
Honors & Awards
Ira Herskowitz Award, Genetics Society of America (2016)
Dr. V. Ramalingaswami Chair, Indian National Science Academy (2014-2015)
EMBO Member, European Molecular Biology Organization (2013)
ERC Advanced Investigator, European Research Council (2012, 2017)
Emmy Noether-Program Young Investigator, Deutsche Forschungsgemeinschaft (2004-2010)
Howard Hughes Predoctoral Fellow, Howard Hughes Medical Institute (1998-2001)
Boards, Advisory Committees, Professional Organizations
Advisory Editorial Board Member, Molecular Systems Biology (2014 - Present)
Associate Editor, Genetics (2014 - Present)
EMBO Member, European Molecular Biology Organization (2013 - Present)
Scientific Advisory Board Member, Saccharomyces Genome Database (SGD, www.yeastgenome.org) (2013 - Present)
Associate Editor, G3: Genes | Genomes | Genetics (2011 - Present)
President, Scientific Advisory Board, Sophia Genetics (2011 - Present)
Ph.D., Stanford University, Genetics (2001)
Bachelor of Science, Yale University, Molecular Biophysics and Biochemistry (1997)
Current Research and Scholarly Interests
One of the most daunting challenges in medicine is the complex nature of most diseases (including cancer, diabetes, heart disease and some forms of rare disease) due to interactions between multiple genetic variants and environmental influences. Our research is directed at understanding such complex traits; to do so, we develop novel genomic approaches to investigate the molecular processes that link genotype to phenotype, identify the causal underlying factors, and quantify their contributions. We investigate variation at the level of the genome, transcriptome and proteome, which we integrate with higher-level phenotypes. We also use the resulting molecular networks to predict and evaluate intervention points that enable modulation of phenotype. In particular, our projects are in the following areas:
1) Functions and mechanisms of gene regulation
We have developed several technologies to characterize and quantify pervasive transcription at the genome-wide level as well as its functional impact. In particular, we are interested in the function and regulation of non-coding RNAs, antisense transcription, and the molecular phenotypes that arise from transcriptome complexity.
2) Systems genetics
We have also piloted new technologies to dissect the genetic basis of complex, multifactorial phenotypes. We are interested in studying how genetic variation is inherited through recombination, the consequences of genetic variation, learning to predict phenotype from genotype, and integrating multiple layers of molecular data in order to define intervention points that can be targeted to modulate phenotypes of interest.
3) Disease models
We have used multiple model systems, primarily yeast and human cells, to characterize the genetic and cellular systems affected in particular diseases and assess potential therapeutic strategies. We are studying diseases in patient-derived cells using systematic, followed by mechanistic approaches, to unravel mechanisms and discover novel treatments.
Notably, we place a strong emphasis on the development of new technologies, leveraging the Stanford Genome Technology Center's strengths in this area. Some engineering applications we aim to develop include novel biosensors for detection of minute quantities of biological material and single-cell approaches to investigate genome and transcriptome complexity. Through developing such technologies, we will work with the Center to reduce the cost and increase the accuracy and throughput of biomedical research and health care. Ultimately, we strive to develop approaches that will enable personalized, preventive medicine.
Graduate and Fellowship Programs
- Single-cell transcriptome analysis reveals coordinated ectopic gene-expression patterns in medullary thymic epithelial cells NATURE IMMUNOLOGY 2015; 16 (9): 933-941
Widespread Co-translational RNA Decay Reveals Ribosome Dynamics.
2015; 161 (6): 1400-1412
It is generally assumed that mRNAs undergoing translation are protected from decay. Here, we show that mRNAs are, in fact, co-translationally degraded. This is a widespread and conserved process affecting most genes, where 5'-3' transcript degradation follows the last translating ribosome, producing an in vivo ribosomal footprint. By sequencing the ends of 5' phosphorylated mRNA degradation intermediates, we obtain a genome-wide drug-free measurement of ribosome dynamics. We identify general translation termination pauses in both normal and stress conditions. In addition, we describe novel codon-specific ribosomal pausing sites in response to oxidative stress that are dependent on the RNase Rny1. Our approach is simple and straightforward and does not require the use of translational inhibitors or in vitro RNA footprinting that can alter ribosome protection patterns.
View details for DOI 10.1016/j.cell.2015.05.008
View details for PubMedID 26046441
Mitochondrial protein sorting as a therapeutic target for ATP synthase disorders
Mitochondrial diseases are systemic, prevalent and often fatal; yet treatments remain scarce. Identifying molecular intervention points that can be therapeutically targeted remains a major challenge, which we confronted via a screening assay we developed. Using yeast models of mitochondrial ATP synthase disorders, we screened a drug repurposing library, and applied genomic and biochemical techniques to identify pathways of interest. Here we demonstrate that modulating the sorting of nuclear-encoded proteins into mitochondria, mediated by the TIM23 complex, proves therapeutic in both yeast and patient-derived cells exhibiting ATP synthase deficiency. Targeting TIM23-dependent protein sorting improves an array of phenotypes associated with ATP synthase disorders, including biogenesis and activity of the oxidative phosphorylation machinery. Our study establishes mitochondrial protein sorting as an intervention point for ATP synthase disorders, and because of the central role of this pathway in mitochondrial biogenesis, it holds broad value for the treatment of mitochondrial diseases.
View details for DOI 10.1038/ncomms6585
View details for Web of Science ID 000347224100001
View details for PubMedID 25519239
Genotype-environment interactions reveal causal pathways that mediate genetic effects on phenotype.
2013; 9 (9)
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.
View details for DOI 10.1371/journal.pgen.1003803
View details for PubMedID 24068968
View details for PubMedCentralID PMC3778020
The Genomic and Transcriptomic Landscape of a HeLa Cell Line
G3-GENES GENOMES GENETICS
2013; 3 (8): 1213-1224
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.
View details for DOI 10.1534/g3.113.005777
View details for Web of Science ID 000322822300005
View details for PubMedID 23550136
Extensive transcriptional heterogeneity revealed by isoform profiling.
2013; 497 (7447): 127-131
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.
View details for DOI 10.1038/nature12121
View details for PubMedID 23615609
Gene Loops Enhance Transcriptional Directionality
2012; 338 (6107): 671-675
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.
View details for DOI 10.1126/science.1224350
View details for Web of Science ID 000310516000054
View details for PubMedID 23019609
Dissecting the Genetic Basis of Resistance to Malaria Parasites in Anopheles gambiae
2009; 326 (5949): 147-150
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.
View details for DOI 10.1126/science.1175241
View details for Web of Science ID 000270355600057
View details for PubMedID 19797663
Bidirectional promoters generate pervasive transcription in yeast
2009; 457 (7232): 1033-U7
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.
View details for DOI 10.1038/nature07728
View details for Web of Science ID 000263425400047
View details for PubMedID 19169243
High-resolution mapping of meiotic crossovers and non-crossovers in yeast
2008; 454 (7203): 479-U1
Meiotic recombination has a central role in the evolution of sexually reproducing organisms. The two recombination outcomes, crossover and non-crossover, increase genetic diversity, but have the potential to homogenize alleles by gene conversion. Whereas crossover rates vary considerably across the genome, non-crossovers and gene conversions have only been identified in a handful of loci. To examine recombination genome wide and at high spatial resolution, we generated maps of crossovers, crossover-associated gene conversion and non-crossover gene conversion using dense genetic marker data collected from all four products of fifty-six yeast (Saccharomyces cerevisiae) meioses. Our maps reveal differences in the distributions of crossovers and non-crossovers, showing more regions where either crossovers or non-crossovers are favoured than expected by chance. Furthermore, we detect evidence for interference between crossovers and non-crossovers, a phenomenon previously only known to occur between crossovers. Up to 1% of the genome of each meiotic product is subject to gene conversion in a single meiosis, with detectable bias towards GC nucleotides. To our knowledge the maps represent the first high-resolution, genome-wide characterization of the multiple outcomes of recombination in any organism. In addition, because non-crossover hotspots create holes of reduced linkage within haplotype blocks, our results stress the need to incorporate non-crossovers into genetic linkage analysis.
View details for DOI 10.1038/nature07135
View details for Web of Science ID 000257860300043
View details for PubMedID 18615017
Genome sequencing and comparative analysis of Saccharomyces cerevisiae strain YJM789
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2007; 104 (31): 12825-12830
We sequenced the genome of Saccharomyces cerevisiae strain YJM789, which was derived from a yeast isolated from the lung of an AIDS patient with pneumonia. The strain is used for studies of fungal infections and quantitative genetics because of its extensive phenotypic differences to the laboratory reference strain, including growth at high temperature and deadly virulence in mouse models. Here we show that the approximately 12-Mb genome of YJM789 contains approximately 60,000 SNPs and approximately 6,000 indels with respect to the reference S288c genome, leading to protein polymorphisms with a few known cases of phenotypic changes. Several ORFs are found to be unique to YJM789, some of which might have been acquired through horizontal transfer. Localized regions of high polymorphism density are scattered over the genome, in some cases spanning multiple ORFs and in others concentrated within single genes. The sequence of YJM789 contains clues to pathogenicity and spurs the development of more powerful approaches to dissecting the genetic basis of complex hereditary traits.
View details for DOI 10.1073/pnas.0701291104
View details for Web of Science ID 000248603900043
View details for PubMedID 17652520
View details for PubMedCentralID PMC1933262
A high-resolution map of transcription in the yeast genome
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2006; 103 (14): 5320-5325
There is abundant transcription from eukaryotic genomes unaccounted for by protein coding genes. A high-resolution genome-wide survey of transcription in a well annotated genome will help relate transcriptional complexity to function. By quantifying RNA expression on both strands of the complete genome of Saccharomyces cerevisiae using a high-density oligonucleotide tiling array, this study identifies the boundary, structure, and level of coding and noncoding transcripts. A total of 85% of the genome is expressed in rich media. Apart from expected transcripts, we found operon-like transcripts, transcripts from neighboring genes not separated by intergenic regions, and genes with complex transcriptional architecture where different parts of the same gene are expressed at different levels. We mapped the positions of 3' and 5' UTRs of coding genes and identified hundreds of RNA transcripts distinct from annotated genes. These nonannotated transcripts, on average, have lower sequence conservation and lower rates of deletion phenotype than protein coding genes. Many other transcripts overlap known genes in antisense orientation, and for these pairs global correlations were discovered: UTR lengths correlated with gene function, localization, and requirements for regulation; antisense transcripts overlapped 3' UTRs more than 5' UTRs; UTRs with overlapping antisense tended to be longer; and the presence of antisense associated with gene function. These findings may suggest a regulatory role of antisense transcription in S. cerevisiae. Moreover, the data show that even this well studied genome has transcriptional complexity far beyond current annotation.
View details for DOI 10.1073/pnas.0601091103
View details for Web of Science ID 000236636400017
View details for PubMedID 16569694
View details for PubMedCentralID PMC1414796
Assembly of functionally integrated human forebrain spheroids
2017; 545 (7652): 54-?
The development of the nervous system involves a coordinated succession of events including the migration of GABAergic (γ-aminobutyric-acid-releasing) neurons from ventral to dorsal forebrain and their integration into cortical circuits. However, these interregional interactions have not yet been modelled with human cells. Here we generate three-dimensional spheroids from human pluripotent stem cells that resemble either the dorsal or ventral forebrain and contain cortical glutamatergic or GABAergic neurons. These subdomain-specific forebrain spheroids can be assembled in vitro to recapitulate the saltatory migration of interneurons observed in the fetal forebrain. Using this system, we find that in Timothy syndrome-a neurodevelopmental disorder that is caused by mutations in the CaV1.2 calcium channel-interneurons display abnormal migratory saltations. We also show that after migration, interneurons functionally integrate with glutamatergic neurons to form a microphysiological system. We anticipate that this approach will be useful for studying neural development and disease, and for deriving spheroids that resemble other brain regions to assemble circuits in vitro.
View details for DOI 10.1038/nature22330
View details for Web of Science ID 000400480400029
View details for PubMedID 28445465
- Modulating Crossover Frequency and Interference for Obligate Crossovers in Saccharomyces cerevisiae Meiosis G3-GENES GENOMES GENETICS 2017; 7 (5): 1511-1524
Myt1l safeguards neuronal identity by actively repressing many non-neuronal fates
2017; 544 (7649): 245-?
Normal differentiation and induced reprogramming require the activation of target cell programs and silencing of donor cell programs. In reprogramming, the same factors are often used to reprogram many different donor cell types. As most developmental repressors, such as RE1-silencing transcription factor (REST) and Groucho (also known as TLE), are considered lineage-specific repressors, it remains unclear how identical combinations of transcription factors can silence so many different donor programs. Distinct lineage repressors would have to be induced in different donor cell types. Here, by studying the reprogramming of mouse fibroblasts to neurons, we found that the pan neuron-specific transcription factor Myt1-like (Myt1l) exerts its pro-neuronal function by direct repression of many different somatic lineage programs except the neuronal program. The repressive function of Myt1l is mediated via recruitment of a complex containing Sin3b by binding to a previously uncharacterized N-terminal domain. In agreement with its repressive function, the genomic binding sites of Myt1l are similar in neurons and fibroblasts and are preferentially in an open chromatin configuration. The Notch signalling pathway is repressed by Myt1l through silencing of several members, including Hes1. Acute knockdown of Myt1l in the developing mouse brain mimicked a Notch gain-of-function phenotype, suggesting that Myt1l allows newborn neurons to escape Notch activation during normal development. Depletion of Myt1l in primary postmitotic neurons de-repressed non-neuronal programs and impaired neuronal gene expression and function, indicating that many somatic lineage programs are actively and persistently repressed by Myt1l to maintain neuronal identity. It is now tempting to speculate that similar 'many-but-one' lineage repressors exist for other cell fates; such repressors, in combination with lineage-specific activators, would be prime candidates for use in reprogramming additional cell types.
View details for DOI 10.1038/nature21722
View details for Web of Science ID 000398897900040
View details for PubMedID 28379941
Human haematopoietic stem cell lineage commitment is a continuous process
NATURE CELL BIOLOGY
2017; 19 (4): 271-?
Blood formation is believed to occur through stepwise progression of haematopoietic stem cells (HSCs) following a tree-like hierarchy of oligo-, bi- and unipotent progenitors. However, this model is based on the analysis of predefined flow-sorted cell populations. Here we integrated flow cytometric, transcriptomic and functional data at single-cell resolution to quantitatively map early differentiation of human HSCs towards lineage commitment. During homeostasis, individual HSCs gradually acquire lineage biases along multiple directions without passing through discrete hierarchically organized progenitor populations. Instead, unilineage-restricted cells emerge directly from a 'continuum of low-primed undifferentiated haematopoietic stem and progenitor cells' (CLOUD-HSPCs). Distinct gene expression modules operate in a combinatorial manner to control stemness, early lineage priming and the subsequent progression into all major branches of haematopoiesis. These data reveal a continuous landscape of human steady-state haematopoiesis downstream of HSCs and provide a basis for the understanding of haematopoietic malignancies.
View details for DOI 10.1038/ncb3493
View details for Web of Science ID 000397917000006
View details for PubMedID 28319093
Modulation of mRNA and lncRNA expression dynamics by the Set2-Rpd3S pathway
H3K36 methylation by Set2 targets Rpd3S histone deacetylase to transcribed regions of mRNA genes, repressing internal cryptic promoters and slowing elongation. Here we explore the function of this pathway by analysing transcription in yeast undergoing a series of carbon source shifts. Approximately 80 mRNA genes show increased induction upon SET2 deletion. A majority of these promoters have overlapping lncRNA transcription that targets H3K36me3 and deacetylation by Rpd3S to the mRNA promoter. We previously reported a similar mechanism for H3K4me2-mediated repression via recruitment of the Set3C histone deacetylase. Here we show that the distance between an mRNA and overlapping lncRNA promoter determines whether Set2-Rpd3S or Set3C represses. This analysis also reveals many previously unreported cryptic ncRNAs induced by specific carbon sources, showing that cryptic promoters can be environmentally regulated. Therefore, in addition to repression of cryptic transcription and modulation of elongation, H3K36 methylation maintains optimal expression dynamics of many mRNAs and ncRNAs.
View details for DOI 10.1038/ncomms13534
View details for Web of Science ID 000388574100001
View details for PubMedID 27892458
View details for PubMedCentralID PMC5133700
A privacy-preserving solution for compressed storage and selective retrieval of genomic data.
In clinical genomics, the continuous evolution of bioinformatic algorithms and sequencing platforms makes it beneficial to store patients' complete aligned genomic data in addition to variant calls relative to a reference sequence. Due to the large size of human genome sequence data files (varying from 30 GB to 200 GB depending on coverage), two major challenges facing genomics laboratories are the costs of storage and the efficiency of the initial data processing. In addition, privacy of genomic data is becoming an increasingly serious concern, yet no standard data storage solutions exist that enable compression, encryption, and selective retrieval. Here we present a privacy-preserving solution named SECRAM (Selective retrieval on Encrypted and Compressed Reference-oriented Alignment Map) for the secure storage of compressed aligned genomic data. Our solution enables selective retrieval of encrypted data and improves the efficiency of downstream analysis (e.g., variant calling). Compared with BAM, the de facto standard for storing aligned genomic data, SECRAM uses 18% less storage. Compared with CRAM, one of the most compressed nonencrypted formats (using 34% less storage than BAM), SECRAM maintains efficient compression and downstream data processing, while allowing for unprecedented levels of security in genomic data storage. Compared with previous work, the distinguishing features of SECRAM are that (1) it is position-based instead of read-based, and (2) it allows random querying of a subregion from a BAM-like file in an encrypted form. Our method thus offers a space-saving, privacy-preserving, and effective solution for the storage of clinical genomic data.
View details for PubMedID 27789525
A global genetic interaction network maps a wiring diagram of cellular function
2016; 353 (6306): 1381-?
We generated a global genetic interaction network for Saccharomyces cerevisiae, constructing more than 23 million double mutants, identifying about 550,000 negative and about 350,000 positive genetic interactions. This comprehensive network maps genetic interactions for essential gene pairs, highlighting essential genes as densely connected hubs. Genetic interaction profiles enabled assembly of a hierarchical model of cell function, including modules corresponding to protein complexes and pathways, biological processes, and cellular compartments. Negative interactions connected functionally related genes, mapped core bioprocesses, and identified pleiotropic genes, whereas positive interactions often mapped general regulatory connections among gene pairs, rather than shared functionality. The global network illustrates how coherent sets of genetic interactions connect protein complex and pathway modules to map a functional wiring diagram of the cell.
View details for DOI 10.1126/science.aaf1420
View details for Web of Science ID 000383708700030
View details for PubMedID 27708008
Principles for RNA metabolism and alternative transcription initiation within closely spaced promoters
2016; 48 (9): 984-?
Mammalian transcriptomes are complex and formed by extensive promoter activity. In addition, gene promoters are largely divergent and initiate transcription of reverse-oriented promoter upstream transcripts (PROMPTs). Although PROMPTs are commonly terminated early, influenced by polyadenylation sites, promoters often cluster so that the divergent activity of one might impact another. Here we found that the distance between promoters strongly correlates with the expression, stability and length of their associated PROMPTs. Adjacent promoters driving divergent mRNA transcription support PROMPT formation, but owing to polyadenylation site constraints, these transcripts tend to spread into the neighboring mRNA on the same strand. This mechanism to derive new alternative mRNA transcription start sites (TSSs) is also evident at closely spaced promoters supporting convergent mRNA transcription. We suggest that basic building blocks of divergently transcribed core promoter pairs, in combination with the wealth of TSSs in mammalian genomes, provide a framework with which evolution shapes transcriptomes.
View details for DOI 10.1038/ng.3616
View details for Web of Science ID 000382398800008
View details for PubMedID 27455346
Meiotic Interactors of a Mitotic Gene TAO3 Revealed by Functional Analysis of its Rare Variant
G3-GENES GENOMES GENETICS
2016; 6 (8): 2255-2263
Studying the molecular consequences of rare genetic variants has the potential to identify novel and hitherto uncharacterized pathways causally contributing to phenotypic variation. Here, we characterize the functional consequences of a rare coding variant of TAO3, previously reported to contribute significantly to sporulation efficiency variation in Saccharomyces cerevisiae During mitosis, the common TAO3 allele interacts with CBK1-a conserved NDR kinase. Both TAO3 and CBK1 are components of the RAM signaling network that regulates cell separation and polarization during mitosis. We demonstrate that the role of the rare allele TAO3(4477C) in meiosis is distinct from its role in mitosis by being independent of ACE2-a RAM network target gene. By quantitatively measuring cell morphological dynamics, and expressing the TAO3(4477C) allele conditionally during sporulation, we show that TAO3 has an early role in meiosis. This early role of TAO3 coincides with entry of cells into meiotic division. Time-resolved transcriptome analyses during early sporulation identified regulators of carbon and lipid metabolic pathways as candidate mediators. We show experimentally that, during sporulation, the TAO3(4477C) allele interacts genetically with ERT1 and PIP2, regulators of the tricarboxylic acid cycle and gluconeogenesis metabolic pathways, respectively. We thus uncover a meiotic functional role for TAO3, and identify ERT1 and PIP2 as novel regulators of sporulation efficiency. Our results demonstrate that studying the causal effects of genetic variation on the underlying molecular network has the potential to provide a more extensive understanding of the pathways driving a complex trait.
View details for DOI 10.1534/g3.116.029900/-/DC1
View details for Web of Science ID 000381282300002
View details for PubMedID 27317780
Protein Abundance Control by Non-coding Antisense Transcription
2016; 15 (12): 2625-2636
Stable unannotated transcripts (SUTs), some of which overlap protein-coding genes in antisense direction, are a class of non-coding RNAs. While case studies have reported important regulatory roles for several of such RNAs, their general impact on protein abundance regulation of the overlapping gene is not known. To test this, we employed seamless gene manipulation to repress antisense SUTs of 162 yeast genes by using a unidirectional transcriptional terminator and a GFP tag. We found that the mere presence of antisense SUTs was not sufficient to influence protein abundance, that observed effects of antisense SUTs correlated with sense transcript start site overlap, and that the effects were generally weak and led to reduced protein levels. Antisense regulated genes showed increased H3K4 di- and trimethylation and had slightly lower than expected noise levels. Our results suggest that the functionality of antisense RNAs has gene and condition-specific components.
View details for DOI 10.1016/j.celrep.2016.05.043
View details for Web of Science ID 000378255900007
View details for PubMedID 27292640
Functional interplay between MSL1 and CDK7 controls RNA polymerase II Ser5 phosphorylation
NATURE STRUCTURAL & MOLECULAR BIOLOGY
2016; 23 (6): 580-589
Proper gene expression requires coordinated interplay among transcriptional coactivators, transcription factors and the general transcription machinery. We report here that MSL1, a central component of the dosage compensation complex in Drosophila melanogaster and Drosophila virilis, displays evolutionarily conserved sex-independent binding to promoters. Genetic and biochemical analyses reveal a functional interaction of MSL1 with CDK7, a subunit of the Cdk-activating kinase (CAK) complex of the general transcription factor TFIIH. Importantly, MSL1 depletion leads to decreased phosphorylation of Ser5 of RNA polymerase II. In addition, we demonstrate that MSL1 is a phosphoprotein, and transgenic flies expressing MSL1 phosphomutants show mislocalization of the histone acetyltransferase MOF and histone H4 K16 acetylation, thus ultimately causing male lethality due to a failure of dosage compensation. We propose that, by virtue of its interaction with components of the general transcription machinery, MSL1 exists in different phosphorylation states, thereby modulating transcription in flies.
View details for DOI 10.1038/nsmb.3233
View details for Web of Science ID 000377421800021
View details for PubMedID 27183194
Translational Capacity of a Cell Is Determined during Transcription Elongation via the Ccr4-Not Complex
2016; 15 (8): 1782-1794
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.
View details for DOI 10.1016/j.celrep.2016.04.055
View details for Web of Science ID 000376654500014
View details for PubMedID 27184853
The cellular growth rate controls overall mRNA turnover, and modulates either transcription or degradation rates of particular gene regulons
NUCLEIC ACIDS RESEARCH
2016; 44 (8): 3643-3658
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.
View details for DOI 10.1093/nar/gkv1512
View details for Web of Science ID 000376389000023
View details for PubMedID 26717982
DChIPRep, an R/Bioconductor package for differential enrichment analysis in chromatin studies
The genome-wide study of epigenetic states requires the integrative analysis of histone modification ChIP-seq data. Here, we introduce an easy-to-use analytic framework to compare profiles of enrichment in histone modifications around classes of genomic elements, e.g. transcription start sites (TSS). Our framework is available via the user-friendly R/Bioconductor package DChIPRep. DChIPRep uses biological replicate information as well as chromatin Input data to allow for a rigorous assessment of differential enrichment. DChIPRep is available for download through the Bioconductor project at http://bioconductor.org/packages/DChIPRep. Contact. DChIPRep@gmail.com.
View details for DOI 10.7717/peerj.1981
View details for Web of Science ID 000375186000008
View details for PubMedID 27168989
Landscape and Dynamics of Transcription Initiation in the Malaria Parasite Plasmodium falciparum
2016; 14 (10): 2463-2475
A comprehensive map of transcription start sites (TSSs) across the highly AT-rich genome of P. falciparum would aid progress toward deciphering the molecular mechanisms that underlie the timely regulation of gene expression in this malaria parasite. Using high-throughput sequencing technologies, we generated a comprehensive atlas of transcription initiation events at single-nucleotide resolution during the parasite intra-erythrocytic developmental cycle. This detailed analysis of TSS usage enabled us to define architectural features of plasmodial promoters. We demonstrate that TSS selection and strength are constrained by local nucleotide composition. Furthermore, we provide evidence for coordinate and stage-specific TSS usage from distinct sites within the same transcription unit, thereby producing transcript isoforms, a subset of which are developmentally regulated. This work offers a framework for further investigations into the interactions between genomic sequences and regulatory factors governing the complex transcriptional program of this major human pathogen.
View details for DOI 10.1016/j.celrep.2016.02.025
View details for Web of Science ID 000371998700017
View details for PubMedID 26947071
Integrating Cell Phone Imaging with Magnetic Levitation (i-LEV) for Label-Free Blood Analysis at the Point-of-Living.
2016; 12 (9): 1222-1229
There is an emerging need for portable, robust, inexpensive, and easy-to-use disease diagnosis and prognosis monitoring platforms to share health information at the point-of-living, including clinical and home settings. Recent advances in digital health technologies have improved early diagnosis, drug treatment, and personalized medicine. Smartphones with high-resolution cameras and high data processing power enable intriguing biomedical applications when integrated with diagnostic devices. Further, these devices have immense potential to contribute to public health in resource-limited settings where there is a particular need for portable, rapid, label-free, easy-to-use, and affordable biomedical devices to diagnose and continuously monitor patients for precision medicine, especially those suffering from rare diseases, such as sickle cell anemia, thalassemia, and chronic fatigue syndrome. Here, a magnetic levitation-based diagnosis system is presented in which different cell types (i.e., white and red blood cells) are levitated in a magnetic gradient and separated due to their unique densities. Moreover, an easy-to-use, smartphone incorporated levitation system for cell analysis is introduced. Using our portable imaging magnetic levitation (i-LEV) system, it is shown that white and red blood cells can be identified and cell numbers can be quantified without using any labels. In addition, cells levitated in i-LEV can be distinguished at single-cell resolution, potentially enabling diagnosis and monitoring, as well as clinical and research applications.
View details for DOI 10.1002/smll.201501845
View details for PubMedID 26523938
Genome-wide quantification of 5'-phosphorylated mRNA degradation intermediates for analysis of ribosome dynamics.
2016; 11 (2): 359-376
Co-translational mRNA degradation is a widespread process in which 5'-3' exonucleolytic degradation follows the last translating ribosome, thus producing an in vivo ribosomal footprint that delimits the 5' position of the mRNA molecule within the ribosome. To study this degradation process and ribosome dynamics, we developed 5PSeq, which is a method that profiles the genome-wide abundance of mRNA degradation intermediates by virtue of their 5'-phosphorylated (5'P) ends. The approach involves targeted ligation of an oligonucleotide to the 5'P end of mRNA degradation intermediates, followed by depletion of rRNA molecules, reverse transcription of 5'P mRNAs and Illumina high-throughput sequencing. 5PSeq can identify translational pauses at rare codons that are often masked when using alternative methods. This approach can be applied to previously extracted RNA samples, and it is straightforward and does not require polyribosome purification or in vitro RNA footprinting. The protocol we describe here can be applied to Saccharomyces cerevisiae and potentially to other eukaryotic organisms. Three days are required to generate 5PSeq libraries.
View details for DOI 10.1038/nprot.2016.026
View details for PubMedID 26820793
Comprehensive Identification of RNA-Binding Proteins by RNA Interactome Capture.
Methods in molecular biology (Clifton, N.J.)
2016; 1358: 131-139
RNA associates with RNA-binding proteins (RBPs) from synthesis to decay, forming dynamic ribonucleoproteins (RNPs). In spite of the preeminent role of RBPs regulating RNA fate, the scope of cellular RBPs has remained largely unknown. We have recently developed a novel and comprehensive method to identify the repertoire of active RBPs of cultured cells, called RNA interactome capture. Using in vivo UV cross-linking on cultured cells, proteins are covalently bound to RNA if the contact between the two is direct ("zero distance"). Protein-RNA complexes are purified by poly(A) tail-dependent oligo(dT) capture and analyzed by quantitative mass spectrometry. Because UV irradiation is applied to living cells and purification is performed using highly stringent washes, RNA interactome capture identifies physiologic and direct protein-RNA interactions. Applied to HeLa cells, this protocol revealed the near-complete repertoire of RBPs, including hundreds of novel RNA binders. Apart from its RBP discovery capacity, quantitative and comparative RNA interactome capture can also be used to study the responses of the RBP repertoire to different physiological cues and processes, including metabolic stress, differentiation, development, or the response to drugs.
View details for DOI 10.1007/978-1-4939-3067-8_8
View details for PubMedID 26463381
Quantitative CRISPR interference screens in yeast identify chemical-genetic interactions and new rules for guide RNA design.
2016; 17 (1): 45-?
Genome-scale CRISPR interference (CRISPRi) has been used in human cell lines; however, the features of effective guide RNAs (gRNAs) in different organisms have not been well characterized. Here, we define rules that determine gRNA effectiveness for transcriptional repression in Saccharomyces cerevisiae.We create an inducible single plasmid CRISPRi system for gene repression in yeast, and use it to analyze fitness effects of gRNAs under 18 small molecule treatments. Our approach correctly identifies previously described chemical-genetic interactions, as well as a new mechanism of suppressing fluconazole toxicity by repression of the ERG25 gene. Assessment of multiple target loci across treatments using gRNA libraries allows us to determine generalizable features associated with gRNA efficacy. Guides that target regions with low nucleosome occupancy and high chromatin accessibility are clearly more effective. We also find that the best region to target gRNAs is between the transcription start site (TSS) and 200 bp upstream of the TSS. Finally, unlike nuclease-proficient Cas9 in human cells, the specificity of truncated gRNAs (18 nt of complementarity to the target) is not clearly superior to full-length gRNAs (20 nt of complementarity), as truncated gRNAs are generally less potent against both mismatched and perfectly matched targets.Our results establish a powerful functional and chemical genomics screening method and provide guidelines for designing effective gRNAs, which consider chromatin state and position relative to the target gene TSS. These findings will enable effective library design and genome-wide programmable gene repression in many genetic backgrounds.
View details for DOI 10.1186/s13059-016-0900-9
View details for PubMedID 26956608
Genome-Wide Identification of Alternative Polyadenylation Events Using 3'T-Fill.
Methods in molecular biology (Clifton, N.J.)
2016; 1358: 295-302
Due to the increasing appreciation of the impact of alternative polyadenylation on cellular biology, our straightforward, scalable method is of interest to any researcher studying eukaryotic transcription. In addition to high quality gene expression measurements, it precisely maps poly(A) sites and thereby permits the distinction between differential 3'UTR isoforms. As sequencing through long homopolymer stretches is not possible on the Illumina platform, we developed a method that fills up the poly(A) stretch with dTTPs before the sequencing reaction starts.
View details for DOI 10.1007/978-1-4939-3067-8_18
View details for PubMedID 26463391
- Loss of the Yeast SR Protein Npl3 Alters Gene Expression Due to Transcription Readthrough PLOS GENETICS 2015; 11 (12)
- Chromatin Dynamics and the RNA Exosome Function in Concert to Regulate Transcriptional Homeostasis CELL REPORTS 2015; 13 (8): 1610-1622
- Inflammation-Induced Emergency Megakaryopoiesis Driven by Hematopoietic Stem Cell-like Megakaryocyte Progenitors CELL STEM CELL 2015; 17 (4): 422-434
- Genetic Control of Chromatin States in Humans Involves Local and Distal Chromosomal Interactions CELL 2015; 162 (5): 1051-1065
Multitarget, quantitative nanoplasmonic electrical field-enhanced resonating device (NE2RD) for diagnostics.
Proceedings of the National Academy of Sciences of the United States of America
2015; 112 (32): E4354-63
Recent advances in biosensing technologies present great potential for medical diagnostics, thus improving clinical decisions. However, creating a label-free general sensing platform capable of detecting multiple biotargets in various clinical specimens over a wide dynamic range, without lengthy sample-processing steps, remains a considerable challenge. In practice, these barriers prevent broad applications in clinics and at patients' homes. Here, we demonstrate the nanoplasmonic electrical field-enhanced resonating device (NE(2)RD), which addresses all these impediments on a single platform. The NE(2)RD employs an immunodetection assay to capture biotargets, and precisely measures spectral color changes by their wavelength and extinction intensity shifts in nanoparticles without prior sample labeling or preprocessing. We present through multiple examples, a label-free, quantitative, portable, multitarget platform by rapidly detecting various protein biomarkers, drugs, protein allergens, bacteria, eukaryotic cells, and distinct viruses. The linear dynamic range of NE(2)RD is five orders of magnitude broader than ELISA, with a sensitivity down to 400 fg/mL This range and sensitivity are achieved by self-assembling gold nanoparticles to generate hot spots on a 3D-oriented substrate for ultrasensitive measurements. We demonstrate that this precise platform handles multiple clinical samples such as whole blood, serum, and saliva without sample preprocessing under diverse conditions of temperature, pH, and ionic strength. The NE(2)RD's broad dynamic range, detection limit, and portability integrated with a disposable fluidic chip have broad applications, potentially enabling the transition toward precision medicine at the point-of-care or primary care settings and at patients' homes.
View details for DOI 10.1073/pnas.1510824112
View details for PubMedID 26195743
- Multitarget, quantitative nanoplasmonic electrical field-enhanced resonating device ((NERD)-R-2) for diagnostics PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA 2015; 112 (32): E4354-E4363
Magnetic levitation of single cells.
Proceedings of the National Academy of Sciences of the United States of America
2015; 112 (28): E3661-8
Several cellular events cause permanent or transient changes in inherent magnetic and density properties of cells. Characterizing these changes in cell populations is crucial to understand cellular heterogeneity in cancer, immune response, infectious diseases, drug resistance, and evolution. Although magnetic levitation has previously been used for macroscale objects, its use in life sciences has been hindered by the inability to levitate microscale objects and by the toxicity of metal salts previously applied for levitation. Here, we use magnetic levitation principles for biological characterization and monitoring of cells and cellular events. We demonstrate that each cell type (i.e., cancer, blood, bacteria, and yeast) has a characteristic levitation profile, which we distinguish at an unprecedented resolution of 1 × 10(-4) g⋅mL(-1). We have identified unique differences in levitation and density blueprints between breast, esophageal, colorectal, and nonsmall cell lung cancer cell lines, as well as heterogeneity within these seemingly homogenous cell populations. Furthermore, we demonstrate that changes in cellular density and levitation profiles can be monitored in real time at single-cell resolution, allowing quantification of heterogeneous temporal responses of each cell to environmental stressors. These data establish density as a powerful biomarker for investigating living systems and their responses. Thereby, our method enables rapid, density-based imaging and profiling of single cells with intriguing applications, such as label-free identification and monitoring of heterogeneous biological changes under various physiological conditions, including antibiotic or cancer treatment in personalized medicine.
View details for DOI 10.1073/pnas.1509250112
View details for PubMedID 26124131
The Nuclear PolyA-Binding Protein Nab2p Is Essential for mRNA Production
2015; 12 (1): 128-139
Polyadenylation of mRNA is a key step in eukaryotic gene expression. However, despite the major impact of poly(A) tails on mRNA metabolism, the precise roles of poly(A)-binding proteins (PABPs) in nuclear mRNA biogenesis remain elusive. Here, we demonstrate that rapid nuclear depletion of the S. cerevisiae PABP Nab2p leads to a global loss of cellular mRNA, but not of RNA lacking poly(A) tails. Disappearance of mRNA is a nuclear event, but not due to decreased transcription. Instead, the absence of Nab2p results in robust nuclear mRNA decay by the ribonucleolytic RNA exosome in a polyadenylation-dependent process. We conclude that Nab2p is required to protect early mRNA and therefore constitutes a crucial nuclear mRNA biogenesis factor.
View details for DOI 10.1016/j.celrep.2015.06.008
View details for Web of Science ID 000357673300014
View details for PubMedID 26119729
Yeast as a system for modeling mitochondrial disease mechanisms and discovering therapies.
Disease models & mechanisms
2015; 8 (6): 509-526
Mitochondrial diseases are severe and largely untreatable. Owing to the many essential processes carried out by mitochondria and the complex cellular systems that support these processes, these diseases are diverse, pleiotropic, and challenging to study. Much of our current understanding of mitochondrial function and dysfunction comes from studies in the baker's yeast Saccharomyces cerevisiae. Because of its good fermenting capacity, S. cerevisiae can survive mutations that inactivate oxidative phosphorylation, has the ability to tolerate the complete loss of mitochondrial DNA (a property referred to as 'petite-positivity'), and is amenable to mitochondrial and nuclear genome manipulation. These attributes make it an excellent model system for studying and resolving the molecular basis of numerous mitochondrial diseases. Here, we review the invaluable insights this model organism has yielded about diseases caused by mitochondrial dysfunction, which ranges from primary defects in oxidative phosphorylation to metabolic disorders, as well as dysfunctions in maintaining the genome or in the dynamics of mitochondria. Owing to the high level of functional conservation between yeast and human mitochondrial genes, several yeast species have been instrumental in revealing the molecular mechanisms of pathogenic human mitochondrial gene mutations. Importantly, such insights have pointed to potential therapeutic targets, as have genetic and chemical screens using yeast.
View details for DOI 10.1242/dmm.020438
View details for PubMedID 26035862
- Yeast as a system for modeling mitochondrial disease mechanisms and discovering therapies DISEASE MODELS & MECHANISMS 2015; 8 (6): 509-526
Temporal expression profiling identifies pathways mediating effect of causal variant on phenotype.
2015; 11 (6)
Even with identification of multiple causal genetic variants for common human diseases, understanding the molecular processes mediating the causal variants' effect on the disease remains a challenge. This understanding is crucial for the development of therapeutic strategies to prevent and treat disease. While static profiling of gene expression is primarily used to get insights into the biological bases of diseases, it makes differentiating the causative from the correlative effects difficult, as the dynamics of the underlying biological processes are not monitored. Using yeast as a model, we studied genome-wide gene expression dynamics in the presence of a causal variant as the sole genetic determinant, and performed allele-specific functional validation to delineate the causal effects of the genetic variant on the phenotype. Here, we characterized the precise genetic effects of a functional MKT1 allelic variant in sporulation efficiency variation. A mathematical model describing meiotic landmark events and conditional activation of MKT1 expression during sporulation specified an early meiotic role of this variant. By analyzing the early meiotic genome-wide transcriptional response, we demonstrate an MKT1-dependent role of novel modulators, namely, RTG1/3, regulators of mitochondrial retrograde signaling, and DAL82, regulator of nitrogen starvation, in additively effecting sporulation efficiency. In the presence of functional MKT1 allele, better respiration during early sporulation was observed, which was dependent on the mitochondrial retrograde regulator, RTG3. Furthermore, our approach showed that MKT1 contributes to sporulation independent of Puf3, an RNA-binding protein that steady-state transcription profiling studies have suggested to mediate MKT1-pleiotropic effects during mitotic growth. These results uncover interesting regulatory links between meiosis and mitochondrial retrograde signaling. In this study, we highlight the advantage of analyzing allele-specific transcriptional dynamics of mediating genes. Applications in higher eukaryotes can be valuable for inferring causal molecular pathways underlying complex dynamic processes, such as development, physiology and disease progression.
View details for DOI 10.1371/journal.pgen.1005195
View details for PubMedID 26039065
Variation in Crossover Frequencies Perturb Crossover Assurance Without Affecting Meiotic Chromosome Segregation in Saccharomyces cerevisiae
2015; 199 (2): 399-412
The segregation of homologous chromosomes during the Meiosis I division requires an obligate crossover per homolog pair (crossover assurance). In Saccharomyces cerevisiae and mammals, Msh4 and Msh5 proteins stabilize Holliday junctions and its progenitors to facilitate crossing over. S. cerevisiae msh4/5 hypomorphs that reduce crossover levels up to twofold at specific loci on chromosomes VII, VIII, and XV without affecting homolog segregation were identified recently. We use the msh4-R676W hypomorph to ask if the obligate crossover is insulated from variation in crossover frequencies, using a S. cerevisiae S288c/YJM789 hybrid to map recombination genome-wide. The msh4-R676W hypomorph made on average 64 crossovers per meiosis compared to 94 made in wild type and 49 in the msh4Δ mutant confirming the defect seen at individual loci on a genome-wide scale. Crossover reductions in msh4-R676W and msh4Δ were significant across chromosomes regardless of size, unlike previous observations made at specific loci. The msh4-R676W hypomorph showed reduced crossover interference. Although crossover reduction in msh4-R676W is modest, 42% of the four viable spore tetrads showed nonexchange chromosomes. These results, along with modeling of crossover distribution, suggest the significant reduction in crossovers across chromosomes and the loss of interference compromises the obligate crossover in the msh4 hypomorph. The high spore viability of the msh4 hypomorph is maintained by efficient segregation of the natural nonexchange chromosomes. Our results suggest that variation in crossover frequencies can compromise the obligate crossover and also support a mechanistic role for interference in obligate crossover formation.
View details for DOI 10.1534/genetics.114.172320
View details for Web of Science ID 000349459400009
View details for PubMedID 25467183
Chromatin-dependent regulation of RNA polymerases II and III activity throughout the transcription cycle
NUCLEIC ACIDS RESEARCH
2015; 43 (2): 787-802
The particular behaviour of eukaryotic RNA polymerases along different gene regions and amongst distinct gene functional groups is not totally understood. To cast light onto the alternative active or backtracking states of RNA polymerase II, we have quantitatively mapped active RNA polymerases at a high resolution following a new biotin-based genomic run-on (BioGRO) technique. Compared with conventional profiling with chromatin immunoprecipitation, the analysis of the BioGRO profiles in Saccharomyces cerevisiae shows that RNA polymerase II has unique activity profiles at both gene ends, which are highly dependent on positioned nucleosomes. This is the first demonstration of the in vivo influence of positioned nucleosomes on transcription elongation. The particular features at the 5' end and around the polyadenylation site indicate that this polymerase undergoes extensive specific-activity regulation in the initial and final transcription elongation phases. The genes encoding for ribosomal proteins show distinctive features at both ends. BioGRO also provides the first nascentome analysis for RNA polymerase III, which indicates that transcription of tRNA genes is poorly regulated at the individual copy level. The present study provides a novel perspective of the transcription cycle that incorporates inactivation/reactivation as an important aspect of RNA polymerase dynamics.
View details for DOI 10.1093/nar/gku1349
View details for Web of Science ID 000350209000017
View details for PubMedID 25550430
The conserved histone deacetylase Rpd3 and its DNA binding subunit Ume6 control dynamic transcript architecture during mitotic growth and meiotic development
NUCLEIC ACIDS RESEARCH
2015; 43 (1): 115-128
It was recently reported that the sizes of many mRNAs change when budding yeast cells exit mitosis and enter the meiotic differentiation pathway. These differences were attributed to length variations of their untranslated regions. The function of UTRs in protein translation is well established. However, the mechanism controlling the expression of distinct transcript isoforms during mitotic growth and meiotic development is unknown. In this study, we order developmentally regulated transcript isoforms according to their expression at specific stages during meiosis and gametogenesis, as compared to vegetative growth and starvation. We employ regulatory motif prediction, in vivo protein-DNA binding assays, genetic analyses and monitoring of epigenetic amino acid modification patterns to identify a novel role for Rpd3 and Ume6, two components of a histone deacetylase complex already known to repress early meiosis-specific genes in dividing cells, in mitotic repression of meiosis-specific transcript isoforms. Our findings classify developmental stage-specific early, middle and late meiotic transcript isoforms, and they point to a novel HDAC-dependent control mechanism for flexible transcript architecture during cell growth and differentiation. Since Rpd3 is highly conserved and ubiquitously expressed in many tissues, our results are likely relevant for development and disease in higher eukaryotes.
View details for DOI 10.1093/nar/gku1185
View details for Web of Science ID 000350207100017
View details for PubMedID 25477386
Single-cell polyadenylation site mapping reveals 3' isoform choice variability.
Molecular systems biology
2015; 11 (6): 812-?
Cell-to-cell variability in gene expression is important for many processes in biology, including embryonic development and stem cell homeostasis. While heterogeneity of gene expression levels has been extensively studied, less attention has been paid to mRNA polyadenylation isoform choice. 3' untranslated regions regulate mRNA fate, and their choice is tightly controlled during development, but how 3' isoform usage varies within genetically and developmentally homogeneous cell populations has not been explored. Here, we perform genome-wide quantification of polyadenylation site usage in single mouse embryonic and neural stem cells using a novel single-cell transcriptomic method, BATSeq. By applying BATBayes, a statistical framework for analyzing single-cell isoform data, we find that while the developmental state of the cell globally determines isoform usage, single cells from the same state differ in the choice of isoforms. Notably this variation exceeds random selection with equal preference in all cells, a finding that was confirmed by RNA FISH data. Variability in 3' isoform choice has potential implications on functional cell-to-cell heterogeneity as well as utility in resolving cell populations.
View details for DOI 10.15252/msb.20156198
View details for PubMedID 26040288
Negative feedback buffers effects of regulatory variants
MOLECULAR SYSTEMS BIOLOGY
2015; 11 (1)
Mechanisms conferring robustness against regulatory variants have been controversial. Previous studies suggested widespread buffering of RNA misexpression on protein levels during translation. We do not find evidence that translational buffering is common. Instead, we find extensive buffering at the level of RNA expression, exerted through negative feedback regulation acting in trans, which reduces the effect of regulatory variants on gene expression. Our approach is based on a novel experimental design in which allelic differential expression in a yeast hybrid strain is compared to allelic differential expression in a pool of its spores. Allelic differential expression in the hybrid is due to cis-regulatory differences only. Instead, in the pool of spores allelic differential expression is not only due to cis-regulatory differences but also due to local trans effects that include negative feedback. We found that buffering through such local trans regulation is widespread, typically compensating for about 15% of cis-regulatory effects on individual genes. Negative feedback is stronger not only for essential genes, indicating its functional relevance, but also for genes with low to middle levels of expression, for which tight regulation matters most. We suggest that negative feedback is one mechanism of Waddington's canalization, facilitating the accumulation of genetic variants that might give selective advantage in different environments.
View details for Web of Science ID 000349105100009
View details for PubMedID 25634765
A high-throughput ChIP-Seq for large-scale chromatin studies
MOLECULAR SYSTEMS BIOLOGY
2015; 11 (1)
We present a modified approach of chromatin immuno-precipitation followed by sequencing (ChIP-Seq), which relies on the direct ligation of molecular barcodes to chromatin fragments, thereby permitting experimental scale-up. With Bar-ChIP now enabling the concurrent profiling of multiple DNA-protein interactions, we report the simultaneous generation of 90 ChIP-Seq datasets without any robotic instrumentation. We demonstrate that application of Bar-ChIP to a panel of Saccharomyces cerevisiae chromatin-associated mutants provides a rapid and accurate genome-wide overview of their chromatin status. Additionally, we validate the utility of this technology to derive novel biological insights by identifying a role for the Rpd3S complex in maintaining H3K14 hypo-acetylation in gene bodies. We also report an association between the presence of intragenic H3K4 tri-methylation and the emergence of cryptic transcription in a Set2 mutant. Finally, we uncover a crosstalk between H3K14 acetylation and H3K4 methylation in this mutant. These results show that Bar-ChIP enables biological discovery through rapid chromatin profiling at single-nucleosome resolution for various conditions and protein modifications at once.
View details for Web of Science ID 000349105100007
View details for PubMedID 25583149
Expression of nuclear and mitochondrial genes encoding ATP synthase is synchronized by disassembly of a multisynthetase complex.
2014; 56 (6): 763-776
In eukaryotic cells, oxidative phosphorylation involves multisubunit complexes of mixed genetic origin. Assembling these complexes requires an organelle-independent synchronizing system for the proper expression of nuclear and mitochondrial genes. Here we show that proper expression of the F1FO ATP synthase (complex V) depends on a cytosolic complex (AME) made of two aminoacyl-tRNA synthetases (cERS and cMRS) attached to an anchor protein, Arc1p. When yeast cells adapt to respiration the Snf1/4 glucose-sensing pathway inhibits ARC1 expression triggering simultaneous release of cERS and cMRS. Free cMRS and cERS relocate to the nucleus and mitochondria, respectively, to synchronize nuclear transcription and mitochondrial translation of ATP synthase genes. Strains releasing asynchronously the two aminoacyl-tRNA synthetases display aberrant expression of nuclear and mitochondrial genes encoding subunits of complex V resulting in severe defects of the oxidative phosphorylation mechanism. This work shows that the AME complex coordinates expression of enzymes that require intergenomic control.
View details for DOI 10.1016/j.molcel.2014.10.015
View details for PubMedID 25453761
Roadblock Termination by Reb1p Restricts Cryptic and Readthrough Transcription
2014; 56 (5): 667-680
Widely transcribed compact genomes must cope with the major challenge of frequent overlapping or concurrent transcription events. Efficient and timely transcription termination is crucial to control pervasive transcription and prevent transcriptional interference. In yeast, transcription termination of RNA polymerase II (RNAPII) occurs via two possible pathways that both require recognition of termination signals on nascent RNA by specific factors. We describe here an additional mechanism of transcription termination for RNAPII and demonstrate its biological significance. We show that the transcriptional activator Reb1p bound to DNA is a roadblock for RNAPII, which pauses and is ubiquitinated, thus triggering termination. Reb1p-dependent termination generates a class of cryptic transcripts that are degraded in the nucleus by the exosome. We also observed transcriptional interference between neighboring genes in the absence of Reb1p. This work demonstrates the importance of roadblock termination for controlling pervasive transcription and preventing transcription through gene regulatory regions.
View details for DOI 10.1016/j.molcel.2014.10.026
View details for Web of Science ID 000346653300008
View details for PubMedID 25479637
Transcription mediated insulation and interference direct gene cluster expression switches
In yeast, many tandemly arranged genes show peak expression in different phases of the metabolic cycle (YMC) or in different carbon sources, indicative of regulation by a bi-modal switch, but it is not clear how these switches are controlled. Using native elongating transcript analysis (NET-seq), we show that transcription itself is a component of bi-modal switches, facilitating reciprocal expression in gene clusters. HMS2, encoding a growth-regulated transcription factor, switches between sense- or antisense-dominant states that also coordinate up- and down-regulation of transcription at neighbouring genes. Engineering HMS2 reveals alternative mono-, di- or tri-cistronic and antisense transcription units (TUs), using different promoter and terminator combinations, that underlie state-switching. Promoters or terminators are excluded from functional TUs by read-through transcriptional interference, while antisense TUs insulate downstream genes from interference. We propose that the balance of transcriptional insulation and interference at gene clusters facilitates gene expression switches during intracellular and extracellular environmental change.
View details for DOI 10.7554/eLife.03635
View details for Web of Science ID 000345639000007
View details for PubMedID 25407679
Induced Mutations in Yeast Cell Populations Adapting to an Unforeseen Challenge
2014; 9 (10)
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.
View details for DOI 10.1371/journal.pone.0111133
View details for Web of Science ID 000343662800087
View details for PubMedID 25340744
A Genome-Wide Map of Mitochondrial DNA Recombination in Yeast
2014; 198 (2): 755-U428
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.
View details for DOI 10.1534/genetics.114.166637
View details for Web of Science ID 000343885300033
View details for PubMedID 25081569
Ultrasensitive proteome analysis using paramagnetic bead technology
MOLECULAR SYSTEMS BIOLOGY
2014; 10 (10)
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.
View details for DOI 10.15252/msb.20145625
View details for Web of Science ID 000344595300003
View details for PubMedID 25358341
The Not5 Subunit of the Ccr4-Not Complex Connects Transcription and Translation
2014; 10 (10)
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.
View details for DOI 10.1371/journal.pgen.1004569
View details for Web of Science ID 000344650700002
View details for PubMedID 25340856
Heritability and genetic basis of protein level variation in an outbred population
2014; 24 (8): 1363-1370
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.
View details for DOI 10.1101/gr.170506.113
View details for Web of Science ID 000339860200011
View details for PubMedID 24823668
- Genome-wide identification of transcript start and end sites by transcript isoform sequencing NATURE PROTOCOLS 2014; 9 (7): 1740-1759
Yeast Growth Plasticity Is Regulated by Environment-Specific Multi-QTL Interactions
G3-GENES GENOMES GENETICS
2014; 4 (5): 769-777
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.
View details for DOI 10.1534/g3.113.009142
View details for Web of Science ID 000336483900001
View details for PubMedID 24474169
View details for PubMedCentralID PMC4025475
Role of histone modifications and early termination in pervasive transcription and antisense-mediated gene silencing in yeast.
Nucleic acids research
2014; 42 (7): 4348-4362
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.
View details for DOI 10.1093/nar/gku100
View details for PubMedID 24497191
An Evaluation of High-Throughput Approaches to QTL Mapping in Saccharomyces cerevisiae.
2014; 196 (3): 853-865
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 or of polymorphisms between parental genomes. We address these challenges by applying three independent high-throughput approaches for QTL mapping to map the genetic variants underlying 11 phenotypes in two genetically distant Saccharomyces cerevisiae strains, namely (1) individual analysis of >700 meiotic segregants, (2) bulk segregant analysis, and (3) reciprocal hemizygosity scanning, a new genome-wide method that we developed. We reveal differences in the performance of each approach and, by combining them, identify eight polymorphic genes that affect eight different phenotypes: colony shape, flocculation, growth on two nonfermentable carbon sources, and resistance to two drugs, salt, and high temperature. Our results demonstrate the power of individual segregant analysis to dissect QTL 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.
View details for DOI 10.1534/genetics.113.160291
View details for PubMedID 24374355
- Control of Cdc28 CDK1 by a Stress-Induced lncRNA MOLECULAR CELL 2014; 53 (4): 549-561
- Alternative polyadenylation diversifies post-transcriptional regulation by selective RNA-protein interactions MOLECULAR SYSTEMS BIOLOGY 2014; 10 (2)
High-density tiling microarray analysis of the full transcriptional activity of yeast.
Methods in molecular biology (Clifton, N.J.)
2014; 1205: 257-273
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.
View details for DOI 10.1007/978-1-4939-1363-3_16
View details for PubMedID 25213250
The Role of Ctk1 Kinase in Termination of Small Non-Coding RNAs
2013; 8 (12)
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.
View details for DOI 10.1371/journal.pone.0080495
View details for Web of Science ID 000327949300022
View details for PubMedID 24324601
Gene regulation by antisense transcription.
Nature reviews. Genetics
2013; 14 (12): 880-893
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.
View details for DOI 10.1038/nrg3594
View details for PubMedID 24217315
Extensive Variation in Chromatin States Across Humans
2013; 342 (6159): 750-752
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.
View details for DOI 10.1126/science.1242510
View details for Web of Science ID 000326647600047
View details for PubMedID 24136358
Drift and conservation of differential exon usage across tissues in primate species.
Proceedings of the National Academy of Sciences of the United States of America
2013; 110 (38): 15377-15382
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.
View details for DOI 10.1073/pnas.1307202110
View details for PubMedID 24003148
View details for PubMedCentralID PMC3780897
Polyadenylation site-induced decay of upstream transcripts enforces promoter directionality
NATURE STRUCTURAL & MOLECULAR BIOLOGY
2013; 20 (8): 923-?
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.
View details for DOI 10.1038/nsmb.2640
View details for Web of Science ID 000322715300004
View details for PubMedID 23851456
Multiple Genomic Changes Associated with Reorganization of Gene Regulation and Adaptation in Yeast
MOLECULAR BIOLOGY AND EVOLUTION
2013; 30 (7): 1514-1526
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.
View details for DOI 10.1093/molbev/mst071
View details for Web of Science ID 000321056200004
View details for PubMedID 23589456
System-wide identification of RNA-binding proteins by interactome capture
2013; 8 (3): 491-500
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.
View details for DOI 10.1038/nprot.2013.020
View details for Web of Science ID 000317110600005
View details for PubMedID 23411631
An efficient method for genome-wide polyadenylation site mapping and RNA quantification.
Nucleic acids research
2013; 41 (5)
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 3'T-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.
View details for DOI 10.1093/nar/gks1249
View details for PubMedID 23295673
Genotyping 1000 yeast strains by next-generation sequencing
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.Up to 384 whole-genome libraries can be prepared from yeast cells in one week using this method, for less than 15 euros per sample. We demonstrate the robustness of this protocol by sequencing over 1000 yeast genomes at ~30x coverage. The sequence information from 768 yeast segregants derived from two divergent S. cerevisiae strains was used to generate a meiotic recombination map at unprecedented resolution. Comparisons to other datasets indicate a high conservation of recombination at a chromosome-wide scale, but differences at the local scale. Additionally, we detected a high degree of aneuploidy (3.6%) by examining the sequencing coverage in these segregants. Differences in allele frequency allowed us to attribute instances of aneuploidy to gains of chromosomes during meiosis or mitosis, both of which showed a strong tendency to missegregate specific chromosomes.Here we present a high throughput workflow to sequence genomes of large number of yeast strains at a low price. We have used this workflow to obtain recombination and aneuploidy data from hundreds of segregants, which can serve as a foundation for future studies of linkage, recombination, and chromosomal aberrations in yeast and higher eukaryotes.
View details for DOI 10.1186/1471-2164-14-90
View details for Web of Science ID 000315034700001
View details for PubMedID 23394869
Silencing of genes and alleles by RNAi in Anopheles gambiae.
Methods in molecular biology (Clifton, N.J.)
2013; 923: 161-176
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.
View details for PubMedID 22990777
Natural sequence variants of yeast environmental sensors confer cell-to-cell expression variability.
Molecular systems biology
2013; 9: 695-?
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 Pmet17-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 Pmet17-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.
View details for DOI 10.1038/msb.2013.53
View details for PubMedID 24104478
RNA Polymerase II Collision Interrupts Convergent Transcription
2012; 48 (3): 365-374
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.
View details for DOI 10.1016/j.molce1.2012.08.027
View details for Web of Science ID 000311260900006
View details for PubMedID 23041286
Extensive Degradation of RNA Precursors by the Exosome in Wild-Type Cells
2012; 48 (3): 409-421
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.
View details for DOI 10.1016/j.molce1.2012.08.018
View details for Web of Science ID 000311260900010
View details for PubMedID 23000176
easyRNASeq: a bioconductor package for processing RNA-Seq data
2012; 28 (19): 2532-2533
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.We developed 'easyRNASeq', an R package that simplifies the processing of RNA sequencing data, hiding the complex interplay of the required packages behind a single functionality.The package is implemented in R (as of version 2.15) and is available from Bioconductor (as of version 2.10) at the URL: http://bioconductor.org/packages/release/bioc/html/easyRNASeq.html, where installation and usage instructions can be email@example.com.
View details for DOI 10.1093/bioinformatics/bts477
View details for Web of Science ID 000309687500022
View details for PubMedID 22847932
Set3 HDAC Mediates Effects of Overlapping Noncoding Transcription on Gene Induction Kinetics
2012; 150 (6): 1158-1169
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.
View details for DOI 10.1016/j.cell.2012.08.016
View details for Web of Science ID 000308928100010
View details for PubMedID 22959268
Genetic Modifiers of Chromatin Acetylation Antagonize the Reprogramming of Epi-Polymorphisms
2012; 8 (9)
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.
View details for DOI 10.1371/journal.pgen.1002958
View details for Web of Science ID 000309817900033
View details for PubMedID 23028365
- Experimental Relocation of the Mitochondrial ATP9 Gene to the Nucleus Reveals Forces Underlying Mitochondrial Genome Evolution PLOS GENETICS 2012; 8 (8)
Rrp6p Controls mRNA Poly(A) Tail Length and Its Decoration with Poly(A) Binding Proteins
2012; 47 (2): 267-280
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.
View details for DOI 10.1016/j.molcel.2012.05.005
View details for Web of Science ID 000307084000013
View details for PubMedID 22683267
Insights into RNA Biology from an Atlas of Mammalian mRNA-Binding Proteins
2012; 149 (6): 1393-1406
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.
View details for DOI 10.1016/j.cell.2012.04.031
View details for Web of Science ID 000305119600021
View details for PubMedID 22658674
Experimental relocation of the mitochondrial ATP9 gene to the nucleus reveals forces underlying mitochondrial genome evolution.
2012; 8 (8)
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.
View details for DOI 10.1371/journal.pgen.1002876
View details for PubMedID 22916027
Genome-wide H4 K16 acetylation by SAS-I is deposited independently of transcription and histone exchange
NUCLEIC ACIDS RESEARCH
2012; 40 (1): 65-74
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.
View details for DOI 10.1093/nar/gkr649
View details for Web of Science ID 000298733500015
View details for PubMedID 21908408
Genome-wide polyadenylation site mapping.
Methods in enzymology
2012; 513: 271-296
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 3'UTR 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.
View details for DOI 10.1016/B978-0-12-391938-0.00012-4
View details for PubMedID 22929774
Accumulation of noncoding RNA due to an RNase P defect in Saccharomyces cerevisiae.
RNA (New York, N.Y.)
2011; 17 (8): 1441-1450
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.
View details for DOI 10.1261/rna.2737511
View details for PubMedID 21665995
A yeast-based assay identifies drugs active against human mitochondrial disorders
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2011; 108 (29): 11989-11994
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.
View details for DOI 10.1073/pnas.1101478108
View details for Web of Science ID 000292876900057
View details for PubMedID 21715656
Functional consequences of bidirectional promoters
TRENDS IN GENETICS
2011; 27 (7): 267-276
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.
View details for DOI 10.1016/j.tig.2011.04.002
View details for Web of Science ID 000292413300003
View details for PubMedID 21601935
Antisense expression increases gene expression variability and locus interdependency.
Molecular systems biology
2011; 7: 468-?
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.
View details for DOI 10.1038/msb.2011.1
View details for PubMedID 21326235
- Antisense expression increases gene expression variability and locus interdependency MOLECULAR SYSTEMS BIOLOGY 2011; 7
Execution of the meiotic noncoding RNA expression program and the onset of gametogenesis in yeast require the conserved exosome subunit Rrp6
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2011; 108 (3): 1058-1063
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.
View details for DOI 10.1073/pnas.1016459108
View details for Web of Science ID 000286310300034
View details for PubMedID 21149693
Yeast Sen1 Helicase Protects the Genome from Transcription-Associated Instability
2011; 41 (1): 21-32
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.
View details for DOI 10.1016/j.molcel.2010.12.007
View details for Web of Science ID 000286692400006
View details for PubMedID 21211720
Genome-wide transcriptome analysis in yeast using high-density tiling arrays.
Methods in molecular biology (Clifton, N.J.)
2011; 759: 107-123
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.
View details for DOI 10.1007/978-1-61779-173-4_7
View details for PubMedID 21863484
Genome-wide survey of post-meiotic segregation during yeast recombination
2011; 12 (4)
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.By genotyping tens of thousands of genetic markers in yeast segregants and their clonal progeny, we analyzed post-meiotic segregation at a genome-wide scale. We show that post-meiotic segregation occurs in close to 10% of recombination events. Although the overall number of markers affected in a single meiosis is small, the rate of post-meiotic segregation is more than five orders of magnitude larger than the base substitution mutation rate. Post-meiotic segregation took place with equal relative frequency in crossovers and non-crossovers, and usually at the edges of gene conversion tracts. Furthermore, post-meiotic segregation tended to occur in markers that are isolated from other heterozygosities and preferentially at polymorphism types that are relatively uncommon in the yeast species.Overall, our survey reveals the genome-wide characteristics of post-meiotic segregation. The results show that post-meiotic segregation is widespread in meiotic recombination and could be a significant determinant of allelic inheritance and allele frequencies at the population level.
View details for DOI 10.1186/gb-2011-12-4-r36
View details for Web of Science ID 000292871100004
View details for PubMedID 21481229
The Baker's Yeast Diploid Genome Is Remarkably Stable in Vegetative Growth and Meiosis
2010; 6 (9)
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.
View details for DOI 10.1371/journal.pgen.1001109
View details for Web of Science ID 000282369200047
View details for PubMedID 20838597
Antagonistic Changes in Sensitivity to Antifungal Drugs by Mutations of an Important ABC Transporter Gene in a Fungal Pathogen
2010; 5 (6)
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.
View details for DOI 10.1371/journal.pone.0011309
View details for Web of Science ID 000279140800010
View details for PubMedID 20593017
Genetic analysis of variation in transcription factor binding in yeast
2010; 464 (7292): 1187-U106
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.
View details for DOI 10.1038/nature08934
View details for Web of Science ID 000276891100036
View details for PubMedID 20237471
Natural Single-Nucleosome Epi-Polymorphisms in Yeast
2010; 6 (4)
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.
View details for DOI 10.1371/journal.pgen.1000913
View details for Web of Science ID 000277354200030
View details for PubMedID 20421933
High-resolution transcription atlas of the mitotic cell cycle in budding yeast
2010; 11 (3)
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.We discovered 523 antisense transcripts, of which 80 cycle or are located opposite periodically expressed mRNAs, 135 unannotated intergenic non-coding RNAs, of which 11 cycle, and 109 cell-cycle-regulated protein-coding genes that had not previously been shown to cycle. We detected periodic expression coupling of sense and antisense transcript pairs, including antisense transcripts opposite of key cell-cycle regulators, like FAR1 and TAF2.Our dataset presents the most comprehensive resource to date on gene expression during the budding yeast cell cycle. It reveals periodic expression of both protein-coding and non-coding RNA and profiles the expression of non-annotated RNAs throughout the cell cycle for the first time. This data enables hypothesis-driven mechanistic studies concerning the functions of non-coding RNAs.
View details for DOI 10.1186/gb-2010-11-3-r24
View details for Web of Science ID 000277309100007
View details for PubMedID 20193063
Trans-acting antisense RNAs mediate transcriptional gene cosuppression in S. cerevisiae
GENES & DEVELOPMENT
2009; 23 (13): 1534-1545
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.
View details for DOI 10.1101/gad.522509
View details for Web of Science ID 000267954800008
View details for PubMedID 19571181
Genome-wide allele- and strand-specific expression profiling
MOLECULAR SYSTEMS BIOLOGY
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.
View details for DOI 10.1038/msb.2009.31
View details for Web of Science ID 000267629300008
View details for PubMedID 19536197
Array-based genotyping in S. cerevisiae using semi-supervised clustering
2009; 25 (8): 1056-1062
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.We developed ssGenotyping (ssG), a multivariate, semi-supervised approach for using microarrays to genotype haploid individuals at thousands of polymorphic sites. Using a meiotic recombination dataset, we show that ssG is more accurate than existing supervised classification methods, and that it produces denser marker coverage. The ssG algorithm is able to fit probe-specific affinity differences and to detect and filter spurious signal, permitting high-confidence genotyping at nucleotide resolution. We also demonstrate that oligonucleotide probe response depends significantly on genomic background, even when the probe's specific target sequence is unchanged. As a result, supervised classifiers trained on reference strains may not generalize well to diverged strains; ssG's semi-supervised approach, on the other hand, adapts automatically.
View details for DOI 10.1093/bioinformatics/btp104
View details for Web of Science ID 000265094400011
View details for PubMedID 19237444
Widespread bidirectional promoters are the major source of cryptic transcripts in yeast
2009; 457 (7232): 1038-U8
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.
View details for DOI 10.1038/nature07747
View details for Web of Science ID 000263425400048
View details for PubMedID 19169244
Identification of mitochondrial disease genes through integrative analysis of multiple datasets
2008; 46 (4): 248-255
Determining the genetic factors in a disease is crucial to elucidating its molecular basis. This task is challenging due to a lack of information on gene function. The integration of large-scale functional genomics data has proven to be an effective strategy to prioritize candidate disease genes. Mitochondrial disorders are a prevalent and heterogeneous class of diseases that are particularly amenable to this approach. Here we explain the application of integrative approaches to the identification of mitochondrial disease genes. We first examine various datasets that can be used to evaluate the involvement of each gene in mitochondrial function. The data integration methodology is then described, accompanied by examples of common implementations. Finally, we discuss how gene networks are constructed using integrative techniques and applied to candidate gene prioritization. Relevant public data resources are indicated. This report highlights the success and potential of data integration as well as its applicability to the search for mitochondrial disease genes.
View details for DOI 10.1016/j.ymeth.2008.10.002
View details for Web of Science ID 000261647200003
View details for PubMedID 18930150
Sequential Elimination of Major-Effect Contributors Identifies Additional Quantitative Trait Loci Conditioning High-Temperature Growth in Yeast
2008; 180 (3): 1661-1670
Several quantitative trait loci (QTL) mapping strategies can successfully identify major-effect loci, but often have poor success detecting loci with minor effects, potentially due to the confounding effects of major loci, epistasis, and limited sample sizes. To overcome such difficulties, we used a targeted backcross mapping strategy that genetically eliminated the effect of a previously identified major QTL underlying high-temperature growth (Htg) in yeast. This strategy facilitated the mapping of three novel QTL contributing to Htg of a clinically derived yeast strain. One QTL, which is linked to the previously identified major-effect QTL, was dissected, and NCS2 was identified as the causative gene. The interaction of the NCS2 QTL with the first major-effect QTL was background dependent, revealing a complex QTL architecture spanning these two linked loci. Such complex architecture suggests that more genes than can be predicted are likely to contribute to quantitative traits. The targeted backcrossing approach overcomes the difficulties posed by sample size, genetic linkage, and epistatic effects and facilitates identification of additional alleles with smaller contributions to complex traits.
View details for DOI 10.1534/genetics.108.092932
View details for Web of Science ID 000261036200033
View details for PubMedID 18780730
Systematic screens for human disease genes, from yeast to human and back
2008; 4 (1): 18-29
Systematic screens for human disease genes have emerged in recent years, due to the wealth of information provided by genome sequences and large scale datasets. Here we review how integration of genomic data in yeast and human is helping to elucidate the genetic basis of mitochondrial diseases. The identification of nearly all yeast mitochondrial proteins and many of their functional interactions provides insight into the role of mitochondria in cellular processes. This information enables prioritization of the candidate genes underlying mitochondrial disorders. In an iterative fashion, the link between predicted human candidate genes and their disease phenotypes can be experimentally tested back in yeast.
View details for DOI 10.1039/b709494a
View details for Web of Science ID 000251673300002
View details for PubMedID 18075670
Mosaic Genome Architecture of the Anopheles gambiae Species Complex
2007; 2 (11)
Attempts over the last three decades to reconstruct the phylogenetic history of the Anopheles gambiae species complex have been important for developing better strategies to control malaria transmission.We used fingerprint genotyping data from 414 field-collected female mosquitoes at 42 microsatellite loci to infer the evolutionary relationships of four species in the A. gambiae complex, the two major malaria vectors A. gambiae sensu stricto (A. gambiae s.s.) and A. arabiensis, as well as two minor vectors, A. merus and A. melas.We identify six taxonomic units, including a clear separation of West and East Africa A. gambiae s.s. S molecular forms. We show that the phylogenetic relationships vary widely between different genomic regions, thus demonstrating the mosaic nature of the genome of these species. The two major malaria vectors are closely related and closer to A. merus than to A. melas at the genome-wide level, which is also true if only autosomes are considered. However, within the Xag inversion region of the X chromosome, the M and two S molecular forms are most similar to A. merus. Near the X centromere, outside the Xag region, the two S forms are highly dissimilar to the other taxa. Furthermore, our data suggest that the centromeric region of chromosome 3 is a strong discriminator between the major and minor malaria vectors.Although further studies are needed to elucidate the basis of the phylogenetic variation among the different regions of the genome, the preponderance of sympatric admixtures among taxa strongly favor introgression of different genomic regions between species, rather than lineage sorting of ancestral polymorphism, as a possible mechanism.
View details for DOI 10.1371/journal.pone.0001249
View details for Web of Science ID 000207459300030
View details for PubMedID 18043756
Antisense artifacts in transcriptome microarray experiments are resolved by actinomycin D
NUCLEIC ACIDS RESEARCH
2007; 35 (19)
Recent transcription profiling studies have revealed an unanticipatedly large proportion of antisense transcription across eukaryotic and bacterial genomes. However, the extent and significance of antisense transcripts is controversial partly because experimental artifacts are suspected. Here, we present a method to generate clean genome-wide transcriptome profiles, using actinomycin D (ActD) during reverse transcription. We show that antisense artifacts appear to be triggered by spurious synthesis of second-strand cDNA during reverse transcription reactions. Strand-specific hybridization signals obtained from Saccharomyces cerevisiae tiling arrays were compared between samples prepared with and without ActD. Use of ActD removed about half of the detectable antisense transcripts, consistent with their being artifacts, while sense expression levels and about 200 antisense transcripts were not affected. Our findings thus facilitate a more accurate assessment of the extent and position of antisense transcription, towards a better understanding of its role in cells.
View details for DOI 10.1093/nar/gkm683
View details for Web of Science ID 000250827000034
View details for PubMedID 17897965
Assessing systems properties of yeast mitochondria through an interaction map of the organelle
2006; 2 (10): 1612-1624
Mitochondria carry out specialized functions; compartmentalized, yet integrated into the metabolic and signaling processes of the cell. Although many mitochondrial proteins have been identified, understanding their functional interrelationships has been a challenge. Here we construct a comprehensive network of the mitochondrial system. We integrated genome-wide datasets to generate an accurate and inclusive mitochondrial parts list. Together with benchmarked measures of protein interactions, a network of mitochondria was constructed in their cellular context, including extra-mitochondrial proteins. This network also integrates data from different organisms to expand the known mitochondrial biology beyond the information in the existing databases. Our network brings together annotated and predicted functions into a single framework. This enabled, for the entire system, a survey of mutant phenotypes, gene regulation, evolution, and disease susceptibility. Furthermore, we experimentally validated the localization of several candidate proteins and derived novel functional contexts for hundreds of uncharacterized proteins. Our network thus advances the understanding of the mitochondrial system in yeast and identifies properties of genes underlying human mitochondrial disorders.
View details for DOI 10.1371/journal.pgen.0020170
View details for Web of Science ID 000241922200011
View details for PubMedID 17054397
Transcript mapping with high-density oligonucleotide tiling arrays
2006; 22 (16): 1963-1970
High-density DNA tiling microarrays are a powerful tool for the characterization of complete transcriptomes. The two major analytical challenges are the segmentation of the hybridization signal along genomic coordinates to accurately determine transcript boundaries and the adjustment of the sequence-dependent response of the oligonucleotide probes to achieve quantitative comparability of the signal between different probes.We describe a dynamic programming algorithm for finding a globally optimal fit of a piecewise constant expression profile along genomic coordinates. We developed a probe-specific background correction and scaling method that employs empirical probe response parameters determined from reference hybridizations with no need for paired mismatch probes. This combined analysis approach allows the accurate determination of dynamical changes in transcription architectures from hybridization data and will help to study the biological significance of complex transcriptional phenomena in eukaryotic genomes.R package tilingArray at http://www.bioconductor.org.
View details for DOI 10.1093/bioinformatics/btl289
View details for Web of Science ID 000239900200006
View details for PubMedID 16787969
Capturing cellular machines by systematic screens of protein complexes
TRENDS IN MICROBIOLOGY
2006; 14 (8): 336-339
Two recent studies have provided the most complete screen for protein complexes in yeast to date, in which partners were identified for approximately half of the proteome. A comparison shows that these two datasets are complementary. In addition, one of the analyses points to a modular organization of the cellular protein network. These data will prove useful in defining principles and trends that arise when combining large-scale datasets of different natures, and in deriving properties of protein machines in cellular systems.
View details for DOI 10.1016/j.tim.2006.06.002
View details for Web of Science ID 000239886600003
View details for PubMedID 16782340
Complex genetic interactions in a quantitative trait locus
2006; 2 (2): 140-147
Whether in natural populations or between two unrelated members of a species, most phenotypic variation is quantitative. To analyze such quantitative traits, one must first map the underlying quantitative trait loci. Next, and far more difficult, one must identify the quantitative trait genes (QTGs), characterize QTG interactions, and identify the phenotypically relevant polymorphisms to determine how QTGs contribute to phenotype. In this work, we analyzed three Saccharomyces cerevisiae high-temperature growth (Htg) QTGs (MKT1, END3, and RHO2). We observed a high level of genetic interactions among QTGs and strain background. Interestingly, while the MKT1 and END3 coding polymorphisms contribute to phenotype, it is the RHO2 3'UTR polymorphisms that are phenotypically relevant. Reciprocal hemizygosity analysis of the Htg QTGs in hybrids between S288c and ten unrelated S. cerevisiae strains reveals that the contributions of the Htg QTGs are not conserved in nine other hybrids, which has implications for QTG identification by marker-trait association. Our findings demonstrate the variety and complexity of QTG contributions to phenotype, the impact of genetic background, and the value of quantitative genetic studies in S. cerevisiae.
View details for DOI 10.1371/journal.pgen.0020013
View details for Web of Science ID 000239494300006
View details for PubMedID 16462944
Re-analysis of data and its integration
2005; 579 (8): 1802-1807
To understand a biological process it is clear that a single approach will not be sufficient, just like a single measurement on a protein--such as its expression level--does not describe protein function. Using reference sets of proteins as benchmarks different approaches can be scaled and integrated. Here, we demonstrate the power of data re-analysis and integration by applying it in a case study to data from deletion phenotype screens and mRNA expression profiling.
View details for DOI 10.1016/j.febslet.2005.02.006
View details for Web of Science ID 000227657800007
View details for PubMedID 15763555
Elevated evolutionary rates in the laboratory strain of Saccharomyces cerevisiae
PROCEEDINGS OF THE NATIONAL ACADEMY OF SCIENCES OF THE UNITED STATES OF AMERICA
2005; 102 (4): 1092-1097
By using the maximum likelihood method, we made a genome-wide comparison of the evolutionary rates in the lineages leading to the laboratory strain (S288c) and a wild strain (YJM789) of Saccharomyces cerevisiae and found that genes in the laboratory strain tend to evolve faster than in the wild strain. The pattern of elevated evolution suggests that relaxation of selection intensity is the dominant underlying reason, which is consistent with recurrent bottlenecks in the S. cerevisiae laboratory strain population. Supporting this conclusion are the following observations: (i) the increases in nonsynonymous evolutionary rate occur for genes in all functional categories; (ii) most of the synonymous evolutionary rate increases in S288c occur in genes with strong codon usage bias; (iii) genes under stronger negative selection have a larger increase in nonsynonymous evolutionary rate; and (iv) more genes with adaptive evolution were detected in the laboratory strain, but they do not account for the majority of the increased evolution. The present discoveries suggest that experimental and possible industrial manipulations of the laboratory strain of yeast could have had a strong effect on the genetic makeup of this model organism. Furthermore, they imply an evolution of laboratory model organisms away from their wild counterparts, questioning the relevancy of the models especially when extensive laboratory cultivation has occurred. In addition, these results shed light on the evolution of livestock and crop species that have been under human domestication for years.
View details for DOI 10.1073/pnas.0409159102
View details for Web of Science ID 000226617900026
View details for PubMedID 15647350
Integrative analysis of the mitochondrial proteome in yeast
2004; 2 (6): 795-804
In this study yeast mitochondria were used as a model system to apply, evaluate, and integrate different genomic approaches to define the proteins of an organelle. Liquid chromatography mass spectrometry applied to purified mitochondria identified 546 proteins. By expression analysis and comparison to other proteome studies, we demonstrate that the proteomic approach identifies primarily highly abundant proteins. By expanding our evaluation to other types of genomic approaches, including systematic deletion phenotype screening, expression profiling, subcellular localization studies, protein interaction analyses, and computational predictions, we show that an integration of approaches moves beyond the limitations of any single approach. We report the success of each approach by benchmarking it against a reference set of known mitochondrial proteins, and predict approximately 700 proteins associated with the mitochondrial organelle from the integration of 22 datasets. We show that a combination of complementary approaches like deletion phenotype screening and mass spectrometry can identify over 75% of the known mitochondrial proteome. These findings have implications for choosing optimal genome-wide approaches for the study of other cellular systems, including organelles and pathways in various species. Furthermore, our systematic identification of genes involved in mitochondrial function and biogenesis in yeast expands the candidate genes available for mapping Mendelian and complex mitochondrial disorders in humans.
View details for DOI 10.1371/journal.pbio.0020160
View details for Web of Science ID 000222380400019
View details for PubMedID 15208715
- Maximizing the potential of functional genomics NATURE REVIEWS GENETICS 2004; 5 (3): 190-201
Role of duplicate genes in genetic robustness against null mutations
2003; 421 (6918): 63-66
Deleting a gene in an organism often has little phenotypic effect, owing to two mechanisms of compensation. The first is the existence of duplicate genes: that is, the loss of function in one copy can be compensated by the other copy or copies. The second mechanism of compensation stems from alternative metabolic pathways, regulatory networks, and so on. The relative importance of the two mechanisms has not been investigated except for a limited study, which suggested that the role of duplicate genes in compensation is negligible. The availability of fitness data for a nearly complete set of single-gene-deletion mutants of the Saccharomyces cerevisiae genome has enabled us to carry out a genome-wide evaluation of the role of duplicate genes in genetic robustness against null mutations. Here we show that there is a significantly higher probability of functional compensation for a duplicate gene than for a singleton, a high correlation between the frequency of compensation and the sequence similarity of two duplicates, and a higher probability of a severe fitness effect when the duplicate copy that is more highly expressed is deleted. We estimate that in S. cerevisiae at least a quarter of those gene deletions that have no phenotype are compensated by duplicate genes.
View details for DOI 10.1038/nature01198
View details for Web of Science ID 000180165500037
View details for PubMedID 12511954
Gene function on a genomic scale
JOURNAL OF CHROMATOGRAPHY B-ANALYTICAL TECHNOLOGIES IN THE BIOMEDICAL AND LIFE SCIENCES
2002; 782 (1-2): 151-163
The ability to obtain experimental measurements for thousands of genes has revolutionized our view of biological systems. While traditional studies of gene function evaluated many different properties for a single gene, genomic approaches can measure a single property for thousands of genes. Over the last years, genomic approaches have been developed to measure many different properties, including gene expression, deletion phenotype, and protein characteristics. The promise of integrating these datasets has made it attractive to test whether genomic approaches can determine gene function with accuracy comparable to single gene approaches.
View details for Web of Science ID 000179785100013
View details for PubMedID 12458004
Systematic screen for human disease genes in yeast
2002; 31 (4): 400-404
High similarity between yeast and human mitochondria allows functional genomic study of Saccharomyces cerevisiae to be used to identify human genes involved in disease. So far, 102 heritable disorders have been attributed to defects in a quarter of the known nuclear-encoded mitochondrial proteins in humans. Many mitochondrial diseases remain unexplained, however, in part because only 40-60% of the presumed 700-1,000 proteins involved in mitochondrial function and biogenesis have been identified. Here we apply a systematic functional screen using the pre-existing whole-genome pool of yeast deletion mutants to identify mitochondrial proteins. Three million measurements of strain fitness identified 466 genes whose deletions impaired mitochondrial respiration, of which 265 were new. Our approach gave higher selection than other systematic approaches, including fivefold greater selection than gene expression analysis. To apply these advantages to human disorders involving mitochondria, human orthologs were identified and linked to heritable diseases using genomic map positions.
View details for DOI 10.1038/ng929
View details for Web of Science ID 000177147100016
View details for PubMedID 12134146
Evolutionary rate in the protein interaction network
2002; 296 (5568): 750-752
High-throughput screens have begun to reveal the protein interaction network that underpins most cellular functions in the yeast Saccharomyces cerevisiae. How the organization of this network affects the evolution of the proteins that compose it is a fundamental question in molecular evolution. We show that the connectivity of well-conserved proteins in the network is negatively correlated with their rate of evolution. Proteins with more interactors evolve more slowly not because they are more important to the organism, but because a greater proportion of the protein is directly involved in its function. At sites important for interaction between proteins, evolutionary changes may occur largely by coevolution, in which substitutions in one protein result in selection pressure for reciprocal changes in interacting partners. We confirm one predicted outcome of this process-namely, that interacting proteins evolve at similar rates.
View details for Web of Science ID 000175281700060
View details for PubMedID 11976460
Dissecting the architecture of a quantitative trait locus in yeast
2002; 416 (6878): 326-330
Most phenotypic diversity in natural populations is characterized by differences in degree rather than in kind. Identification of the actual genes underlying these quantitative traits has proved difficult. As a result, little is known about their genetic architecture. The failures are thought to be due to the different contributions of many underlying genes to the phenotype and the ability of different combinations of genes and environmental factors to produce similar phenotypes. This study combined genome-wide mapping and a new genetic technique named reciprocal-hemizygosity analysis to achieve the complete dissection of a quantitative trait locus (QTL) in Saccharomyces cerevisiae. A QTL architecture was uncovered that was more complex than expected. Functional linkages both in cis and in trans were found between three tightly linked quantitative trait genes that are neither necessary nor sufficient in isolation. This arrangement of alleles explains heterosis (hybrid vigour), the increased fitness of the heterozygote compared with homozygotes. It also demonstrates a deficiency in current approaches to QTL dissection with implications extending to traits in other organisms, including human genetic diseases.
View details for Web of Science ID 000174482200044
View details for PubMedID 11907579
Transcriptional regulation and function during the human cell cycle
2001; 27 (1): 48-54
We report here the transcriptional profiling of the cell cycle on a genome-wide scale in human fibroblasts. We identified approximately 700 genes that display transcriptional fluctuation with a periodicity consistent with that of the cell cycle. Systematic analysis of these genes revealed functional organization within groups of coregulated transcripts. A diverse set of cytoskeletal reorganization genes exhibit cell-cycle-dependent regulation, indicating that biological pathways are redirected for the execution of cell division. Many genes involved in cell motility and remodeling of the extracellular matrix are expressed predominantly in M phase, indicating a mechanism for balancing proliferative and invasive cellular behavior. Transcripts upregulated during S phase displayed extensive overlap with genes induced by DNA damage; cell-cycle-regulated transcripts may therefore constitute coherent programs used in response to external stimuli. Our data also provide clues to biological function for hundreds of previously uncharacterized human genes.
View details for Web of Science ID 000166187900014
View details for PubMedID 11137997
- Combining genome sequences and new technologies for dissecting the genetics of complex phenotypes TRENDS IN PLANT SCIENCE 2000; 5 (9): 397-401
- High-density arrays and insights into genome function BIOTECHNOLOGY & GENETIC ENGINEERING REVIEWS, VOL 17 2000; 17: 109-146
A genome-wide transcriptional analysis of the mitotic cell cycle
1998; 2 (1): 65-73
Progression through the eukaryotic cell cycle is known to be both regulated and accompanied by periodic fluctuation in the expression levels of numerous genes. We report here the genome-wide characterization of mRNA transcript levels during the cell cycle of the budding yeast S. cerevisiae. Cell cycle-dependent periodicity was found for 416 of the 6220 monitored transcripts. More than 25% of the 416 genes were found directly adjacent to other genes in the genome that displayed induction in the same cell cycle phase, suggesting a mechanism for local chromosomal organization in global mRNA regulation. More than 60% of the characterized genes that displayed mRNA fluctuation have already been implicated in cell cycle period-specific biological roles. Because more than 20% of human proteins display significant homology to yeast proteins, these results also link a range of human genes to cell cycle period-specific biological functions.
View details for Web of Science ID 000075174500007
View details for PubMedID 9702192