Genomic Analysis of QTLs and Genes Altering Natural Variation in Stochastic Noise

Quantitative genetic analysis has long been used to study how natural variation of genotype can influence an organism's phenotype. While most studies have focused on genetic determinants of phenotypic average, it is rapidly becoming understood that stochastic noise is genetically determined. However, it is not known how many traits display genetic control of stochastic noise nor how broadly these stochastic loci are distributed within the genome. Understanding these questions is critical to our understanding of quantitative traits and how they relate to the underlying causal loci, especially since stochastic noise may be directly influenced by underlying changes in the wiring of regulatory networks. We identified QTLs controlling natural variation in stochastic noise of glucosinolates, plant defense metabolites, as well as QTLs for stochastic noise of related transcripts. These loci included stochastic noise QTLs unique for either transcript or metabolite variation. Validation of these loci showed that genetic polymorphism within the regulatory network alters stochastic noise independent of effects on corresponding average levels. We examined this phenomenon more globally, using transcriptomic datasets, and found that the Arabidopsis transcriptome exhibits significant, heritable differences in stochastic noise. Further analysis allowed us to identify QTLs that control genomic stochastic noise. Some genomic QTL were in common with those altering average transcript abundance, while others were unique to stochastic noise. Using a single isogenic population, we confirmed that natural variation at ELF3 alters stochastic noise in the circadian clock and metabolism. Since polymorphisms controlling stochastic noise in genomic phenotypes exist within wild germplasm for naturally selected phenotypes, this suggests that analysis of Arabidopsis evolution should account for genetic control of stochastic variance and average phenotypes. It remains to be determined if natural genetic variation controlling stochasticity is equally distributed across the genomes of other multi-cellular eukaryotes.

[1]  Albert-László Barabási,et al.  Statistical mechanics of complex networks , 2001, ArXiv.

[2]  S. Bell,et al.  Large-Scale Reverse Genetics in Arabidopsis: Case Studies from the Chloroplast 2010 Project1[C][W][OA] , 2009, Plant Physiology.

[3]  W. Ewens.  Genetics and analysis of quantitative traits , 1999 .

[4]  M. Kreitman,et al.  Fitness costs of R-gene-mediated resistance in Arabidopsis thaliana , 2003, Nature.

[5]  S. Lindquist,et al.  HSP90-buffered genetic variation is common in Arabidopsis thaliana , 2008, Proceedings of the National Academy of Sciences.

[6]  C. Lister,et al.  Recombinant inbred lines for mapping RFLP and phenotypic markers in Arabidopsis thaliana , 1993 .

[7]  Richard M. Clark,et al.  Common Sequence Polymorphisms Shaping Genetic Diversity in Arabidopsis thaliana , 2007, Science.

[8]  Rachel E. Kerwin,et al.  A Novel 2-Oxoacid-Dependent Dioxygenase Involved in the Formation of the Goiterogenic 2-Hydroxybut-3-enyl Glucosinolate and Generalist Insect Resistance in Arabidopsis[C][W][OA] , 2008, Plant Physiology.

[9]  R. Albert,et al.  The large-scale organization of metabolic networks , 2000, Nature.

[10]  M. Kreitman,et al.  A Genome-Wide Survey of R Gene Polymorphisms in Arabidopsis[W] , 2006, The Plant Cell Online.

[11]  Thomas Mitchell-Olds,et al.  Glucosinolate and Trichome Defenses in a Natural Arabidopsis lyrata Population , 2006, Journal of Chemical Ecology.

[12]  D. Kliebenstein,et al.  Differential Levels of Insect Herbivory in the Field Associated with Genotypic Variation in Glucosinolates in Arabidopsis thaliana , 2008, Journal of Chemical Ecology.

[13]  L. Kruglyak,et al.  Genetic Dissection of Transcriptional Regulation in Budding Yeast , 2002, Science.

[14]  M. Koornneef,et al.  The phenotype of some late-flowering mutants is enhanced by a locus on chromosome 5 that is not effective in the Landsberg erecta wild-type , 1994 .

[15]  Joachim Selbig,et al.  Starch as a major integrator in the regulation of plant growth , 2009, Proceedings of the National Academy of Sciences.

[16]  M. Schmid,et al.  Genome-Wide Insertional Mutagenesis of Arabidopsis thaliana , 2003, Science.

[17]  Brook T. Moyers,et al.  Effects of Genetic Perturbation on Seasonal Life History Plasticity , 2009, Science.

