- Research article
- Open Access
Large-scale functional RNAi screen in C. elegans identifies genes that regulate the dysfunction of mutant polyglutamine neurons
BMC Genomics volume 13, Article number: 91 (2012)
A central goal in Huntington's disease (HD) research is to identify and prioritize candidate targets for neuroprotective intervention, which requires genome-scale information on the modifiers of early-stage neuron injury in HD.
Here, we performed a large-scale RNA interference screen in C. elegans strains that express N-terminal huntingtin (htt) in touch receptor neurons. These neurons control the response to light touch. Their function is strongly impaired by expanded polyglutamines (128Q) as shown by the nearly complete loss of touch response in adult animals, providing an in vivo model in which to manipulate the early phases of expanded-polyQ neurotoxicity. In total, 6034 genes were examined, revealing 662 gene inactivations that either reduce or aggravate defective touch response in 128Q animals. Several genes were previously implicated in HD or neurodegenerative disease, suggesting that this screen has effectively identified candidate targets for HD. Network-based analysis emphasized a subset of high-confidence modifier genes in pathways of interest in HD including metabolic, neurodevelopmental and pro-survival pathways. Finally, 49 modifiers of 128Q-neuron dysfunction that are dysregulated in the striatum of either R/2 or CHL2 HD mice, or both, were identified.
Collectively, these results highlight the relevance to HD pathogenesis, providing novel information on the potential therapeutic targets for neuroprotection in HD.
Huntington's disease (HD) is a dominantly-inherited disorder caused by expanded polyglutamine (polyQ) tracts in the N-terminal portion of huntingtin (htt) and characterized by striatal and cortical degeneration . While HD pathogenesis might involve a gain of toxic properties by mutant htt as well as a loss of normal htt function, several studies have emphasized a critical role of misfolded N-terminal fragments of mutant htt [2, 3] that are natural products of htt processing . Huntingtin is thought to have a large number of partner proteins involved in a variety of biological processes [5–8], suggesting that polyglutamine expansion in htt may alter several biological processes that are essential to cellular homeostasis and neuron survival. Consistent with this possibility, transcriptomic analyses have revealed that mutant/polyQ-expanded htt expression may change the expression of a large number of genes in the caudate nucleus of HD subjects  and in the brain of fragment, knock-in and full-length mouse models of HD . Given the complexity of alterations in cells and tissues expressing expanded polyQ/htt, gene perturbation screens are expected to shed light on the mechanisms that may be critical for cell response to mutant polyQ expression. Additionally, the integration of gene perturbation data with other types of 'omics data such as for example gene expression profiles in HD mice  might fulfill the need for target profiling and prioritisation in HD. This may be particularly important when comparing target gene activity for the regulation of neuronal dysfunction versus that of neuronal death. Before the occurrence of neuronal cell loss, neuron dysfunction is indeed thought to be an early-stage component of the pathogenic process in HD [3, 11]. However, neuronal cell dysfunction is not easily tractable in vivo. Interestingly, the genetic regulation of response to light touch by mechanosensory neurons in C. elegans has been the subject of detailed analysis . By targeting polyQ expression to C. elegans touch cells, we previously generated a model that recapitulates an early phase of expanded polyQ neurotoxicity, namely neuronal dysfunction in the absence of cell death [13, 14]. In these animals, exon 1 like htt fused to CFP was expressed in touch receptor neurons under the control of the promoter for the mec-3 gene (mec-3p) . This promoter is active from the early stages of embryonic development through adulthood . Young adults expressing normal polyQs (19Q) have a moderate loss of touch response with no or little change in touch cell morphology, whereas young adults expressing expanded polyQs (128Q) show a dramatic loss of touch response, a behavioral phenotype that is accompanied by a strong level of axonal dystrophy . Axonal dystrophy is followed by late-onset axonal degeneration as observed in older adult animals (unpublished data). The loss of touch response in expanded-polyQ nematodes can be manipulated by genetic means, as exemplified by the protective effect of sirtuin overexpression and aggravating effect of loss-of-function in the sir-2.1/SIRT1--daf-16/ FOXO cell survival pathway . These animals thus provide an in vivo model in which to screen for modifiers of neuron dysfunction induced by expanded polyQs. Here, we used 128Q nematodes in a large-scale functional RNA interference (RNAi) screen for modulation of response to touch. Additionally, we used 19Q nematodes to test a subset of the 128Q-confirmed hits for modulation of response to touch in 19Q nematodes. To this end, we generated polyQ strains with enhanced response of touch receptor neurons to double-standed RNA delivered by feeding. We then undertook a semi-directed screening approach to test a total of 6034 genes including 5128 genes tested at random and 906 genes that belong to major signaling pathways (16-19). This led us to identify 662 genes that either reduce or aggravate the loss of response to light touch when knocked-down by RNAi in 128Q nematodes. Screening in 19Q nematodes suggested that greater than 90% of these genes specifically modify 128Q cytotoxicity. In-depth data analysis emphasized a role for cell differentiation/survival, metabolic and neurodegenerative disease pathways in the modulation of expanded-polyQ cytotoxicity. Additionally, data analysis and comparisons to gene expression data from HD mice emphasized the potential value of this RNAi dataset for target discovery in HD.
C. elegans strains
Nematode strains were handled using standard methods . Integrated polyQ strains used in this study were previously described [17, 18]. Briefly, ID245 (igIs245[mec-3p::htt57-19Q::CFP; lin-15(+); mec-7p::YFP]) and ID1 (igIs1[mec-3p::htt57-128Q::CFP; lin-15(+); mec-7p::YFP]) stably co-express the first 57 amino acids of human huntingtin (htt) fused to CFP under the control of the promoter of the mec-3 gene and YFP under the control of the promoter of the mec-7 gene. We previously reported that ID245 animals show a moderate loss of response to light touch at the tail whereas ID1 animals show a strong loss of response to light touch at the tail, the latter phenotype accompanied by axonal swelling . The enhanced RNAi strain YY13 which carries the rrf-3(mg373) mutation  was provided by the Caenorhabditis Genetics Center, which is funded by the NIH National Center for Research Resources (NCRR). To construct polyQ strains with enhanced RNAi sensitivity in neurons, we crossed YY13 animals with either ID245 or ID1 animals to generate strain ID448 (igIs245[mec-3p::htt57-19Q::CFP;lin-15(+);mec-7p::YFP];rrf-3(pk1426)) and strain ID447 (igIs1[mec-3p::htt57-128Q::CFP;lin-15(+);mec-7p::YFP];rrf-3(pk1426)), respectively. The ID448 and ID447 strains were tested for sensitivity of touch receptor neurons to RNAi and subsequently used to screen the C. elegans ORF-RNAi library (Open Biosystems, AL) that includes clones from the C. elegans ORFeome library v1.1  as detailed below.
Bacterial clones from the C. elegans ORF-RNAi library were grown in 500 μl of LB culture medium supplemented with ampicillin (100 μg/ml) and tetracyclin (12.5 μg/ml), and they were induced by adding 10 μl of IPTG (12.5 mg/ml) 4 hours prior to incubation with nematodes. Synchronized L1 larvae were obtained by hypochlorite extraction using standard methods  and they were grown at 20°C on Petri dishes until they reach the young adult stage. The day of incubation in 96-well plates (D0), two young adult animals were incubated per well. Each of the wells contained 33.6 μl of bacterial culture and 16.6 μl of M9 media supplemented with cholesterol, ampicillin, tetracyclin, fungizone and IPTG. Incubation was performed at 20°C during 4 days until the animals reached adulthood and produce late L4 larvae. To ensure appropriate feeding of the nematodes in the well, 16.6 μl of bacterial cultures expressing RNAi clones were added at day 1. At day 4, animals were transferred to agar plates and allowed to recover for 15 minutes. RNAi clones producing lethality, larval arrest, egg-laying defect, developmental delay or any other abnormalities such as morphological abnormalities (2243 RNAi clones) were excluded from further analysis. For RNAi clones that did not alter nematode growth and morphology, the two P0 animals (old and dark nematodes) initially incubated in the well were manually removed from the plate and young adults (F1 progeny) were assayed for light touch response as previously described . Briefly, touch tests involved scoring for the response to light touch at the tail by using a fine hair. Touch test were performed by scoring 5 touches at the tail of the animal for 70-90 animals/RNAi clone/test. Ordinarily, wild-type animals will respond by backing away from the touch. The responses were recorded for every animal such that, for example, 4 responses out of 5 at the tail is given as 80% responsiveness, and the mean values for responsiveness were retained for further analysis of the data. A score for change in touch response (S) was calculated as [(Percent response - mean baseline)/mean baseline]. Using this formula, the maximally-achievable score is 4.55 (maximal suppression of 128Q cytotoxicity, i.e. 100% responsiveness) and the minimally-achievable score is -1 (maximal aggravation of 128Q cytotoxicity). The baseline in the assay was monitored across the entire time of the screening from the incubation of the animals treated with the empty vector L4440 in all of the 96-well plates assayed. The cumulated baseline over the entire time of the screening was observed to be stable (20.6 ± 3%) and the cumulated baseline/month was used to calculate S scores. Statistical analysis was performed using one-way ANOVA and the Tukey's Multiple Comparison Test. To decrease the risk for false positives, RNAi clones were scored as positive hits (suppressors or enhancers) if the S ± SD value was outside the baseline ± (SD*2.6). Initial hits (suppressors, enhancers) were tested for confirmation or invalidation of effect on touch response and the resulting data were retained for biological content analysis. All of the screening assays were performed by investigators that were trained to score touch response of 128Q animals in a similar way using blind assays.
Prior to performing the RNAi screen, the efficiency of gene inactivation in touch receptor neurons was tested by incubating ID447 L1 larvae in 96-well plates (1 larvae/well) with bacteria expressing either RNAi against GFP (plasmid pL4417) or empty vector (plasmid pL4440). Incubation was performed as described above until F1 young adults are produced. At the end of incubation, 70-90 animals/genotype were moved to agar plates and mounted on slides, and YFP levels in touch cells (mec-7p target cells) were measured at high (63X) magnification using a fluorescence microscope (Zeiss Axioplan Imaging II) and image analysis software (Image J). Additionally, YFP mRNA levels were measured by qRT-PCR from whole animals using standard methods. All of the nematode experiments were performed at least three times. Student's t-tests were used for statistical analysis of YFP mRNA levels and YFP expression levels.
The final screening results were initially classified using Gene Ontology (GO) annotations for Biological Process, Molecular Function and Cellular Component http://www.geneontology.org/. GO enrichment tests were performed using Ontologizer v 2.0 . Additionally, we used hypergeometric tests to test for gene set enrichment [23, 24]. The gene sets used for enrichment analysis were established using core components of (i) canonical pathways based on Wormbook http://www.wormbook.org, KEGG http://www.genome.ad.jp/kegg/ and Panther http://www.pantherdb.org/pathway/.
To analyze further the biological content of RNAi clones confirmed to modulate 128Q toxicity, we performed network-boosted data analysis. To this end, we subjected the RNAi data to spectral decomposition analysis using the integrated C. elegans network Wormnet v.2 . This large-coverage network contains gene interaction data that were established according to several types of evidences derived from multiple organisms, including several types of gene interactions based on gene co-expression (microarray measurements of the expression of C. elegans mRNAs), protein-protein interactions including genome-wide yeast two-hybrid interactions between C. elegans proteins and protein-protein interactions from the Worm Interactome database , genetic interactions from WormBase , human and fly 'interologs' i.e. human and fly protein interactions [29, 30], Yeast 'associalogs' i.e. functional interactions transferred from the yeast functional gene network , phylogenic profiles (comparative genomics linkages derived from the analysis of bacterial and archaeal genomes), gene neighbors and co-citation (literature-mined C. elegans gene associations). To increase the power of spectral decomposition analysis in detecting pathways, we also used a more homogenous network containing high confidence interactions, namely Core Wormnet  in which we omitted co-expression and co-citation links. To perform spectral decomposition analysis, we developed an enhanced version of the previously described  Fourier-like analysis. This approach is based on the spectral decomposition of gene activity profiles with respect to the connectivity of the reference network. S scores in the RNAi screen are considered as a 'signal' and the Fourier transformation decomposes the signal as a superposition of signals at different frequencies, with 40 successive attenuations performed . The aim of these transformations is to keep the low-frequency signals and discard the high frequency signals, considered as noise, in the context of gene connectivity. The logic behind this procedure is that highly inter-connected genes may be implicated in the same pathway and may fluctuate in the same direction . Starting from the raw RNAi screen data, we iteratively removed the high-frequency components 2.5% at a time, resulting in 40 increasingly smoothed versions of the data. Given a smoothed gene activity data, we considered modules as groups of connected genes with a signal > mean S-score + SD or < mean S-score-SD. Additionally, and to generate an even more relevant signal and robust information, we selected the modules that were present after the 4th attenuation and that showed a stable gene content (≥ 40% of the genes with a value over the threshold) for at least 2 successive attenuations. The resulting modules were filtered for redundant gene content and functional content of the remaining modules analyzed using Ontologizer v 2.0 . Additionally, we used hypergeometric tests to annotate for signaling pathways (gene lists extracted from Wormbook, KEGG and Panther), druggable genes , htt partners [7, 34] and genes involved in autophagy, mitochondrial function and synaptic activity (gene lists kindly provided by the CHDI Foundation). Information on conservation in humans was derived from Inparanoid clusters . In addition to tables included herein, the results of network-boosted data analysis are displayed in a freely-accessible database (see http://www.broca.inserm.fr/EHDN2/RNAiscreen) for the vizualisation of modules. This database allows the users to browse the results using GO/KEGG/Panther annotations and other annotations of interest, and it uses VisANT 3.5 for graph visualization . We used confidence graphs for display of modules by normalizing the statistical weights for functional interactions as available in Wormnet v 2.0 .
The touch receptor neurons of F1 128Q;rrf-3 animals are sensitive to RNAi
To construct polyQ strains with enhanced RNAi sensitivity in neurons, we crossed the enhanced RNAi strain YY13 which carries the rrf-3(mg373)II mutation  with integrated strains that co-express N-terminal (amino acids 1-57) htt under the control of mec3-p and YFP under the control of mec7-p. This resulted in strain ID448 (19Q;rrf-3) and strain ID447(128Q;rrf-3). The efficiency of gene inactivation in touch receptor neurons was tested by incubating strain ID447 L1 larvae in 96-well plates with bacteria expressing either RNAi against GFP (plasmid L4417) or empty vector (plasmid L4440) until F1 young adults are produced. Treatment with RNAi against GFP induced an overall reduction of YFP expression levels in mec7-p target cells as illustrated in Figure 1A/B. The level of YFP mRNA expression was strongly reduced (about 80%) by RNAi against GFP compared to empty vector (Figure 1C). Additionally, the level of YFP expression in the mechanosensory neurons that control the response to touch at the tail (PLM neurons) was reduced by RNAi against GFP compared to empty vector (Figure 1D). Collectively, these results indicated that the touch receptor neurons of F1 young adults are sensitive to RNAi. Consistent with previous reports , no or little reduction of YFP levels was detected in F0 young adults (data not shown).
A large-scale functional RNAi screen identifies 662 genes that modulate neuron dysfunction in expanded-polyglutamine nematodes
Using 128Q nematodes in the rrf-3 sensitized background, we then conducted a primary RNAi screen for the modulation of neuronal dysfunction. This screen was based on the manual quantification of response to light touch at the tail of F1 young adult animals. We screened 6260 RNAi clones equaling 6034 genes in the C. elegans genome (Figure 2). To decrease the risk for false positive, we applied a stringent threshold where the modification of touch response was considered significant if the score for touch response was outside the mean baseline ± SD*2.6 (Additional file 1: Figure S1). This screen identified 2243 RNAi clones (2211 genes) that caused toxicity phenotypes such as lethality, larval arrest, egg-laying defect, developmental delay and other abnormalities such as morphological abnormalities, 3101 RNAi clones (3054 genes) that showed no effect and 916 RNAi clones (769 genes) that modulated neuronal dysfunction (Figure 2, Additional file 2: Table S1; Additional file 3: Table S2). The latter class of RNAi clones was subjected to a second screen in 128Q nematodes, which led to the identification of 664 RNAi clones confirmed to modulate neuronal dysfunction. These clones corresponded to 662 genes, consisting of 399 genes that reduced neuronal dysfunction when knocked-down by RNAi and 263 genes that enhanced neuronal dysfunction when knocked-down by RNAi (Additional file 4: Table S3). As previously reported, 128Q nematodes show at least 80% decreased touch response at the tail with a strong level of axonal dystrophy, whereas 19Q nematodes show a 45-50% reduction of touch response at the tail without significant axonal dystrophy . Since our screen was based on the manual quantification of touch response in 70-90 animals/RNAi clone, which is considerably more labour-intensive and time-consuming compared to faster RNAi screens such as those looking for modifiers of mutant htt aggregation [38, 39], we evaluated whether the modifiers of 128Q-neuron dysfunction might also modify neuronal dysfunction in 19Q;rrf-3 nematodes by testing a subset of the 662 genes (245 genes choosen at random). Only 15/245 genes (6.1%) were found to modulate neuronal dysfunction in 19Q;rrf-3 nematodes when knocked-down by RNAi (Additional file 5: Table S4), suggesting that most of the 662 genes emphasized by screens in 128Q;rrf-3 nematodes specifically modify 128Q cytotoxicity. Genome-wide RNAi screens previously aimed at identifying modifiers of polyQ aggregation in C. elegans[40, 41], Consistent with the notion that aggregation and cytotoxicity can be genetically uncoupled , the overlap with our RNAi screen was small with 12 modifiers of 128Q-neuron dysfunction previously found to enhance the aggregation of 35Q-peptides when knocked-down by RNAi  and three other modifiers of 128Q-neuron dysfunction previously found to suppress the aggregation of 35Q-peptides when knocked-down by RNAi  (Additional file 6: Table S5).
In-depth data analysis emphasizes modifiers of 128Q-neuron dysfunction in pathways associated with metabolism, development, survival and neurodegenerative disease
The 662 genes associated with the modification of 128Q-neuron dysfunction appeared to encompass a variety of biological processes (cell death, protein folding, intracellular transport, metabolic processes, response to stress, stress-activated pathways) that may have a role in neurodegenerative disease pathogenesis as suggested by their functional classification using GO annotations (Figure 3, Additional file 7: Tables S6; Additional file 8: Table S7). To further assess the biological content corresponding to modifiers of 128Q-neuron dysfunction, we performed two types of analysis. First, we used gene lists describing major signaling pathways (see Materials and Methods) to perform gene set enrichment analyses. This analysis was performed for genes that either aggravate (herein referred to as 'protection-associated' or PA genes) or suppress (herein referred to as 'pathogenic process-associated' or PPA genes) 128Q-neuron dysfunction when knocked-down by RNAi. Pathway analysis for PA genes identified eight metabolic pathways such as for example those involved in glycolysis and amino acid metabolism (Table 1). Consistent with the notion that mutant htt expression may dysregulate a large number of genes and alter a variety of biological processes [3, 10], pathway analysis for PPA genes identified three main pathway categories, including pathways primarily involved in metabolism (glycolysis, amino acid metabolism), cytokine signaling/immune system (interferon-gamma signaling, interleukin signaling, Toll receptor signaling) and cellular differentiation and homeostasis (IGF pathway, angiogenesis, integrin signaling, p53 pathway by glucose deprivation, FGF signaling, EGF signaling, Toll receptor signaling) (Table 2). Interestingly, the latter category included the Alzheimer's disease (AD)-amyloid secretase pathway, pathways previously associated to neuroprotection in HD such as IGF signaling  and FGF signaling , pathways that could play a role in neurodegenerative disease pathogenesis such as EGF signaling , and pathways having a dual role in neurodegeneration and neuroprotection such as Toll receptor signaling . These observations illustrated the revelance of our RNAi dataset to neurodegeneration. Given that PPA genes might constitute preferred targets for therapeutic intervention, these results suggested that several of the targets for early-stage neuroprotection in HD may lie in the pathways that are important to neuronal cell differentiation and survival.
In contrast to metabolic pathways which may show a rather simple structure, pathways that respond to external stimuli to regulate gene expression may be highly inter-connected, and they may have highly-branched structures and several biological outcomes. Additionally, RNAi-based screens may generate noise due to off-target effects and variations in the degree of knockdown of target genes. Therefore, to map RNAi activity to functional clusters and to identify a more reliable set of genes from the RNAi hits, we used a network-based approach. To this end, we used a Fourier-like analysis  based on the reference network Wormnet . To enhance the discriminative power and reliability of this analysis, we developed and applied an optimized version of the Fourier-like analysis (see Materials and Methods). We analyzed the gene set for which a score for activity (S score) was available (genes that either modulate 128Q-neuron dysfunction or produce no effect when knocked-down by RNAi). The modules resulting from this analysis were annotated by using Gene Ontology, KEGG and Panther pathways and gene lists of interest (conserved genes, druggable genes, htt partners, genes involved in autophagy, mitochondrial function and synaptic activity). These analyses identified a more reliable set of 286 modifier genes as distributed in 137 modules for suppression of 128Q-neuron dysfunction by RNAi and 105 modules for aggravation of 128Q-neuron dysfunction by RNAi (Additional file 9: Table S8; Additional file 10: Table S9). Modules of variable sizes were obtained (2-50 genes) for either suppression or aggravation of 128Q-neuron dysfunction. Of note, small modules are equally significant compared to large modules as the size of the resulting modules was always driven by the stable and best balance between gene connectivity and gene activity during the spectral decomposition process (see Materials and Methods). The modules and their content can be visualized using a freely-accessible database (see http://www.broca.inserm.fr/EHDN2/RNAiscreen) that contains a query page and interactive graphical interface . Our network-based approach confirmed the trends observed in gene set enrichment analysis, namely the detection of pathways involved in metabolism, immune response and development. Additionally, it provided insights into the possible mechanisms of how RNAi hits potentially interact to modify 128Q-neuron dysfunction. We notably detected modules seeded by RNAi hits that were annotated for 'Huntington's disease', 'Parkinson's disease' and 'panthotenate kinase activity', the latter previously associated with neurodegeneration . More largely, we detected modules annotated for biological processes established or strongly suspected (see  for review of HD mechanisms) to have a role in neurodegenerative disease pathogenesis such as 'oxidative phosphorylation', 'cholesterol biosynthesis' and 'cell cycle'. Interestingly, several modules were annotated for cell differentiation and survival pathways such as 'Wnt signaling', a pathway suspected to have a role in HD pathogenesis , 'apoptosis signaling' and 'integrin signaling', suggesting that such pathways might be targeted for therapeutic intervention in HD. Annotations for metabolic pathways such as 'sucrose metabolism', 'pentose phosphate pathway' and 'cholesterol biosynthesis' were also detected, consistent with the notion that metabolic pathways may be targeted for therapeutic intervention in HD [3, 47]. Additionally, several modules were annotated for genes of interest in HD such as druggable genes (as notably observed in the modules for suppression of 128Q-neuron dysfunction by RNAi; see Additional file 9: Table S8), htt partners, and genes involved in autophagy, mitochondrial function and synaptic activity, which further illustrated the relevance of our dataset to neuroprotective target discovery in HD.
Modifiers of 128Q-neuron dysfunction in C. elegans are dysregulated in the striatum of either CHL2 knock-in mice or R6/2 transgenic mice
To further assess the relevance of our RNAi data to candidate target profiling and discovery in HD, we next compared our dataset to striatal gene expression data collected from mouse models of HD. Given that 128Q nematodes recapitulate progressive neuron dysfunction, which may occur in HD before cell loss, we compared our RNAi data to striatal gene expression data collected from the CHL2(Q150/Q150) knock-in mice, a model that may recapitulate progressive neuronal injury similar to that found early in the course of HD . To reach a more specific level of comparison, we also compared our data to striatal gene expression data collected from the R6/2 transgenic mice, a model that exhibits rapid phenotypic changes as induced by a short N-terminal fragment of mutant htt . Among the 662 modifiers of 128Q-neuron dysfunction in C. elegans, 239 genes were conserved in the mouse. Forty-nine genes were found at the intersection of the C. elegans RNAi data and gene dysregulation data in HD mice (Figure 4). Six of these 49 genes including Acy1, Sh3gl3, Snw1, Aars, Ppp5 and Gpi1 were oppositely dysregulated in the CHL2 mice compared to R6/2 mice (Figure 4), illustrating the differences that may exist between these two models at the molecular level. In contrast, genes that are involved in a variety of biological processes and that modulate 128Q-neuron dysfunction when knocked-down by RNAi were dysregulated in either the CHL2 or R6/2 mice (Figure 4), illustrating how data comparison using substractive analysis may be useful for target prioritisation. Notably, 12 genes that suppress 128Q-neuron dysfunction when knocked-down by RNAi in C. elegans were observed to be up-regulated in either the CHL2 or R6/2 mice, suggesting that inhibiting these genes might elicit neuroprotection in HD. One of them, the phenylalanine-4-hydroxylase Pah, appeared to be particularly interesting since it is a druggable gene that acts upstream of tyrosine-dopamine biosynthesis and since increased dopamine levels were previously associated with early HD phenotypes such as hyperkinesia .
The aim of our study was to collect information on the genes and pathways that may act as modifiers of the early stage of mutant htt toxicity in HD such as neuron dysfunction. To this end, we used an in vivo model to screen for gene inactivations that modulate the neuronal dysfunction induced by expanded-polyQ expression in C. elegans, and we used network-boosted analysis to enhance the identification of candidate pathways that may be further evaluated in mammalian systems. We identified 662 hits from the evaluation of 6034 genes. Compared to RNAi screens based on biochemical processes such as protein aggregation [38–40], the relatively high number of hits in our study is likely to reflect the multigenic nature of behavioral phenotypes such as response to touch rather than changes in polyQ-expanded htt expression levels. Previously-identified modulators of neuron dysfunction in expanded-polyQ nematodes do not change transgene expression levels [14, 18], further suggesting most RNAi effects are not due to changes in htt expression levels.
In contrast to 128Q-neuron dysfunction, neuronal dysfunction produced by normal polyQ expression in our trangenic nematodes is not directly relevant to HD pathogenesis. Microarray analysis of gene expression in touch receptor cells isolated by fluorescence-activated cell sorting indeed indicated that only 18 out of 2000 genes dysregulated by 128Q expression are also dysregulated by 19Q expression (submitted elsewhere), and several of the expanded-polyQ cytotoxicity that we previously identified showed no effect on normal polyQ cytotoxicity [14, 18]. Additionally, we observed that greater than 90% of the modifiers of 128Q-neuron dysfunction showed no effect in 19Q nematodes using a subset of 243 genes identified by RNAi screening in 128Q nematodes. Collectively, these observations suggest that the mechanisms of neuron dysfunction in 19Q nematodes  are different from those underlying the strong level of neuron dysfunction in 128Q nematodes. While specific modifiers of expanded-polyQ cytotoxicity might constitute preferred HD targets, one cannot exclude the possibility that non-specific modifiers of neuron dysfunction in polyQ nematodes may belong to the core machinery for neuron differentiation and survival, thus being of potential interest, at least from a screening perspective. This was for example illustrated by two N-acetylgalactosaminyltransferases (gly-9 and gly-3), a family of enzymes involved in neuron differentiation and survival .
A major feature of network-boosted analysis is to reduce data complexicity by highlighting individual pathways and possibly inter-connected pathways that are embedded into modules seeded with gene activity. In this respect, the network-boosted data analysis of our RNAi dataset was more instructive compared to the sole use of GO annotations or gene set enrichment analysis. While a role for a variety of pathways was suggested by gene set enrichment analysis, network-boosted analysis indeed emphasized a role for these pathways in a more comprehensive fashion, providing new insights of how the genes associated to these pathways may interact within modules. Furthermore, network-boosted analysis identified pathways not detected by gene set enrichment analysis. As indicated by the pathways embedded in the modules resulting from network-boosted analysis (see the on-line database at http://www.broca.inserm.fr/EHDN2/RNAiscreen), a role for neurodegenerative disease pathways and neurodevelopmental/cell survival pathways in the modulation of expanded-polyQ cytotoxicity was emphasized, which may provide guidance to the selection of candidate targets for neuroprotection in HD. The annotation of the modules resulting from network-boosted analysis identified clusters enriched in genes of interest such as druggable genes, htt partners and genes involved in autophagy, mitochondrial function and synaptic activity (all gene classes which can be queried using the on-line database), providing additional guidance to neuroprotective target discovery in HD. Here, it is important to note that the hypothesis underlying network-based approaches is that highly-interconnected genes might show similar biological activities. In this respect, the network-boosted analysis of our data is predictive of the genes that might decrease or increase 128Q-neuron dysfunction when knocked-down by RNAi, even though they showed no effect or borderline effects in the RNAi screen. While this may increase the number of candidate targets to be considered for neuroprotective intervention, it decreases the risk for noise associated with RNAi screens and favors the inclusion of potentially-interesting information for multi-model target discovery.
Multi-model target discovery is strongly expected to enhance the drug discovery process for neuroprotection in HD. A standard method for the validation of modifiers of neuronal dysfunction as instructed by simple models of HD pathogenesis relies on target validation in the mouse. Because it is not yet feasible to test in the mouse all genes supported as modifiers of neuronal decline in simple models, the outcome of our functional screen and other screens will be analyzed using data integration procedures in which the selection of the best targets for validation will rely on target expression levels and activity in multiple models of HD. One of the objectives of data integration is to assess the relevance of simple-model data to mouse models, and more largely to human neurodegenerative disease. In this respect, it is notable that several modifiers of 128Q-neuron dysfunction found in our screen were previously associated with the pathogenesis of neurodegenerative diseases such as HD or AD, suggesting that our C. elegans dataset is relevant to HD pathogenesis in higher order models of the disease. This is for example the case for four genes retained by network-boosted analysis including (i) dnj-1/ DNAJ, a gene that encodes a DNAJ chaperone , aggravates neuron dysfunction when knocked-down by RNAi and is down-regulated in the CHL2 mice, (ii) sel-8, a C. elegans gene that encodes a nuclear protein required for lin-12/NOTCH signaling, which may have neuron survival activities , and that aggravates 128Q-neuron dysfunction when knocked-down by RNAi, (iii) jnk-1/JNK, a gene that promote HD pathology in a rat model of the disease  and suppresses 128Q-neuron dysfunction when knocked-down by RNAi, and (iv) pnk-1/PanK4, a gene that encodes a panthotenate kinase and suppresses 128Q-neuron dysfunction when knocked-down by RNAi, which is consistent with the pathological role of panthotenate kinases in neurodegeneration . Additionally, the comparison of our data to gene dysregulation in the striatum of CHL2 knock-in mice versus that of R6/2 transgenic mice further indicated that our dataset contains information relevant to HD target selection as illustrated by the proposed inhibition of phenylalanine-4-hydroxylase as a candidate strategy to reduce early HD phenotypes that may be associated to a transient excess of dopamine .
In conclusion, our results provide new insights into the genes and pathways that might regulate the dysfunction and death of neurons expressing mutant polyQs, emphasizing a role for specific metabolic, developmental and pro-survival pathways. Several models of HD are being used to screen for modifiers of mutant polyQ toxicity and for genes dysregulated by mutant htt. In this context, the use of cross-species data integration for the analysis of our RNAi dataset and other datasets is anticipated to boost the prediction and prioritisation of neuroprotective targets in HD, and this will be investigated in additional studies.
Gil JM, Rego AC: Mechanisms of neurodegeneration in Huntington's disease. Eur J Neurosci. 2008, 27 (11): 2803-2820. 10.1111/j.1460-9568.2008.06310.x.
Graham RK, Deng Y, Slow EJ, Haigh B, Bissada N, Lu G, Pearson J, Shehadeh J, Bertram L, Murphy Z, et al: Cleavage at the caspase-6 site is required for neuronal dysfunction and degeneration due to mutant huntingtin. Cell. 2006, 125 (6): 1179-1191. 10.1016/j.cell.2006.04.026.
Zuccato C, Valenza M, Cattaneo E: Molecular mechanisms and potential therapeutical targets in Huntington's disease. Physiol Rev. 2010, 90 (3): 905-981. 10.1152/physrev.00041.2009.
Sun B, Fan W, Balciunas A, Cooper JK, Bitan G, Steavenson S, Denis PE, Young Y, Adler B, Daugherty L, et al: Polyglutamine repeat length-dependent proteolysis of huntingtin. Neurobiol Dis. 2002, 11 (1): 111-122. 10.1006/nbdi.2002.0539.
Holbert S, Denghien I, Kiechle T, Rosenblatt A, Wellington C, Hayden MR, Margolis RL, Ross CA, Dausset J, Ferrante RJ, et al: The Gln-Ala repeat transcriptional activator CA150 interacts with huntingtin: Neuropathologic and genetic evidence for a role in Huntington's disease pathogenesis. Proc Natl Acad Sci USA. 2001, 98 (4): 1811-1816. 10.1073/pnas.041566798.
Goehler H, Lalowski M, Stelzl U, Waelter S, Stroedicke M, Worm U, Droege A, Lindenberg KS, Knoblich M, Haenig C, et al: A protein interaction network links GIT1, an enhancer of huntingtin aggregation, to Huntington's disease. Mol Cell. 2004, 15 (6): 853-865. 10.1016/j.molcel.2004.09.016.
Kaltenbach LS, Romero E, Becklin RR, Chettier R, Bell R, Phansalkar A, Strand A, Torcassi C, Savage J, Hurlburt A, et al: Huntingtin interacting proteins are genetic modifiers of neurodegeneration. PLoS Genet. 2007, 3 (5): e82-10.1371/journal.pgen.0030082.
Holbert S, Dedeoglu A, Humbert S, Saudou F, Ferrante RJ, Neri C: Cdc42-interacting protein 4 binds to huntingtin: neuropathologic and biological evidence for a role in Huntington's disease. Proc Natl Acad Sci USA. 2003, 100 (5): 2712-2717. 10.1073/pnas.0437967100.
Hodges A, Strand AD, Aragaki AK, Kuhn A, Sengstag T, Hughes G, Elliston LA, Hartog C, Goldstein DR, Thu D, et al: Regional and cellular gene expression changes in human Huntington's disease brain. Hum Mol Genet. 2006, 15 (6): 965-977. 10.1093/hmg/ddl013.
Kuhn A, Goldstein DR, Hodges A, Strand AD, Sengstag T, Kooperberg C, Becanovic K, Pouladi MA, Sathasivam K, Cha JH, et al: Mutant huntingtin's effects on striatal gene expression in mice recapitulate changes observed in human Huntington's disease brain and do not differ with mutant huntingtin length or wild-type huntingtin dosage. Hum Mol Genet. 2007, 16 (15): 1845-1861. 10.1093/hmg/ddm133.
Smith R, Brundin P, Li JY: Synaptic dysfunction in Huntington's disease: a new perspective. Cell Mol Life Sci. 2005, 62 (17): 1901-1912. 10.1007/s00018-005-5084-5.
Chalfie M: Neurosensory mechanotransduction. Nat Rev Mol Cell Biol. 2009, 10 (1): 44-52. 10.1038/nrm2595.
Parker JA, Connolly JB, Wellington C, Hayden M, Dausset J, Neri C: Expanded polyglutamines in Caenorhabditis elegans cause axonal abnormalities and severe dysfunction of PLM mechanosensory neurons without cell death. Proc Natl Acad Sci USA. 2001, 98 (23): 13318-13323. 10.1073/pnas.231476398.
Parker JA, Arango M, Abderrahmane S, Lambert E, Tourette C, Catoire H, Neri C: Resveratrol rescues mutant polyglutamine cytotoxicity in nematode and mammalian neurons. Nat Genet. 2005, 37 (4): 349-350. 10.1038/ng1534.
Du H, Chalfie M: Genes regulating touch cell development in Caenorhabditis elegans. Genetics. 2001, 158 (1): 197-207.
Brenner S: The genetics of Caenorhabditis elegans. Genetics. 1974, 77 (1): 71-94.
Parker JA, Connolly JB, Wellington C, Hayden M, Dausset J, Neri C: Expanded polyglutamines in Caenorhabditis elegans cause axonal abnormalities and severe dysfunction of PLM mechanosensory neurons without cell death. Proc Natl Acad Sci USA. 2001, 98 (23): 13318-13323. 10.1073/pnas.231476398.
Parker JA, Metzler M, Georgiou J, Mage M, Roder JC, Rose AM, Hayden MR, Neri C: Huntingtin-interacting protein 1 influences worm and mouse presynaptic function and protects Caenorhabditis elegans neurons against mutant polyglutamine toxicity. J Neurosci. 2007, 27 (41): 11056-11064. 10.1523/JNEUROSCI.1941-07.2007.
Pavelec DM, Lachowiec J, Duchaine TF, Smith HE, Kennedy S: Requirement for the ERI/DICER complex in endogenous RNA interference and sperm development in Caenorhabditis elegans. Genetics. 2009, 183 (4): 1283-1295. 10.1534/genetics.109.108134.
Reboul J, Vaglio P, Rual JF, Lamesch P, Martinez M, Armstrong CM, Li S, Jacotot L, Bertin N, Janky R, et al: C. elegans ORFeome version 1.1: experimental verification of the genome annotation and resource for proteome-scale protein expression. Nat Genet. 2003, 34 (1): 35-41. 10.1038/ng1140.
Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, Davis AP, Dolinski K, Dwight SS, Eppig JT, et al: Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000, 25 (1): 25-29. 10.1038/75556.
Bauer S, Grossmann S, Vingron M, Robinson PN: Ontologizer 2.0--a multifunctional tool for GO term enrichment analysis and data exploration. Bioinformatics. 2008, 24 (14): 1650-1651. 10.1093/bioinformatics/btn250.
Mootha VK, Lindgren CM, Eriksson KF, Subramanian A, Sihag S, Lehar J, Puigserver P, Carlsson E, Ridderstrale M, Laurila E, et al: PGC-1alpha-responsive genes involved in oxidative phosphorylation are coordinately downregulated in human diabetes. Nat Genet. 2003, 34 (3): 267-273. 10.1038/ng1180.
Subramanian A, Tamayo P, Mootha VK, Mukherjee S, Ebert BL, Gillette MA, Paulovich A, Pomeroy SL, Golub TR, Lander ES, et al: Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc Natl Acad Sci USA. 2005, 102 (43): 15545-15550. 10.1073/pnas.0506580102.
Kanehisa M, Goto S: KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 2000, 28 (1): 27-30. 10.1093/nar/28.1.27.
Lee I, Lehner B, Crombie C, Wong W, Fraser AG, Marcotte EM: A single gene network accurately predicts phenotypic effects of gene perturbation in Caenorhabditis elegans. Nat Genet. 2008, 40 (2): 181-188. 10.1038/ng.2007.70.
Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, Vidalain PO, Han JD, Chesneau A, Hao T, et al: A map of the interactome network of the metazoan C. elegans. Science. 2004, 303 (5657): 540-543. 10.1126/science.1091403.
Chen N, Harris TW, Antoshechkin I, Bastiani C, Bieri T, Blasiar D, Bradnam K, Canaran P, Chan J, Chen CK: WormBase: a comprehensive data resource for Caenorhabditis biology and genomics. Nucleic Acids Res. 2005, D383-D389. 33 Database issue
Rual JF, Venkatesan K, Hao T, Hirozane-Kishikawa T, Dricot A, Li N, Berriz GF, Gibbons FD, Dreze M, Ayivi-Guedehoussou N, et al: Towards a proteome-scale map of the human protein-protein interaction network. Nature. 2005, 437 (7062): 1173-1178. 10.1038/nature04209.
Giot L, Bader JS, Brouwer C, Chaudhuri A, Kuang B, Li Y, Hao YL, Ooi CE, Godwin B, Vitols E, et al: A protein interaction map of Drosophila melanogaster. Science. 2003, 302 (5651): 1727-1736. 10.1126/science.1090289.
Lee I, Li Z, Marcotte EM: An improved, bias-reduced probabilistic functional gene network of baker's yeast, Saccharomyces cerevisiae. PLoS One. 2007, 2 (10): e988-10.1371/journal.pone.0000988.
Rapaport F, Zinovyev A, Dutreix M, Barillot E, Vert JP: Classification of microarray data using gene networks. BMC Bioinforma. 2007, 8: 35-10.1186/1471-2105-8-35.
Russ AP, Lampel S: The druggable genome: an update. Drug Discov Today. 2005, 10 (23-24): 1607-1610. 10.1016/S1359-6446(05)03666-4.
Stelzl U, Worm U, Lalowski M, Haenig C, Brembeck FH, Goehler H, Stroedicke M, Zenkner M, Schoenherr A, Koeppen S, et al: A human protein-protein interaction network: a resource for annotating the proteome. Cell. 2005, 122 (6): 957-968. 10.1016/j.cell.2005.08.029.
Remm M, Storm CE, Sonnhammer EL: Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Biol. 2001, 314 (5): 1041-1052. 10.1006/jmbi.2000.5197.
Hu Z, Hung JH, Wang Y, Chang YC, Huang CL, Huyck M, DeLisi C: VisANT 3.5: multi-scale network visualization, analysis and inference based on the gene ontology. Nucleic Acids Res. 2009, W115-W121. 37 Web Server
Calixto A, Chelur D, Topalidou I, Chen X, Chalfie M: Enhanced neuronal RNAi in C. elegans using SID-1. Nat Methods. 2010, 7 (7): 554-559. 10.1038/nmeth.1463.
Doumanis J, Wada K, Kino Y, Moore AW, Nukina N: RNAi screening in Drosophila cells identifies new modifiers of mutant huntingtin aggregation. PLoS One. 2009, 4 (9): e7275-10.1371/journal.pone.0007275.
Zhang S, Binari R, Zhou R, Perrimon N: A genome-wide RNAi screen for modifiers of aggregates formation by mutant huntingtin in Drosophila. Genetics. 2010, 184 (4): 1165-1179. 10.1534/genetics.109.112516.
Nollen EA, Garcia SM, van Haaften G, Kim S, Chavez A, Morimoto RI, Plasterk RH: Genome-wide RNA interference screen identifies previously undescribed regulators of polyglutamine aggregation. Proc Natl Acad Sci USA. 2004, 101 (17): 6403-6408. 10.1073/pnas.0307697101.
Silva MC, Fox S, Beam M, Thakkar H, Amaral MD, Morimoto RI: A genetic screening strategy identifies novel regulators of the proteostasis network. PLoS Genet. 2011, 7 (12): e1002438-10.1371/journal.pgen.1002438.
Jin K, LaFevre-Bernt M, Sun Y, Chen S, Gafni J, Crippen D, Logvinova A, Ross CA, Greenberg DA, Ellerby LM: FGF-2 promotes neurogenesis and neuroprotection and prolongs survival in a transgenic mouse model of Huntington's disease. Proc Natl Acad Sci USA. 2005, 102 (50): 18189-18194. 10.1073/pnas.0506375102.
Lievens JC, Rival T, Iche M, Chneiweiss H, Birman S: Expanded polyglutamine peptides disrupt EGF receptor signaling and glutamate transporter expression in Drosophila. Hum Mol Genet. 2005, 14 (5): 713-724. 10.1093/hmg/ddi067.
Nguyen MD, Julien JP, Rivest S: Innate immunity: the missing link in neuroprotection and neurodegeneration?. Nat Rev Neurosci. 2002, 3 (3): 216-227. 10.1038/nrn752.
Kruer MC, Hiken M, Gregory A, Malandrini A, Clark D, Hogarth P, Grafe M, Hayflick SJ, Woltjer RL: Novel histopathologic findings in molecularly-confirmed pantothenate kinase-associated neurodegeneration. Brain. 2011, 134 (Pt 4): 947-958.
Gines S, Ivanova E, Seong IS, Saura CA, MacDonald ME: Enhanced Akt signaling is an early pro-survival response that reflects N-methyl-D-aspartate receptor activation in Huntington's disease knock-in striatal cells. J Biol Chem. 2003, 278 (50): 50514-50522. 10.1074/jbc.M309348200.
Zwilling D, Huang SY, Sathyasaikumar KV, Notarangelo FM, Guidetti P, Wu HQ, Lee J, Truong J, Andrews-Zwilling Y, Hsieh EW, et al: Kynurenine 3-monooxygenase inhibition in blood ameliorates neurodegeneration. Cell. 2011, 145 (6): 863-874. 10.1016/j.cell.2011.05.020.
Lin CH, Tallaksen-Greene S, Chien WM, Cearley JA, Jackson WS, Crouse AB, Ren S, Li XJ, Albin RL, Detloff PJ: Neurological abnormalities in a knock-in mouse model of Huntington's disease. Hum Mol Genet. 2001, 10 (2): 137-144. 10.1093/hmg/10.2.137.
Raymond LA, Andre VM, Cepeda C, Gladding CM, Milnerwood AJ, Levine MS: Pathophysiology of Huntington's disease: time-dependent alterations in synaptic and receptor function. Neuroscience. 2011, 198: 252-273.
Yamashita T, Wu YP, Sandhoff R, Werth N, Mizukami H, Ellis JM, Dupree JL, Geyer R, Sandhoff K, Proia RL: Interruption of ganglioside synthesis produces central nervous system degeneration and altered axon-glial interactions. Proc Natl Acad Sci USA. 2005, 102 (8): 2725-2730. 10.1073/pnas.0407785102.
Hageman J, van Waarde MA, Zylicz A, Walerych D, Kampinga HH: The diverse members of the mammalian HSP70 machine show distinct chaperone-like activities. Biochem J. 2011, 435 (1): 127-142. 10.1042/BJ20101247.
Ables JL, Breunig JJ, Eisch AJ, Rakic P: Not(ch) just development: Notch signalling in the adult brain. Nat Rev Neurosci. 2011, 12 (5): 269-283. 10.1038/nrn3024.
Perrin V, Dufour N, Raoul C, Hassig R, Brouillet E, Aebischer P, Luthi-Carter R, Deglon N: Implication of the JNK pathway in a rat model of Huntington's disease. Exp Neurol. 2009, 215 (1): 191-200. 10.1016/j.expneurol.2008.10.008.
We are grateful to Cheryl Wellington and Michael Hayden (University of British Columbia, Canada) for providing constructs encoding htt with 128 Glns and Jean Philippe Vert for help with bioinformatics. We thank Hazare Nouari, David Bragado, Salima Abderrahmane, Sabine Cleofax, and Sidonie Vivequin for technical assistance, and Jim Wang and George Yohrling (CHDI Foundation) for lists of genes involved in autophagy, mitochondrial activity, and synaptic activity. J.A.P. was supported by a Young Researcher Award from Inserm. R.V. is supported by a Poste Vert fellowship from Inserm. This work was supported by Inserm, the Agence Nationale de la Recherche (ANR), the Fondation pour la Recherche Médicale (FRM), Paris, France, the Hereditary Disease Foundation (USA) and the European Huntington Disease Network (Euro-HD, Germany).
The authors declare that they have no competing interests.
FXL, LM and FP designed and performed data analysis and developed database tools. CB contributed to database developments. JPV contributed essential tools for data analysis. RVM, JAP and CT designed and performed C. elegans experiments. CN conceived and supervised experiments. FXL and CN wrote the manuscript. All authors read and approved the final manuscript.
François-Xavier Lejeune, Lilia Mesrob, Frédéric Parmentier contributed equally to this work.
Electronic supplementary material
Additional file 1: . Distribution of S scores for the primary screen in 128Q;rrf-3 nematodes. The graph shows the S scores for 4017 RNAi clones that did not elicit lethality or developmental abnormalities. The S score was calculated as [(Percent response - mean baseline)/mean baseline]. The maximally-achievable score is 4.55 (100% response to touch) and the smallest score is -1 (complete loss of touch response). The blue line indicate the mean baseline (S = 0). Dooted blue lines show the mean baseline ± SD*2.6. Clones with S scores falling outside the mean baseline ± SD*2.6 interval were retained for the secondary screen. (JPEG 200 KB)
About this article
Cite this article
Lejeune, F., Mesrob, L., Parmentier, F. et al. Large-scale functional RNAi screen in C. elegans identifies genes that regulate the dysfunction of mutant polyglutamine neurons. BMC Genomics 13, 91 (2012). https://doi.org/10.1186/1471-2164-13-91
- Neuronal Dysfunction
- Light Touch
- RNAi Screen
- RNAi Clone
- Axonal Dystrophy