Temporal Expression Profiling Identifies Pathways Mediating Effect of Causal Variant on Phenotype

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.

[1]  C. Clayton,et al.  Trypanosome MKT1 and the RNA-binding protein ZC3H11: interactions and potential roles in post-transcriptional regulatory networks , 2014, Nucleic acids research.

[2]  Ben Lehner,et al.  The effects of genetic variation on gene expression dynamics during development , 2013, Nature.

[3]  B. Kaufman,et al.  Aconitase Couples Metabolic Regulation to Mitochondrial DNA Maintenance , 2005, Science.

[4]  B. Conne,et al.  Dcp1-bodies in mouse oocytes. , 2009, Molecular biology of the cell.

[5]  Trey Ideker,et al.  Genotype to phenotype via network analysis. , 2013, Current opinion in genetics & development.

[6]  B. Cohen,et al.  Causal Variation in Yeast Sporulation Tends to Reside in a Pathway Bottleneck , 2014, PLoS genetics.

[7]  A. Battle,et al.  Determining causality and consequence of expression quantitative trait loci , 2014, Human Genetics.

[8]  W. Saunders,et al.  Large-scale functional genomic analysis of sporulation and meiosis in Saccharomyces cerevisiae. , 2003, Genetics.

[9]  Wolfgang Huber,et al.  Antisense expression increases gene expression variability and locus interdependency , 2011, Molecular systems biology.

[10]  R. A. Butow,et al.  Mitochondrial retrograde signaling. , 2006, Annual review of genetics.

[11]  J R Johnston,et al.  Genealogy of principal strains of the yeast genetic stock center. , 1986, Genetics.

[12]  A. Neiman,et al.  Ascospore Formation in the Yeast Saccharomyces cerevisiae , 2005, Microbiology and Molecular Biology Reviews.

[13]  Robert P. Davey,et al.  Population genomics of domestic and wild yeasts , 2008, Nature.

[14]  R. D. Gietz,et al.  Transformation of yeast by lithium acetate/single-stranded carrier DNA/polyethylene glycol method. , 2002, Methods in enzymology.

[15]  A. Amon,et al.  Regulation of entry into gametogenesis , 2011, Philosophical Transactions of the Royal Society B: Biological Sciences.

[16]  Ben Lehner,et al.  Predicting phenotypic variation from genotypes, phenotypes and a combination of the two. , 2013, Current opinion in biotechnology.

[17]  Julien Gagneur,et al.  Genotype-Environment Interactions Reveal Causal Pathways That Mediate Genetic Effects on Phenotype , 2013, PLoS genetics.

[18]  Lude Franke,et al.  From genome to function by studying eQTLs. , 2014, Biochimica et biophysica acta.

[19]  Michal Linial,et al.  Novel Unsupervised Feature Filtering of Biological Data , 2006, ISMB.

[20]  L. Liang,et al.  Mapping complex disease traits with global gene expression , 2009, Nature Reviews Genetics.

[21]  C. Jacq,et al.  Posttranscriptional control of mitochondrial biogenesis: Spatio‐temporal regulation of the protein import process , 2010, FEBS letters.

[22]  N. Turner,et al.  Mitochondrial Stress Signaling Promotes Cellular Adaptations , 2014, International journal of cell biology.

[23]  Ronald A. Butow,et al.  RTG1 and RTG2: Two yeast genes required for a novel path of communication from mitochondria to the nucleus , 1993, Cell.

[24]  I. Simon,et al.  Studying and modelling dynamic biological processes using time-series gene expression data , 2012, Nature Reviews Genetics.

[25]  J. Fay,et al.  Dissecting the pleiotropic consequences of a quantitative trait nucleotide. , 2009, FEMS yeast research.

[26]  J. Fay The molecular basis of phenotypic variation in yeast. , 2013, Current opinion in genetics & development.

[27]  Rachel B. Brem,et al.  Integrating large-scale functional genomic data to dissect the complexity of yeast regulatory networks , 2008, Nature Genetics.

[28]  E. Dermitzakis,et al.  From expression QTLs to personalized transcriptomics , 2011, Nature Reviews Genetics.

