Natural Genetic Variation Influences Protein Abundances in C. elegans Developmental Signalling Pathways

The impact of genomic variation on protein phosphorylation states and regulatory networks

Genomic variation impacts on cellular networks by affecting the abundance (e.g., protein levels) and the functional states (e.g., protein phosphorylation) of their components. Previous work has focused on the former, while in this context, the functional states of proteins have largely remained neglected. Here, we generated high‐quality transcriptome, proteome, and phosphoproteome data for a panel of 112 genomically well‐defined yeast strains. Genetic effects on transcripts were generally transmitted to the protein layer, but specific gene groups, such as ribosomal proteins, showed diverging effects on protein levels compared with RNA levels. Phosphorylation states proved crucial to unravel genetic effects on signaling networks. Correspondingly, genetic variants that cause phosphorylation changes were mostly different from those causing abundance changes in the respective proteins. Underscoring their relevance for cell physiology, phosphorylation traits were more strongly correlated with cell physiological traits such as chemical compound resistance or cell morphology, compared with transcript or protein abundance. This study demonstrates how molecular networks mediate the effects of genomic variants to cellular traits and highlights the particular importance of protein phosphorylation.

Multi-Omics Approach Reveals Dysregulation of Protein Phosphorylation Correlated with Lipid Metabolism in Mouse Non-Alcoholic Fatty Liver

Obesity caused by overnutrition is a major risk factor for non-alcoholic fatty liver disease (NAFLD). Several lipid intermediates such as fatty acids, glycerophospholipids and sphingolipids are implicated in NAFLD, but detailed characterization of lipids and their functional links to proteome and phosphoproteome remain to be elucidated. To characterize this complex molecular relationship, we used a multi-omics approach by conducting comparative proteomic, phoshoproteomic and lipidomic analyses of high fat (HFD) and low fat (LFD) diet fed mice livers. We quantified 2447 proteins and 1339 phosphoproteins containing 1650 class I phosphosites, of which 669 phosphosites were significantly different between HFD and LFD mice livers. We detected alterations of proteins associated with cellular metabolic processes such as small molecule catabolic process, monocarboxylic acid, long- and medium-chain fatty acid, and ketone body metabolic processes, and peroxisome organization. We observed a significant downregulation of protein phosphorylation in HFD fed mice liver in general. Untargeted lipidomics identified upregulation of triacylglycerols, glycerolipids and ether glycerophosphocholines and downregulation of glycerophospholipids, such as lysoglycerophospholipids, as well as ceramides and acylcarnitines. Analysis of differentially regulated phosphosites revealed phosphorylation dependent deregulation of insulin signaling as well as lipogenic and lipolytic pathways during HFD induced obesity. Thus, this study reveals a molecular connection between decreased protein phosphorylation and lipolysis, as well as lipid-mediated signaling in diet-induced obesity.

Proteome-wide systems genetics identifies UFMylation as a regulator of skeletal muscle function

