- Research article
- Open Access
Genetic background modifies phenotypic and transcriptional responses in a C. elegans model of α-synuclein toxicity
BMC Genomicsvolume 20, Article number: 232 (2019)
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.
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.
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.
Dementia is a growing global problem affecting large numbers of people. The most common causes of dementia are associated with neurodegeneration, such as that resulting from Alzheimer’s disease (AD) and Parkinson’s disease (PD) [1,2,3]. These neurodegenerative dementias share neuroanatomical and biochemical similarities and result from protein misfolding. A major difficulty in determining the mechanisms that produce neurodegenerative dementia is that the underlying cellular-level pathology varies greatly among patients (reviewed in [4, 5]). These differences in pathology allow for genome-wide association studies (GWAS) of PD patients, identifying genetic risk-variants and demonstrate the complex genetic architecture of neurodegenerative dementia [6,7,8].
Model organisms, such as the nematode Caenorhabditis elegans, are of great value for studying neurodegenerative diseases [9,10,11,12,13]. This is due to their experimental tractability and to the broad and general conservation of genetic pathways across species. Importantly, C. elegans allows for sophisticated genome-wide genetic screens, which are much less complex than similar studies on mammals. In C. elegans, most of these screens have sought to identify loci that represent potential candidate diagnostic and therapeutic targets, or to understand the underlying pathological processes. For example, analysis of transgenic C. elegans expressing the human amyloid-β peptide – the main component of the amyloid plaques found in AD – linked toxicity to insulin/insulin-like growth factor (IGF), dietary restriction and the heat shock response [9, 10, 12, 14,15,16]. Similarly for PD, analysis of worms expressing α-synuclein identified associations with ageing and insulin-like signalling  and with autophagy and lysosomal function .
Most C. elegans studies are however limited to the canonical Bristol N2 genetic background, and therefore do not provide insight into how genetic variation between individuals might affect disease-associated traits. This is a major issue as it is clear, both specifically in C. elegans (e.g. [19,20,21,22]) and more generally in other systems [23,24,25], that the phenotype of any given mutation, transgene or allele can vary depending on the genetic background . For C. elegans, natural genetic variation is known to result in extensive phenotypic variation (e.g. [27, 28], and see [29, 30] for reviews for older studies) and differentially affects both the proteome [31, 32] and transcriptome [22, 31, 33,34,35,36,37]. Therefore, only studying mutational effects in a single genetic background biases our understanding of disease phenotypes, and represents a missed opportunity to elucidate disease mechanisms.
Exemplifying this pattern of a reliance on a single genetic background, only one study has, to date, looked at protein misfolding disease in multiple wild isolate genetic backgrounds of C. elegans . Crucially, this showed that natural genetic variation can uncouple different phenotypic effects of polyglutamine (polyQ40) expression . This strongly suggests that similar important variation between genetic backgrounds will be found for other protein misfolding diseases. Given that for protein misfolding diseases the nature of the modifying alleles segregating within human populations remain largely elusive , the experimental tractability of C. elegans in combination with genetic variation makes the species an excellent system in which to address this issue.
In C. elegans, expression of an α-synuclein and yellow fluorescent protein (YFP) fusion in the body wall muscle results in an age-dependent accumulation of inclusions . These inclusions of α-synuclein form aggregates in aging worms that are similar to the pathological inclusions seen in humans with PD . This therefore represents an appropriate model of α-synuclein toxicity . To investigate the effect of genetic background on the consequences of α-synuclein expression, we have created introgression lines (ILs) containing this α-synuclein and YFP transgene in the background of four wild isolates of C. elegans. Our analyses of these new ILs, and of α-synuclein in an N2 genetic background, identify both general and genotype-specific changes in gene expression. These changes predict a range of phenotypic effects that we then experimentally confirm. Importantly, given the reliance on N2 in C. elegans research, we show both that some effects are N2-specific and that for other phenotypes the analysis of other genetic backgrounds uncovers substantial variation not seen in N2. Our approach therefore evaluates how genetic background conditions transgene effect(s), and specifically illustrates how the consequences of ectopic overexpression of α-synuclein depends on genetic background.
Introgression line construction and validation
We introgressed the pkIs2386 transgene [unc-54p::α-synuclein::YFP + unc-119(+)] from NL5901 , which has an N2 genetic background, into the genetically divergent wild isolates JU1511, JU1926, JU1931, and JU1941 (see  for information on genetic distance between these lines). The pkIs2386 transgene results in α-synuclein expression in the body wall musculature and the vulval muscles. After back-crossing and selfing, four new ILs were obtained, SCH1511, SCH1926, SCH1931, and SCH1941, with for example, SCH1511 containing the transgene in a JU1511 background. In combination with N2 and NL5901 as controls, we were therefore able to investigate the phenotypic and genomic effects of α-synuclein in five genetic backgrounds (N2, JU1511, JU1926, JU1931, and JU1941). An additional IL, SCH4856, was subsequently created in which the transgene has been introgressed into a CB4856 background. SCH4856 was used to test for the effects of N2 alleles (see below).
We sought to identify the site of the introgression and to determine how much of the N2 genome surrounding the transgene had also been introgressed into the wild isolates. PCR-based genotyping located the transgene in chromosome IV, indicating that the homozygous introgressions in the new ILs spanned between 4.2 and 13.2Mb of chromosome IV (Additional file 1). By using the set of genetic markers from Volkers et al. , we identified significantly differential expressed genes on chromosome IV of the α-synuclein lines, and also identified an additional introgression on chromosome V in SCH1931 (Additional file 2). PCR-based genotyping together with the transcriptomic analysis indicated a consistent genomic location for the α-synuclein transgene in each of the genetic backgrounds, but did not allow detection of the precise introgression boundaries.
Transcriptome effects of the α-synuclein introgression depend on the genetic background
The effects of α-synuclein expression on gene expression was measured in five different genetic backgrounds, with three biological replicates per genotype (see Additional file 3 and ArrayExpress accession E-MTAB-6960 for details). Analysis of genome-wide transcriptional changes in the α-synuclein expressing worms indicated that there were differences in developmental rate between the ILs. Variation between ILs in developmental rate was also observed during IL construction. We therefore estimated the age of the samples by their gene expression profile (as in [41, 42]) and used principal component analysis (PCA) on genome-wide expression levels to investigate differences between isolates and the effects of α-synuclein. This analysis revealed separation between the wild isolates and their corresponding transgenic ILs (Fig. 1a). Although genotypes were more scattered on these first two PCA axes, this indicated that developmental delay and reduced lifespan were associated with the introgressed α-synuclein transgene (Fig. 1a). A large part of the global gene expression differences are therefore explained by slower development in the α-synuclein lines. Hence, many of the genes identified as responding to α-synuclein in the different backgrounds were likely to be associated with their differential age. Using a threshold of -log10(p) > 3.4 (FDR < 0.05) we found 3521 genes affected by Age, 1509 affected by Genotype, 646 by α-synuclein and 428 by an interaction between Genotype and α-synuclein (Fig. 1b, Additional file 4). Genes affected by α-synuclein or the interaction between α-synuclein and age were most often also affected by genotype. The overlap among these four factors yielded 78 genes highly specific for α-synuclein relative to both Age and Genotype (Additional file 4). This showed that both genotype-specific and universal (genotype independent) changes in gene expression were induced by α-synuclein expression, i.e. we identified a core set of genes of which the expression was altered in all genetic backgrounds and others that were genotype-specific. Strikingly, many genes were not expressed and/or regulated in the same way in wild isolate genetic backgrounds compared with N2. These data therefore indicate that the transcriptome effects of the α-synuclein introgression depend on the genetic background.
Gene expression enrichment analysis
Gene ontology analysis (GO) of differentially expressed genes revealed the molecular, cellular, and biological processes affected by development, α-synuclein and genetic background (Table 1, with full results in Additional file 4). Genes that change expression in response to genetic background were enriched for genes involved in the innate immune response and oxidation-reduction process (Additional file 4) as was previously found by . Widespread changes in genes related to muscle function were observed (Table 1), an expected response given the changes in cellular environment induced by expression of α-synuclein in the body wall muscle. This analysis also identified changes in genes involved in pharyngeal pumping (Table 1).
As expected there were also changes in various pathways associated with cellular stress responses (Table 1), and protein homeostasis (Additional file 4). Taking the effect of both α-synuclein and age into consideration, genes associated with metabolic processes, transporters of ions and lipids, and kinase activities for ATP were enriched (Additional file 4).
Phenotypic effects of α-synuclein expression vary among genetic backgrounds
Expression of α-synuclein slows development in some genetic backgrounds
Gene expression analysis indicated that α-synuclein lines were developmentally delayed. To test this directly, we scored development time, i.e. the time to the first appearance of eggs. Analysis of these data indicated that development was affected by α-synuclein (aS, p < 2e-8), genetic background (Genotype, p < 2e-8) and the interaction between α-synuclein and genetic background (aS x Genotype, p < 0.0002). In most α-synuclein expressing lines development was delayed compared to the corresponding wild isolate (Fig. 2a), with SCH1941 and SCH1926 showing significantly longer larval development time periods than the corresponding wild type controls (p < 1e-6, and p < 0.008 for SCH1941 and SCH1926, respectively) (Fig. 2a).
As found in the phenotypic assay, the transcriptomic samples of the α-synuclein lines were estimated to be younger than the corresponding wild isolate and so showed delayed development (Fig. 2b, ANOVA model Estimated_Age ~ aS*Genotype; aS p < 6e-5), which was independent of the genetic background (Genotype, p = 0.72) and of the interaction between α-synuclein and genetic background (aS x Genotype, p = 0.84). Comparing the WT and α-synuclein lines for each individual line showed that SCH1931 and JU1931 again displayed the smallest developmental difference as found in the phenotypic assay. Overall, we therefore conclude that the α-synuclein introgression had a genotype-specific impact, differentially decreasing the developmental rate in each genetic background.
Expression of α-synuclein decreases pharyngeal pumping rate in some genetic backgrounds
Given the observation of differential expression of genes with a function in pharyngeal pumping (Table 1), we measured pharyngeal pumping rate in all lines 48 and 72 h after recovery from L1 arrest (Fig. 3). No significant differences were seen at 48 h, but 72 h after recovery from L1 arrest pharyngeal pumping was affected by α-synuclein (aS, p < 3e-16), by the genetic background (Genotype, p < 3e-10) and by the interaction of α-synuclein with genetic background (aS x Genotype, p < 0.004; Fig. 3b). At this point, the pharyngeal pumping rate at 72 h in the α-synuclein expressing lines had slowed down compared to their corresponding wild isolate (Fig. 3b; N2 p = 0.93; JU1511 p < 5e-4; JU1926 p < 2e-5; JU1931 p = 0.54; JU1941 p < 1e-8). These data therefore indicated that the presence of an α-synuclein introgression had a noticeable genotype-specific impact on pharyngeal pumping rate. As pumping rate primarily reflects pharyngeal muscle activity, a tissue in which the α-synuclein transgene is not expressed, these changes reflect an indirect systemic effect of α-synuclein expression.
Expression of α-synuclein decreases movement in some genetic backgrounds
Given the site of α-synuclein expression and the detection of changes in the expression of genes involved in muscle structure and function (Table 1) we tested for differences between lines in movement. Here, we assayed movement of worms in liquid (Fig. 3c). These data indicate that N2 and NL5901 show similar movement rates and that the α-synuclein expressing lines expression show reduced activity in comparison to their corresponding wild isolate. These data therefore indicate that the presence of an α-synuclein introgression produces a genotype-specific impact, likely a direct consequence of α-synuclein expression in the body wall muscles.
Expression of α-synuclein decreases lifespan in some genetic backgrounds
Given the age-related changes in pharyngeal pumping rate and enrichment of genes associated with aging and stress response pathways in response to α-synuclein expression, we hypothesised that the introgressed α-synuclein transgene could affect longevity. Therefore, we measured lifespan in all lines (Additional file 5). Comparison of N2 and NL5901 (α-synuclein in an N2 background) showed that lifespan was not affected by α-synuclein expression in the N2 genetic background (Fig. 4; N2 log-rank p = 0.14, mean age p = 0.99). However, in the other genetic backgrounds, lines containing the α-synuclein introgression displayed significantly accelerated death (Fig. 4a) and a shortened lifespan (Fig. 4b), a global effect of α-synuclein expression on animal physiology and function. These data also indicated that the α-synuclein introgression in the wild isolate genetic backgrounds resulted in increased rates of maternal hatching (bagging) (Additional file 6, Fig. 5). Given that the α-synuclein transgene is expressed in the vulval muscles this result mirrors the direct effect of α-synuclein expression on muscle function seen for thrashing (Fig. 3c). These observations indicate that α-synuclein expression lowers lifespan in some genotypes, but not others.
Genetic background affects the outcome of α-synuclein expression
Comparison of phenotypic effects across the genetic backgrounds tested suggests that the phenotypes associated with the α-synuclein transgene introgression vary between genotypes in a consistent manner (Fig. 5). Multiple phenotypes were affected by α-synuclein in the JU1511, JU1926, and JU1941 background, whereas they were less affected in the JU1931 and N2 backgrounds. This suggests a shared genetic component across many of the phenotypes we have observed, but does also indicate that there is genotype-specific variation, e.g. no pair of lines shows the same pattern of phenotypic effects (Fig. 5). These data also indicate that there are both local effects – e.g. bagging and trashing, where effects can be directly attributed to expression of α-synuclein in the muscles – and global effects on the whole animal.
As alleles from different C. elegans backgrounds can produce a range of synthetic deleterious effects resulting in full or partially genomic incompatibilities between their genomes [43, 44], we sought to test if the N2 region alone replicated the phenotypes we observed here. We therefore introgressed the pkIs2386 transgene into a CB4856 genetic background, generating the SCH4856 line and undertook comparisons of this line with CB4856 and CBN93, a line with an introgression of the N2 genome spanning the 3.3–12.8Mbp region on chromosome IV in an CB4856 background (Additional file 1). This control was undertaken as deleterious interactions between alleles from N2 and those from CB4856 are well characterized [43, 44]. Comparisons of pumping rate and of development time between these lines indicated that the expression of α-synuclein in SCH4856 produced effects not seen in CBN93 (Fig. 6) and that these effects mirror those seen in the other genetic backgrounds. This provides strong support for the view that the phenotypic effects we observe were a consequence of α-synuclein expression and not the introgression of the N2 region.
We have introgressed a transgene that results in the expression of α-synuclein in the body wall and vulval muscles into different genetic backgrounds of C. elegans. Analysis of these newly created α-synuclein ILs indicates that genetic background effects both the response of the animals to α-synuclein expression at the level of gene expression and in terms of the phenotypic consequences. Further, this can be seen at both a local and global level, with effects seen in traits associated with muscles in which α-synuclein is expressed (bagging and thrashing) and traits that represent a global, whole organism, effects (pumping, development rate and lifespan). Our gene expression analysis identifies a range of changes associated with α-synuclein expression. Many of these – particularly those involved in muscle function and in various stress responses – matched expectations given the site of expression and the known effects of α-synuclein aggregation on cellular function. We also identify changes in lipid metabolism (Table 1 and Additional file 4), an important result given that this has been previously linked to α-synuclein pathology in a yeast model  and that there is evidence of a direct association between α-synuclein and lipid droplets . Changes in lipid levels of C. elegans expressing α-synuclein in the muscles have previously been observed , our results therefore provide a set of candidate genes to investigate how α-synuclein expression may mediate lipid metabolism. We also show that such effects can be expected to vary across C. elegans genetic backgrounds.
On a phenotypic level we found that, in comparison to the corresponding wild isolates, the α-synuclein expressing lines: developed more slowly; suffered an increased rate of decline in pharyngeal pumping; exhibited a shorter lifespan; and died from matricidal hatching of eggs at an increased rate. Similar phenotypic effects are observed in other C. elegans models of protein misfolding disease. For example, transgenic worms expressing polyQ proteins and Aβ both show a significantly shortened lifespan [38, 48]. A reduced lifespan could therefore be a general toxicity phenotype, providing an indirect measure of organismal dysfunction caused by misfolded protein aggregations in the body wall.
Many of the phenotypic changes we see are either only found in wild isolate genetic backgrounds or are more pronounced in these backgrounds (Fig. 5). This demonstrates that there is variation between isolates that appears to differently affect local and global consequences of α-synuclein expression. Previously, only a single C. elegans study has directly investigated the effect of genetic background in the context of protein misfolding disease. That study found complex variation in polyQ40 aggregation and toxicity in three wild isolate backgrounds and in a panel of 21 recombinant inbred lines (RILs) . The RIL analysis also showed that the various effects of polyQ40 expression could be uncoupled . In this study, we also show some evidence that natural genetic variation in C. elegans can uncouple different effects of α-synuclein expression (Fig. 5), suggesting that this may be a general pattern.
A caveat with this interpretation is however that a range of incompatibilities between alleles from different C. elegans genomes have been identified [43, 44] and that the region surrounding the site of the α-synuclein transgene contains a number of mapped quantitative trait loci affecting various life history traits. For example, ILs with introgressed regions of chromosome IV from CB4856 in an N2 genetic background have identified QTLs affecting lifespan and pumping rate that partially overlap our introgression . Similarly, a complex interaction between alleles from the N2 and CB4856 genomes has been found on chromosome IV . However, little is known about the frequency of such synthetic deleterious effects between other genotypes of C. elegans. Proteomic analysis of age-dependent changes in protein solubility by Reis-Rodrigues et al.  also indicates that the genes encoding proteins that become insoluble with age are enriched for modifiers of lifespan. Additionally, a heat-stress specific QTL for recovery  as well as an eQTL-hotspot have been identified on the left arm of chromosome IV (IV: 1.0–2.5 Mb) , while a QTL for maternal hatching rate is at the position of the introgression (chromosome IV ~6M) . That such effects might alter disease-related processes complicates the interpretation of both our results here and the previous work of Gidalevitz et al. . To experimentally address this caveat, we introgressed the α-synuclein transgene into a CB4856 genetic background and undertook comparisons with a line that contains a comparable region of the N2 genome, but no transgene, in a CB4856 background. These comparisons indicated that the N2 introgression alone does not recapitulate the phenotypic effects seen in the α-synuclein expressing line (Fig. 6). As these comparisons again show that the expression of α-synuclein in a CB4856 genetic background results in different and more severe effects than those seen in an N2 background, this strongly supports the view that genetic background needs to be considered more generally when C. elegans is used as a disease model. A move to the use of more defined modifications in a range of genetic backgrounds that limited the introduction of other alleles – perhaps via the use of CRISPR to introduce the relevant transgenes – would therefore be ideal. It would also be of interest to determine the extent to which genotype-specific were observed in other models, for example where α-synuclein is expressed in neurons . As there is extensive variation between isolates of C. elegans in lifespan and in various stress responses and in lifespan, and that variation has now been identified for two protein misfolding diseases such a development would facilitate the systematic analysis of the role of natural variation in such diseases.
C. elegans maintenance and growth conditions
Worms were maintained at 20°C on nematode growth medium (NGM) plates seeded with Escherichia coli OP50  and all assays were undertaken at 20°C. Assays were initiated using eggs isolated from gravid adults treated with sodium hypochlorite and NaOH . For the microarray analysis of gene expression and for the pharyngeal pumping assays, these synchronized eggs were allowed to hatch on NGM plates seeded with E. coli overnight at 20°C, i.e. these worms were not arrested at the L1 stage. After 48 h, worms in all strains had reached the L4 stage. For the lifespan, development, and thrashing activity assays, the synchronized eggs were hatched on NGM plates without E.coli, in order to obtain synchronized L1 larvae. When they reach L4 stage, these assay were then measured or setup will all lines.
C. elegans lines
The C. elegans wild types JU1511, JU1926, JU1931, JU1941, CB4856 (Hawaii), canonical strain N2 (Bristol), and the transgenic strain NL5901, which contains the pkIs2386 transgene [unc-54p::α-synuclein::YFP + unc-119(+)], were obtained from the Caenorhabditis Genetics Center. To introgress the α-synuclein::YFP transgene into the wild backgrounds, NL5901 hermaphrodites were mated with wild type males, then five to ten F1 fluorescent males were mated with wild type hermaphrodites. Fluorescent males were then isolated in the F2 progeny of this cross and used to backcross to the relevant wild type. This backcrossing to the wild type background was repeated until generation F7 was reached. At this point, worms were allowed to self-fertilise to obtain homozygotes lines and then cryopreserved. This crossing design was used to produce the strains SCH1511, SCH1926, SCH1931, SCH1941, and SCH4856, with for example, SCH1511 containing the transgene in a JU1511 background.
Five genetic backgrounds – N2, JU1511, JU1926, JU1931, and JU1941 – were used to investigate the phenotypic and genomic effects of α-synuclein expression. SCH4856 was used in separate assays to test the effect of N2 alleles surrounding the transgene integration site. Here, comparisons were made to the CBN093 strain which contains a comparable region of the N2 genome introgressed into CB4856, but no transgene. CBN093 was constructed by back-crossing WN071  with CB4856. The strain was back-crossed followed by segregation for 11 generations until a single homozygous region on chromosome IV was obtained. Thereafter, the strain was inbred for 11 generations and the genotype confirmed by sequencing (SRP154243; https://www.ncbi.nlm.nih.gov/sra) using the approach [56, 57].
DNA was isolated from the lysates of 5 individual adults from each line, and then was used for genotyping PCRs. Genotyping primers utilized insertions/deletions between the CB4856 and N2 genomes, with 41 primer pairs covering the genome used (see Additional file 1 for details of primers and marker locations) . PCR was carried out with the GoTaq DNA polymerase kit (Promega) according to the manufacturer’s recommendations. The sizes of amplified products were assessed by electrophoresis in 1.5% agarose gels stained with Ethidium Bromide. This allowed the size of the introgressions to be determined, with these regions shown in Additional file 1.
Marker genes found by DNA hybridizations in Volkers et al. 2013  were used to find the genomic position of the α-synuclein introgression in each genetic background. We tested which genomic region was missing markers for the wildtype background, this region then must be from the NL5901 background.
Sample preparations and RNA microarray analysis
In total three independent replicates of the four JU strains and their corresponding transgenic ILs, as well as N2 and NL5901, were analysed. Worms were assayed at 48 h from egg isolation (see details above) and all strains were at the L4 stage. Worms were therefore the same chronological age, and at the same larval stage, but this design would allow the detection of effects of genetic background, α-synuclein expression and the interaction of these effects on development. At this point worms were generated by hatching alkaline hypochlorite-purified eggs and then harvested by centrifugation, washed with M9 buffer, frozen in liquid N2, and stored at − 80 °C until use. The Maxwell® 16 Tissue LEV Total RNA Purification Kit was used for mRNA isolation, following the manufacturer’s protocol with a modified lysis step. In addition to the lysis buffer, proteinase K was added and the sample was incubated at 65°C while shaking at 1000 rpm for 10 min. Thereafter the standard protocol was followed. PolyA RNA was used to generated Cy3 and Cy5-labeled cRNA samples, which were then hybridized to 4X44K slides V2 (Agilent) C. elegans whole genome GeneChips, processed, and scanned (full microarray data in ArrayExpress with accession E-MTAB-6960). RNA microarray statistical analysis and data processing were performed using the Limma package for the R software environment . To find the genes affected by α-synuclein, genotype and age, these terms were used as explanatory factors in a linear model (gene expression ~ age + aS * genotype). Significance thresholds were determined by permutations of all spots on the array. In the permutations, the RNA hybridization intensities were randomly distributed over the genotypes and batches. Therefore, the p-value that gave a min ratio of false positives/true positives of 0.05 (= − log10(p) > 3.4) was set for convenience, i.e. an FDR at 0.0186 for α-synuclein effect, 0.0249 for genotype effect and 0.073 for the interaction between α-synuclein and genotype.
Enrichments were done on genes groups divined by their significance from the linear model a -log10(p) > 2 was used, then a hyper geometric test in R was used to test each GO term for enrichment.
To partition the variation in gene expression a PCA was done on the transcription profiles (the log2 ratios with the mean) of all samples. The first two axis were used for visualisation.
Age estimation of the worms sampled for transcriptomics was done by comparing the class I age responsive genes from Snoek et al. , as done by van der Bent et al.  and Jovic et al. . These genes show a linear increase during development between 46 and 54 h and can therefore be used to estimate the relative differences in development between samples/populations by their transcriptomes. The α-synuclein expressing SCH lines developed more slowly than the JU lines and none of the lines had started moulting at 48 h, so there were no adults on the plates.
Nematodes were cultured on OP50 bacteria without fluorouracil deoxyribose (FUdR) from synchronised L1 juveniles until young adults. When setting up the experiment, in total 100 worms per strain were randomly selected from the population and transferred onto 10 plates (i.e. each plate contains 10 worms). Then, they were transferred away from their progeny every day during the vigorous reproductive period and every other day during the reduced reproductive period. From the young adult stage, worms were examined for signs of life daily. Individuals that were not moving or twitching after gentle stimulation, followed by vigorous stimulation, or that did not exhibit pharyngeal pumping for 30 s, were considered dead. Individual worms that had died from internal hatching of progeny (bagging), or that had crawled off the plates, were censored from lifespan result but were used for maternal hatching analysis. Two biological replicates were done for lifespan assay. Data analysis was performed in R using the ‘survival’ package. The rates of bagging from these assays were also compared to determine if this was affected by α-synuclein expression.
Synchronised L1 juveniles were obtained by allowing eggs to hatch in M9 buffer after three washes in M9 after hypochlorite treatment. Arrested L1s were then transferred to NGM plates, which fed with E. coli OP50 and were incubated at 20°C. Tracking observations and inspections were done at regular time intervals. Development time was defined as the period between worm inoculation and the moment at which the first appearance of eggs and the period until the reproductions reach peak level. Three biological replicates were done for this assay.
To avoid the effect of short term starvation after inoculation on pharyngeal pumping rate, synchronized eggs, isolated as described above, were allowed to hatch on NGM plates with E.coli (t = 0 h). Individual worms were then isolated and observed at the L4 stage (t = 48 h, day 2) and as young adults (t = 72 h, day 3). The number of contractions in the terminal bulb of pharynx was counted for 30 s (n = 19 per genotype) for the L4s, and for 60 s for 3 day-old worms (n = 30 per genotype). Three biological replicates were done for this assay.
Synchronized starved L1 juveniles were inoculated on NGM plates with E.coli at 20°C. When N2 worms were observed to have reached the L4 stage, worms were transferred into the wells of a 96-well-plate, where each well contained 50μl of M9 buffer . Three biological replicates were undertaken, with 10 worms per well for 3–5 wells per genotype in each of the replicates. Activity was then measured in an 30 min period using a WMicroTracker (PhylumTech). Here, worm movements in the wells are detected as they interfere with an array of micro beams of infrared light and the number of interference events is counted – hence a higher score represents more activity.
Phenotypic differences between the lines were tested using an ANOVA (model: phenotype~ aS * GB) to test for overall effects and Tukey HSD to test for lines specific differences. Survival curves were tested by log-rank test from the “survival” package in R.
Analysis of variance
Caenorhabditis Genetics Center
Clustered regularly interspaced short palindromic repeats
expression quantitative trait loci
- F2 (number i.):
Second (i. th) filial generation
False discovery rate
Gene ontology analysis
Genome-wide association studies
Insulin/insulin-like growth factor
Larval stage 1/4
A megabase pair
Nematode growth medium
Principal component analysis
Polymerase chain reaction
Quantitative trait locus
Recombinant inbred lines
- Tukey HSD:
Tukey’s range test honestly significant difference
Yellow fluorescent protein
Kivipelto M, Ngandu T, Fratiglioni L, et al. Obesity and vascular risk factors at midlife and the risk of dementia and Alzheimer disease. Arch Neurol. 2005;62(10):1556–60.
Aarsland D, Kurz MW. The epidemiology of dementia associated with Parkinson disease. J Neurol Sci. 2010;289(1):18–22.
Aybek S, Lazeyras F, Gronchi-Perrin A, Burkhard PR, Villemure J-G, Vingerhoets FJG. Hippocampal atrophy predicts conversion to dementia after STN-DBS in Parkinson’s disease. Parkinsonism Relat Disord. 2009;15(7):521–4.
Seelaar H, Kamphorst W, Rosso SM, Azmani A, Masdjedi R, de Koning I, et al. Distinct genetic forms of frontotemporal dementia. Neurology. 2008;71(16):1220–6.
Gratwicke J, Jahanshahi M, Foltynie T. Parkinson ’ s disease dementia : a neural networks perspective. Brain. 2015;138(6):1454–76.
Satake W, Nakabayashi Y, Mizuta I, Hirota Y, Ito C, Kubo M, et al. Genome-wide association study identifies common variants at four loci as genetic risk factors for Parkinson's disease. Nat Genet. 2009;41(12):1303.
Nalls MA, Pankratz N, Lill CM, Do CB, Hernandez DG, Saad M, et al. Large-scale meta-analysis of genome-wide association data identifies six new risk loci for Parkinson’s disease. Nat Genet. 2014;46(9):989.
Chen Y, Cao B, Gu X, Ou R, Wei Q, Liu H, et al. No association between 5 new GWAS-linked loci in Parkinson’s disease and multiple system atrophy in a Chinese population. Neurobiol Aging. 2018;67:202–e7.
Teschendorf D, Link CD. What have worm models told us about the mechanisms of neuronal dysfunction in human neurodegenerative diseases? Mol Neurodegener. 2009;4(1):38.
Calahorro F, Ruiz-Rubio M. Caenorhabditis elegans as an experimental tool for the study of complex neurological diseases: Parkinson’s disease, Alzheimer’s disease and autism spectrum disorder. Invertebr Neurosci. 2011;11(2):73–83.
Rodriguez M, Snoek LB, De Bono M, Kammenga JE. Worms under stress: C. elegans stress response and its relevance to complex human disease and aging. Trends Genet. 2013;29(6):367–74.
Sin O, Michels H, Nollen EAA. Genetic screens in Caenorhabditis elegans models for neurodegenerative diseases. Biochim. Biophys. Acta (BBA) - Mol. Basis Dis. 2014;1842(10):1951–9.
Wang YA, Kammenga JE, Harvey SC. Genetic variation in neurodegenerative diseases and its accessibility in the model organism Caenorhabditis elegans. Human Genomics. 2017;11(1):12.
Fonte V, Kipp DR, Yerg J, Merin D, Forrestal M, Wagner E, et al. Suppression of in vivo β-amyloid peptide toxicity by overexpression of the HSP-16.2 small chaperone protein. J Biol Chem. 2008;283(2):784–91.
Steinkraus KA, Smith ED, Davis C, Carr D, Pendergrass WR, Sutphin GL, et al. Dietary restriction suppresses proteotoxicity and enhances longevity by an hsf-1-dependent mechanism in Caenorhabditis elegans. Aging Cell. 2008;7(3):394–404.
McColl G, Roberts BR, Pukala TL, Kenche VB, Roberts CM, Link CD, et al. Utility of an improved model of amyloid-beta (Aβ1-42) toxicity in Caenorhabditis elegans for drug screening for Alzheimer’s disease. Mol Neurodegener. 2012;7(1):57.
van Ham TJ, Thijssen KL, Breitling R, Hofstra RMW, Plasterk RHA, Nollen EAA. C. elegans model identifies genetic modifiers of α-synuclein inclusion formation during aging. PLoS Genet. 2008;4(3):e1000027.
Hamamichi S, Rivas RN, Knight AL, Cao S, Caldwell KA, Caldwell GA. Hypothesis-based RNAi screening identifies neuroprotective genes in a Parkinson’s disease model. PNSA. 2008;105(2):728–33.
Duveau F, Félix MA. Role of pleiotropy in the evolution of a cryptic developmental variation in Caenorhabditis elegans. PLoS Biol. 2012;10(1):e1001230.
Elvin M, Snoek L, Frejno M, Klemstein U, Kammenga J, Poulin G. A fitness assay for comparing RNAi effects across multiple C. elegans genotypes. BMC Genomics. 2011;12(1):510.
Schmid T, Snoek LB, Fröhli E, van der Bent ML, Kammenga J, Hajnal A. Systemic regulation of RAS/MAPK signaling by the serotonin metabolite 5-HIAA. PLoS Genet. 2015;11(5):e1005236.
Sterken MG, van der Plaat LVB, Riksen JAG, Rodriguez M, Schmid T, Hajnal A, et al. Ras/MAPK modifier loci revealed by eQTL in Caenorhabditis elegans. G3: Genes, Genomes, Genetics. 2017;7(9):3185–93.
Wade CM, Daly MJ. Genetic variation in laboratory mice. Nat Genet. 2005;37(11):1175.
Coe TS, Hamilton PB, Griffiths AM, Hodgson DJ, Wahab MA, Tyler CR. Genetic variation in strains of zebrafish (Danio rerio) and the implications for ecotoxicology studies. Ecotoxicology. 2009;18(1):144–50.
Justice MJ, Dhillon P. Using the mouse to model human disease: increasing validity and reproducibility. Dis Model Mech. 2016;9:101–3.
Kammenga JE. The background puzzle : how identical mutations in the same gene lead to different disease symptoms. FEBS J. 2017;284:3362–73.
Reynolds RM, Phillips PC. Natural variation for lifespan and stress response in the nematode Caenorhabditis remanei. PLoS One. 2013;8(4):e58212.
Martin N, Singh J, Aballay A. Natural genetic variation in the Caenorhabditis elegans response to Pseudomonas aeruginosa. G3: Genes, Genomes, Genetics. 2017;7(4):1137–47.
Kammenga JE, Phillips PC, De Bono M, Doroszuk A. Beyond induced mutants : using worms to study natural variation in genetic pathways. Trends Genet. 2008;24(4):178–85.
Gaertner BE, Phillips PC. Caenorhabditis elegans as a platform for molecular quantitative genetics and the systems biology of natural variation. Genet Res. 2010;92(5–6):331–48.
Kamkina P, Snoek LB, Grossmann J, Volkers RJM, Sterken MG, Daube M, et al. Natural genetic variation differentially affects the proteome and transcriptome in C. elegans. Mol Cell Proteomics. 2016;15(5):1670–80.
Singh KD, Roschitzki B, Snoek LB, Grossmann J, Zheng X, Elvin M, et al. Natural genetic variation influences protein abundances in C. elegans developmental signalling pathways. PLoS One. 2016;11(3):e0149418.
Rockman MV, Skrovanek SS, Kruglyak L. Selection at linked sites shapes heritable phenotypic variation in C. elegans. Science. 2010;330(6002):372–6.
Viñuela A, Snoek LB, Riksen JAG, Kammenga JE. Genome-wide gene expression regulation as a function of genotype and age in C. elegans. Genome Res. 2010;20(7):929–37.
Volkers RJM, Snoek LB, Hubar CJ van H, Coopman R, Chen W, Yang W, et al. Gene-environment and protein-degradation signatures characterize genomic and phenotypic diversity in wild Caenorhabditis elegans populations. BMC Biol. 2013;11(1):93.
Snoek BL, Sterken MG, Bevers RPJ, Volkers RJM, van’t Hof A, Brenchley R, et al. Contribution of trans regulatory eQTL to cryptic genetic variation in C. elegans. BMC Genomics. 2017;18(1):500.
Li Y, Breitling R, Snoek LB, van der Velde KJ, Swertz MA, Riksen J, Jansen RC, Kammenga JE. Global genetic robustness of the alternative splicing machinery in Caenorhabditis elegans. Genetics. 2010;186(1):405–10.
Gidalevitz T, Wang N, Deravaj T, Alexander-Floyd J, Morimoto RI. Natural genetic variation determines susceptibility to aggregation or toxicity in a C. elegans model for polyglutamine disease. BMC Biol. 2013;11(1):100.
Kearney JA. Genetic modifiers of neurological disease. Curr Opin Genet Dev. 2011;21(3):349–53.
Jansen IE, Ye H, Heetveld S, Lechler MC, Michels H, Seinstra RI, et al. Discovery and functional prioritization of Parkinson’s disease candidate genes from large-scale whole exome sequencing. Genome Biol. 2017;18(1):22.
Snoek LB, Sterken MG, Volkers RJM, Klatter M, Bosman KJ, Bevers RPJ, et al. A rapid and massive gene expression shift marking adolescent transition in C. elegans. Sci Rep. 2014a;4:3912.
van der Bent ML, Sterken MG, Volkers RJM, Riksen JAG, Schmid T, Hajnal A, et al. Loss-of-function of β-catenin bar-1 slows development and activates the Wnt pathway in Caenorhabditis elegans. Sci Rep. 2014;4:4926.
Seidel HS, Rockman MV, Kruglyak L. Widespread genetic incompatibility in C. elegans maintained by balancing selection. Science. 2008;319(5863):589–94.
Snoek LB, Orbidans HE, Stastna JJ, Aartse A, Rodriguez M, Riksen JAG, et al. Widespread genomic incompatibilities in Caenorhabditis elegans. G3: genes, genomes. Genetics. 2014b;4(10):1813–23.
Cooper AA, Gitler AD, Cashikar A, Haynes CM, Hill KJ, Bhullar B, et al. α-Synuclein blocks ER-Golgi traffic and Rab1 rescues neuron loss in Parkinson's models. Science. 2006;313(5785):324–8.
Cole NB, Murphy DD, Grider T, Rueter S, Brasaemle D, Nussbaum RL. Lipid droplet binding and oligomerization properties of the Parkinson's disease protein α-synuclein. J Biol Chem. 2002;277(8):6344–52.
Jadiya P, Khan A, Sammi SR, Kaur S, Mir SS, Nazir A. Anti-parkinsonian effects of Bacopa monnieri: insights from transgenic and pharmacological Caenorhabditis elegans models of Parkinson’s disease. Biochem Biophys Res Commun. 2011;413(4):605–10.
Fong S, Teo E, Ng LF, Chen C, Lakshmanan LN, Tsoi SY, et al. Energy crisis precedes global metabolic failure in a novel Caenorhabditis elegans Alzheimer disease model. Sci Rep. 2016;6:33781.
Doroszuk A, Snoek LB, Fradin E, Riksen J, Kammenga J. A genome-wide library of CB4856/N2 introgression lines of CB4856/N2 introgression lines of Caenorhabditis elegans. Nucleic Acids Res. 2009;37(16):e110.
Reis-Rodrigues P, Czerwieniec G, Peters TW, Evani US, Alavez S, Gaman EA, et al. Proteomic analysis of age-dependent changes in protein solubility identifies genes that modulate lifespan. Aging Cell. 2012;11(1):120–7.
Rodriguez M, Snoek LB, Riksen JAG, Bevers RP, Kammenga JE. Genetic variation for stress-response hormesis in C. elegans lifespan. Exp Gerontol. 2012;47(8):581–7.
Stastna JJ, Snoek LB, Kammenga JE, Harvey SC. Genotype-dependent lifespan effects in peptone deprived Caenorhabditis elegans. Sci Rep. 2015;5:16259.
Brenner S. The genetics of Caenorhabditis elegans. Genetics. 1974;77(1):71–94.
Stiernagle T. Maintenance of C. elegans. WormBook. The C. elegans research community; 2006.
Li Y, Álvarez OA, Gutteling EW, Tijsterman M, Fu J, Riksen JAG, et al. Mapping determinants of gene expression plasticity by genetical genomics in C. elegans. PLoS Genet. 2006;2(12):e222.
Gao AW, Sterken MG, uit de Bos J, van Creij J, Kamble R. Natural genetic variation in C . elegans reveals genomic loci controlling metabolite levels. Genome Res. 2018;28(9):1296–308.
Sterken MG. Building towards a multi-dimensional genetic architecture in Caenorhabditis elegans: Doctoral dissertation, Wageningen University; 2016.
Ritchie ME, Phipson B, Wu D, Hu Y, Law CW, Shi W, Smyth GK. “limma powers differential expression analyses for RNA-sequencing and microarray studies.” Nucleic Acids Res. 2015;43(7):e47. Limma package available at https://www.bioconductor.org/packages/release/bioc/html/limma.html.
Jovic K, Sterken MG, Grilli J, Bevers RPJ, Rodriguez M, Riksen JAG, et al. Temporal dynamics of gene expression in heat-stressed Caenorhabditis elegans. PLoS One. 2017;12(12):e0189445.
We thank Yu Nie for comments and assistance with assays. We also thank the CGC for worm strains.
This work is supported by a Leverhulme Trust Research Grant (RPG-2016-040) awarded to SH, and by BLS NWO (project no. 823.01.001). Some strains were provided by the CGC, which is funded by NIH Office of Research Infrastructure Programs (P40 OD010440). These funding bodies had no role in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.
Availability of data and materials
All data generated or analysed during this study are available, either in this published article and its supplementary information files or ArrayExpress (accession E-MTAB-6960).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Markers used in the genotyping of the α-synuclein ILs. (XLSX 22 kb)
A figure of gene expression markers across the genome. Determination of the α-synuclein introgression in the four wild isolate backgrounds. Marker genes with different expression between N2 and the wild isolates (All) were used to detect the α-synuclein introgression and N2 border regions. The marker genes missing in the α-synuclein lines indicate the N2 border regions and position of the α-synuclein introgression (aS Missing). Different genetic backgrounds are indicated by the different colours. The position(s) where all lines have missing markers show the likely α-synuclein locus, the extra missing markers on chromosome V show a possible extra introgression in SCH1931. (DOCX 201 kb)
A matrix containing the microarray data and significance values from the linear model (TXT 3344 kb)
Gene ontology analysis (GO) of genes identified as affected by Age, Genotype (GB), α-synuclein (aS) and the interaction between genotype and aS, and the details of genes shown in Fig. 1b including the p values and average log2 ratios with the mean per genotype for the 78 genes highly specific for aS relative to both Age and GB. (XLSX 246 kb)
Lifespan assays in α-synuclein ILs and relevant wild isolates. (TXT 2 kb)
Rates of matricidal hatching in α-synuclein ILs and relevant wild isolates. (TXT 574 bytes)