1
|
Huang Y, Mao Z, Zhang Y, Zhao J, Luan X, Wu K, Yun L, Yu J, Shi Z, Liao X, Ma H. Omics data analysis reveals the system-level constraint on cellular amino acid composition. Synth Syst Biotechnol 2024; 9:304-311. [PMID: 38510205 PMCID: PMC10951587 DOI: 10.1016/j.synbio.2024.03.001] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/12/2023] [Revised: 03/01/2024] [Accepted: 03/01/2024] [Indexed: 03/22/2024] Open
Abstract
Proteins play a pivotal role in coordinating the functions of organisms, essentially governing their traits, as the dynamic arrangement of diverse amino acids leads to a multitude of folded configurations within peptide chains. Despite dynamic changes in amino acid composition of an individual protein (referred to as AAP) and great variance in protein expression levels under different conditions, our study, utilizing transcriptomics data from four model organisms uncovers surprising stability in the overall amino acid composition of the total cellular proteins (referred to as AACell). Although this value may vary between different species, we observed no significant differences among distinct strains of the same species. This indicates that organisms enforce system-level constraints to maintain a consistent AACell, even amid fluctuations in AAP and protein expression. Further exploration of this phenomenon promises insights into the intricate mechanisms orchestrating cellular protein expression and adaptation to varying environmental challenges.
Collapse
Affiliation(s)
- Yuanyuan Huang
- College of Biotechnology, Tianjin University of Science and Technology, Tianjin, 300457, China
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Zhitao Mao
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Yue Zhang
- College of Biotechnology, Tianjin University of Science and Technology, Tianjin, 300457, China
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Jianxiao Zhao
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
- Frontier Science Center for Synthetic Biology and Key Laboratory of Systems Bioengineering (Ministry of Education), Tianjin University, Tianjin, 300072, China
| | - Xiaodi Luan
- College of Biotechnology, Tianjin University of Science and Technology, Tianjin, 300457, China
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Ke Wu
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Lili Yun
- Tianjin Medical Laboratory, BGI-Tianjin, BGI-Shenzhen, Tianjin, 300308, China
| | - Jing Yu
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Zhenkun Shi
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Xiaoping Liao
- Haihe Laboratory of Synthetic Biology, Tianjin, 300308, China
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| | - Hongwu Ma
- Biodesign Center, Key Laboratory of Engineering Biology for Low-carbon Manufacturing, Tianjin Institute of Industrial Biotechnology, Chinese Academy of Sciences, Tianjin, 300308, China
- National Center of Technology Innovation for Synthetic Biology, Tianjin, 300308, China
| |
Collapse
|
2
|
Jones EF, Haldar A, Oza VH, Lasseigne BN. Quantifying transcriptome diversity: a review. Brief Funct Genomics 2024; 23:83-94. [PMID: 37225889 DOI: 10.1093/bfgp/elad019] [Citation(s) in RCA: 3] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/18/2023] [Revised: 04/14/2023] [Accepted: 05/05/2023] [Indexed: 05/26/2023] Open
Abstract
Following the central dogma of molecular biology, gene expression heterogeneity can aid in predicting and explaining the wide variety of protein products, functions and, ultimately, heterogeneity in phenotypes. There is currently overlapping terminology used to describe the types of diversity in gene expression profiles, and overlooking these nuances can misrepresent important biological information. Here, we describe transcriptome diversity as a measure of the heterogeneity in (1) the expression of all genes within a sample or a single gene across samples in a population (gene-level diversity) or (2) the isoform-specific expression of a given gene (isoform-level diversity). We first overview modulators and quantification of transcriptome diversity at the gene level. Then, we discuss the role alternative splicing plays in driving transcript isoform-level diversity and how it can be quantified. Additionally, we overview computational resources for calculating gene-level and isoform-level diversity for high-throughput sequencing data. Finally, we discuss future applications of transcriptome diversity. This review provides a comprehensive overview of how gene expression diversity arises, and how measuring it determines a more complete picture of heterogeneity across proteins, cells, tissues, organisms and species.
Collapse
Affiliation(s)
- Emma F Jones
- The Department of Cell, Developmental and Integrative Biology, Heersink School of Medicine, The University of Alabama at Birmingham, Birmingham, AL, USA
| | - Anisha Haldar
- The Department of Cell, Developmental and Integrative Biology, Heersink School of Medicine, The University of Alabama at Birmingham, Birmingham, AL, USA
| | - Vishal H Oza
- The Department of Cell, Developmental and Integrative Biology, Heersink School of Medicine, The University of Alabama at Birmingham, Birmingham, AL, USA
| | - Brittany N Lasseigne
- The Department of Cell, Developmental and Integrative Biology, Heersink School of Medicine, The University of Alabama at Birmingham, Birmingham, AL, USA
| |
Collapse
|
3
|
van Eijnatten AL, Sterken MG, Kammenga JE, Nijveen H, Snoek BL. The effect of developmental variation on expression QTLs in a multi parental Caenorhabditis elegans population. G3 (BETHESDA, MD.) 2024; 14:jkad273. [PMID: 38015660 PMCID: PMC10849341 DOI: 10.1093/g3journal/jkad273] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 09/21/2023] [Revised: 09/21/2023] [Accepted: 10/27/2023] [Indexed: 11/30/2023]
Abstract
Regulation of gene expression plays a crucial role in developmental processes and adaptation to changing environments. expression quantitative trait locus (eQTL) mapping is a technique used to study the genetic regulation of gene expression using the transcriptomes of recombinant inbred lines (RILs). Typically, the age of the inbred lines at the time of RNA sampling is carefully controlled. This is necessary because the developmental process causes changes in gene expression, complicating the interpretation of eQTL mapping experiments. However, due to genetics and variation in ambient micro-environments, organisms can differ in their "developmental age," even if they are of the same chronological age. As a result, eQTL patterns are affected by developmental variation in gene expression. The model organism Caenorhabditis elegans is particularly suited for studying the effect of developmental variation on eQTL mapping patterns. In a span of days, C. elegans transitions from embryo through 4 larval stages to adult while undergoing massive changes to its transcriptome. Here, we use C. elegans to investigate the effect of developmental age variation on eQTL patterns and present a normalization procedure. We used dynamical eQTL mapping, which includes the developmental age as a cofactor, to separate the variation in development from genotypic variation and explain variation in gene expression levels. We compare classical single marker eQTL mapping and dynamical eQTL mapping using RNA-seq data of ∼200 multi-parental RILs of C. elegans. The results show that (1) many eQTLs are caused by developmental variation, (2) most trans-bands are developmental QTLs, and (3) dynamical eQTL mapping detects additional eQTLs not found with classical eQTL mapping. We recommend that correction for variation in developmental age should be strongly considered in eQTL mapping studies given the large impact of processes like development on the transcriptome.
Collapse
Affiliation(s)
- Abraham L van Eijnatten
- Theoretical Biology and Bioinformatics, Utrecht University, Padualaan 8,3584 CH Utrecht, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, 6708 PB Wageningen, The Netherlands
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, 6708 PB Wageningen, The Netherlands
| | - Harm Nijveen
- Laboratory of Bioinformatics, Wageningen University, Droevendaalsesteeg 1, Radix West, Building 107, 6708 PB Wageningen, The Netherlands
| | - Basten L Snoek
- Theoretical Biology and Bioinformatics, Utrecht University, Padualaan 8,3584 CH Utrecht, The Netherlands
| |
Collapse
|
4
|
Gao AW, Alam GE, Zhu Y, Li W, Katsyuba E, Sulc J, Li TY, Li X, Overmyer KA, Lalou A, Mouchiroud L, Sleiman MB, Cornaglia M, Morel JD, Houtkooper RH, Coon JJ, Auwerx J. High-content phenotypic analysis of a C. elegans recombinant inbred population identifies genetic and molecular regulators of lifespan. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2024:2024.01.15.575638. [PMID: 38293129 PMCID: PMC10827074 DOI: 10.1101/2024.01.15.575638] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/01/2024]
Abstract
Lifespan is influenced by complex interactions between genetic and environmental factors. Studying those factors in model organisms of a single genetic background limits their translational value for humans. Here, we mapped lifespan determinants in 85 genetically diverse C. elegans recombinant intercross advanced inbred lines (RIAILs). We assessed molecular profiles - transcriptome, proteome, and lipidome - and life-history traits, including lifespan, development, growth dynamics, and reproduction. RIAILs exhibited large variations in lifespan, which positively correlated with developmental time. Among the top candidates obtained from multi-omics data integration and QTL mapping, we validated known and novel longevity modulators, including rict-1, gfm-1 and mltn-1. We translated their relevance to humans using UK Biobank data and showed that variants in RICTOR and GFM1 are associated with an elevated risk of age-related heart disease, dementia, diabetes, kidney, and liver diseases. We organized our dataset as a resource (https://lisp-lms.shinyapps.io/RIAILs/) that allows interactive explorations for new longevity targets.
Collapse
Affiliation(s)
- Arwen W. Gao
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Laboratory Genetic Metabolic Diseases, Amsterdam Gastroenterology, Endocrinology, and Metabolism, Amsterdam UMC, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
| | - Gaby El Alam
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Yunyun Zhu
- Department of Biomolecular Chemistry, University of Wisconsin, Madison, WI 53506, USA
| | - Weisha Li
- Laboratory Genetic Metabolic Diseases, Amsterdam Gastroenterology, Endocrinology, and Metabolism, Amsterdam UMC, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
| | - Elena Katsyuba
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Nagi Bioscience SA, EPFL Innovation Park, CH-1025 Saint-Sulpice, Switzerland
| | - Jonathan Sulc
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Terytty Y. Li
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Present address: State Key Laboratory of Genetic Engineering, Shanghai Key Laboratory of Metabolic Remodeling and Health, Laboratory of Longevity and Metabolic Adaptations, Institute of Metabolism and Integrative Biology, Fudan University, Shanghai, China
| | - Xiaoxu Li
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Katherine A. Overmyer
- Department of Biomolecular Chemistry, University of Wisconsin, Madison, WI 53506, USA
- National Center for Quantitative Biology of Complex Systems, Madison, WI 53706, USA
- Morgridge Institute for Research, Madison, WI 53515, USA
| | - Amelia Lalou
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Laurent Mouchiroud
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
- Nagi Bioscience SA, EPFL Innovation Park, CH-1025 Saint-Sulpice, Switzerland
| | - Maroun Bou Sleiman
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Matteo Cornaglia
- Nagi Bioscience SA, EPFL Innovation Park, CH-1025 Saint-Sulpice, Switzerland
| | - Jean-David Morel
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| | - Riekelt H. Houtkooper
- Laboratory Genetic Metabolic Diseases, Amsterdam Gastroenterology, Endocrinology, and Metabolism, Amsterdam UMC, University of Amsterdam, Meibergdreef 9, 1105 AZ Amsterdam, The Netherlands
| | - Joshua J. Coon
- Department of Biomolecular Chemistry, University of Wisconsin, Madison, WI 53506, USA
- National Center for Quantitative Biology of Complex Systems, Madison, WI 53706, USA
- Morgridge Institute for Research, Madison, WI 53515, USA
- Department of Chemistry, University of Wisconsin, Madison, WI 53506, USA
| | - Johan Auwerx
- Laboratory of Integrative Systems Physiology, Interfaculty Institute of Bioengineering, École Polytechnique Fédérale de Lausanne, CH-1015 Lausanne, Switzerland
| |
Collapse
|
5
|
Bell AD, Chou HT, Valencia F, Paaby AB. Beyond the reference: gene expression variation and transcriptional response to RNA interference in Caenorhabditis elegans. G3 (BETHESDA, MD.) 2023; 13:jkad112. [PMID: 37221008 PMCID: PMC10411595 DOI: 10.1093/g3journal/jkad112] [Citation(s) in RCA: 1] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Subscribe] [Scholar Register] [Received: 03/27/2023] [Revised: 05/11/2023] [Accepted: 05/15/2023] [Indexed: 05/25/2023]
Abstract
Though natural systems harbor genetic and phenotypic variation, research in model organisms is often restricted to a reference strain. Focusing on a reference strain yields a great depth of knowledge but potentially at the cost of breadth of understanding. Furthermore, tools developed in the reference context may introduce bias when applied to other strains, posing challenges to defining the scope of variation within model systems. Here, we evaluate how genetic differences among 5 wild Caenorhabditis elegans strains affect gene expression and its quantification, in general and after induction of the RNA interference (RNAi) response. Across strains, 34% of genes were differentially expressed in the control condition, including 411 genes that were not expressed at all in at least 1 strain; 49 of these were unexpressed in reference strain N2. Reference genome mapping bias caused limited concern: despite hyperdiverse hotspots throughout the genome, 92% of variably expressed genes were robust to mapping issues. The transcriptional response to RNAi was highly strain- and target-gene-specific and did not correlate with RNAi efficiency, as the 2 RNAi-insensitive strains showed more differentially expressed genes following RNAi treatment than the RNAi-sensitive reference strain. We conclude that gene expression, generally and in response to RNAi, differs across C. elegans strains such that the choice of strain may meaningfully influence scientific inferences. Finally, we introduce a resource for querying gene expression variation in this dataset at https://wildworm.biosci.gatech.edu/rnai/.
Collapse
Affiliation(s)
- Avery Davis Bell
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr NW, EBB Building, Atlanta, GA 30332, USA
| | - Han Ting Chou
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr NW, EBB Building, Atlanta, GA 30332, USA
| | - Francisco Valencia
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr NW, EBB Building, Atlanta, GA 30332, USA
| | - Annalise B Paaby
- School of Biological Sciences, Georgia Institute of Technology, 950 Atlantic Dr NW, EBB Building, Atlanta, GA 30332, USA
| |
Collapse
|
6
|
Bell AD, Chou HT, Paaby AB. Beyond the reference: gene expression variation and transcriptional response to RNAi in C. elegans. BIORXIV : THE PREPRINT SERVER FOR BIOLOGY 2023:2023.03.24.533964. [PMID: 36993640 PMCID: PMC10055391 DOI: 10.1101/2023.03.24.533964] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Grants] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 06/03/2023]
Abstract
A universal feature of living systems is that natural variation in genotype underpins variation in phenotype. Yet, research in model organisms is often constrained to a single genetic background, the reference strain. Further, genomic studies that do evaluate wild strains typically rely on the reference strain genome for read alignment, leading to the possibility of biased inferences based on incomplete or inaccurate mapping; the extent of reference bias can be difficult to quantify. As an intermediary between genome and organismal traits, gene expression is well positioned to describe natural variability across genotypes generally and in the context of environmental responses, which can represent complex adaptive phenotypes. C. elegans sits at the forefront of investigation into small-RNA gene regulatory mechanisms, or RNA interference (RNAi), and wild strains exhibit natural variation in RNAi competency following environmental triggers. Here, we examine how genetic differences among five wild strains affect the C. elegans transcriptome in general and after inducing RNAi responses to two germline target genes. Approximately 34% of genes were differentially expressed across strains; 411 genes were not expressed at all in at least one strain despite robust expression in others, including 49 genes not expressed in reference strain N2. Despite the presence of hyper-diverse hotspots throughout the C. elegans genome, reference mapping bias was of limited concern: over 92% of variably expressed genes were robust to mapping issues. Overall, the transcriptional response to RNAi was strongly strain-specific and highly specific to the target gene, and the laboratory strain N2 was not representative of the other strains. Moreover, the transcriptional response to RNAi was not correlated with RNAi phenotypic penetrance; the two germline RNAi incompetent strains exhibited substantial differential gene expression following RNAi treatment, indicating an RNAi response despite failure to reduce expression of the target gene. We conclude that gene expression, both generally and in response to RNAi, differs across C. elegans strains such that choice of strain may meaningfully influence scientific conclusions. To provide a public, easily accessible resource for querying gene expression variation in this dataset, we introduce an interactive website at https://wildworm.biosci.gatech.edu/rnai/ .
Collapse
Affiliation(s)
- Avery Davis Bell
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA
| | - Han Ting Chou
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA
| | - Annalise B. Paaby
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA
| |
Collapse
|
7
|
Sterken MG, Nijveen H, van Zanten M, Jiménez-Gómez JM, Geshnizjani N, Willems LAJ, Rienstra J, Hilhorst HWM, Ligterink W, Snoek BL. Plasticity of maternal environment-dependent expression-QTLs of tomato seeds. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2023; 136:28. [PMID: 36810666 PMCID: PMC9944408 DOI: 10.1007/s00122-023-04322-0] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Figures] [Subscribe] [Scholar Register] [Received: 12/27/2021] [Accepted: 10/27/2022] [Indexed: 06/18/2023]
Abstract
Seeds are essential for plant reproduction, survival, and dispersal. Germination ability and successful establishment of young seedlings strongly depend on seed quality and on environmental factors such as nutrient availability. In tomato (Solanum lycopersicum) and many other species, seed quality and seedling establishment characteristics are determined by genetic variation, as well as the maternal environment in which the seeds develop and mature. The genetic contribution to variation in seed and seedling quality traits and environmental responsiveness can be estimated at transcriptome level in the dry seed by mapping genomic loci that affect gene expression (expression QTLs) in contrasting maternal environments. In this study, we applied RNA-sequencing to construct a linkage map and measure gene expression of seeds of a tomato recombinant inbred line (RIL) population derived from a cross between S. lycopersicum (cv. Moneymaker) and S. pimpinellifolium (G1.1554). The seeds matured on plants cultivated under different nutritional environments, i.e., on high phosphorus or low nitrogen. The obtained single-nucleotide polymorphisms (SNPs) were subsequently used to construct a genetic map. We show how the genetic landscape of plasticity in gene regulation in dry seeds is affected by the maternal nutrient environment. The combined information on natural genetic variation mediating (variation in) responsiveness to the environment may contribute to knowledge-based breeding programs aiming to develop crop cultivars that are resilient to stressful environments.
Collapse
Affiliation(s)
- Mark G. Sterken
- Laboratory of Nematology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Harm Nijveen
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
- Laboratory of Bioinformatics, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Martijn van Zanten
- Plant Stress Resilience, Institute of Environmental Biology, Utrecht University, 3584 CH Utrecht, The Netherlands
| | - Jose M. Jiménez-Gómez
- Department of Plant Breeding and Genetics, Max Planck Institute for Plant Breeding Research, Cologne, Germany
- Institut Jean-Pierre Bourgin, INRAE, AgroParisTech, Université Paris-Saclay, 78000 Versailles, France
| | - Nafiseh Geshnizjani
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Leo A. J. Willems
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Juriaan Rienstra
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Henk W. M. Hilhorst
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Wilco Ligterink
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, 6708 PB Wageningen, The Netherlands
| | - Basten L. Snoek
- Laboratory of Nematology, Wageningen University, 6708 PB Wageningen, The Netherlands
- Theoretical Biology and Bioinformatics, Institute of Biodynamics and Biocomplexity, Utrecht University, 3584 CH Utrecht, The Netherlands
| |
Collapse
|
8
|
The impact of species-wide gene expression variation on Caenorhabditis elegans complex traits. Nat Commun 2022; 13:3462. [PMID: 35710766 PMCID: PMC9203580 DOI: 10.1038/s41467-022-31208-4] [Citation(s) in RCA: 13] [Impact Index Per Article: 6.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/09/2022] [Accepted: 06/08/2022] [Indexed: 12/15/2022] Open
Abstract
Phenotypic variation in organism-level traits has been studied in Caenorhabditis elegans wild strains, but the impacts of differences in gene expression and the underlying regulatory mechanisms are largely unknown. Here, we use natural variation in gene expression to connect genetic variants to differences in organismal-level traits, including drug and toxicant responses. We perform transcriptomic analyses on 207 genetically distinct C. elegans wild strains to study natural regulatory variation of gene expression. Using this massive dataset, we perform genome-wide association mappings to investigate the genetic basis underlying gene expression variation and reveal complex genetic architectures. We find a large collection of hotspots enriched for expression quantitative trait loci across the genome. We further use mediation analysis to understand how gene expression variation could underlie organism-level phenotypic variation for a variety of complex traits. These results reveal the natural diversity in gene expression and possible regulatory mechanisms in this keystone model organism, highlighting the promise of using gene expression variation to understand how phenotypic diversity is generated.
Collapse
|
9
|
The genetic architecture underlying body-size traits plasticity over different temperatures and developmental stages in Caenorhabditis elegans. Heredity (Edinb) 2022; 128:313-324. [PMID: 35383317 PMCID: PMC9076863 DOI: 10.1038/s41437-022-00528-y] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/24/2021] [Revised: 03/15/2022] [Accepted: 03/17/2022] [Indexed: 01/25/2023] Open
Abstract
Most ectotherms obey the temperature-size rule, meaning they grow larger in a colder environment. This raises the question of how the interplay between genes and temperature affects the body size of ectotherms. Despite the growing body of literature on the physiological life-history and molecular genetic mechanism underlying the temperature-size rule, the overall genetic architecture orchestrating this complex phenotype is not yet fully understood. One approach to identify genetic regulators of complex phenotypes is quantitative trait locus (QTL) mapping. Here, we explore the genetic architecture of body-size phenotypes, and plasticity of body-size phenotypes at different temperatures using Caenorhabditis elegans as a model ectotherm. We used 40 recombinant inbred lines (RILs) derived from N2 and CB4856, which were reared at four different temperatures (16, 20, 24, and 26 °C) and measured at two developmental stages (L4 and adult). The animals were measured for body length, width at vulva, body volume, length/width ratio, and seven other body-size traits. The genetically diverse RILs varied in their body-size phenotypes with heritabilities ranging from 0.0 to 0.99. We detected 18 QTL underlying the body-size traits across all treatment combinations, with the majority clustering on Chromosome X. We hypothesize that the Chromosome X QTL could result from a known pleiotropic regulator-npr-1-known to affect the body size of C. elegans through behavioral changes. We also found five plasticity QTL of body-size traits where three colocalized with body-size QTL. In conclusion, our findings shed more light on multiple loci affecting body-size plasticity and the possibility of co-regulation of traits and traits plasticity by the same loci under different environments.
Collapse
|
10
|
van Sluijs L, Bosman KJ, Pankok F, Blokhina T, Wilten JIHA, te Molder DM, Riksen JAG, Snoek BL, Pijlman GP, Kammenga JE, Sterken MG. Balancing Selection of the Intracellular Pathogen Response in Natural Caenorhabditis elegans Populations. Front Cell Infect Microbiol 2022; 11:758331. [PMID: 35174100 PMCID: PMC8841876 DOI: 10.3389/fcimb.2021.758331] [Citation(s) in RCA: 6] [Impact Index Per Article: 3.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2021] [Accepted: 12/21/2021] [Indexed: 12/17/2022] Open
Abstract
Genetic variation in host populations may lead to differential viral susceptibilities. Here, we investigate the role of natural genetic variation in the Intracellular Pathogen Response (IPR), an important antiviral pathway in the model organism Caenorhabditis elegans against Orsay virus (OrV). The IPR involves transcriptional activity of 80 genes including the pals-genes. We examine the genetic variation in the pals-family for traces of selection and explore the molecular and phenotypic effects of having distinct pals-gene alleles. Genetic analysis of 330 global C. elegans strains reveals that genetic diversity within the IPR-related pals-genes can be categorized in a few haplotypes worldwide. Importantly, two key IPR regulators, pals-22 and pals-25, are in a genomic region carrying signatures of balancing selection, suggesting that different evolutionary strategies exist in IPR regulation. We infected eleven C. elegans strains that represent three distinct pals-22 pals-25 haplotypes with Orsay virus to determine their susceptibility. For two of these strains, N2 and CB4856, the transcriptional response to infection was also measured. The results indicate that pals-22 pals-25 haplotype shapes the defense against OrV and host genetic variation can result in constitutive activation of IPR genes. Our work presents evidence for balancing genetic selection of immunity genes in C. elegans and provides a novel perspective on the functional diversity that can develop within a main antiviral response in natural host populations.
Collapse
Affiliation(s)
- Lisa van Sluijs
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
- Laboratory of Virology, Wageningen University and Research, Wageningen, Netherlands
| | - Kobus J. Bosman
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Frederik Pankok
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Tatiana Blokhina
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Jop I. H. A. Wilten
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Dennie M. te Molder
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Joost A. G. Riksen
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Basten L. Snoek
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Gorben P. Pijlman
- Laboratory of Virology, Wageningen University and Research, Wageningen, Netherlands
| | - Jan E. Kammenga
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
| | - Mark G. Sterken
- Laboratory of Nematology, Wageningen University and Research, Wageningen, Netherlands
- Laboratory of Virology, Wageningen University and Research, Wageningen, Netherlands
- *Correspondence: Mark G. Sterken,
| |
Collapse
|
11
|
Andersen EC, Rockman MV. Natural genetic variation as a tool for discovery in Caenorhabditis nematodes. Genetics 2022; 220:iyab156. [PMID: 35134197 PMCID: PMC8733454 DOI: 10.1093/genetics/iyab156] [Citation(s) in RCA: 20] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/20/2021] [Accepted: 09/11/2021] [Indexed: 11/12/2022] Open
Abstract
Over the last 20 years, studies of Caenorhabditis elegans natural diversity have demonstrated the power of quantitative genetic approaches to reveal the evolutionary, ecological, and genetic factors that shape traits. These studies complement the use of the laboratory-adapted strain N2 and enable additional discoveries not possible using only one genetic background. In this chapter, we describe how to perform quantitative genetic studies in Caenorhabditis, with an emphasis on C. elegans. These approaches use correlations between genotype and phenotype across populations of genetically diverse individuals to discover the genetic causes of phenotypic variation. We present methods that use linkage, near-isogenic lines, association, and bulk-segregant mapping, and we describe the advantages and disadvantages of each approach. The power of C. elegans quantitative genetic mapping is best shown in the ability to connect phenotypic differences to specific genes and variants. We will present methods to narrow genomic regions to candidate genes and then tests to identify the gene or variant involved in a quantitative trait. The same features that make C. elegans a preeminent experimental model animal contribute to its exceptional value as a tool to understand natural phenotypic variation.
Collapse
Affiliation(s)
- Erik C Andersen
- Department of Molecular Biosciences, Northwestern University, Evanston, IL 60201, USA
| | - Matthew V Rockman
- Department of Biology and Center for Genomics & Systems Biology, New York University, New York, NY 10003, USA
| |
Collapse
|
12
|
Snoek BL, Sterken MG, Nijveen H, Volkers RJM, Riksen J, Rosenstiel PC, Schulenburg H, Kammenga JE. The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population. G3-GENES GENOMES GENETICS 2021; 11:6347583. [PMID: 34568931 PMCID: PMC8496280 DOI: 10.1093/g3journal/jkab258] [Citation(s) in RCA: 2] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Received: 06/15/2021] [Accepted: 07/17/2021] [Indexed: 11/29/2022]
Abstract
Studying genetic variation of gene expression provides a powerful way to unravel the molecular components underlying complex traits. Expression quantitative trait locus (eQTL) studies have been performed in several different model species, yet most of these linkage studies have been based on the genetic segregation of two parental alleles. Recently, we developed a multiparental segregating population of 200 recombinant inbred lines (mpRILs) derived from four wild isolates (JU1511, JU1926, JU1931, and JU1941) in the nematode Caenorhabditis elegans. We used RNA-seq to investigate how multiple alleles affect gene expression in these mpRILs. We found 1789 genes differentially expressed between the parental lines. Transgression, expression beyond any of the parental lines in the mpRILs, was found for 7896 genes. For expression QTL mapping almost 9000 SNPs were available. By combining these SNPs and the RNA-seq profiles of the mpRILs, we detected almost 6800 eQTLs. Most trans-eQTLs (63%) co-locate in six newly identified trans-bands. The trans-eQTLs found in previous two-parental allele eQTL experiments and this study showed some overlap (17.5–46.8%), highlighting on the one hand that a large group of genes is affected by polymorphic regulators across populations and conditions, on the other hand, it shows that the mpRIL population allows identification of novel gene expression regulatory loci. Taken together, the analysis of our mpRIL population provides a more refined insight into C. elegans complex trait genetics and eQTLs in general, as well as a starting point to further test and develop advanced statistical models for detection of multiallelic eQTLs and systems genetics studying the genotype–phenotype relationship.
Collapse
Affiliation(s)
- Basten L Snoek
- Laboratory of Nematology, Wageningen University, NL-6708 PB Wageningen, The Netherlands.,Theoretical Biology and Bioinformatics, Utrecht University, 3584 CH Utrecht, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, NL-6708 PB Wageningen, The Netherlands
| | - Harm Nijveen
- Bioinformatics Group, Wageningen University, NL-6708 PB Wageningen, The Netherlands
| | - Rita J M Volkers
- Laboratory of Nematology, Wageningen University, NL-6708 PB Wageningen, The Netherlands
| | - Joost Riksen
- Laboratory of Nematology, Wageningen University, NL-6708 PB Wageningen, The Netherlands
| | - Philip C Rosenstiel
- Institute for Clinical Molecular Biology, University of Kiel, 24098 Kiel, Germany.,Competence Centre for Genomic Analysis (CCGA) Kiel, University of Kiel, 24098 Kiel, Germany
| | - Hinrich Schulenburg
- Zoological Institute, University of Kiel, 24098 Kiel, Germany.,Max Planck Institute for Evolutionary Biology, 24306 Ploen, Germany
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, NL-6708 PB Wageningen, The Netherlands
| |
Collapse
|
13
|
Evans KS, van Wijk MH, McGrath PT, Andersen EC, Sterken MG. From QTL to gene: C. elegans facilitates discoveries of the genetic mechanisms underlying natural variation. Trends Genet 2021; 37:933-947. [PMID: 34229867 DOI: 10.1016/j.tig.2021.06.005] [Citation(s) in RCA: 27] [Impact Index Per Article: 9.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/19/2021] [Revised: 06/01/2021] [Accepted: 06/03/2021] [Indexed: 11/15/2022]
Abstract
Although many studies have examined quantitative trait variation across many species, only a small number of genes and thereby molecular mechanisms have been discovered. Without these data, we can only speculate about evolutionary processes that underlie trait variation. Here, we review how quantitative and molecular genetics in the nematode Caenorhabditis elegans led to the discovery and validation of 37 quantitative trait genes over the past 15 years. Using these data, we can start to make inferences about evolution from these quantitative trait genes, including the roles that coding versus noncoding variation, gene family expansion, common versus rare variants, pleiotropy, and epistasis play in trait variation across this species.
Collapse
Affiliation(s)
- Kathryn S Evans
- Molecular Biosciences, Northwestern University, Evanston, IL 60208, USA; Interdisciplinary Biological Sciences Program, Northwestern University, Evanston, IL 60208, USA
| | - Marijke H van Wijk
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands
| | - Patrick T McGrath
- School of Biological Sciences, Georgia Institute of Technology, Atlanta, GA 30332, USA
| | - Erik C Andersen
- Molecular Biosciences, Northwestern University, Evanston, IL 60208, USA.
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands.
| |
Collapse
|
14
|
Sterken MG, van Sluijs L, Wang YA, Ritmahan W, Gultom ML, Riksen JAG, Volkers RJM, Snoek LB, Pijlman GP, Kammenga JE. Punctuated Loci on Chromosome IV Determine Natural Variation in Orsay Virus Susceptibility of Caenorhabditis elegans Strains Bristol N2 and Hawaiian CB4856. J Virol 2021; 95:e02430-20. [PMID: 33827942 PMCID: PMC8315983 DOI: 10.1128/jvi.02430-20] [Citation(s) in RCA: 11] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/21/2020] [Accepted: 03/29/2021] [Indexed: 01/06/2023] Open
Abstract
Host-pathogen interactions play a major role in evolutionary selection and shape natural genetic variation. The genetically distinct Caenorhabditis elegans strains, Bristol N2 and Hawaiian CB4856, are differentially susceptible to the Orsay virus (OrV). Here, we report the dissection of the genetic architecture of susceptibility to OrV infection. We compare OrV infection in the relatively resistant wild-type CB4856 strain to the more susceptible canonical N2 strain. To gain insight into the genetic architecture of viral susceptibility, 52 fully sequenced recombinant inbred lines (CB4856 × N2 RILs) were exposed to OrV. This led to the identification of two loci on chromosome IV associated with OrV resistance. To verify the two loci and gain additional insight into the genetic architecture controlling virus infection, introgression lines (ILs) that together cover chromosome IV, were exposed to OrV. Of the 27 ILs used, 17 had an CB4856 introgression in an N2 background, and 10 had an N2 introgression in a CB4856 background. Infection of the ILs confirmed and fine-mapped the locus underlying variation in OrV susceptibility, and we found that a single nucleotide polymorphism in cul-6 may contribute to the difference in OrV susceptibility between N2 and CB4856. An allele swap experiment showed the strain CB4856 became as susceptible as the N2 strain by having an N2 cul-6 allele, although having the CB4856 cul-6 allele did not increase resistance in N2. In addition, we found that multiple strains with nonoverlapping introgressions showed a distinct infection phenotype from the parental strain, indicating that there are punctuated locations on chromosome IV determining OrV susceptibility. Thus, our findings reveal the genetic complexity of OrV susceptibility in C. elegans and suggest that viral susceptibility is governed by multiple genes.IMPORTANCE Genetic variation determines the viral susceptibility of hosts. Yet, pinpointing which genetic variants determine viral susceptibility remains challenging. Here, we have exploited the genetic tractability of the model organism Caenorhabditis elegans to dissect the genetic architecture of Orsay virus infection. Our results provide novel insight into natural determinants of Orsay virus infection.
Collapse
Affiliation(s)
- Mark G Sterken
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
- Laboratory of Virology, Wageningen University, Wageningen, The Netherlands
| | - Lisa van Sluijs
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
- Laboratory of Virology, Wageningen University, Wageningen, The Netherlands
| | - Yiru A Wang
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Wannisa Ritmahan
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Mitra L Gultom
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Joost A G Riksen
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Rita J M Volkers
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - L Basten Snoek
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
- Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, The Netherlands
| | - Gorben P Pijlman
- Laboratory of Virology, Wageningen University, Wageningen, The Netherlands
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| |
Collapse
|
15
|
Brion C, Lutz SM, Albert FW. Simultaneous quantification of mRNA and protein in single cells reveals post-transcriptional effects of genetic variation. eLife 2020; 9:60645. [PMID: 33191917 PMCID: PMC7707838 DOI: 10.7554/elife.60645] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2020] [Accepted: 11/14/2020] [Indexed: 01/27/2023] Open
Abstract
Trans-acting DNA variants may specifically affect mRNA or protein levels of genes located throughout the genome. However, prior work compared trans-acting loci mapped in separate studies, many of which had limited statistical power. Here, we developed a CRISPR-based system for simultaneous quantification of mRNA and protein of a given gene via dual fluorescent reporters in single, live cells of the yeast Saccharomyces cerevisiae. In large populations of recombinant cells from a cross between two genetically divergent strains, we mapped 86 trans-acting loci affecting the expression of ten genes. Less than 20% of these loci had concordant effects on mRNA and protein of the same gene. Most loci influenced protein but not mRNA of a given gene. One locus harbored a premature stop variant in the YAK1 kinase gene that had specific effects on protein or mRNA of dozens of genes. These results demonstrate complex, post-transcriptional genetic effects on gene expression.
Collapse
Affiliation(s)
- Christian Brion
- Department of Genetics, Cell Biology and Development, University of Minnesota, Minneapolis, United States
| | - Sheila M Lutz
- Department of Genetics, Cell Biology and Development, University of Minnesota, Minneapolis, United States
| | - Frank Wolfgang Albert
- Department of Genetics, Cell Biology and Development, University of Minnesota, Minneapolis, United States
| |
Collapse
|
16
|
Evans KS, Zdraljevic S, Stevens L, Collins K, Tanny RE, Andersen EC. Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans. PLoS Genet 2020; 16:e1008986. [PMID: 33175833 PMCID: PMC7682890 DOI: 10.1371/journal.pgen.1008986] [Citation(s) in RCA: 17] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/06/2020] [Revised: 11/23/2020] [Accepted: 09/23/2020] [Indexed: 12/14/2022] Open
Abstract
Zinc is an essential trace element that acts as a co-factor for many enzymes and transcription factors required for cellular growth and development. Altering intracellular zinc levels can produce dramatic effects ranging from cell proliferation to cell death. To avoid such fates, cells have evolved mechanisms to handle both an excess and a deficiency of zinc. Zinc homeostasis is largely maintained via zinc transporters, permeable channels, and other zinc-binding proteins. Variation in these proteins might affect their ability to interact with zinc, leading to either increased sensitivity or resistance to natural zinc fluctuations in the environment. We can leverage the power of the roundworm nematode Caenorhabditis elegans as a tractable metazoan model for quantitative genetics to identify genes that could underlie variation in responses to zinc. We found that the laboratory-adapted strain (N2) is resistant and a natural isolate from Hawaii (CB4856) is sensitive to micromolar amounts of exogenous zinc supplementation. Using a panel of recombinant inbred lines, we identified two large-effect quantitative trait loci (QTL) on the left arm of chromosome III and the center of chromosome V that are associated with zinc responses. We validated and refined both QTL using near-isogenic lines (NILs) and identified a naturally occurring deletion in sqst-5, a sequestosome-related gene, that is associated with resistance to high exogenous zinc. We found that this deletion is relatively common across strains within the species and that variation in sqst-5 is associated with zinc resistance. Our results offer a possible mechanism for how organisms can respond to naturally high levels of zinc in the environment and how zinc homeostasis varies among individuals.
Collapse
Affiliation(s)
- Kathryn S. Evans
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
- Interdisciplinary Biological Sciences Program, Northwestern University, Evanston, Illinois, United States of America
| | - Stefan Zdraljevic
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
- Interdisciplinary Biological Sciences Program, Northwestern University, Evanston, Illinois, United States of America
| | - Lewis Stevens
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
| | - Kimberly Collins
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
| | - Robyn E. Tanny
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
| | - Erik C. Andersen
- Molecular Biosciences, Northwestern University, Evanston, Illinois, United States of America
- Robert H. Lurie Comprehensive Cancer Center, Northwestern University, Chicago, Illinois, United States of America
- * E-mail:
| |
Collapse
|
17
|
Network Analysis Prioritizes DEWAX and ICE1 as the Candidate Genes for Major eQTL Hotspots in Seed Germination of Arabidopsis thaliana. G3-GENES GENOMES GENETICS 2020; 10:4215-4226. [PMID: 32963085 PMCID: PMC7642920 DOI: 10.1534/g3.120.401477] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 01/08/2023]
Abstract
Seed germination is characterized by a constant change of gene expression across different time points. These changes are related to specific processes, which eventually determine the onset of seed germination. To get a better understanding on the regulation of gene expression during seed germination, we performed a quantitative trait locus mapping of gene expression (eQTL) at four important seed germination stages (primary dormant, after-ripened, six-hour after imbibition, and radicle protrusion stage) using Arabidopsis thaliana Bay x Sha recombinant inbred lines (RILs). The mapping displayed the distinctness of the eQTL landscape for each stage. We found several eQTL hotspots across stages associated with the regulation of expression of a large number of genes. Interestingly, an eQTL hotspot on chromosome five collocates with hotspots for phenotypic and metabolic QTL in the same population. Finally, we constructed a gene co-expression network to prioritize the regulatory genes for two major eQTL hotspots. The network analysis prioritizes transcription factors DEWAX and ICE1 as the most likely regulatory genes for the hotspot. Together, we have revealed that the genetic regulation of gene expression is dynamic along the course of seed germination.
Collapse
|
18
|
Sterken MG, Bevers RPJ, Volkers RJM, Riksen JAG, Kammenga JE, Snoek BL. Dissecting the eQTL Micro-Architecture in Caenorhabditis elegans. Front Genet 2020; 11:501376. [PMID: 33240309 PMCID: PMC7670075 DOI: 10.3389/fgene.2020.501376] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2019] [Accepted: 10/13/2020] [Indexed: 01/11/2023] Open
Abstract
The study of expression quantitative trait loci (eQTL) using natural variation in inbred populations has yielded detailed information about the transcriptional regulation of complex traits. Studies on eQTL using recombinant inbred lines (RILs) led to insights on cis and trans regulatory loci of transcript abundance. However, determining the underlying causal polymorphic genes or variants is difficult, but ultimately essential for the understanding of regulatory networks of complex traits. This requires insight into whether associated loci are single eQTL or a combination of closely linked eQTL, and how this QTL micro-architecture depends on the environment. We addressed these questions by testing for independent replication of previously mapped eQTL in Caenorhabditis elegans using new data from introgression lines (ILs). Both populations indicate that the overall heritability of gene expression, number, and position of eQTL differed among environments. Across environments we were able to replicate 70% of the cis- and 40% of the trans-eQTL using the ILs. Testing eight different simulation models, we suggest that additive effects explain up to 60-93% of RIL/IL heritability for all three environments. Closely linked eQTL explained up to 40% of RIL/IL heritability in the control environment whereas only 7% in the heat-stress and recovery environments. In conclusion, we show that reproducibility of eQTL was higher for cis vs. trans eQTL and that the environment affects the eQTL micro-architecture.
Collapse
Affiliation(s)
- Mark G. Sterken
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
| | - Roel P. J. Bevers
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
| | - Rita J. M. Volkers
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
| | - Joost A. G. Riksen
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
| | - Jan E. Kammenga
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
| | - Basten L. Snoek
- Laboratory of Nematology, Wageningen University & Research, Wageningen, Netherlands
- Theoretical Biology & Bioinformatics, Utrecht University, Utrecht, Netherlands
| |
Collapse
|
19
|
Delaney DM, Hoekstra LA, Janzen FJ. Becoming creatures of habit: Among- and within-individual variation in nesting behaviour shift with age. J Evol Biol 2020; 33:1614-1624. [PMID: 32897610 DOI: 10.1111/jeb.13701] [Citation(s) in RCA: 4] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/30/2019] [Revised: 08/24/2020] [Accepted: 08/31/2020] [Indexed: 11/29/2022]
Abstract
The quantification of repeatability has enabled behavioural and evolutionary ecologists to assess the heritable potential of traits. For behavioural traits that vary across life, age-related variation should be accounted for to prevent biasing the microevolutionary estimate of interest. Moreover, to gain a mechanistic understanding of ontogenetic variation in behaviour, among- and within-individual variance should be quantified across life. We leveraged a 30-year study of painted turtles (Chrysemys picta) to assess how age contributes to variation in the repeatability of nesting behaviours. We found that four components of nesting behaviour were repeatable and that accounting for age increased the repeatability estimate for maternal choice of canopy cover over nests. We detected canalization (diminished within-individual variance with age) of canopy cover choice in a reduced data set despite no shift in repeatability. Additionally, random regression analysis revealed that females became more divergent from each other in their choice of canopy cover with age. Thus, properly modelling age-related variance should more precisely estimate heritable potential, and assessing among- and within-individual variance components in addition to repeatability will offer a more mechanistic understanding of behavioural variation across age.
Collapse
Affiliation(s)
- David M Delaney
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, USA.,Department of Ecology and Evolutionary Biology, University of Colorado Boulder, Boulder, CO, USA
| | - Luke A Hoekstra
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, USA.,Department of Integrative Biology, Oklahoma State University, Stillwater, OK, USA
| | - Fredric J Janzen
- Department of Ecology, Evolution, and Organismal Biology, Iowa State University, Ames, IA, USA.,Kellogg Biological Station, Michigan State University, Hickory Corners, MI, USA
| |
Collapse
|
20
|
The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses. G3-GENES GENOMES GENETICS 2020; 10:2353-2364. [PMID: 32385045 PMCID: PMC7341127 DOI: 10.1534/g3.120.401310] [Citation(s) in RCA: 25] [Impact Index Per Article: 6.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Subscribe] [Scholar Register] [Indexed: 12/14/2022]
Abstract
Pleiotropy, the concept that a single gene controls multiple distinct traits, is prevalent in most organisms and has broad implications for medicine and agriculture. The identification of the molecular mechanisms underlying pleiotropy has the power to reveal previously unknown biological connections between seemingly unrelated traits. Additionally, the discovery of pleiotropic genes increases our understanding of both genetic and phenotypic complexity by characterizing novel gene functions. Quantitative trait locus (QTL) mapping has been used to identify several pleiotropic regions in many organisms. However, gene knockout studies are needed to eliminate the possibility of tightly linked, non-pleiotropic loci. Here, we use a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to identify a single large-effect QTL on the center of chromosome V associated with variation in responses to eight chemotherapeutics. We validate this QTL with near-isogenic lines and pair genome-wide gene expression data with drug response traits to perform mediation analysis, leading to the identification of a pleiotropic candidate gene, scb-1, for some of the eight chemotherapeutics. Using deletion strains created by genome editing, we show that scb-1, which was previously implicated in response to bleomycin, also underlies responses to other double-strand DNA break-inducing chemotherapeutics. This finding provides new evidence for the role of scb-1 in the nematode drug response and highlights the power of mediation analysis to identify causal genes.
Collapse
|
21
|
Song Y, Wang Y, Li Y, Wang L, Zhang W, Cheng J, Zhu Y, Zhang H, Zhang Q, Niu H, zheng Y, Liang M, Deng M, Shi H, Wang H, Zhang F, Zhu Z. The whole transcriptome regulation as a function of mitochondrial polymorphisms and aging in Caenorhabditis elegans. Aging (Albany NY) 2020; 12:2453-2470. [PMID: 32019902 PMCID: PMC7041728 DOI: 10.18632/aging.102754] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/02/2019] [Accepted: 01/07/2020] [Indexed: 06/10/2023]
Abstract
Recently, mitochondrial-nuclear interaction in aging has been widely studied. However, the nuclear genome controlled by natural mitochondrial variations that influence aging has not been comprehensively understood so far. We hypothesized that mitochondrial polymorphisms could play critical roles in the aging process, probably by regulation of the whole-transcriptome expression. Our results showed that mitochondria polymorphisms not only decreased the mitochondrial mass but also miRNA, lncRNA, mRNA, circRNA and metabolite profiles. Furthermore, most genes that are associated with mitochondria show age-related expression features (P = 3.58E-35). We also constructed a differentially expressed circRNA-lncRNA-miRNA-mRNA regulatory network and a ceRNA network affected by the mitochondrial variations. In addition, Kyoto Encyclopedia of Genes and Genomes pathway analyses showed that the genes affected by the mitochondrial variation were enriched in metabolic activity. We finally constructed a multi-level regulatory network with aging which affected by the mitochondrial variation in Caenorhabditis elegans. The interactions between these genes and metabolites have great values for further aging research. In sum, our findings provide new evidence for understanding the molecular mechanisms of how mitochondria influence aging.
Collapse
Affiliation(s)
- Yuanjian Song
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| | - Yuechen Wang
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| | - Ying Li
- Medical Technology School of Xuzhou Medical University, Xuzhou, China
| | - Liang Wang
- Department of Bioinformatics, School of Medical Informatics and Engineering, Xuzhou Medical University, Xuzhou, China
| | - WenDa Zhang
- Department of Urology, Xuzhou Central Hospital, Xuzhou, China
| | - Jing Cheng
- Medical Technology School of Xuzhou Medical University, Xuzhou, China
| | - Yao Zhu
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| | - Haoyu Zhang
- School of Marine Sciences, Nanjing University of Information Science and Technology, Nanjing, China
| | - Qiang Zhang
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| | - Haichen Niu
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| | - Yingwei zheng
- Department of Biochemistry, Xuzhou Medical University, Xuzhou, China
| | - Mengyu Liang
- Clinical College of Xuzhou Medical University, Xuzhou, China
| | - Mengqiong Deng
- Clinical College of Xuzhou Medical University, Xuzhou, China
| | - Hao Shi
- Clinical College of Xuzhou Medical University, Xuzhou, China
| | - Hao Wang
- Clinical College of Xuzhou Medical University, Xuzhou, China
| | - Fang Zhang
- Research Facility Center for Morphology, Xuzhou Medical University, Xuzhou, China
| | - Zuobin Zhu
- Department of Genetics, Xuzhou Medical University, Xuzhou, China
| |
Collapse
|
22
|
Snoek BL, Sterken MG, Hartanto M, van Zuilichem AJ, Kammenga JE, de Ridder D, Nijveen H. WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans. Database (Oxford) 2020; 2020:baz149. [PMID: 31960906 PMCID: PMC6971878 DOI: 10.1093/database/baz149] [Citation(s) in RCA: 11] [Impact Index Per Article: 2.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/14/2019] [Revised: 11/30/2019] [Accepted: 12/13/2019] [Indexed: 12/19/2022]
Abstract
Quantitative genetics provides the tools for linking polymorphic loci to trait variation. Linkage analysis of gene expression is an established and widely applied method, leading to the identification of expression quantitative trait loci (eQTLs). (e)QTL detection facilitates the identification and understanding of the underlying molecular components and pathways, yet (e)QTL data access and mining often is a bottleneck. Here, we present WormQTL2, a database and platform for comparative investigations and meta-analyses of published (e)QTL data sets in the model nematode worm C. elegans. WormQTL2 integrates six eQTL studies spanning 11 conditions as well as over 1000 traits from 32 studies and allows experimental results to be compared, reused and extended upon to guide further experiments and conduct systems-genetic analyses. For example, one can easily screen a locus for specific cis-eQTLs that could be linked to variation in other traits, detect gene-by-environment interactions by comparing eQTLs under different conditions, or find correlations between QTL profiles of classical traits and gene expression. WormQTL2 makes data on natural variation in C. elegans and the identified QTLs interactively accessible, allowing studies beyond the original publications. Database URL: www.bioinformatics.nl/WormQTL2/.
Collapse
Affiliation(s)
- Basten L Snoek
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
- Theoretical Biology and Bioinformatics, Utrecht University, Padualaan 8, 3584 CH Utrecht, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| | - Margi Hartanto
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
- Bioinformatics Group, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| | - Albert-Jan van Zuilichem
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| | - Dick de Ridder
- Bioinformatics Group, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| | - Harm Nijveen
- Bioinformatics Group, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB Wageningen, The Netherlands
| |
Collapse
|
23
|
Jovic K, Grilli J, Sterken MG, Snoek BL, Riksen JAG, Allesina S, Kammenga JE. Transcriptome resilience predicts thermotolerance in Caenorhabditis elegans. BMC Biol 2019; 17:102. [PMID: 31822273 PMCID: PMC6905072 DOI: 10.1186/s12915-019-0725-6] [Citation(s) in RCA: 24] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/13/2019] [Accepted: 11/18/2019] [Indexed: 12/11/2022] Open
Abstract
BACKGROUND The detrimental effects of a short bout of stress can persist and potentially turn lethal, long after the return to normal conditions. Thermotolerance, which is the capacity of an organism to withstand relatively extreme temperatures, is influenced by the response during stress exposure, as well as the recovery process afterwards. While heat-shock response mechanisms have been studied intensively, predicting thermal tolerance remains a challenge. RESULTS Here, we use the nematode Caenorhabditis elegans to measure transcriptional resilience to heat stress and predict thermotolerance. Using principal component analysis in combination with genome-wide gene expression profiles collected in three high-resolution time series during control, heat stress, and recovery conditions, we infer a quantitative scale capturing the extent of stress-induced transcriptome dynamics in a single value. This scale provides a basis for evaluating transcriptome resilience, defined here as the ability to depart from stress-expression dynamics during recovery. Independent replication across multiple highly divergent genotypes reveals that the transcriptional resilience parameter measured after a spike in temperature is quantitatively linked to long-term survival after heat stress. CONCLUSION Our findings imply that thermotolerance is an intrinsic property that pre-determines long-term outcome of stress and can be predicted by the transcriptional resilience parameter. Inferring the transcriptional resilience parameters of higher organisms could aid in evaluating rehabilitation strategies after stresses such as disease and trauma.
Collapse
Affiliation(s)
- Katharina Jovic
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, Wageningen, 6708 PB, The Netherlands
| | - Jacopo Grilli
- Department of Ecology and Evolution, University of Chicago, 1101 E 57th St, Chicago, IL, 60637, USA
- Santa Fe Institute, 1399 Hyde Park Rd, Santa Fe, NM, 87501, USA
- The Abdus Salam International Center for Theoretical Physics (ICTP), Strada Costiera 11, I-34014, Trieste, Italy
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, Wageningen, 6708 PB, The Netherlands
| | - Basten L Snoek
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, Wageningen, 6708 PB, The Netherlands
- Theoretical Biology and Bioinformatics, Utrecht University, Padualaan 8, Utrecht, 3584 CH, The Netherlands
| | - Joost A G Riksen
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, Wageningen, 6708 PB, The Netherlands
| | - Stefano Allesina
- Department of Ecology and Evolution, University of Chicago, 1101 E 57th St, Chicago, IL, 60637, USA.
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, Wageningen, 6708 PB, The Netherlands.
| |
Collapse
|
24
|
Balliu B, Durrant M, Goede OD, Abell N, Li X, Liu B, Gloudemans MJ, Cook NL, Smith KS, Knowles DA, Pala M, Cucca F, Schlessinger D, Jaiswal S, Sabatti C, Lind L, Ingelsson E, Montgomery SB. Genetic regulation of gene expression and splicing during a 10-year period of human aging. Genome Biol 2019; 20:230. [PMID: 31684996 PMCID: PMC6827221 DOI: 10.1186/s13059-019-1840-y] [Citation(s) in RCA: 26] [Impact Index Per Article: 5.2] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 02/11/2019] [Accepted: 09/27/2019] [Indexed: 12/12/2022] Open
Abstract
BACKGROUND Molecular and cellular changes are intrinsic to aging and age-related diseases. Prior cross-sectional studies have investigated the combined effects of age and genetics on gene expression and alternative splicing; however, there has been no long-term, longitudinal characterization of these molecular changes, especially in older age. RESULTS We perform RNA sequencing in whole blood from the same individuals at ages 70 and 80 to quantify how gene expression, alternative splicing, and their genetic regulation are altered during this 10-year period of advanced aging at a population and individual level. We observe that individuals are more similar to their own expression profiles later in life than profiles of other individuals their own age. We identify 1291 and 294 genes differentially expressed and alternatively spliced with age, as well as 529 genes with outlying individual trajectories. Further, we observe a strong correlation of genetic effects on expression and splicing between the two ages, with a small subset of tested genes showing a reduction in genetic associations with expression and splicing in older age. CONCLUSIONS These findings demonstrate that, although the transcriptome and its genetic regulation is mostly stable late in life, a small subset of genes is dynamic and is characterized by a reduction in genetic regulation, most likely due to increasing environmental variance with age.
Collapse
Affiliation(s)
- Brunilda Balliu
- Department of Pathology, Stanford University School of Medicine, Stanford, USA.
| | - Matthew Durrant
- Department of Genetics, Stanford University School of Medicine, Stanford, USA
| | - Olivia de Goede
- Department of Genetics, Stanford University School of Medicine, Stanford, USA
| | - Nathan Abell
- Department of Genetics, Stanford University School of Medicine, Stanford, USA
| | - Xin Li
- Department of Pathology, Stanford University School of Medicine, Stanford, USA
| | - Boxiang Liu
- Department of Biology, Stanford University School of Medicine, Stanford, USA
| | | | - Naomi L Cook
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Kevin S Smith
- Department of Pathology, Stanford University School of Medicine, Stanford, USA
| | | | - Mauro Pala
- Dipartimento di Scienze Biomediche, Universita di Sassari, Sassari, Italy
| | - Francesco Cucca
- Dipartimento di Scienze Biomediche, Universita di Sassari, Sassari, Italy
| | | | - Siddhartha Jaiswal
- Department of Pathology, Stanford University School of Medicine, Stanford, USA
| | - Chiara Sabatti
- Department of Biomedical Data Science, Stanford University School of Medicine, Stanford, USA
| | - Lars Lind
- Department of Medical Sciences, Uppsala University, Uppsala, Sweden
| | - Erik Ingelsson
- Department of Medicine, Division of Cardiovascular Medicine, Stanford University School of Medicine, Stanford, USA.
- Stanford Cardiovascular Institute, Stanford University, Stanford, USA.
- Stanford Diabetes Research Center, Stanford University, Stanford, USA.
| | - Stephen B Montgomery
- Department of Pathology, Stanford University School of Medicine, Stanford, USA.
- Department of Genetics, Stanford University School of Medicine, Stanford, USA.
| |
Collapse
|
25
|
Tarkhov AE, Alla R, Ayyadevara S, Pyatnitskiy M, Menshikov LI, Shmookler Reis RJ, Fedichev PO. A universal transcriptomic signature of age reveals the temporal scaling of Caenorhabditis elegans aging trajectories. Sci Rep 2019; 9:7368. [PMID: 31089188 PMCID: PMC6517414 DOI: 10.1038/s41598-019-43075-z] [Citation(s) in RCA: 27] [Impact Index Per Article: 5.4] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/26/2018] [Accepted: 04/15/2019] [Indexed: 12/13/2022] Open
Abstract
We collected 60 age-dependent transcriptomes for C. elegans strains including four exceptionally long-lived mutants (mean adult lifespan extended 2.2- to 9.4-fold) and three examples of lifespan-increasing RNAi treatments. Principal Component Analysis (PCA) reveals aging as a transcriptomic drift along a single direction, consistent across the vastly diverse biological conditions and coinciding with the first principal component, a hallmark of the criticality of the underlying gene regulatory network. We therefore expected that the organism's aging state could be characterized by a single number closely related to vitality deficit or biological age. The "aging trajectory", i.e. the dependence of the biological age on chronological age, is then a universal stochastic function modulated by the network stiffness; a macroscopic parameter reflecting the network topology and associated with the rate of aging. To corroborate this view, we used publicly available datasets to define a transcriptomic biomarker of age and observed that the rescaling of age by lifespan simultaneously brings together aging trajectories of transcription and survival curves. In accordance with the theoretical prediction, the limiting mortality value at the plateau agrees closely with the mortality rate doubling exponent estimated at the cross-over age near the average lifespan. Finally, we used the transcriptomic signature of age to identify possible life-extending drug compounds and successfully tested a handful of the top-ranking molecules in C. elegans survival assays and achieved up to a +30% extension of mean lifespan.
Collapse
Affiliation(s)
- Andrei E Tarkhov
- Gero LLC, Nizhny Susalny per. 5/4, Moscow, 105064, Russia.
- Skolkovo Institute of Science and Technology, Skolkovo Innovation Center, Bolshoy Boulevard 30, bld. 1, Moscow, 121205, Russia.
| | - Ramani Alla
- Central Arkansas Veterans Healthcare System, Research Service, Little Rock, Arkansas, USA
- Department of Geriatrics, Reynolds Institute on Aging, University of Arkansas for Medical Sciences, Little Rock, Arkansas, USA
| | - Srinivas Ayyadevara
- Central Arkansas Veterans Healthcare System, Research Service, Little Rock, Arkansas, USA
- Department of Geriatrics, Reynolds Institute on Aging, University of Arkansas for Medical Sciences, Little Rock, Arkansas, USA
| | - Mikhail Pyatnitskiy
- Gero LLC, Nizhny Susalny per. 5/4, Moscow, 105064, Russia
- Institute of Biomedical Chemistry, 119121, Moscow, Russia
| | - Leonid I Menshikov
- Gero LLC, Nizhny Susalny per. 5/4, Moscow, 105064, Russia
- National Research Center "Kurchatov Institute", 1, Akademika Kurchatova pl., Moscow, 123182, Russia
| | - Robert J Shmookler Reis
- Central Arkansas Veterans Healthcare System, Research Service, Little Rock, Arkansas, USA
- Department of Geriatrics, Reynolds Institute on Aging, University of Arkansas for Medical Sciences, Little Rock, Arkansas, USA
- Bioinformatics Program, University of Arkansas for Medical Sciences, and University of Arkansas at Little Rock, Little Rock, Arkansas, USA
| | - Peter O Fedichev
- Gero LLC, Nizhny Susalny per. 5/4, Moscow, 105064, Russia.
- Moscow Institute of Physics and Technology, 141700, Institutskii per. 9, Dolgoprudny, Moscow Region, Russia.
| |
Collapse
|
26
|
Lai RW, Lu R, Danthi PS, Bravo JI, Goumba A, Sampathkumar NK, Benayoun BA. Multi-level remodeling of transcriptional landscapes in aging and longevity. BMB Rep 2019. [PMID: 30526773 PMCID: PMC6386224 DOI: 10.5483/bmbrep.2019.52.1.296] [Citation(s) in RCA: 33] [Impact Index Per Article: 6.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/09/2023] Open
Abstract
In multi-cellular organisms, the control of gene expression is key not only for development, but also for adult cellular homeostasis, and gene expression has been observed to be deregulated with aging. In this review, we discuss the current knowledge on the transcriptional alterations that have been described to occur with age in metazoans. First, we discuss age-related transcriptional changes in protein-coding genes, the expected functional impact of such changes, and how known pro-longevity interventions impact these changes. Second, we discuss the changes and impact of emerging aspects of transcription in aging, including age-related changes in splicing, lncRNAs and circRNAs. Third, we discuss the changes and potential impact of transcription of transposable elements with aging. Fourth, we highlight small ncRNAs and their potential impact on the regulation of aging phenotypes. Understanding the aging transcriptome will be key to identify important regulatory targets, and ultimately slow-down or reverse aging and extend healthy lifespan in humans.
Collapse
Affiliation(s)
- Rochelle W Lai
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089, USA
| | - Ryan Lu
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089, USA
| | - Prakroothi S Danthi
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089, USA
| | - Juan I Bravo
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089; Graduate program in the Biology of Aging, University of Southern California, Los Angeles, CA 90089, USA
| | - Alexandre Goumba
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089, USA
| | | | - Bérénice A Benayoun
- Leonard Davis School of Gerontology, University of Southern California, Los Angeles, CA 90089; USC Norris Comprehensive Cancer Center, Epigenetics and Gene Regulation, Los Angeles, CA 90089; USC Stem Cell Initiative, Los Angeles, CA 90089, USA
| |
Collapse
|
27
|
Dose dependent gene expression is dynamically modulated by the history, physiology and age of yeast cells. BIOCHIMICA ET BIOPHYSICA ACTA-GENE REGULATORY MECHANISMS 2019; 1862:457-471. [DOI: 10.1016/j.bbagrm.2019.02.009] [Citation(s) in RCA: 5] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/10/2018] [Revised: 02/21/2019] [Accepted: 02/23/2019] [Indexed: 12/14/2022]
|
28
|
Wang YA, Snoek BL, Sterken MG, Riksen JAG, Stastna JJ, Kammenga JE, Harvey SC. Genetic background modifies phenotypic and transcriptional responses in a C. elegans model of α-synuclein toxicity. BMC Genomics 2019; 20:232. [PMID: 30894116 PMCID: PMC6427842 DOI: 10.1186/s12864-019-5597-1] [Citation(s) in RCA: 8] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/10/2018] [Accepted: 03/11/2019] [Indexed: 11/13/2022] Open
Abstract
Background Accumulation of protein aggregates are a major hallmark of progressive neurodegenerative disorders such as Parkinson’s disease and Alzheimer’s disease. Transgenic Caenorhabditis elegans nematodes expressing the human synaptic protein α-synuclein in body wall muscle show inclusions of aggregated protein, which affects similar genetic pathways as in humans. It is not however known how the effects of α-synuclein expression in C. elegans differs among genetic backgrounds. Here, we compared gene expression patterns and investigated the phenotypic consequences of transgenic α-synuclein expression in five different C. elegans genetic backgrounds. Results Transcriptome analysis indicates that α-synuclein expression effects pathways associated with nutrient storage, lipid transportation and ion exchange and that effects vary depending on the genetic background. These gene expression changes predict that a range of phenotypes will be affected by α-synuclein expression. We confirm this, showing that α-synuclein expression delayed development, reduced lifespan, increased rate of matricidal hatching, and slows pharyngeal pumping. Critically, these phenotypic effects depend on the genetic background and coincide with the core changes in gene expression. Conclusions Together, our results show genotype-specific effects and core alterations in both gene expression and in phenotype in response to α-synuclein expression. We conclude that the effects of α-synuclein expression are substantially modified by the genetic background, illustrating that genetic background needs to be considered in C. elegans models of neurodegenerative disease. Electronic supplementary material The online version of this article (10.1186/s12864-019-5597-1) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Yiru A Wang
- Biomolecular Research Group, School of Human and Life Sciences, Canterbury Christ Church University, North Holmes Road, Canterbury, CT1 1QU, UK.,Laboratory of Nematology, Wageningen University, 6708, PB, Wageningen, The Netherlands
| | - Basten L Snoek
- Laboratory of Nematology, Wageningen University, 6708, PB, Wageningen, The Netherlands.,Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, 6708, PB, Wageningen, The Netherlands
| | - Joost A G Riksen
- Laboratory of Nematology, Wageningen University, 6708, PB, Wageningen, The Netherlands
| | - Jana J Stastna
- Biomolecular Research Group, School of Human and Life Sciences, Canterbury Christ Church University, North Holmes Road, Canterbury, CT1 1QU, UK
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, 6708, PB, Wageningen, The Netherlands
| | - Simon C Harvey
- Biomolecular Research Group, School of Human and Life Sciences, Canterbury Christ Church University, North Holmes Road, Canterbury, CT1 1QU, UK.
| |
Collapse
|
29
|
Snoek BL, Volkers RJM, Nijveen H, Petersen C, Dirksen P, Sterken MG, Nakad R, Riksen JAG, Rosenstiel P, Stastna JJ, Braeckman BP, Harvey SC, Schulenburg H, Kammenga JE. A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits. BMC Biol 2019; 17:24. [PMID: 30866929 PMCID: PMC6417139 DOI: 10.1186/s12915-019-0642-8] [Citation(s) in RCA: 17] [Impact Index Per Article: 3.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/17/2018] [Accepted: 02/26/2019] [Indexed: 01/29/2023] Open
Abstract
BACKGROUND The nematode Caenorhabditis elegans has been extensively used to explore the relationships between complex traits, genotypes, and environments. Complex traits can vary across different genotypes of a species, and the genetic regulators of trait variation can be mapped on the genome using quantitative trait locus (QTL) analysis of recombinant inbred lines (RILs) derived from genetically and phenotypically divergent parents. Most RILs have been derived from crossing two parents from globally distant locations. However, the genetic diversity between local C. elegans populations can be as diverse as between global populations and could thus provide means of identifying genetic variation associated with complex traits relevant on a broader scale. RESULTS To investigate the effect of local genetic variation on heritable traits, we developed a new RIL population derived from 4 parental wild isolates collected from 2 closely located sites in France: Orsay and Santeuil. We crossed these 4 genetically diverse parental isolates to generate a population of 200 multi-parental RILs and used RNA-seq to obtain sequence polymorphisms identifying almost 9000 SNPs variable between the 4 genotypes with an average spacing of 11 kb, doubling the mapping resolution relative to currently available RIL panels for many loci. The SNPs were used to construct a genetic map to facilitate QTL analysis. We measured life history traits such as lifespan, stress resistance, developmental speed, and population growth in different environments, and found substantial variation for most traits. We detected multiple QTLs for most traits, including novel QTLs not found in previous QTL analysis, including those for lifespan and pathogen responses. This shows that recombining genetic variation across C. elegans populations that are in geographical close proximity provides ample variation for QTL mapping. CONCLUSION Taken together, we show that using more parents than the classical two parental genotypes to construct a RIL population facilitates the detection of QTLs and that the use of wild isolates facilitates the detection of QTLs. The use of multi-parent RIL populations can further enhance our understanding of local adaptation and life history trade-offs.
Collapse
Affiliation(s)
- Basten L Snoek
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands. .,Theoretical Biology and Bioinformatics, Utrecht University, Padualaan 8, 3584 CH, Utrecht, The Netherlands.
| | - Rita J M Volkers
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands
| | - Harm Nijveen
- Bioinformatics Group, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands
| | - Carola Petersen
- Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Philipp Dirksen
- Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands
| | - Rania Nakad
- Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Joost A G Riksen
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands
| | - Philip Rosenstiel
- Institute for Clinical Molecular Biology, University of Kiel, 24098, Kiel, Germany
| | - Jana J Stastna
- Biomolecular Research Group, School of Human and Life Sciences, Canterbury Christ Church University, North Holmes Road, Canterbury, CT1 1QU, UK
| | - Bart P Braeckman
- Department of Biology, Ghent University, K. L. Ledeganckstraat 35, B-9000, Ghent, Belgium
| | - Simon C Harvey
- Biomolecular Research Group, School of Human and Life Sciences, Canterbury Christ Church University, North Holmes Road, Canterbury, CT1 1QU, UK
| | - Hinrich Schulenburg
- Zoological Institute, University of Kiel, 24098, Kiel, Germany. .,Max Planck Institute for Evolutionary Biology, August-Thienemann-Str. 2, 24306, Plön, Germany.
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Droevendaalsesteeg 1, NL-6708 PB, Wageningen, The Netherlands.
| |
Collapse
|
30
|
Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses. Genetics 2018; 210:1509-1525. [PMID: 30341085 PMCID: PMC6283156 DOI: 10.1534/genetics.118.301311] [Citation(s) in RCA: 33] [Impact Index Per Article: 5.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/28/2018] [Accepted: 10/16/2018] [Indexed: 01/25/2023] Open
Abstract
Phenotypic complexity is caused by the contributions of environmental factors and multiple genetic loci, interacting or acting independently. Studies of yeast and Arabidopsis often find that the majority of natural variation across phenotypes is attributable to independent additive quantitative trait loci (QTL). Detected loci in these organisms explain most of the estimated heritable variation. By contrast, many heritable components underlying phenotypic variation in metazoan models remain undetected. Before the relative impacts of additive and interactive variance components on metazoan phenotypic variation can be dissected, high replication and precise phenotypic measurements are required to obtain sufficient statistical power to detect loci contributing to this missing heritability. Here, we used a panel of 296 recombinant inbred advanced intercross lines of Caenorhabditis elegans and a high-throughput fitness assay to detect loci underlying responses to 16 different toxins, including heavy metals, chemotherapeutic drugs, pesticides, and neuropharmaceuticals. Using linkage mapping, we identified 82 QTL that underlie variation in responses to these toxins, and predicted the relative contributions of additive loci and genetic interactions across various growth parameters. Additionally, we identified three genomic regions that impact responses to multiple classes of toxins. These QTL hotspots could represent common factors impacting toxin responses. We went further to generate near-isogenic lines and chromosome substitution strains, and then experimentally validated these QTL hotspots, implicating additive and interactive loci that underlie toxin-response variation.
Collapse
|
31
|
Gao AW, Sterken MG, Uit de Bos J, van Creij J, Kamble R, Snoek BL, Kammenga JE, Houtkooper RH. Natural genetic variation in C. elegans identified genomic loci controlling metabolite levels. Genome Res 2018; 28:1296-1308. [PMID: 30108180 PMCID: PMC6120624 DOI: 10.1101/gr.232322.117] [Citation(s) in RCA: 20] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/10/2017] [Accepted: 07/20/2018] [Indexed: 12/31/2022]
Abstract
Metabolic homeostasis is sustained by complex biological networks that respond to nutrient availability. Genetic and environmental factors may disrupt this equilibrium, leading to metabolic disorders, including obesity and type 2 diabetes. To identify the genetic factors controlling metabolism, we performed quantitative genetic analysis using a population of 199 recombinant inbred lines (RILs) in the nematode Caenorhabditis elegans We focused on the genomic regions that control metabolite levels by measuring fatty acid (FA) and amino acid (AA) composition in the RILs using targeted metabolomics. The genetically diverse RILs showed a large variation in their FA and AA levels with a heritability ranging from 32% to 82%. We detected strongly co-correlated metabolite clusters and 36 significant metabolite quantitative trait loci (mQTL). We focused on mQTL displaying highly significant linkage and heritability, including an mQTL for the FA C14:1 on Chromosome I, and another mQTL for the FA C18:2 on Chromosome IV. Using introgression lines (ILs), we were able to narrow down both mQTL to a 1.4-Mbp and a 3.6-Mbp region, respectively. RNAi-based screening focusing on the Chromosome I mQTL identified several candidate genes for the C14:1 mQTL, including lagr-1, Y87G2A.2, nhr-265, nhr-276, and nhr-81 Overall, this systems approach provides us with a powerful platform to study the genetic basis of C. elegans metabolism. Furthermore, it allows us to investigate interventions such as nutrients and stresses that maintain or disturb the regulatory network controlling metabolic homeostasis, and identify gene-by-environment interactions.
Collapse
Affiliation(s)
- Arwen W Gao
- Laboratory Genetic Metabolic Diseases, Amsterdam UMC, University of Amsterdam, Amsterdam Gastroenterology and Metabolism, 1105 AZ Amsterdam, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands
| | - Jelmi Uit de Bos
- Laboratory Genetic Metabolic Diseases, Amsterdam UMC, University of Amsterdam, Amsterdam Gastroenterology and Metabolism, 1105 AZ Amsterdam, The Netherlands
| | - Jelle van Creij
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands
| | - Rashmi Kamble
- Laboratory Genetic Metabolic Diseases, Amsterdam UMC, University of Amsterdam, Amsterdam Gastroenterology and Metabolism, 1105 AZ Amsterdam, The Netherlands
| | - Basten L Snoek
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University and Research, 6708 PB, Wageningen, The Netherlands
| | - Riekelt H Houtkooper
- Laboratory Genetic Metabolic Diseases, Amsterdam UMC, University of Amsterdam, Amsterdam Gastroenterology and Metabolism, 1105 AZ Amsterdam, The Netherlands
| |
Collapse
|
32
|
Viñuela A, Brown AA, Buil A, Tsai PC, Davies MN, Bell JT, Dermitzakis ET, Spector TD, Small KS. Age-dependent changes in mean and variance of gene expression across tissues in a twin cohort. Hum Mol Genet 2018; 27:732-741. [PMID: 29228364 PMCID: PMC5886097 DOI: 10.1093/hmg/ddx424] [Citation(s) in RCA: 60] [Impact Index Per Article: 10.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/31/2017] [Revised: 11/10/2017] [Accepted: 11/29/2017] [Indexed: 12/13/2022] Open
Abstract
Changes in the mean and variance of gene expression with age have consequences for healthy aging and disease development. Age-dependent changes in phenotypic variance have been associated with a decline in regulatory functions leading to increase in disease risk. Here, we investigate age-related mean and variance changes in gene expression measured by RNA-seq of fat, skin, whole blood and derived lymphoblastoid cell lines (LCLs) expression from 855 adult female twins. We see evidence of up to 60% of age effects on transcription levels shared across tissues, and 47% of those on splicing. Using gene expression variance and discordance between genetically identical MZ twin pairs, we identify 137 genes with age-related changes in variance and 42 genes with age-related discordance between co-twins; implying the latter are driven by environmental effects. We identify four eQTLs whose effect on expression is age-dependent (FDR 5%). Combined, these results show a complicated mix of environmental and genetically driven changes in expression with age. Using the twin structure in our data, we show that additive genetic effects explain considerably more of the variance in gene expression than aging, but less that other environmental factors, potentially explaining why reliable expression-derived biomarkers for healthy-aging have proved elusive compared with those derived from methylation.
Collapse
Affiliation(s)
- Ana Viñuela
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
- Department of Genetic Medicine and Development, University of Geneva Medical School, 1211 Geneva, Switzerland
- Institute for Genetics and Genomics in Geneva (iGE3), University of Geneva, 1211 Geneva, Switzerland
- Swiss Institute of Bioinformatics, 1211 Geneva, Switzerland
| | - Andrew A Brown
- Wellcome Trust Sanger Institute, Hinxton CB10 1SA, Cambridge, UK
- Division of Mental Health and Addiction, NORMENT, KG Jebsen Centre for Psychosis Research, Oslo University Hospital, Oslo 0450, Norway
| | - Alfonso Buil
- Department of Genetic Medicine and Development, University of Geneva Medical School, 1211 Geneva, Switzerland
- Institute for Genetics and Genomics in Geneva (iGE3), University of Geneva, 1211 Geneva, Switzerland
- Swiss Institute of Bioinformatics, 1211 Geneva, Switzerland
| | - Pei-Chien Tsai
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
| | - Matthew N Davies
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
| | - Jordana T Bell
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
| | - Emmanouil T Dermitzakis
- Department of Genetic Medicine and Development, University of Geneva Medical School, 1211 Geneva, Switzerland
- Institute for Genetics and Genomics in Geneva (iGE3), University of Geneva, 1211 Geneva, Switzerland
- Swiss Institute of Bioinformatics, 1211 Geneva, Switzerland
| | - Timothy D Spector
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
| | - Kerrin S Small
- Department of Twin Research and Genetic Epidemiology, King's College London, St Thomas' Campus, SE1 7EH London, UK
| |
Collapse
|
33
|
Distinctive roles of age, sex, and genetics in shaping transcriptional variation of human immune responses to microbial challenges. Proc Natl Acad Sci U S A 2017; 115:E488-E497. [PMID: 29282317 PMCID: PMC5776984 DOI: 10.1073/pnas.1714765115] [Citation(s) in RCA: 135] [Impact Index Per Article: 19.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/14/2022] Open
Abstract
Identifying the drivers of the interindividual diversity of the human immune system is crucial to understand their consequences on immune-mediated diseases. By examining the transcriptional responses of 1,000 individuals to various microbial challenges, we show that age and sex influence the expression of many immune-related genes, but their effects are overall moderate, whereas genetic factors affect a smaller gene set but with a stronger effect. We identify numerous genetic variants that affect transcriptional variation on infection, many of which are associated with autoimmune or inflammatory disorders. These results enable additional exploration of the role of regulatory variants in the pathogenesis of immune-related diseases and improve our understanding of the respective effects of age, sex, and genetics on immune response variation. The contribution of host genetic and nongenetic factors to immunological differences in humans remains largely undefined. Here, we generated bacterial-, fungal-, and viral-induced immune transcriptional profiles in an age- and sex-balanced cohort of 1,000 healthy individuals and searched for the determinants of immune response variation. We found that age and sex affected the transcriptional response of most immune-related genes, with age effects being more stimulus-specific relative to sex effects, which were largely shared across conditions. Although specific cell populations mediated the effects of age and sex on gene expression, including CD8+ T cells for age and CD4+ T cells and monocytes for sex, we detected a direct effect of these intrinsic factors for the majority of immune genes. The mapping of expression quantitative trait loci (eQTLs) revealed that genetic factors had a stronger effect on immune gene regulation than age and sex, yet they affected a smaller number of genes. Importantly, we identified numerous genetic variants that manifested their regulatory effects exclusively on immune stimulation, including a Candida albicans-specific master regulator at the CR1 locus. These response eQTLs were enriched in disease-associated variants, particularly for autoimmune and inflammatory disorders, indicating that differences in disease risk may result from regulatory variants exerting their effects only in the presence of immune stress. Together, this study quantifies the respective effects of age, sex, genetics, and cellular heterogeneity on the interindividual variability of immune responses and constitutes a valuable resource for further exploration in the context of different infection risks or disease outcomes.
Collapse
|
34
|
Serin EAR, Snoek LB, Nijveen H, Willems LAJ, Jiménez-Gómez JM, Hilhorst HWM, Ligterink W. Construction of a High-Density Genetic Map from RNA-Seq Data for an Arabidopsis Bay-0 × Shahdara RIL Population. Front Genet 2017; 8:201. [PMID: 29259624 PMCID: PMC5723289 DOI: 10.3389/fgene.2017.00201] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2017] [Accepted: 11/21/2017] [Indexed: 12/17/2022] Open
Abstract
High-density genetic maps are essential for high resolution mapping of quantitative traits. Here, we present a new genetic map for an Arabidopsis Bayreuth × Shahdara recombinant inbred line (RIL) population, built on RNA-seq data. RNA-seq analysis on 160 RILs of this population identified 30,049 single-nucleotide polymorphisms (SNPs) covering the whole genome. Based on a 100-kbp window SNP binning method, 1059 bin-markers were identified, physically anchored on the genome. The total length of the RNA-seq genetic map spans 471.70 centimorgans (cM) with an average marker distance of 0.45 cM and a maximum marker distance of 4.81 cM. This high resolution genotyping revealed new recombination breakpoints in the population. To highlight the advantages of such high-density map, we compared it to two publicly available genetic maps for the same population, comprising 69 PCR-based markers and 497 gene expression markers derived from microarray data, respectively. In this study, we show that SNP markers can effectively be derived from RNA-seq data. The new RNA-seq map closes many existing gaps in marker coverage, saturating the previously available genetic maps. Quantitative trait locus (QTL) analysis for published phenotypes using the available genetic maps showed increased QTL mapping resolution and reduced QTL confidence interval using the RNA-seq map. The new high-density map is a valuable resource that facilitates the identification of candidate genes and map-based cloning approaches.
Collapse
Affiliation(s)
- Elise A R Serin
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, Wageningen, Netherlands
| | - L B Snoek
- Laboratory of Nematology, Wageningen University, Wageningen, Netherlands.,Theoretical Biology and Bioinformatics, Utrecht University, Utrecht, Netherlands
| | - Harm Nijveen
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, Wageningen, Netherlands.,Laboratory of Bioinformatics, Wageningen University, Wageningen, Netherlands
| | - Leo A J Willems
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, Wageningen, Netherlands
| | - Jose M Jiménez-Gómez
- Department of Plant Breeding and Genetics, Max Planck Institute for Plant Breeding Research, Cologne, Germany.,Institut Jean-Pierre Bourgin, Institut National de la Recherche Agronomique, AgroParisTech, Centre National de la Recherche Scientifique, Université Paris-Saclay, Versailles Cedex, France
| | - Henk W M Hilhorst
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, Wageningen, Netherlands
| | - Wilco Ligterink
- Wageningen Seed Lab, Laboratory of Plant Physiology, Wageningen University, Wageningen, Netherlands
| |
Collapse
|
35
|
Gosik K, Kong L, Chinchilli VM, Wu R. iFORM/eQTL: an ultrahigh-dimensional platform for inferring the global genetic architecture of gene transcripts. Brief Bioinform 2017; 18:250-259. [PMID: 26944084 DOI: 10.1093/bib/bbw014] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/13/2015] [Indexed: 01/03/2023] Open
Abstract
Knowledge about how changes in gene expression are encoded by expression quantitative trait loci (eQTLs) is a key to construct the genotype-phenotype map for complex traits or diseases. Traditional eQTL mapping is to associate one transcript with a single marker at a time, thereby limiting our inference about a complete picture of the genetic architecture of gene expression. Here, we implemented an ultrahigh-dimensional variable selection model to build a computing platform that can systematically scan main effects and interaction effects among all possible loci and identify a set of significant eQTLs modulating differentiation and function of gene expression. This platform, named iFORM/eQTL, was assembled by forward-selection-based procedures to tackle complex covariance structures of gene-gene interactions. iFORM/eQTL can particularly discern the role of cis-QTLs, trans-QTLs and their epistatic interactions in gene expression. Results from the reanalysis of a published genetic and genomic data set through iFORM/eQTL gain new discoveries on the genetic origin of gene expression differentiation in Caenorhabditis elegans, which could not be detected by a traditional one-locus/one-transcript analysis approach.
Collapse
Affiliation(s)
- Kirk Gosik
- Department of Statistics, The Pennsylvania State University, University Park, PA, USA
| | - Lan Kong
- Department of Public Health Sciences Pennsylvania State University College of Medicine, Hershey, PA, USA
| | - Vernon M Chinchilli
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA
| | - Rongling Wu
- Division of Biostatistics and Bioinformatics, Department of Public Health Sciences, Penn State College of Medicine, Hershey, PA, USA
| |
Collapse
|
36
|
Gao AW, Uit de Bos J, Sterken MG, Kammenga JE, Smith RL, Houtkooper RH. Forward and reverse genetics approaches to uncover metabolic aging pathways in Caenorhabditis elegans. Biochim Biophys Acta Mol Basis Dis 2017; 1864:2697-2706. [PMID: 28919364 DOI: 10.1016/j.bbadis.2017.09.006] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2017] [Revised: 09/05/2017] [Accepted: 09/07/2017] [Indexed: 01/08/2023]
Abstract
The biological mechanisms of aging have been studied in depth and prominent findings in this field promote the development of new therapies for age-associated disorders. Various model organisms are used for research on aging; among these, the nematode Caenorhabditis elegans has been widely used and has provided valuable knowledge in determining the regulatory mechanisms driving the aging process. Many genes involved in lifespan regulation are associated with metabolic pathways and are influenced by genetic and environmental factors. In line with this, C. elegans provides a promising platform to study such gene by environment interactions, in either a reverse or forward genetics approach. In this review, we discuss longevity mechanisms related to metabolic networks that have been discovered in C. elegans. We also highlight the use of wild populations to study the complex genetic basis of natural variation for quantitative traits that mediate longevity.
Collapse
Affiliation(s)
- Arwen W Gao
- Laboratory Genetic Metabolic Diseases, Academic Medical Center of Amsterdam, 1105 AZ Amsterdam, The Netherlands
| | - Jelmi Uit de Bos
- Laboratory Genetic Metabolic Diseases, Academic Medical Center of Amsterdam, 1105 AZ Amsterdam, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University and Research, 6708 PB Wageningen, The Netherlands
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University and Research, 6708 PB Wageningen, The Netherlands
| | - Reuben L Smith
- Laboratory Genetic Metabolic Diseases, Academic Medical Center of Amsterdam, 1105 AZ Amsterdam, The Netherlands
| | - Riekelt H Houtkooper
- Laboratory Genetic Metabolic Diseases, Academic Medical Center of Amsterdam, 1105 AZ Amsterdam, The Netherlands.
| |
Collapse
|
37
|
Abstract
The oncogenic Ras/MAPK pathway is evolutionarily conserved across metazoans. Yet, almost all our knowledge on this pathway comes from studies using single genetic backgrounds, whereas mutational effects can be highly background dependent. Therefore, we lack insight in the interplay between genetic backgrounds and the Ras/MAPK-signaling pathway. Here, we used a Caenorhabditis elegans RIL population containing a gain-of-function mutation in the Ras/MAPK-pathway gene let-60 and measured how gene expression regulation is affected by this mutation. We mapped eQTL and found that the majority (∼73%) of the 1516 detected cis-eQTL were not specific for the let-60 mutation, whereas most (∼76%) of the 898 detected trans-eQTL were associated with the let-60 mutation. We detected six eQTL trans-bands specific for the interaction between the genetic background and the mutation, one of which colocalized with the polymorphic Ras/MAPK modifier amx-2. Comparison between transgenic lines expressing allelic variants of amx-2 showed the involvement of amx-2 in 79% of the trans-eQTL for genes mapping to this trans-band. Together, our results have revealed hidden loci affecting Ras/MAPK signaling using sensitized backgrounds in C. elegans. These loci harbor putative polymorphic modifier genes that would not have been detected using mutant screens in single genetic backgrounds.
Collapse
|
38
|
Westhues M, Schrag TA, Heuer C, Thaller G, Utz HF, Schipprack W, Thiemann A, Seifert F, Ehret A, Schlereth A, Stitt M, Nikoloski Z, Willmitzer L, Schön CC, Scholten S, Melchinger AE. Omics-based hybrid prediction in maize. TAG. THEORETICAL AND APPLIED GENETICS. THEORETISCHE UND ANGEWANDTE GENETIK 2017. [PMID: 28647896 DOI: 10.1007/s00122-017-2934-0] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 05/02/2023]
Abstract
Complementing genomic data with other "omics" predictors can increase the probability of success for predicting the best hybrid combinations using complex agronomic traits. Accurate prediction of traits with complex genetic architecture is crucial for selecting superior candidates in animal and plant breeding and for guiding decisions in personalized medicine. Whole-genome prediction has revolutionized these areas but has inherent limitations in incorporating intricate epistatic interactions. Downstream "omics" data are expected to integrate interactions within and between different biological strata and provide the opportunity to improve trait prediction. Yet, predicting traits from parents to progeny has not been addressed by a combination of "omics" data. Here, we evaluate several "omics" predictors-genomic, transcriptomic and metabolic data-measured on parent lines at early developmental stages and demonstrate that the integration of transcriptomic with genomic data leads to higher success rates in the correct prediction of untested hybrid combinations in maize. Despite the high predictive ability of genomic data, transcriptomic data alone outperformed them and other predictors for the most complex heterotic trait, dry matter yield. An eQTL analysis revealed that transcriptomic data integrate genomic information from both, adjacent and distant sites relative to the expressed genes. Together, these findings suggest that downstream predictors capture physiological epistasis that is transmitted from parents to their hybrid offspring. We conclude that the use of downstream "omics" data in prediction can exploit important information beyond structural genomics for leveraging the efficiency of hybrid breeding.
Collapse
Affiliation(s)
- Matthias Westhues
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599, Stuttgart, Germany
| | - Tobias A Schrag
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599, Stuttgart, Germany
| | - Claas Heuer
- Institute of Animal Breeding and Husbandry, Christian-Albrechts-University Kiel, 24098, Kiel, Germany
- Inguran LLC dba STGenetics, 22575 SH6 South, Navasota, TX, 77868, USA
| | - Georg Thaller
- Institute of Animal Breeding and Husbandry, Christian-Albrechts-University Kiel, 24098, Kiel, Germany
| | - H Friedrich Utz
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599, Stuttgart, Germany
| | - Wolfgang Schipprack
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599, Stuttgart, Germany
| | - Alexander Thiemann
- Biocenter Klein Flottbek, Developmental Biology and Biotechnology, University of Hamburg, 22609, Hamburg, Germany
| | - Felix Seifert
- Biocenter Klein Flottbek, Developmental Biology and Biotechnology, University of Hamburg, 22609, Hamburg, Germany
| | - Anita Ehret
- Institute of Animal Breeding and Husbandry, Christian-Albrechts-University Kiel, 24098, Kiel, Germany
| | - Armin Schlereth
- Max-Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Mark Stitt
- Max-Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Zoran Nikoloski
- Max-Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Lothar Willmitzer
- Max-Planck Institute of Molecular Plant Physiology, 14476, Potsdam, Germany
| | - Chris C Schön
- Plant Breeding, Technische Universität München, 85354, Freising, Germany
| | - Stefan Scholten
- Biocenter Klein Flottbek, Developmental Biology and Biotechnology, University of Hamburg, 22609, Hamburg, Germany.
| | - Albrecht E Melchinger
- Institute of Plant Breeding, Seed Science and Population Genetics, University of Hohenheim, 70599, Stuttgart, Germany.
| |
Collapse
|
39
|
Snoek BL, Sterken MG, Bevers RPJ, Volkers RJM, Van't Hof A, Brenchley R, Riksen JAG, Cossins A, Kammenga JE. Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans. BMC Genomics 2017; 18:500. [PMID: 28662696 PMCID: PMC5492678 DOI: 10.1186/s12864-017-3899-8] [Citation(s) in RCA: 54] [Impact Index Per Article: 7.7] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/29/2017] [Accepted: 06/22/2017] [Indexed: 11/10/2022] Open
Abstract
Background Cryptic genetic variation (CGV) is the hidden genetic variation that can be unlocked by perturbing normal conditions. CGV can drive the emergence of novel complex phenotypes through changes in gene expression. Although our theoretical understanding of CGV has thoroughly increased over the past decade, insight into polymorphic gene expression regulation underlying CGV is scarce. Here we investigated the transcriptional architecture of CGV in response to rapid temperature changes in the nematode Caenorhabditis elegans. We analyzed regulatory variation in gene expression (and mapped eQTL) across the course of a heat stress and recovery response in a recombinant inbred population. Results We measured gene expression over three temperature treatments: i) control, ii) heat stress, and iii) recovery from heat stress. Compared to control, exposure to heat stress affected the transcription of 3305 genes, whereas 942 were affected in recovering animals. These affected genes were mainly involved in metabolism and reproduction. The gene expression pattern in recovering animals resembled both the control and the heat-stress treatment. We mapped eQTL using the genetic variation of the recombinant inbred population and detected 2626 genes with an eQTL in the heat-stress treatment, 1797 in the control, and 1880 in the recovery. The cis-eQTL were highly shared across treatments. A considerable fraction of the trans-eQTL (40–57%) mapped to 19 treatment specific trans-bands. In contrast to cis-eQTL, trans-eQTL were highly environment specific and thus cryptic. Approximately 67% of the trans-eQTL were only induced in a single treatment, with heat-stress showing the most unique trans-eQTL. Conclusions These results illustrate the highly dynamic pattern of CGV across three different environmental conditions that can be evoked by a stress response over a relatively short time-span (2 h) and that CGV is mainly determined by response related trans regulatory eQTL. Electronic supplementary material The online version of this article (doi:10.1186/s12864-017-3899-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Basten L Snoek
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands
| | - Mark G Sterken
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands
| | - Roel P J Bevers
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands
| | - Rita J M Volkers
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands
| | - Arjen Van't Hof
- Centre for Genome research, Institute of Integrative Biology, Biosciences Building, University of Liverpool, L69 7ZB, Liverpool, UK
| | - Rachel Brenchley
- Centre for Genome research, Institute of Integrative Biology, Biosciences Building, University of Liverpool, L69 7ZB, Liverpool, UK
| | - Joost A G Riksen
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands
| | - Andrew Cossins
- Centre for Genome research, Institute of Integrative Biology, Biosciences Building, University of Liverpool, L69 7ZB, Liverpool, UK
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University and Research, Droevendaalsesteeg 1, 6708PB, Wageningen, The Netherlands.
| |
Collapse
|
40
|
Variation in gene expression within clones of the earthworm Dendrobaena octaedra. PLoS One 2017; 12:e0174960. [PMID: 28384196 PMCID: PMC5383104 DOI: 10.1371/journal.pone.0174960] [Citation(s) in RCA: 6] [Impact Index Per Article: 0.9] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/19/2016] [Accepted: 03/19/2017] [Indexed: 12/23/2022] Open
Abstract
Gene expression is highly plastic, which can help organisms to both acclimate and adapt to changing environments. Possible variation in gene expression among individuals with the same genotype (among clones) is not widely considered, even though it could impact the results of studies that focus on gene expression phenotypes, for example studies using clonal lines. We examined the extent of within and between clone variation in gene expression in the earthworm Dendrobaena octaedra, which reproduces through apomictic parthenogenesis. Five microsatellite markers were developed and used to confirm that offspring are genetic clones of their parent. After that, expression of 12 genes was measured from five individuals each from six clonal lines after exposure to copper contaminated soil. Variation in gene expression was higher over all genotypes than within genotypes, as initially assumed. A subset of the genes was also examined in the offspring of exposed individuals in two of the clonal lines. In this case, variation in gene expression within genotypes was as high as that observed over all genotypes. One gene in particular (chymotrypsin inhibitor) also showed significant differences in the expression levels among genetically identical individuals. Gene expression can vary considerably, and the extent of variation may depend on the genotypes and genes studied. Ensuring a large sample, with many different genotypes, is critical in studies comparing gene expression phenotypes. Researchers should be especially cautious inferring gene expression phenotypes when using only a single clonal or inbred line, since the results might be specific to only certain genotypes.
Collapse
|
41
|
Time-dependent genetic effects on gene expression implicate aging processes. Genome Res 2017; 27:545-552. [PMID: 28302734 PMCID: PMC5378173 DOI: 10.1101/gr.207688.116] [Citation(s) in RCA: 25] [Impact Index Per Article: 3.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/30/2016] [Accepted: 01/23/2017] [Indexed: 01/04/2023]
Abstract
Gene expression is dependent on genetic and environmental factors. In the last decade, a large body of research has significantly improved our understanding of the genetic architecture of gene expression. However, it remains unclear whether genetic effects on gene expression remain stable over time. Here, we show, using longitudinal whole-blood gene expression data from a twin cohort, that the genetic architecture of a subset of genes is unstable over time. In addition, we identified 2213 genes differentially expressed across time points that we linked with aging within and across studies. Interestingly, we discovered that most differentially expressed genes were affected by a subset of 77 putative causal genes. Finally, we observed that putative causal genes and down-regulated genes were affected by a loss of genetic control between time points. Taken together, our data suggest that instability in the genetic architecture of a subset of genes could lead to widespread effects on the transcriptome with an aging signature.
Collapse
|
42
|
Zych K, Snoek BL, Elvin M, Rodriguez M, Van der Velde KJ, Arends D, Westra HJ, Swertz MA, Poulin G, Kammenga JE, Breitling R, Jansen RC, Li Y. reGenotyper: Detecting mislabeled samples in genetic data. PLoS One 2017; 12:e0171324. [PMID: 28192439 PMCID: PMC5305221 DOI: 10.1371/journal.pone.0171324] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/22/2016] [Accepted: 01/19/2017] [Indexed: 12/11/2022] Open
Abstract
In high-throughput molecular profiling studies, genotype labels can be wrongly assigned at various experimental steps; the resulting mislabeled samples seriously reduce the power to detect the genetic basis of phenotypic variation. We have developed an approach to detect potential mislabeling, recover the “ideal” genotype and identify “best-matched” labels for mislabeled samples. On average, we identified 4% of samples as mislabeled in eight published datasets, highlighting the necessity of applying a “data cleaning” step before standard data analysis.
Collapse
Affiliation(s)
- Konrad Zych
- Groningen Bioinformatics Centre, University of Groningen, Groningen, The Netherlands
| | - Basten L. Snoek
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Mark Elvin
- Faculty of Life Sciences, University of Manchester, Manchester, United Kingdom
| | - Miriam Rodriguez
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - K. Joeri Van der Velde
- Genomics Coordination Center, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
| | - Danny Arends
- Groningen Bioinformatics Centre, University of Groningen, Groningen, The Netherlands
| | - Harm-Jan Westra
- Divisions of Genetics and Rheumatology, Department of Medicine, Brigham and Women’s Hospital and Harvard Medical School, Boston, Massachusetts, United States of America
- Partners Center for Personalized Genetic Medicine, Boston, Massachusetts, United States of America
- Program in Medical and Population Genetics, Broad Institute of MIT and Harvard, Cambridge, Massachusetts, United States of America
| | - Morris A. Swertz
- Genomics Coordination Center, University of Groningen, University Medical Center Groningen, Groningen, The Netherlands
| | - Gino Poulin
- Faculty of Life Sciences, University of Manchester, Manchester, United Kingdom
| | - Jan E. Kammenga
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Rainer Breitling
- Manchester Institute of Biotechnology, Faculty of Life Sciences, University of Manchester, Manchester, United Kingdom
| | - Ritsert C. Jansen
- Groningen Bioinformatics Centre, University of Groningen, Groningen, The Netherlands
| | - Yang Li
- Groningen Bioinformatics Centre, University of Groningen, Groningen, The Netherlands
- University of Groningen, University Medical Center Groningen, Department of Genetics, Groningen, The Netherlands
- * E-mail:
| |
Collapse
|
43
|
Woo HR, Koo HJ, Kim J, Jeong H, Yang JO, Lee IH, Jun JH, Choi SH, Park SJ, Kang B, Kim YW, Phee BK, Kim JH, Seo C, Park C, Kim SC, Park S, Lee B, Lee S, Hwang D, Nam HG, Lim PO. Programming of Plant Leaf Senescence with Temporal and Inter-Organellar Coordination of Transcriptome in Arabidopsis. PLANT PHYSIOLOGY 2016; 171:452-67. [PMID: 26966169 PMCID: PMC4854694 DOI: 10.1104/pp.15.01929] [Citation(s) in RCA: 92] [Impact Index Per Article: 11.5] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Subscribe] [Scholar Register] [Received: 12/11/2015] [Accepted: 03/07/2016] [Indexed: 05/20/2023]
Abstract
Plant leaves, harvesting light energy and fixing CO2, are a major source of foods on the earth. Leaves undergo developmental and physiological shifts during their lifespan, ending with senescence and death. We characterized the key regulatory features of the leaf transcriptome during aging by analyzing total- and small-RNA transcriptomes throughout the lifespan of Arabidopsis (Arabidopsis thaliana) leaves at multidimensions, including age, RNA-type, and organelle. Intriguingly, senescing leaves showed more coordinated temporal changes in transcriptomes than growing leaves, with sophisticated regulatory networks comprising transcription factors and diverse small regulatory RNAs. The chloroplast transcriptome, but not the mitochondrial transcriptome, showed major changes during leaf aging, with a strongly shared expression pattern of nuclear transcripts encoding chloroplast-targeted proteins. Thus, unlike animal aging, leaf senescence proceeds with tight temporal and distinct interorganellar coordination of various transcriptomes that would be critical for the highly regulated degeneration and nutrient recycling contributing to plant fitness and productivity.
Collapse
Affiliation(s)
- Hye Ryun Woo
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Hee Jung Koo
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Jeongsik Kim
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Hyobin Jeong
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Jin Ok Yang
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Il Hwan Lee
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Ji Hyung Jun
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Seung Hee Choi
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Su Jin Park
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Byeongsoo Kang
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - You Wang Kim
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Bong-Kwan Phee
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Jin Hee Kim
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Chaehwa Seo
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Charny Park
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Sang Cheol Kim
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Seongjin Park
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Byungwook Lee
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Sanghyuk Lee
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Daehee Hwang
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Hong Gil Nam
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| | - Pyung Ok Lim
- Department of New Biology, DGIST, Daegu, Republic of Korea (H.R.W., D.H., H.G.N., P.O.L.);Center for Plant Aging Research, Institute for Basic Science (IBS), Daegu, Republic of Korea (H.J.K., J.K., H.J., I.H.L., S.H.C., S.J.P., B.K., Y.W.K., B.-K.P., J.H.K., D.H., H.G.N.);School of Interdisciplinary Bioscience and Bioengineering, POSTECH, Pohang, Republic of Korea (H.J.K.);Korean Bioinformation Center, KRIBB, Daejeon, Republic of Korea (J.O.Y., S.C.K., S.P., B.L.);Division of Molecular Life Sciences, POSTECH, Pohang, Republic of Korea (I.H.L., J.H.J., S.H.C.);Division of Integrative Biosciences and Biotechnologies, POSTECH, Pohang, Republic of Korea (S.J.P.);DNA Link Inc., Seoul, Republic of Korea (C.S.); andEwha Research Center for Systems Biology and Department of Life Science, Ewha Womans University, Seoul, Republic of Korea (C.P., S.L.)
| |
Collapse
|
44
|
Nakad R, Snoek LB, Yang W, Ellendt S, Schneider F, Mohr TG, Rösingh L, Masche AC, Rosenstiel PC, Dierking K, Kammenga JE, Schulenburg H. Contrasting invertebrate immune defense behaviors caused by a single gene, the Caenorhabditis elegans neuropeptide receptor gene npr-1. BMC Genomics 2016; 17:280. [PMID: 27066825 PMCID: PMC4827197 DOI: 10.1186/s12864-016-2603-8] [Citation(s) in RCA: 34] [Impact Index Per Article: 4.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 11/07/2015] [Accepted: 03/25/2016] [Indexed: 01/22/2023] Open
Abstract
Background The invertebrate immune system comprises physiological mechanisms, physical barriers and also behavioral responses. It is generally related to the vertebrate innate immune system and widely believed to provide nonspecific defense against pathogens, whereby the response to different pathogen types is usually mediated by distinct signalling cascades. Recent work suggests that invertebrate immune defense can be more specific at least at the phenotypic level. The underlying genetic mechanisms are as yet poorly understood. Results We demonstrate in the model invertebrate Caenorhabditis elegans that a single gene, a homolog of the mammalian neuropeptide Y receptor gene, npr-1, mediates contrasting defense phenotypes towards two distinct pathogens, the Gram-positive Bacillus thuringiensis and the Gram-negative Pseudomonas aeruginosa. Our findings are based on combining quantitative trait loci (QTLs) analysis with functional genetic analysis and RNAseq-based transcriptomics. The QTL analysis focused on behavioral immune defense against B. thuringiensis, using recombinant inbred lines (RILs) and introgression lines (ILs). It revealed several defense QTLs, including one on chromosome X comprising the npr-1 gene. The wildtype N2 allele for the latter QTL was associated with reduced defense against B. thuringiensis and thus produced an opposite phenotype to that previously reported for the N2 npr-1 allele against P. aeruginosa. Analysis of npr-1 mutants confirmed these contrasting immune phenotypes for both avoidance behavior and nematode survival. Subsequent transcriptional profiling of C. elegans wildtype and npr-1 mutant suggested that npr-1 mediates defense against both pathogens through p38 MAPK signaling, insulin-like signaling, and C-type lectins. Importantly, increased defense towards P. aeruginosa seems to be additionally influenced through the induction of oxidative stress genes and activation of GATA transcription factors, while the repression of oxidative stress genes combined with activation of Ebox transcription factors appears to enhance susceptibility to B. thuringiensis. Conclusions Our findings highlight the role of a single gene, npr-1, in fine-tuning nematode immune defense, showing the ability of the invertebrate immune system to produce highly specialized and potentially opposing immune responses via single regulatory genes. Electronic supplementary material The online version of this article (doi:10.1186/s12864-016-2603-8) contains supplementary material, which is available to authorized users.
Collapse
Affiliation(s)
- Rania Nakad
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany.,Cologne Excellence Cluster for Cellular Stress Responses in Ageing-Associated Diseases (CECAD) and Systems Biology of Ageing, University of Cologne, Joseph-Stelzmann-Str. 26, 50931, Cologne, Germany
| | - L Basten Snoek
- Laboratory of Nematology, Wageningen University, Wageningen, 6708 PB, The Netherlands
| | - Wentao Yang
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Sunna Ellendt
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Franziska Schneider
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Timm G Mohr
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Lone Rösingh
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Anna C Masche
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Philip C Rosenstiel
- Institute for Clinical Molecular Biology, University of Kiel, 24098, Kiel, Germany
| | - Katja Dierking
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany
| | - Jan E Kammenga
- Laboratory of Nematology, Wageningen University, Wageningen, 6708 PB, The Netherlands
| | - Hinrich Schulenburg
- Department of Evolutionary Ecology and Genetics, Zoological Institute, University of Kiel, 24098, Kiel, Germany.
| |
Collapse
|
45
|
Singh KD, Roschitzki B, Snoek LB, Grossmann J, Zheng X, Elvin M, Kamkina P, Schrimpf SP, Poulin GB, Kammenga JE, Hengartner MO. Natural Genetic Variation Influences Protein Abundances in C. elegans Developmental Signalling Pathways. PLoS One 2016; 11:e0149418. [PMID: 26985669 PMCID: PMC4795773 DOI: 10.1371/journal.pone.0149418] [Citation(s) in RCA: 26] [Impact Index Per Article: 3.3] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/18/2015] [Accepted: 01/30/2016] [Indexed: 12/11/2022] Open
Abstract
Complex traits, including common disease-related traits, are affected by many different genes that function in multiple pathways and networks. The apoptosis, MAPK, Notch, and Wnt signalling pathways play important roles in development and disease progression. At the moment we have a poor understanding of how allelic variation affects gene expression in these pathways at the level of translation. Here we report the effect of natural genetic variation on transcript and protein abundance involved in developmental signalling pathways in Caenorhabditis elegans. We used selected reaction monitoring to analyse proteins from the abovementioned four pathways in a set of recombinant inbred lines (RILs) generated from the wild-type strains N2 (Bristol) and CB4856 (Hawaii) to enable quantitative trait locus (QTL) mapping. About half of the cases from the 44 genes tested showed a statistically significant change in protein abundance between various strains, most of these were however very weak (below 1.3-fold change). We detected a distant QTL on the left arm of chromosome II that affected protein abundance of the phosphatidylserine receptor protein PSR-1, and two separate QTLs that influenced embryonic and ionizing radiation-induced apoptosis on chromosome IV. Our results demonstrate that natural variation in C. elegans is sufficient to cause significant changes in signalling pathways both at the gene expression (transcript and protein abundance) and phenotypic levels.
Collapse
Affiliation(s)
- Kapil Dev Singh
- Institute of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
| | - Bernd Roschitzki
- Functional Genomics Center Zurich, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - L. Basten Snoek
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Jonas Grossmann
- Functional Genomics Center Zurich, University of Zurich and ETH Zurich, Zurich, Switzerland
| | - Xue Zheng
- Institute of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
| | - Mark Elvin
- Faculty of Life Sciences, The University of Manchester, Manchester, United Kingdom
| | - Polina Kamkina
- Institute of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
| | - Sabine P. Schrimpf
- Institute of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
| | - Gino B. Poulin
- Faculty of Life Sciences, The University of Manchester, Manchester, United Kingdom
| | - Jan E. Kammenga
- Laboratory of Nematology, Wageningen University, Wageningen, The Netherlands
| | - Michael O. Hengartner
- Institute of Molecular Life Sciences, University of Zurich, Zurich, Switzerland
- * E-mail:
| |
Collapse
|
46
|
Kamkina P, Snoek LB, Grossmann J, Volkers RJM, Sterken MG, Daube M, Roschitzki B, Fortes C, Schlapbach R, Roth A, von Mering C, Hengartner MO, Schrimpf SP, Kammenga JE. Natural Genetic Variation Differentially Affects the Proteome and Transcriptome in Caenorhabditis elegans. Mol Cell Proteomics 2016; 15:1670-80. [PMID: 26944343 DOI: 10.1074/mcp.m115.052548] [Citation(s) in RCA: 16] [Impact Index Per Article: 2.0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/10/2015] [Indexed: 11/06/2022] Open
Abstract
Natural genetic variation is the raw material of evolution and influences disease development and progression. An important question is how this genetic variation translates into variation in protein abundance. To analyze the effects of the genetic background on gene and protein expression in the nematode Caenorhabditis elegans, we quantitatively compared the two genetically highly divergent wild-type strains N2 and CB4856. Gene expression was analyzed by microarray assays, and proteins were quantified using stable isotope labeling by amino acids in cell culture. Among all transcribed genes, we found 1,532 genes to be differentially transcribed between the two wild types. Of the total 3,238 quantified proteins, 129 proteins were significantly differentially expressed between N2 and CB4856. The differentially expressed proteins were enriched for genes that function in insulin-signaling and stress-response pathways, underlining strong divergence of these pathways in nematodes. The protein abundance of the two wild-type strains correlates more strongly than protein abundance versus transcript abundance within each wild type. Our findings indicate that in C. elegans only a fraction of the changes in protein abundance can be explained by the changes in mRNA abundance. These findings corroborate with the observations made across species.
Collapse
Affiliation(s)
- Polina Kamkina
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland; §Ph.D. Program in Molecular Life Sciences Zurich, 8057 Zurich, Switzerland
| | - L Basten Snoek
- ‖Laboratory of Nematology, Wageningen University, Wageningen 6708 PB, The Netherlands
| | - Jonas Grossmann
- **Functional Genomics Center Zurich, University of Zurich and Swiss Federal Institute of Technology Zurich, 8057 Zurich, Switzerland
| | - Rita J M Volkers
- ‖Laboratory of Nematology, Wageningen University, Wageningen 6708 PB, The Netherlands
| | - Mark G Sterken
- ‖Laboratory of Nematology, Wageningen University, Wageningen 6708 PB, The Netherlands
| | - Michael Daube
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland
| | - Bernd Roschitzki
- **Functional Genomics Center Zurich, University of Zurich and Swiss Federal Institute of Technology Zurich, 8057 Zurich, Switzerland
| | - Claudia Fortes
- **Functional Genomics Center Zurich, University of Zurich and Swiss Federal Institute of Technology Zurich, 8057 Zurich, Switzerland
| | - Ralph Schlapbach
- **Functional Genomics Center Zurich, University of Zurich and Swiss Federal Institute of Technology Zurich, 8057 Zurich, Switzerland
| | - Alexander Roth
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland
| | - Christian von Mering
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland
| | - Michael O Hengartner
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland
| | - Sabine P Schrimpf
- From the ‡Institute of Molecular Life Sciences, University of Zurich, 8057 Zurich, Switzerland;
| | - Jan E Kammenga
- ‖Laboratory of Nematology, Wageningen University, Wageningen 6708 PB, The Netherlands;
| |
Collapse
|
47
|
Wang N, Gosik K, Li R, Lindsay B, Wu R. A block mixture model to map eQTLs for gene clustering and networking. Sci Rep 2016; 6:21193. [PMID: 26892775 PMCID: PMC4759821 DOI: 10.1038/srep21193] [Citation(s) in RCA: 1] [Impact Index Per Article: 0.1] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/13/2015] [Accepted: 01/19/2016] [Indexed: 01/13/2023] Open
Abstract
To study how genes function in a cellular and physiological process, a general procedure is to classify gene expression profiles into categories based on their similarity and reconstruct a regulatory network for functional elements. However, this procedure has not been implemented with the genetic mechanisms that underlie the organization of gene clusters and networks, despite much effort made to map expression quantitative trait loci (eQTLs) that affect the expression of individual genes. Here we address this issue by developing a computational approach that integrates gene clustering and network reconstruction with genetic mapping into a unifying framework. The approach can not only identify specific eQTLs that control how genes are clustered and organized toward biological functions, but also enable the investigation of the biological mechanisms that individual eQTLs perturb in a signaling pathway. We applied the new approach to characterize the effects of eQTLs on the structure and organization of gene clusters in Caenorhabditis elegans. This study provides the first characterization, to our knowledge, of the effects of genetic variants on the regulatory network of gene expression. The approach developed can also facilitate the genetic dissection of other dynamic processes, including development, physiology and disease progression in any organisms.
Collapse
Affiliation(s)
- Ningtao Wang
- Department of Biostatistics, University of Texas School of Public Health, Houston, TX 77030, USA.,Department of Public Health Sciences, The Pennsylvania State University, Hershey, PA 17033, USA
| | - Kirk Gosik
- Department of Statistics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Runze Li
- Department of Biostatistics, University of Texas School of Public Health, Houston, TX 77030, USA.,Department of Statistics, The Pennsylvania State University, University Park, PA 16802, USA
| | - Bruce Lindsay
- Department of Biostatistics, University of Texas School of Public Health, Houston, TX 77030, USA
| | - Rongling Wu
- Department of Biostatistics, University of Texas School of Public Health, Houston, TX 77030, USA.,Department of Statistics, The Pennsylvania State University, University Park, PA 16802, USA
| |
Collapse
|
48
|
Rao SV, Muralidhara, Yenisetti SC, Rajini PS. Evidence of neuroprotective effects of saffron and crocin in a Drosophila model of parkinsonism. Neurotoxicology 2016; 52:230-42. [DOI: 10.1016/j.neuro.2015.12.010] [Citation(s) in RCA: 32] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/27/2015] [Revised: 12/05/2015] [Accepted: 12/10/2015] [Indexed: 01/04/2023]
|
49
|
Zhu Z, Lu Q, Zeng F, Wang J, Huang S. Compatibility between mitochondrial and nuclear genomes correlates with the quantitative trait of lifespan in Caenorhabditis elegans. Sci Rep 2015; 5:17303. [PMID: 26601686 PMCID: PMC4658563 DOI: 10.1038/srep17303] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/02/2015] [Accepted: 10/28/2015] [Indexed: 12/12/2022] Open
Abstract
Mutations in mitochondrial genome have epistatic effects on organisms depending on
the nuclear background, but a role for the compatibility of mitochondrial-nuclear
genomes (mit-n) in the quantitative nature of a complex trait remains unexplored. We
studied a panel of recombinant inbred advanced intercrossed lines (RIAILs) of C.
elegans that were established from a cross between the N2 and HW strains. We
determined the HW nuclear genome content and the mitochondrial type (HW or N2) of
each RIAIL strain. We found that the degree of mit-n compatibility was correlated
with the lifespans but not the foraging behaviors of RIAILs. Several known
aging-associated QTLs individually showed no relationship with mitotypes but
collectively a weak trend consistent with a role in mit-n compatibility. By
association mapping, we identified 293 SNPs that showed linkage with lifespan and a
relationship with mitotypes consistent with a role in mit-n compatibility. We
further found an association between mit-n compatibility and several functional
characteristics of mitochondria as well as the expressions of genes involved in the
respiratory oxidation pathway. The results provide the first evidence implicating
mit-n compatibility in the quantitative nature of a complex trait, and may be
informative to certain evolutionary puzzles on hybrids.
Collapse
Affiliation(s)
- Zuobin Zhu
- State Key Laboratory of Medical Genetics, School of Life Sciences, Xiangya Medical School, Central South University, 110 Xiangya Road, Changsha, Hunan, 410078, China
| | - Qing Lu
- State Key Laboratory of Medical Genetics, School of Life Sciences, Xiangya Medical School, Central South University, 110 Xiangya Road, Changsha, Hunan, 410078, China
| | - Fangfang Zeng
- State Key Laboratory of Medical Genetics, School of Life Sciences, Xiangya Medical School, Central South University, 110 Xiangya Road, Changsha, Hunan, 410078, China
| | - Junjing Wang
- State Key Laboratory of Medical Genetics, School of Life Sciences, Xiangya Medical School, Central South University, 110 Xiangya Road, Changsha, Hunan, 410078, China
| | - Shi Huang
- State Key Laboratory of Medical Genetics, School of Life Sciences, Xiangya Medical School, Central South University, 110 Xiangya Road, Changsha, Hunan, 410078, China
| |
Collapse
|
50
|
Collective effects of common SNPs in foraging decisions in Caenorhabditis elegans and an integrative method of identification of candidate genes. Sci Rep 2015; 5:16904. [PMID: 26581252 PMCID: PMC4652280 DOI: 10.1038/srep16904] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/21/2015] [Accepted: 10/22/2015] [Indexed: 01/27/2023] Open
Abstract
Optimal foraging decision is a quantitative flexible behavior, which describes the time at which animals choose to abandon a depleting food supply. The total minor allele content (MAC) in an individual has been shown to correlate with quantitative variations in complex traits. We have studied the role of MAC in the decision to leave a food lawn in recombinant inbred advanced intercross lines (RIAILs) of Caenorhabditis elegans. We found a strong link between MAC and the food lawn leaving rates (Spearman r = 0.4, P = 0.005). We identified 28 genes of unknown functions whose expression levels correlated with both MAC and leaving rates. When examined by RNAi experiments, 8 of 10 tested among the 28 affected leaving rates, whereas only 2 of 9 did among genes that were only associated with leaving rates but not MAC (8/10 vs 2/9, P < 0.05). The results establish a link between MAC and the foraging behavior and identify 8 genes that may play a role in linking MAC with the quantitative nature of the trait. The method of correlations with both MAC and traits may find broad applications in high efficiency identification of target genes for other complex traits in model organisms and humans.
Collapse
|