[29]  L. Kruglyak,et al.  Gene–Environment Interaction in Yeast Gene Expression , 2008, PLoS biology.

[30]  P. Ye,et al.  Characterization of the Metabolic Requirements in Yeast Meiosis , 2013, PloS one.

[31]  Jeffrey T. Leek,et al.  EDGE: extraction and analysis of differential gene expression , 2006, Bioinform..

[32]  E. Génin,et al.  How important are rare variants in common disease? , 2014, Briefings in functional genomics.

[33]  Daniel R. Richards,et al.  Dissecting the architecture of a quantitative trait locus in yeast , 2002, Nature.

[34]  J. Broach Nutritional Control of Growth and Development in Yeast , 2012, Genetics.

[35]  A. Gasch,et al.  Genetic Architecture of Ethanol-Responsive Transcriptome Variation in Saccharomyces cerevisiae Strains , 2014, Genetics.

[36]  Rainer Breitling,et al.  Expression Quantitative Trait Loci Are Highly Sensitive to Cellular Differentiation State , 2009, PLoS genetics.

[37]  S. M. Honigberg,et al.  Signal pathway integration in the switch from the mitotic cell cycle to meiosis in yeast , 2003, Journal of Cell Science.

[38]  T. Cooper,et al.  Roles of the Dal82p Domains in Allophanate/Oxalurate-dependent Gene Expression inSaccharomyces cerevisiae * , 2000, The Journal of Biological Chemistry.

[39]  Ben Lehner Genotype to phenotype: lessons from model organisms for human genetics , 2013, Nature Reviews Genetics.

[40]  Himanshu Sinha,et al.  Complex Genetic Interactions in a Quantitative Trait Locus , 2006, PLoS genetics.

[41]  J. Mariani,et al.  Changes in the Distribution of the α3 Na+/K+ ATPase Subunit in Heterozygous Lurcher Purkinje Cells as a Genetic Model of Chronic Depolarization during Development , 2014, International journal of cell biology.

[42]  C. Loader Local Likelihood Density Estimation , 1996 .

[43]  Ronald W Davis,et al.  Parallel phenotypic analysis of sporulation and postgermination growth in Saccharomyces cerevisiae , 2002, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Martin Vingron,et al.  Variance stabilization applied to microarray data calibration and to the quantification of differential expression , 2002, ISMB.

[45]  P. Reynier,et al.  Low oocyte mitochondrial DNA content in ovarian insufficiency. , 2005, Human reproduction.

[46]  G. Bhanot,et al.  Sporulation Genes Associated with Sporulation Efficiency in Natural Isolates of Yeast , 2013, PloS one.

[47]  Riccardo Bellazzi,et al.  TimeClust: a clustering tool for gene expression time series , 2008, Bioinform..

[48]  J. Mccusker,et al.  Three new dominant drug resistance cassettes for gene disruption in Saccharomyces cerevisiae , 1999, Yeast.

[49]  Leopold Parts,et al.  Population genomics of domestic and wild yeasts , 2008 .

[50]  Roy Parker,et al.  P bodies promote stress granule assembly in Saccharomyces cerevisiae , 2008, The Journal of cell biology.

[51]  J. Thornton,et al.  Mitochondria-to-nuclear signaling is regulated by the subcellular localization of the transcription factors Rtg1p and Rtg3p. , 2000, Molecular biology of the cell.

[52]  A. Lusis,et al.  Systems genetics approaches to understand complex traits , 2013, Nature Reviews Genetics.

[53]  P. Anderson,et al.  Mammalian stress granules and processing bodies. , 2007, Methods in enzymology.

[54]  Anders Blomberg,et al.  Trait Variation in Yeast Is Defined by Population History , 2011, PLoS genetics.

[55]  B. Rothermel,et al.  A basic helix-loop-helix-leucine zipper transcription complex in yeast functions in a signaling pathway from mitochondria to the nucleus , 1997, Molecular and cellular biology.

[56]  Leonid Kruglyak,et al.  Identification and Dissection of a Complex DNA Repair Sensitivity Phenotype in Baker's Yeast , 2008, PLoS genetics.