[18]  Barbara Ann Halkier,et al.  A Complex Interplay of Three R2R3 MYB Transcription Factors Determines the Profile of Aliphatic Glucosinolates in Arabidopsis1[C][W][OA] , 2010, Plant Physiology.

[19]  Barbara A Halkier,et al.  Identification of a flavin-monooxygenase as the S-oxygenating enzyme in aliphatic glucosinolate biosynthesis in Arabidopsis. , 2007, The Plant journal : for cell and molecular biology.

[20]  S. Clough,et al.  Floral dip: a simplified method for Agrobacterium-mediated transformation of Arabidopsis thaliana. , 1998, The Plant journal : for cell and molecular biology.

[21]  R. Last,et al.  Arabidopsis Map-Based Cloning in the Post-Genome Era , 2002, Plant Physiology.

[22]  S. Lindquist,et al.  Hsp90 as a capacitor of phenotypic variation , 2002, Nature.

[23]  John D. Storey,et al.  Mapping the Genetic Architecture of Gene Expression in Human Liver , 2008, PLoS biology.

[24]  Jingyuan Fu,et al.  The genetics of plant metabolism , 2006, Nature Genetics.

[25]  Eric E. Schadt,et al.  The Quantitative Genetics of Phenotypic Robustness , 2010, PloS one.

[26]  R W Doerge,et al.  Genomic Survey of Gene Expression Diversity in Arabidopsis thaliana , 2006, Genetics.

[27]  C. Ballaré,et al.  Chromosomal loci important for cotyledon opening under UV-B in Arabidopsis thaliana , 2010, BMC Plant Biology.

[28]  Richard A. Lankau.  Specialist and generalist herbivores exert opposing selection on a chemical defense. , 2007, The New phytologist.

[29]  J. Lhomme,et al.  Diversity-stability relationships in community ecology: re-examination of the portfolio effect. , 2002, Theoretical population biology.

[30]  Daniel J. Kliebenstein,et al.  Linking Metabolic QTLs with Network and cis-eQTLs Controlling Biosynthetic Pathways , 2007, PLoS genetics.

[31]  M Koornneef,et al.  Development of an AFLP based linkage map of Ler, Col and Cvi Arabidopsis thaliana ecotypes and construction of a Ler/Cvi recombinant inbred line population. , 1998, The Plant journal : for cell and molecular biology.

[32]  D. Kliebenstein,et al.  Systems Biology Uncovers the Foundation of Natural Genetic Diversity1 , 2009, Plant Physiology.

[33]  R. Doerge,et al.  Permutation tests for multiple loci affecting a quantitative character. , 1996, Genetics.

[34]  Albert-László Barabási,et al.  Error and attack tolerance of complex networks , 2000, Nature.

[35]  U. Flügge,et al.  HAG2/MYB76 and HAG3/MYB29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana. , 2008, The New phytologist.

[36]  Jingyuan Fu,et al.  Regulatory network construction in Arabidopsis by using genome-wide gene expression quantitative trait loci , 2007, Proceedings of the National Academy of Sciences.

[37]  O. Loudet,et al.  Bay-0 × Shahdara recombinant inbred line population: a powerful tool for the genetic dissection of complex traits in Arabidopsis , 2002, Theoretical and Applied Genetics.

[38]  Rachel E. Kerwin,et al.  Network Quantitative Trait Loci Mapping of Circadian Clock Outputs Identifies Metabolic Pathway-to-Clock Linkages in Arabidopsis[C][W] , 2011, Plant Cell.

[39]  R. Doerge.  Multifactorial genetics: Mapping and analysis of quantitative trait loci in experimental populations , 2002, Nature Reviews Genetics.

[40]  D. Floreano,et al.  Revealing strengths and weaknesses of methods for gene network inference , 2010, Proceedings of the National Academy of Sciences.

[41]  Pedro de Atauri,et al.  Feedback control of stochastic noise in the yeast galactose utilization pathway , 2006 .

[42]  M. Purugganan,et al.  Epistatic interaction between Arabidopsis FRI and FLC flowering time genes generates a latitudinal cline in a life history trait. , 2004, Proceedings of the National Academy of Sciences of the United States of America.

[43]  D. Kliebenstein,et al.  Understanding the Evolution of Defense Metabolites in Arabidopsis thaliana Using Genome-wide Association Mapping , 2010, Genetics.

[44]  Thomas Mitchell-Olds,et al.  Evolutionary dynamics of an Arabidopsis insect resistance quantitative trait locus , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[45]  Gary D Bader,et al.  Global Mapping of the Yeast Genetic Interaction Network , 2004, Science.