Improving muscle function has great potential to improve the quality of life. To identify novel regulators of skeletal muscle metabolism and function, we performed a proteomic analysis of gastrocnemius muscle from 73 genetically distinct inbred mouse strains, and integrated the data with previously acquired genomics and >300 molecular/phenotypic traits via quantitative trait loci mapping and correlation network analysis. These data identified thousands of associations between protein abundance and phenotypes and can be accessed online (https://muscle.coffeeprot.com/) to identify regulators of muscle function. We used this resource to prioritize targets for a functional genomic screen in human bioengineered skeletal muscle. This identified several negative regulators of muscle function including UFC1, an E2 ligase for protein UFMylation. We show UFMylation is up-regulated in a mouse model of amyotrophic lateral sclerosis, a disease that involves muscle atrophy. Furthermore, in vivo knockdown of UFMylation increased contraction force, implicating its role as a negative regulator of skeletal muscle function.

The genetics of gene expression in a Caenorhabditis elegans multiparental recombinant inbred line population

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.

Proteome-wide Systems Genetics to Identify Functional Regulators of Complex Traits

Proteomic technologies now enable the rapid quantification of thousands of proteins across genetically diverse samples. Integration of these data with systems-genetics analyses is a powerful approach to identify new regulators of economically important or disease-relevant phenotypes in various populations. In this review, we summarize the latest proteomic technologies and discuss technical challenges for their use in population studies. We demonstrate how the analysis of correlation structure and loci mapping can be used to identify genetic factors regulating functional protein networks and complex traits. Finally, we provide an extensive summary of the use of proteome-wide systems genetics throughout fungi, plant, and animal kingdoms and discuss the power of this approach to identify candidate regulators and drug targets in large human consortium studies.

Natural variation in the sequestosome-related gene, sqst-5, underlies zinc homeostasis in Caenorhabditis elegans

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.

The Gene scb-1 Underlies Variation in Caenorhabditis elegans Chemotherapeutic Responses

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.

A Cell Fate Switch in the Caenorhabditis elegans Seam Cell Lineage Occurs Through Modulation of the Wnt Asymmetry Pathway in Response to Temperature Increase

Populations often display consistent developmental phenotypes across individuals despite inevitable biological stochasticity. Nevertheless, developmental robustness has limits, and systems can fail upon change in the environment or the genetic background. We use here the seam cells, a population of epidermal stem cells in Caenorhabditis elegans, to study the influence of temperature change and genetic variation on cell fate. Seam cell development has mostly been studied so far in the laboratory reference strain (N2), grown at 20° temperature. We demonstrate that an increase in culture temperature to 25° introduces variability in the wild-type seam cell lineage, with a proportion of animals showing an increase in seam cell number. We map this increase to lineage-specific symmetrization events of normally asymmetric cell divisions at the fourth larval stage, leading to the retention of seam cell fate in both daughter cells. Using genetics and single-molecule imaging, we demonstrate that this symmetrization occurs via changes in the Wnt asymmetry pathway, leading to aberrant Wnt target activation in anterior cell daughters. We find that intrinsic differences in the Wnt asymmetry pathway already exist between seam cells at 20° and this may sensitize cells toward a cell fate switch at increased temperature. Finally, we demonstrate that wild isolates of C. elegans display variation in seam cell sensitivity to increased culture temperature, although their average seam cell number is comparable at 20°. Our results highlight how temperature can modulate cell fate decisions in an invertebrate model of stem cell patterning.

WormQTL2: an interactive platform for systems genetics in Caenorhabditis elegans

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/.

Pmr-1 gene affects susceptibility of Caenorhabditis elegans to Staphylococcus aureus infection through glycosylation and stress response pathways' alterations

Calcium signaling can elicit different pathways involved in an extreme variety of biological processes. Calcium levels must be tightly regulated in a spatial and temporal manner in order to be efficiently and properly utilized in the host physiology. The Ca²⁺-ATPase, encoded by pmr-1 gene, was first identified in yeast and localized to the Golgi and it appears to be involved in calcium homeostasis. PMR-1 function is evolutionary conserved from yeast to human, where mutations in the orthologous gene ATP2C1 cause Hailey-Hailey disease. In this work, we used the Caenorhabditis elegans model system to gain insight into the downstream response elicited by the loss of pmr-1 gene. We found that pmr-1 knocked down animals not only showed defects in the oligosaccharide structure of glycoproteins at the cell surface but also were characterized by reduced susceptibility to bacterial infection. Although increased resistance to the infection might be related to lack of regular recognition of C. elegans surface glycoproteins by microbial agents, we provide genetic evidence that pmr-1 interfered nematodes mounted a stronger innate immune response to Gram-positive bacterial infection. Thus, our observations indicate pmr-1 as a candidate gene implicated in mediating the worm's innate immune response.

The intrinsic and extrinsic factors that contribute to proteostasis decline and pathological protein misfolding

Proteostasis refers to the ability of cells to maintain the health of the proteome. Highly conserved quality control mechanisms exist to maintain proteostasis. These include the heat shock response, the unfolded protein response, and protein clearance/degradation pathways. Together, these mechanisms and others comprise the proteostasis network. This network is under constant assault and is strikingly sensitive to changes in the protein folding environment, resulting in proteostasis collapse under certain conditions. Here, the intrinsic and extrinsic stresses experienced by the proteostasis network are explored. The intrinsic stresses include genetic background as well as transcriptional and translational fidelity. These cause changes in the abundance or amino acid sequence of cellular proteins. Extrinsic stresses refer to environmental perturbation of the proteome, such as those caused by temperature stress, oxidative stress, air pollution and cigarette smoke. As the stress to the proteome exceeds the capacity of the proteostasis network, progressive neurodegenerative diseases of aging, such as Alzheimer's disease and Huntington's disease are more likely to ensue.

Illuminating Biological Interactions with in Vivo Protein Footprinting

Protein footprinting coupled with mass spectrometry is being increasingly used for the study of protein interactions and conformations. The hydroxyl radical footprinting method, fast photochemical oxidation of proteins (FPOP), utilizes hydroxyl radicals to oxidatively modify solvent accessible amino acids. Here, we describe the further development of FPOP for protein structural analysis in vivo (IV-FPOP) with Caenorhabditis elegans. C. elegans, part of the nematode family, are used as model systems for many human diseases. The ability to perform structural studies in these worms would provide insight into the role of structure in disease pathogenesis. Many parameters were optimized for labeling within the worms including the microfluidic flow system and hydrogen peroxide concentration. IV-FPOP was able to modify several hundred proteins in various organs within the worms. The method successfully probed solvent accessibility similarily to in vitro FPOP, demonstrating its potential for use as a structural technique in a multiorgan system. The coupling of the method with mass spectrometry allows for amino-acid-residue-level structural information, a higher resolution than currently available in vivo methods.

Genetic background modifies phenotypic and transcriptional responses in a C. elegans model of α-synuclein toxicity

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.

A multi-parent recombinant inbred line population of C. elegans allows identification of novel QTLs for complex life history traits

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.

Shared Genomic Regions Underlie Natural Variation in Diverse Toxin Responses

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.

Use of a Sibling Subtraction Method for Identifying Causal Mutations in Caenorhabditis elegans by Whole-Genome Sequencing

Whole-genome sequencing (WGS) is an indispensable tool for identifying causal mutations obtained from genetic screens. To reduce the number of causal mutation candidates typically uncovered by WGS, Caenorhabditis elegans researchers have developed several strategies. One involves crossing N2-background mutants to the polymorphic Hawaiian (HA) strain, which can be used to simultaneously identify mutant strain variants and obtain high-density mapping information. This approach, however, is not well suited for uncovering mutations in complex genetic backgrounds, and HA polymorphisms can alter phenotypes. Other approaches make use of DNA variants present in the initial background or introduced by mutagenesis. This information is used to implicate genomic regions with high densities of DNA lesions that persist after backcrossing, but these methods can provide lower resolution than HA mapping. To identify suppressor mutations using WGS, we developed an approach termed the sibling subtraction method (SSM). This method works by eliminating variants present in both mutants and their nonmutant siblings, thus greatly reducing the number of candidates. We used this method with two members of the C. elegans NimA-related kinase family, nekl-2 and nekl-3. Combining weak aphenotypic alleles of nekl-2 and nekl-3 leads to penetrant molting defects and larval arrest. We isolated ∼50 suppressors of nekl-2; nekl-3 synthetic lethality using F1 clonal screening methods and a peel-1–based counterselection strategy. When applied to five of the suppressors, SSM led to only one to four suppressor candidates per strain. Thus SSM is a powerful approach for identifying causal mutations in any genetic background and provides an alternative to current methods.

Differential regulation of germ line apoptosis and germ cell differentiation by CPEB family members in C. elegans

Cytoplasmic polyadenylation element binding (CPEB) proteins are evolutionary conserved RNA-binding proteins that control mRNA polyadenylation and translation. Orthologs in humans and other vertebrates are mainly involved in oogenesis. This is also the case for the C. elegans CPEB family member CPB-3, whereas two further CPEB proteins (CPB-1 and FOG-1) are involved in spermatogenesis. Here we describe the characterisation of a new missense allele of cpb-3 and show that loss of cpb-3 function leads to an increase in physiological germ cell death. To better understand the interaction and effect of C. elegans CPEB proteins on processes such as physiological apoptosis, germ cell differentiation, and regulation of gene expression, we characterised changes in the transcriptome and proteome of C. elegans CPEB mutants. Our results show that, despite their sequence similarities CPEB family members tend to have distinct overall effects on gene expression (both at the transcript and protein levels). This observation is consistent with the distinct phenotypes observed in the various CPEB family mutants.

Genetic variation in neurodegenerative diseases and its accessibility in the model organism Caenorhabditis elegans

BACKGROUND

Neurodegenerative diseases (NGDs) such as Alzheimer's and Parkinson's are debilitating and largely untreatable conditions strongly linked to age. The clinical, neuropathological, and genetic components of NGDs indicate that neurodegeneration is a complex trait determined by multiple genes and by the environment.

MAIN BODY

The symptoms of NGDs differ among individuals due to their genetic background, and this variation affects the onset and progression of NGD and NGD-like states. Such genetic variation affects the molecular and cellular processes underlying NGDs, leading to differential clinical phenotypes. So far, we have a limited understanding of the mechanisms of individual background variation. Here, we consider how variation between genetic backgrounds affects the mechanisms of aging and proteostasis in NGD phenotypes. We discuss how the nematode Caenorhabditis elegans can be used to identify the role of variation between genetic backgrounds. Additionally, we review advances in C. elegans methods that can facilitate the identification of NGD regulators and/or networks.

CONCLUSION

Genetic variation both in disease genes and in regulatory factors that modulate onset and progression of NGDs are incompletely understood. The nematode C. elegans represents a valuable system in which to address such questions.