[57]  H. Bussemaker,et al.  Harnessing Natural Sequence Variation to Dissect Posttranscriptional Regulatory Networks in Yeast , 2014, G3: Genes, Genomes, Genetics.

[58]  Yang Yang,et al.  Positive Feedback of NDT80 Expression Ensures Irreversible Meiotic Commitment in Budding Yeast , 2014, PLoS genetics.

[59]  Rachel B. Brem,et al.  Polymorphisms in Multiple Genes Contribute to the Spontaneous Mitochondrial Genome Instability of Saccharomyces cerevisiae S288C Strains , 2009, Genetics.

[60]  M. Sorokin,et al.  Mitochondrial signaling in Saccharomyces cerevisiae pseudohyphae formation induced by butanol. , 2013, FEMS yeast research.

[61]  Y. Kassir,et al.  IME1, a positive regulator gene of meiosis in S. cerevisiae , 1988, Cell.

[62]  Michael D. Eastwood,et al.  Developmentally programmed nuclear destruction during yeast gametogenesis. , 2012, Developmental cell.

[63]  D. Thompson,et al.  Determinants of Divergent Adaptation and Dobzhansky-Muller Interaction in Experimental Yeast Populations , 2010, Current Biology.

[64]  P. Brown,et al.  Extensive Association of Functionally and Cytotopically Related mRNAs with Puf Family RNA-Binding Proteins in Yeast , 2004, PLoS biology.

[65]  Ronald W. Davis,et al.  The core meiotic transcriptome in budding yeasts , 2000, Nature Genetics.

[66]  Jingyuan Fu,et al.  Systems genetics: From GWAS to disease pathways. , 2014, Biochimica et biophysica acta.

[67]  Lieven Clement,et al.  Identification of novel causative genes determining the complex trait of high ethanol tolerance in yeast using pooled-segregant whole-genome sequence analysis , 2012, Genome research.

[68]  D. Botstein,et al.  The transcriptional program of sporulation in budding yeast. , 1998, Science.

[69]  Ronald W. Davis,et al.  Quantitative trait loci mapped to single-nucleotide resolution in yeast , 2005, Nature Genetics.

[70]  A. Neiman,et al.  Sporulation in the Budding Yeast Saccharomyces cerevisiae , 2011, Genetics.

[71]  K. Irie,et al.  Posttranscriptional Regulation of HO Expression by the Mkt1-Pbp1 Complex , 2004, Molecular and Cellular Biology.

[72]  Hyunjoong Kim,et al.  Functional Analysis I , 2017 .

[73]  P. Sudarsanam,et al.  Single Nucleotide Variants in Transcription Factors Associate More Tightly with Phenotype than with Gene Expression , 2014, PLoS genetics.

[74]  A. Amon,et al.  Control of Meiosis by Respiration , 2008, Current Biology.

[75]  G. Bellí,et al.  Functional analysis of yeast essential genes using a promoter‐substitution cassette and the tetracycline‐regulatable dual expression system , 1998 .

[76]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[77]  David A. Drubin,et al.  Learning a Prior on Regulatory Potential from eQTL Data , 2009, PLoS genetics.

[78]  Claude Jacq,et al.  Yeast Mitochondrial Biogenesis: A Role for the PUF RNA-Binding Protein Puf3p in mRNA Localization , 2008, PloS one.

[79]  Attila Tóth,et al.  A screen for genes required for meiosis and spore formation based on whole-genome expression , 2001, Current Biology.

[80]  Frédéric Chalmel,et al.  Execution of the meiotic noncoding RNA expression program and the onset of gametogenesis in yeast require the conserved exosome subunit Rrp6 , 2010, Proceedings of the National Academy of Sciences.

[81]  A. Mitchell Control of meiotic gene expression in Saccharomyces cerevisiae. , 1994, Microbiological reviews.

[82]  Catarina Costa,et al.  The YEASTRACT database: an upgraded information system for the analysis of gene and genomic transcription regulation in Saccharomyces cerevisiae , 2013, Nucleic Acids Res..