[46]  Mads Kærn,et al.  A chance at survival: gene expression noise and phenotypic diversification strategies , 2009, Molecular microbiology.

[47]  J. Bergelson,et al.  SAR INCREASES FITNESS OF ARABIDOPSIS THALIANA IN THE PRESENCE OF NATURAL BACTERIAL PATHOGENS , 2007, Evolution; international journal of organic evolution.

[48]  J. Gershenzon,et al.  Comparative quantitative trait loci mapping of aliphatic, indolic and benzylic glucosinolate production in Arabidopsis thaliana leaves and seeds. , 2001, Genetics.

[49]  Richard A. Lankau,et al.  Mutual Feedbacks Maintain Both Genetic and Species Diversity in a Plant Community , 2007, Science.

[50]  P. Raven,et al.  BUTTERFLIES AND PLANTS: A STUDY IN COEVOLUTION , 1964 .

[51]  Barbara Ann Halkier,et al.  Biology and biochemistry of glucosinolates. , 2006, Annual review of plant biology.

[52]  Justin O. Borevitz,et al.  Global Analysis of Genetic, Epigenetic and Transcriptional Polymorphisms in Arabidopsis thaliana Using Whole Genome Tiling Arrays , 2008, PLoS genetics.

[53]  R. Amasino,et al.  FLOWERING LOCUS C Encodes a Novel MADS Domain Protein That Acts as a Repressor of Flowering , 1999, Plant Cell.

[54]  Albert-László Barabási,et al.  Genetic Dissection of Transcriptional Regulation in Budding Yeast , 2002 .

[55]  Hiroaki Kitano,et al.  Biological robustness , 2008, Nature Reviews Genetics.

[56]  P. Schulze-Lefert,et al.  A Glucosinolate Metabolism Pathway in Living Plant Cells Mediates Broad-Spectrum Antifungal Defense , 2009, Science.

[57]  F. V. van Eeuwijk,et al.  Linkage Disequilibrium Mapping of Yield and Yield Stability in Modern Spring Barley Cultivars , 2004, Genetics.

[58]  Ben Hui Liu,et al.  Statistical Genomics: Linkage, Mapping, and QTL Analysis , 1997 .

[59]  Michael F. Covington,et al.  Global transcriptome analysis reveals circadian regulation of key pathways in plant growth and development , 2008, Genome Biology.

[60]  Yunde Zhao,et al.  Auxin biosynthesis by the YUCCA flavin monooxygenases controls the formation of floral organs and vascular tissues in Arabidopsis. , 2006, Genes & development.

[61]  D. Kliebenstein,et al.  Regulatory networks of glucosinolates shape Arabidopsis thaliana fitness. , 2010, Current opinion in plant biology.

[62]  R. Doerge,et al.  Global eQTL Mapping Reveals the Complex Genetic Architecture of Transcript-Level Variation in Arabidopsis , 2007, Genetics.

[63]  M. L. Simpson,et al.  Gene network shaping of inherent noise spectra , 2006, Nature.

[64]  Pedro de Atauri,et al.  Dual feedback loops in the GAL regulon suppress cellular heterogeneity in yeast , 2006, Nature Genetics.

[65]  S. Leibler,et al.  Robustness in simple biochemical networks , 1997, Nature.

[66]  M. Hirai,et al.  Omics-based identification of Arabidopsis Myb transcription factors regulating aliphatic glucosinolate biosynthesis , 2007, Proceedings of the National Academy of Sciences.

[67]  G. Jander,et al.  Arabidopsis myrosinases TGG1 and TGG2 have redundant function in glucosinolate breakdown and insect defense. , 2006, The Plant journal : for cell and molecular biology.

[68]  A. Oudenaarden,et al.  Nature, Nurture, or Chance: Stochastic Gene Expression and Its Consequences , 2008, Cell.

[69]  J. François,et al.  Cell-to-Cell Stochastic Variation in Gene Expression Is a Complex Genetic Trait , 2008, PLoS genetics.

[70]  S. Shannon,et al.  Early-Flowering Mutants of Arabidopsis thaliana , 1992 .

[71]  R. Doerge,et al.  Empirical threshold values for quantitative trait mapping. , 1994, Genetics.

[72]  T. Korves,et al.  A Novel Cost of R Gene Resistance in the Presence of Disease , 2004, The American Naturalist.

[73]  R W Doerge,et al.  High-density haplotyping with microarray-based expression and single feature polymorphism markers in Arabidopsis. , 2006, Genome research.

[74]  B. Halkier,et al.  Piecing together the transport pathway of aliphatic glucosinolates , 2008, Phytochemistry Reviews.

[75]  Adam M. Wentzell,et al.  Genotype, Age, Tissue, and Environment Regulate the Structural Outcome of Glucosinolate Activation1[W][OA] , 2008, Plant Physiology.

[76]  A. Shelton.  Within-Plant Variation In Glucosinolate Concentrations of Raphanus sativus Across Multiple Scales , 2005, Journal of Chemical Ecology.

[77]  Daniel J. Kliebenstein,et al.  Identification of QTLs controlling gene expression networks defined a priori , 2006, BMC Bioinformatics.

[78]  B. Berger,et al.  The R 2 R 3-MYB transcription factor HAG 1 / MYB 28 is a regulator of methionine-derived glucosinolate biosynthesis in Arabidopsis thaliana , 2007 .

[79]  D. Kliebenstein.  Quantitative genomics: analyzing intraspecific variation using global gene expression polymorphisms or eQTLs. , 2009, Annual review of plant biology.

[80]  P. Keightley.  Models of quantitative variation of flux in metabolic pathways. , 1989, Genetics.

[81]  B. Berger,et al.  The R2R3-MYB transcription factor HAG1/MYB28 is a regulator of methionine-derived glucosinolate biosynthesis in Arabidopsis thaliana. , 2007, The Plant journal : for cell and molecular biology.

[82]  Rachel B. Brem,et al.  The landscape of genetic complexity across 5,700 gene expression traits in yeast. , 2005, Proceedings of the National Academy of Sciences of the United States of America.

[83]  M. Kreitman,et al.  Low Levels of Polymorphism in Genes That Control the Activation of Defense Response in Arabidopsis thaliana , 2008, Genetics.

[84]  D. Kliebenstein,et al.  A Systems Biology Approach Identifies a R2R3 MYB Gene Subfamily with Distinct and Overlapping Functions in Regulation of Aliphatic Glucosinolates , 2007, PloS one.

[85]  S. Kay,et al.  Quantitative Analysis of Drosophila period Gene Transcription in Living Animals , 1997, Journal of biological rhythms.

[86]  D. Charlesworth,et al.  Rates and patterns of molecular evolution in inbred and outbred Arabidopsis. , 2002, Molecular biology and evolution.

[87]  Terence P. Speed,et al.  A comparison of normalization methods for high density oligonucleotide array data based on variance and bias , 2003, Bioinform..

[88]  H. Grüneberg,et al.  Introduction to quantitative genetics , 1960 .

[89]  M. Purugganan,et al.  Genetics of microenvironmental canalization in Arabidopsis thaliana , 2007, Proceedings of the National Academy of Sciences.

[90]  Thomas Mitchell-Olds,et al.  Epistasis and balanced polymorphism influencing complex trait variation , 2005, Nature.

[91]  H. Kitano.  Towards a theory of biological robustness , 2007, Molecular systems biology.

[92]  U. Flügge,et al.  HAG 2 / MYB 76 and HAG 3 / MYB 29 exert a specific and coordinated control on the regulation of aliphatic glucosinolate biosynthesis in Arabidopsis thaliana , 2008 .

[93]  J. Maloof,et al.  Network Analysis Identifies ELF3 as a QTL for the Shade Avoidance Response in Arabidopsis , 2010, PLoS genetics.

[94]  Jeffrey W. Smith,et al.  Stochastic Gene Expression in a Single Cell , 2022 .

[95]  Klaus Aktories,et al.  Noise Can Induce Bimodality in Positive Transcriptional Feedback Loops Without Bistability , 2010 .

[96]  Z B Zeng,et al.  Estimating the genetic architecture of quantitative traits. , 1999, Genetical research.

[97]  S. Harmer,et al.  The circadian system in higher plants. , 2009, Annual review of plant biology.

[98]  Rachel B. Brem,et al.  Trans-acting regulatory variation in Saccharomyces cerevisiae and the role of transcription factors , 2003, Nature Genetics.

[99]  D. Kliebenstein.  A Role for Gene Duplication and Natural Variation of Gene Expression in the Evolution of Metabolism , 2008, PloS one.

[100]  D. Kliebenstein,et al.  The Gene Controlling the Quantitative Trait Locus EPITHIOSPECIFIER MODIFIER1 Alters Glucosinolate Hydrolysis and Insect Resistance in Arabidopsis[W] , 2006, The Plant Cell Online.

[101]  A. Shelton.  Variation in chemical defences of plants may improve the effectiveness of defence , 2004 .

[102]  J. Keurentjes.  Genetical metabolomics: closing in on phenotypes. , 2009, Current opinion in plant biology.

[103]  Thomas Mitchell-Olds,et al.  Positive selection driving diversification in plant secondary metabolism. , 2006, Proceedings of the National Academy of Sciences of the United States of America.

[104]  Adam M. Wentzell,et al.  Genetic Networks Controlling Structural Outcome of Glucosinolate Activation across Development , 2008, PLoS genetics.

[105]  J. Raser,et al.  Control of Stochasticity in Eukaryotic Gene Expression , 2004, Science.

[106]  D. Kliebenstein.  Advancing Genetic Theory and Application by Metabolic Quantitative Trait Loci Analysis , 2009, The Plant Cell Online.

[107]  Montgomery Slatkin,et al.  Epigenetic Inheritance and the Missing Heritability Problem , 2009, Genetics.

[108]  G. Church,et al.  Modular epistasis in yeast metabolism , 2005, Nature Genetics.

[109]  Susan Lindquist,et al.  Hsp90 and Environmental Stress Transform the Adaptive Value of Natural Genetic Variation , 2010, Science.

[110]  Richard A. Lankau,et al.  Community Complexity Drives Patterns of Natural Selection on a Chemical Defense of Brassica nigra , 2007, The American Naturalist.

[111]  C. Waddington.  Canalization of Development and the Inheritance of Acquired Characters , 1942, Nature.

[112]  Richard A. Lankau,et al.  Competition, herbivory and genetics interact to determine the accumulation and fitness consequences of a defence metabolite , 2009 .

[113]  Ben Lehner,et al.  Genes Confer Similar Robustness to Environmental, Stochastic, and Genetic Perturbations in Yeast , 2010, PloS one.

[114]  Stig W Omholt,et al.  Statistical Epistasis Is a Generic Feature of Gene Regulatory Networks , 2007, Genetics.

[115]  I. Schmalhausen.  Factors of evolution : the theory of stabilizing selection , 1946 .

[116]  O. Kuipers,et al.  Bistability, epigenetics, and bet-hedging in bacteria. , 2008, Annual review of microbiology.

[117]  I. Sønderby,et al.  Biosynthesis of glucosinolates--gene discovery and beyond. , 2010, Trends in plant science.

[118]  Wenfeng Qian,et al.  Positive selection for elevated gene expression noise in yeast , 2009, Molecular systems biology.

[119]  H. Vogel,et al.  The Gene Controlling the Indole Glucosinolate Modifier1 Quantitative Trait Locus Alters Indole Glucosinolate Structures and Aphid Resistance in Arabidopsis[W] , 2009, The Plant Cell Online.

[120]  Thomas Mitchell-Olds,et al.  Comparative analysis of quantitative trait loci controlling glucosinolates, myrosinase and insect resistance in Arabidopsis thaliana. , 2002, Genetics.

[121]  John D. Storey,et al.  Genetic interactions between polymorphisms that affect gene expression in yeast , 2005, Nature.

[122]  Bjarne Gram Hansen,et al.  Subclade of Flavin-Monooxygenases Involved in Aliphatic Glucosinolate Biosynthesis1[W] , 2008, Plant Physiology.

[123]  Dimitra L. Milioni,et al.  Genomic organization of hsp90 gene family in Arabidopsis , 1997, Plant Molecular Biology.

[124]  Massimo Marchiori,et al.  Error and attacktolerance of complex network s , 2004 .

[125]  Mark G. M. Aarts,et al.  The Impact of the Absence of Aliphatic Glucosinolates on Insect Herbivory in Arabidopsis , 2008, PloS one.

[126]  M. Reichelt,et al.  Gene Duplication in the Diversification of Secondary Metabolism: Tandem 2-Oxoglutarate–Dependent Dioxygenases Control Glucosinolate Biosynthesis in Arabidopsis , 2001, Plant Cell.

[127]  Frederick M. Ausubel,et al.  Glucosinolate Metabolites Required for an Arabidopsis Innate Immune Response , 2009, Science.

[128]  Muhammad Ali Amer,et al.  Genome-wide association study of 107 phenotypes in a common set of Arabidopsis thaliana inbred lines , 2010, Nature.