- Research article
- Open Access
Polyglutamine toxicity in yeast induces metabolic alterations and mitochondrial defects
BMC Genomics volume 16, Article number: 662 (2015)
Protein aggregation and its pathological effects are the major cause of several neurodegenerative diseases. In Huntington’s disease an elongated stretch of polyglutamines within the protein Huntingtin leads to increased aggregation propensity. This induces cellular defects, culminating in neuronal loss, but the connection between aggregation and toxicity remains to be established.
To uncover cellular pathways relevant for intoxication we used genome-wide analyses in a yeast model system and identify fourteen genes that, if deleted, result in higher polyglutamine toxicity. Several of these genes, like UGO1, ATP15 and NFU1 encode mitochondrial proteins, implying that a challenged mitochondrial system may become dysfunctional during polyglutamine intoxication. We further employed microarrays to decipher the transcriptional response upon polyglutamine intoxication, which exposes an upregulation of genes involved in sulfur and iron metabolism and mitochondrial Fe-S cluster formation. Indeed, we find that in vivo iron concentrations are misbalanced and observe a reduction in the activity of the prominent Fe-S cluster containing protein aconitase. Like in other yeast strains with impaired mitochondria, non-fermentative growth is impossible after intoxication with the polyglutamine protein. NMR-based metabolic analyses reveal that mitochondrial metabolism is reduced, leading to accumulation of metabolic intermediates in polyglutamine-intoxicated cells.
These data show that damages to the mitochondrial system occur in polyglutamine intoxicated yeast cells and suggest an intricate connection between polyglutamine-induced toxicity, mitochondrial functionality and iron homeostasis in this model system.
During the last decades, protein misfolding and aggregation of certain proteins were found to play a major role in a variety of diseases, commonly called proteopathies [1–4]. For many of these diseases, age is a major risk factor. Some of them exhibit intracellular protein deposits, like in Parkinson’s disease, amyotrophic lateral sclerosis or Huntington’s disease, while others are diseases of the intra- and extracellular space, like systemic amyloidoses or Alzheimer’s disease . In Huntington’s disease, Huntingtin (Htt), a protein with yet unattributed function, exhibits a tendency to aggregate within cells if mutated . The age of disease onset and the aggregation propensity of Huntingtin are closely related to the length of a stretch of polyglutamine (polyQ) residues in its N-terminal domain . The threshold to transform into a pathogenic protein was shown to be around 35–45 consecutive glutamine residues [8, 9]. Importantly, the phenomenon of polyQ aggregation is not limited to Huntington’s disease, but approximately a dozen other diseases are similarly associated with extended glutamine stretches in specific proteins [9–11].
The origin of toxicity of these aggregation-prone proteins is still under debate. One hypothesis states that the presumed toxic species ties all available chaperone activity to the aggregation process, thereby interfering with the protein quality control system. This would lead to a propagation of folding defects onto other cellular proteins [12–14]. In mammalian cells more than 200 proteins, including several chaperones, were identified in intracellular aggregates . Experiments on model membranes also suggest that soluble oligomers of polyQ proteins are able to compromise the integrity of cellular membranes . In addition, apoptosis and a purely mechanical constriction of neuronal axonal trafficking by aggregates has been suggested to be responsible for the loss of neurons [17–19]. The diversity of observations concerning the mechanism of pathogenesis highlights that potentially different cellular processes are affected in parallel.
The simplest and genetically most accessible eukaryotic model organism is Saccharomyces cerevisiae, in which several systems to study polyQ aggregation have been established [20–22]. Here diverse morphological effects have been described in response to polyQ-expression, including DNA-fragmentation, damage to respiratory chain complexes, apoptosis like effects, spindle formation defects, mislocalization of septin proteins and altered regulation of the prion state [PSI +] [20, 23–26]. In our model system, using Q56-YFP, the polyQ-induced cellular arrest phenotype (pica) is evident from small colony growth, enlarged cell size and incomplete septin assemblies during G1 phase prior to budding . In contrast to other systems it is independent of the prion state of Rnq1 but dependent on the ploidy status of the yeast cell .
Comparing Q0-YFP and Q56-YFP expressing cells in this study we uncover further genetic interactors and analyze the transcriptional state of pica cells. In this process, we define genes, which help to reduce the toxicity of the polyQ protein and identify mitochondrial pathways, which likely are participating during establishment of toxicity.
Q56-YFP toxicity is suppressed by a set of mitochondrial genes
In order to study polyQ induced toxicity we used a yeast model system, which consists of three different constructs fusing either zero, 30 or 56 glutamine residues to YFP . In previous work the 56 amino acid stretch was found to be toxic, while the two other constructs were not harmful. Toxicity is evident from small colony growth after transformation with plasmids containing the constitutively expressed polyQ-encoding gene. To delineate the chain of events responsible for polyQ toxicity, we had performed a genome-wide screen of genomic deletion strains and identified yeast deletion mutants, which showed decreased toxicity in comparison to the wild type (WT) strain . Using the same approach we now focused on deletion strains, which show increased toxicity. From 5160 strains, we retrieved fourteen knock-out strains, some of which had entirely lost the capability to form small colonies even after 15 days of incubation at 30 °C (Fig. 1, Table 1). The presence of these non-essential genes thus is required for the residual growth after pQ56 transformation.
Most of the identified toxicity-suppressors participate in metabolic processes. Four of the 14 genes (ATP15, CEM1, NFU1 and UGO1) are directly localized in mitochondria, one of the suppressors of toxicity (IES1) is required for anaerobic growth and at least 4 others are part of metabolic pathways shared between the cytosol and mitochondria (HIS7, HOM6, MAP1, NCL1). Thus, these deleted genes optimize essential systems, which are functional under normal growth conditions but become dysfunctional at conditions like polyQ intoxication. The increased proportion of mitochondrial genes indicates that in particular the mitochondrial system may be challenged in Q56-YFP producing yeast.
Toxic and non-toxic polyQ stretches disturb the phosphate balance of the cells
To obtain more information on the condition of the pica yeast cells, we investigated the transcriptomic status of Q56-YFP intoxicated yeasts. We determined gene expression differences between intoxicated pQ56 and normally growing pQ0 transformed cells. We used four data sets to approach this question – Q0_3d (pQ0 after 3 days), Q0_2d, Q56_3d and Q56_4d – and obtained average relative expression changes for each gene (Additional file 1). We identified 76 genes, whose expression is reduced in pQ56 transformed cells to less than 33 % of the pQ0 transformed yeasts (Additional file 2).
To define and visualize transcriptional clusters down-regulated in Q56-YFP expressing yeasts, we clustered our hits based on co-regulation patterns from co-expression databases . In this way, most of the 76 genes down-regulated in the microarray experiments can be assembled into an interconnected network (Fig. 2a), as they apparently originate from two to three interconnected expression clusters. Beyond the initial hits we looked at further genes, which usually are part of these clusters: Using the SPELL database we automatically determined several co-regulated candidates with the highest connectivity and included them in the network of down-regulated genes (predicted genes are highlighted by a pink frame in Fig. 2a and listed in Additional file 3). This also helped to integrate further hits into the clusters of the network. Beyond that, these predicted candidates can be used to assess the predictive power of our network. Of the 50 candidates added to the network by the algorithm described in the methods section, 29 were in fact down-regulated to less than 60 % (log2 (Q0/Q56) > 0.5), while only 2 were up-regulated to the same extent (Fig. 2b), providing a sound statistical basis for the clusters formed in the network. It is evident from this network that transcripts related to phosphate metabolism and transport (PHO5, PHO89, PHO84, SPL2, PHM6, VTC1, VTC2, VTC3, VTC4) are down-regulated in Q56-YFP producing cells. This cluster had been observed before in a microarray study of yeasts expressing non-toxic polyQ proteins in liquid cultures . The clear overlap to this study highlights that specific parts of the response to polyQ proteins are remarkably robust even if entirely different toxicity levels and growth conditions are examined. Beyond that, another large cluster of genes, which is usually up-regulated during diauxic shift, is lower expressed in Q56-YFP intoxicated cells when compared to Q0-YFP (Fig. 2a). This suggests that the normally growing pQ0 transformed cells in our study are undergoing diauxic shift whereas the pQ56 transformed cells do not, which might coincide with the smaller colonies formed or other metabolic abnormalities.
To see whether these expression differences coincide with toxicity, we determined differentially expressed clusters in response to the non-toxic Q30-YFP sample (Q30_3d). Noteworthy, yeast cells synthesizing Q30-YFP, show no phenotype and are comparable to Q0-YFP expressing cells in terms of growth rate and colony size. Due to generally smaller differences in expression levels, we utilized a lower threshold – 0.6 log2 scales in both directions – to retrieve a reasonable number of hits for the analysis, while compromising on the p-value of the hits (Additional file 4). Correlations between the hits again were obtained from the co-regulation patterns in the SPELL database (Additional file 5A). Again further co-regulated candidates were derived from the SPELL database and incorporated into the network to evaluate the predictive power of the initial hit set (7 of these candidates were indeed down-regulated and 1 up-regulated of 16 predicted genes, Additional file 3 and Additional file 5B). Genes down-regulated in response to Q30-YFP form mostly one connected cluster. This cluster contains genes of phosphate metabolism and phosphate transport (SPL2, PHM6, PHO5, PHO84, PHO89, VTC1, VTC2, VTC3, VTC4). The consistent presence of this phosphate-cluster in all our experiments and in the previous study  is remarkable (Additional file 6A) and suggests that the presence of polyQ stretches results in alterations to phosphate regulation. Our data imply that this reaction may be a very sensitive marker in the reaction chain leading to the formation of the pica phenotype, but it is barely dependent on the toxicity level.
Polyphosphate levels are elevated in Q56-YFP expressing yeasts
We aimed at determining the metabolic impact of polyQ stretches to understand the cause for the differences in gene expression patterns. We thus analyzed whether dysregulation of the phosphate metabolism was evident in pica cells. To this end, we prepared suspensions of live Q0-YFP, Q30-YFP or Q56-YFP producing yeast cells and recorded 31P NMR spectra of these samples. Notably, distinct phosphate NMR signals could be detected in the cell suspensions (Fig. 3a; Additional file 7). Specifically, Peak 1 at 2 ppm was tentatively assigned to free phosphate, Peak 2 and 3 at −5.5 ppm and −24 ppm were tentatively assigned to polyphosphate [28–30], a storage form of phosphate, usually maintained by yeast cells in the vacuole and cytosol. The 31P NMR spectra of Q0-YFP and Q30-YFP producing yeast were virtually identical, whereas Q56-YFP expressing cells were different (Fig. 3a). Interestingly, the signal intensities for the inorganic phosphate were similar in all samples, implying that free phosphate levels are mostly comparable, whereas the signals originating from polyphosphate were increased by an approximate factor of 2 in pQ56 transformed yeast. This suggests that the phosphate balance is still maintained in Q30-YFP expressing yeasts despite the differences in induction of PHO-genes. In Q56-YFP expressing cells instead the phosphate homeostasis is misbalanced.
We aimed at determining the level of extracellular phosphatase activity by established assays in Q56-YFP, Q30-YFP and Q0-YFP yeast cells to see, whether the down-regulation of the extracellular phosphatases PHO5, PHO11 and PHO12 is observable on the protein level (Fig. 3b). Indeed Q56-YFP expressing yeast show reduced dephosphorylation activity compared to the two control strains, implying that the intoxicated yeast strains contain reduced amounts of phosphatases.
On this basis, it is tempting to speculate that the down-regulation of genes responsible for phosphate uptake and phosphate distribution in the cell could be related to the elevated polyphosphate levels in Q56-YFP synthesizing cells. Apparently in Q30-YFP the down-regulation of these genes enables the yeast cells to stabilize phosphate levels, while in Q56-YFP expressing yeasts the dysregulation manifests. Given that the transcriptional response already is observable in Q30-YFP producing yeasts, it is to be assumed that this part of the response is not related to toxicity.
Expression changes in toxic Q56-YFP show a specific upregulated response
We then looked at the genes upregulated in the response to toxic Q56-YFP. Here we obtained 72 genes, whose expression is more than 2.5 fold higher in pQ56- compared to pQ0-transformed cells in different biological samples (Additional file 8), many of which show p-values with high significance (p < 0.05). For these genes we likewise determined a network of expression clusters utilizing the strategy outlined before (Fig. 4a). This network includes a cluster of sulfur-regulated genes, many of which are part of sulfur uptake and early steps of methionine metabolism (SUL1, SUL2, STR3, MET1, MET16, MET8, HOM3, MET13, MET2, MET22, MET28, MET5, MET14, MET3). Furthermore, a cluster of iron-regulated genes is up-regulated. This includes three of the four iron-transporters encoded in the genome, ENB1, SIT1 and ARN2 and genes related to iron-starvation (FET3, FIT3, FIT2, TIS11 and VMR1). The sulfur and iron clusters are interconnected via several genes co-regulated with both clusters (Fig. 4a). We again determined co-regulated genes from the SPELL-database. Here, the connected and co-regulated candidate genes obtained from the SPELL database (highlighted with a blue frame) are indeed mostly up-regulated (26 of 33 predisctions are up-regulated, 1 down-regulated, Fig. 4b and Additional file 3). For many of them, the change in expression level was just below the 2.5 fold threshold for our initial hit list showing that many of the added genes are among the most strongly up-regulated genes within the data set (Additional file 3). Thus, from this microarray set we obtain upregulated expression networks with remarkable predictive power pointing to iron and methionine starvation in Q56-YFP producing yeast cells. Importantly, other stress responses, like the heat shock response, are not upregulated (see Additional file 9). Also nutrient supply appears sufficiently high despite the intoxication as evident from the down-regulation of clusters related to starvation (see Fig. 2a). Interestingly, the up-regulated genes do not overlap with earlier microarray studies on non-toxic polyQ . They thus may indicate a potential correlation to toxicity.
We again looked at the Q30-YFP control cells to analyze, whether a similar response can be observed here. In these cells, though, only very few genes are up-regulated significantly (Additional file 10, Additional file 11A). Nevertheless we attempted to construct a connected network, but only few of the genes could be connected and only very few further genes could be uncovered from the SPELL database due to a general lack of connectivity. The predictive power of clusters in this Q30-YFP/Q0-YFP data set is much less pronounced (2 up-regulated and 1 down-regulated of 8 predicted genes, Additional file 3, Additional file 11B). Genes involved in iron and sulfur regulation are not retrieved, if Q30-YFP is used to stress the cell. This suggests that the alteration of iron and sulfur homeostasis is specific for pica cells exhibiting the growth phenotype (Additional file 6B).
We tested the up-regulation of some genes from the iron- and sulfur-clusters by investigating genomic GFP-fusions of the respective proteins in the Q56-mCherry intoxicated and Q0-mCherry control background. Indeed for several of them, like Met10p, Met13p, Met22p, Met28p and Tmt1p we find a marked up-regulation observable by fluorescence microscopy of the GFP-fused proteins (Fig. 5), implying that the up-regulation is also evident at the level of protein concentrations within the cell.
Iron-homeostasis is affected in Q56-YFP producing yeasts
Having observed a strong up-regulation of iron and sulfur associated expression clusters in pica yeast we tested, whether a lack of iron ions causes this response. We determined the total intracellular free iron concentration in intact cells by electron paramagnetic resonance, (EPR) (Additional file 12). Membrane-permeable desferroxiamine (DFO) was used to complex free Fe3+ and oxidize intracellular Fe2+ [32, 33]. Surprisingly, Q56-YFP cells contain two-fold higher levels of intracellular iron compared to the control cells. Cells producing Q30-YFP instead show no alteration of the free iron level (Fig. 6a). The strong increase of intracellular iron levels in pQ56 transformed pica cells reflects a clearly disturbed intracellular iron balance. The higher levels of free iron though seem to contradict the observed up-regulation of iron-transporters. However, iron-regulation and its metabolism are both very complex and closely tied to the synthesis of Fe-S cluster proteins. In fact, up-regulation of iron-regulated genes as observed in Q56-YFP yeast is usually induced via the transcription factor Aft1, if Fe-S clusters are not produced to sufficient extent in mitochondria [34, 35].
To evaluate, whether pica yeasts show deficiencies in the generation of Fe-S cluster proteins in general, we tested the activity level of the prominent Fe-S cluster protein aconitase. While Q0-YFP and Q30-YFP producing yeast cells show very similar aconitase activities, we observed a strong reduction in aconitase activity in Q56-YFP yeasts (Fig. 6b). Hence, the higher level of free iron in the cell is apparently not sufficient to ensure the production of the Fe-S cluster containing enzyme aconitase. The reduced production of Fe-S cluster proteins in the mitochondria of pica yeast might thus cause the up-regulation of the iron-regulatory expression cluster observed.
Q56-YFP reduces mitochondrial carbon source utilization
We had found several indications that the mitochondrial system in pica cells is compromised and related to toxicity. To get direct information on the ability to perform non-fermentative metabolism, transformed yeast samples were placed on SD media containing 2 % sorbitol or 3 % glycerol instead of 2 % glucose [36, 37]. Both alternative carbon sources make yeasts more dependent on respiration compared to glucose and require functional mitochondria. As a control, we used galactose, which can be metabolized in a similar manner as glucose. When transformants of control plasmids pQ0 or pQ30 were plated on 3 % glycerol medium, growth was generally slower compared to glucose plates, but the colony patterns after four additional days of incubation resembled those of yeasts grown on glucose (Fig. 7a). Q56-YFP cells instead behaved differently. While few of the large colonies grew, which are resistant to the pica phenotype due to polyploidization or shortened polyQ stretches , the vast amount of small colonies is absent. Also, on 2 % sorbitol the growth suppression by Q56-YFP expression was stronger leading to the absence of small colonies (Fig. 7a). Usage of galactose as carbon source instead did not lead to absence of pica colonies, implying that the residual growth of these cells can be maintained, as long as carbon sources are used that allow sufficient energy production under fermentative conditions. To further study the importance of respiratory activity for growth of pica yeast we inhibited the respiratory chain by cultivating them on oligomycin. Indeed, we observed a loss in growth of pQ56 transformed yeasts whereas the pQ0 cells were still viable at the same oligomycin concentration (Fig. 7a).
To determine whether also differences in mitochondrial activity can be observed on the metabolic level, we tested the capability of pica yeasts to metabolize glucose. To this end we added [U-13C6] glucose to yeast suspensions of Q0-YFP or Q56-YFP cells and analyzed the samples by 13C-NMR spectroscopy at different time points after addition of the substrate (Fig. 7b). Within two hours control yeast produce metabolites of different kinds including ethanol, whereas the signal sets for glucose are diminished (Additional file 13, Additional file 14). The same signal patterns are observed with suspensions of pQ56 transformed yeast cells demonstrating that the conversion of glucose in the fermentative pathway in general is not affected by the presence of extended polyQ stretches (Fig. 7b). Looking into the weaker signals it becomes obvious that throughout the time course of the measurement small peak sets, which based on the chemical shifts and the J-couplings can be attributed to pyruvate, are observable only in pQ56 transformed suspensions (Fig. 7b and Additional file 14F and G). Further, a slight reduction in the carbonate peak at 160.3 ppm is evident in pQ56 transformed yeasts (Fig. 7b). This difference confirms an impaired respiratory activity of mitochondria as apparently the respiratory conversion of pyruvate to CO2 is reduced in the pica cells.
Q56-YFP intoxication results in decreased Cox4p levels
We aimed at determining, whether the altered metabolic pathways induce the growth arrest by energy shortage, or whether other mitochondrial functions, like generation of Fe-S cluster proteins, are more relevant for the pica phenotype. To this end, we investigated ATP-levels of Q56-YFP cells and normal growing Q0-YFP cells. Despite the reduced metabolic activity, ATP-levels themselves seem not to be altered in pica yeasts compared to control yeasts (Fig. 7c). We further investigated NADH levels, as these also could be misbalanced due to their relationship to mitochondrial metabolism. The detected differences in NADH level between pica and control yeasts are minor and not significant (Fig. 7d). Even though these values relate to full cell lysate and not isolated mitochondria, they imply that energy metabolism is still working in these yeasts. We finally aimed at testing by fluorescence microscopy whether proteins within the mitochondria are affected by the polyglutamine intoxication, which could also explain the deficiency in aconitase and the up-regulation of iron-regulated genes. Based on staining with the dye MitoTracker the mitochondrial system seems to be intact and the tubular network is formed (Fig. 7e). Two mitochondrial proteins were analyzed then by fluorescence microscopy, Om45 and Cox4. Both are slightly lower expressed in Q56-YFP yeasts based on our microarray analysis (−0.96 log2 and −0.75 log2 respectively). The mitochondrial membrane protein Om45 does not detectably change its cellular localization or its protein levels (Fig. 7e). Looking at Cox4p, a component of the respiratory chain, instead, the protein levels are strongly diminished in Q56-YFP expressing yeasts, implying that at the level of distinct proteins alterations to the mitochondrial system become observable (Fig. 7e). These data suggest that, while mitochondria generally appear functional in respect to energy metabolism, alterations within the mitochondria of pica cells become observable in connection to toxicity.
Elongated polyQ stretches have been reported as genetic cause for several neurodegenerative diseases, including Huntington’s disease [20, 22, 38, 39]. The reasons for this toxicity and the cellular pathways relevant for this process are under debate and many alterations have been reported, including DNA-fragmentation, apoptosis, spindle disorganization and involvement of the chaperone system or diverse prion proteins [20, 23, 24, 26]. Here we used a toxic polyQ system composed of zero, 30 or 56 glutamine residues to reveal affected cellular pathways based on an un-biased approach starting from a genome-wide screen of the EUROSCARF deletion strain library and a genome-wide assessment of expression changes upon intoxication.
Fourteen genomic deletion strains were uncovered that act synergistically with the polyQ induced phenotype. The knock-out of these genes prevents the residual growth, which in the wild type background is evident from many small colonies on agar plates. Interestingly, the majority of these genes are associated with mitochondrial functions, including the mitochondrial fusion mediator Ugo1p, the subunit of the F1F0-ATPase Atp15p and the protein Nfu1p, which is involved in Fe-S cluster synthesis. Several of these genes, including UGO1 and ATP15, have before been associated with the petite-phenotype, a small-colony formation phenotype, which can be intensified by cultivation on non-fermentative growth media. Given the similar appearance of polyQ-intoxicated pica yeasts, a related cause for these phenotypes could be speculated on, in particular as pica responds to the growth on non-fermentative carbon sources likewise with a general inability to form any colonies . Our screening approach had excluded strains with significant growth defects due to their genomic deletion, in particular when residual growth of small colonies could be observed after pQ56-transformation. We thus omitted some strains, which show the petite-phenotype, even though they were retained in the first round of screening. Two of these strains have deletions in MGM1 and FZO1. Given that UGO1, MGM1 and FZO1 are functionally related in mitochondrial fusion we reanalyzed their behavior in response to pQ56-transformation (Additional file 15). Indeed the differences in colony pattern between Q0-YFP expression and Q56-YFP expression suggest growth defects that are stronger than the additive effects of mutation and Q56-YFP toxicity.
The transcriptional profiling reported here additionally uncovers pathways with mitochondrial association and suggests the iron- and sulfur-metabolism as a potential target. Both pathways converge in the mitochondria during Fe-S cluster synthesis [34, 41]. We observe here the up-regulation of a transcriptional cluster, which is usually under the control of the transcription factor Aft1. This cluster includes the iron importers in the plasma membrane and several other proteins, which regulate the distribution and metabolism of iron ions in the cell. Generally this transcriptional cluster is induced by the lack of or shortage of Fe-S cluster proteins via the Fe-S containing sensor complex Grx [42, 43]. This reduction becomes obvious in our system by the reduced activity of aconitase, a very prominent Fe-S cluster protein in the cytosol. This shortage is not due to down-regulation of aconitase as evident from our genome-wide expression data, which show only slightly up-regulated ACO1 and slightly down-regulated ACO2 expression, but is apparently due to reduced levels of active protein (ACO1: −0.35 log2; ACO2: 0.52 log2 ). The increased levels of iron in the cell do not necessarily have to imply that Fe2+ is sufficiently present in the mitochondria for generation of Fe-S clusters. Fe3+ could also be present in the cytosol or even in complex with polyphosphate in the vacuole, all of which could lead to a shortage of mitochondrial iron in the presence of excessive Fe3+ . Likewise the accumulation of polyphosphate could be caused by the accumulation of iron in the cells given that these two are found in complex with each other on multiple occasions. It is interesting to note that Cox4p, a subunit of the cytochrome C oxidase complex (Complex IV) in the mitochondria and highly sensitive to polyQ intoxication, is part of a protein complex, which complexes iron ions . Thus proteome-wide data may indeed be necessary to establish, whether in particular iron-containing proteins are affected by the presence of polyQ containing aggregates. It is very interesting to further see, that the highly responsive heat shock and stress response network of chaperones does not react in polyQ intoxicated cells. This could indeed indicate that this part of the proteostasis network is not directly involved or at least that shortage of chaperones is not causative to the toxicity experienced in our yeast model system.
Several genome-wide studies had been performed in the past on similar aggregation systems in higher eukaryotes [46–50]. These uncovered a broad range of genetic and physical interactors. As such an RNAi screens in C. elegans uncovered metabolic influences, but in connection with the ascorbate and aldarate system . Additionally the protein homeostasis system was shown to suppress polyglutamine aggregation . Transcriptional changes in huntingtin mice and Huntington's patients point to genes in signal transduction , likewise do studies in mammalian cell culture based on RNA interference [48, 49]. Also studies in Drosophila melanogaster, employing huntingtin and ataxin as model systems find genes in signal transduction and cellular proteostasis .
Beyond that, the damage to the mitochondrial system has been observed in other aggregation model systems repeatedly [52, 53]. Also in subsets of patients suffering from neurodegenerative diseases the involvement of mitochondrial damage or the involvement of iron-metabolism has been encountered [54, 55]. In several cases the reduction in aconitase and the reduction in respiratory chain complexes have been observed [25, 56]. Another striking example linking the iron-metabolism to neurodegenerative disease is the protein frataxin, which itself is the iron-importer into mitochondria. This protein, if its polyQ stretch is expanded, tends to aggregate and result in neurodegenerative ataxias . While it is speculative to assume a common mechanism in generation of neurodegenerative defects in polyQ diseases, it emerges that alterations and damages to the mitochondrial system are highly relevant for the toxicity observable upon expression of polyQ proteins also in yeast.
Protein aggregation of elongated polyglutamine stretches induces cellular defects and death. Analyzing the transcriptional response in yeast we found two clusters connected to the iron and sulfur metabolism to be up-regulated in presence of extended toxic polyglutamine stretches. We report an accumulation of iron in vivo and a reduction in the activity of the prominent Fe-S cluster containing enzyme aconitase. In a genome wide approach, genes related to mitochondrial proteins were found to be crucial for cellular survival. Our data suggest an intricate connection between polyglutamine-induced toxicity, mitochondrial functionality and iron homeostasis in this model system.
Yeast cultivation and storage
S. cerevisiae strains were generally cultivated at 30 °C. Media were chosen according to the selection requirements. WT as well as knockout strains were grown on YPD-plates or in YPD liquid cultures . Strains transformed with plasmids pQ0, pQ30 or pQ56 were grown on appropriate SD media plates. All carbon sources were added at a concentration of 2 % (w/v) except for glycerol, which was added at 3 % (v/v). If not explicitly stated otherwise, the carbon source was glucose. The respiratory chain inhibition was performed by adding 8 mg/ml of an oligomycin isomers mixture (Merck Chemicals GmbH, Schwalbach, Germany) to the SD media plates.
The haploid Saccharomyces Genome Deletion Project library (EUROSCARF, Frankfurt, Germany) in the BY4741 background (MAT a; his3Δ1; leu2Δ0; met15Δ0; ura3Δ0) was used for screening. Screening conditions were as described previously . In the first round of screening, we retained all strains, growing weakly after pQ56-transformation. In subsequent rounds, we discarded strains, which show a strong growth defect by themselves due to the genomic deletion in cases, where residual growth of small pica colonies was observed, as under these conditions it is difficult to clearly distinguish between synergistic and additive effects. For the analysis of cellular structures affected by the presence of polyQ proteins, the respective strains of the Yeast GFP fusion collection in the background of EY0986 (MAT a; his3Δ1; leu2Δ0; met15Δ0; ura3Δ0) were deployed . All GFP-strains used in this manuscript are listed in Additional file 16.
Yeast transformations were performed using an adapted version of the simplified lithium acetate transformation method . The same transformation protocol was used for single transformations as well as for transformations in the 96-well format. 200 μl of a 2 day culture were spun down gently and resuspended in PLATE-solution (40 % PEG4000, 100 mM LiOAc, 10 mM Tris/HCl pH 7.5, 1 mM EDTA, 46 mM DTT). 5.1 μg/ml salmon carrier DNA were added together with 100 ng of plasmid DNA. The mixture was incubated at room temperature for 16 h and subsequently a heat shock was performed for 1 h at 42 °C. The transformation mixture was plated onto appropriate minimal medium and incubated at 30 °C. For 96-well transformations PLATE, salmon carrier DNA, DTT and plasmid were part of a premix to resuspend the yeast cell pellet.
Growth for GeneChip experiments
To analyze gene expression changes in response to polyQ proteins of different lengths, yeast samples were washed directly off the transformation plates. Plates bearing pQ56 transformants contain mostly small colonies. The few large colonies found on these plates were removed before harvesting the cells. As colony growth rates also correlate with carbon source consumption, yeasts were analyzed after three or four days on plates. Several plates had to be harvested for Q56-YFP after 3 days (Q56_3d) or Q56 after 4 days (Q56_4d) due to small colony size to obtain the same amount of biomass, while one plate was sufficient for the samples Q0_2d, Q0_3d and Q30_3d. Cells were pelleted for 10 min at 2500 x g and shock frozen in liquid N2. The isolation and enrichment of mRNA and further preparation including microarray analysis on Affymetrix GeneChip Yeast Genome 2.0  were carried out by the Kompetenzzentrum für Fluoreszente Bioanalytik (Regensburg, Germany). The GeneChip Yeast Genome 2.0 contains 5744 probe sets for genes of S. cerevisiae and 5021 probe stets for genes of Schizosaccharomyces pombe.
Hit selection from microarrays
The microarray data originated from three experimental sets. Q0_3d, Q30_3d and Q56_3d were processed together, likewise Q0_2d and Q56_4d. A further set of Q0_3d and Q30_3d was generated as replicate. The data provided by the Kompetenzzentrum für Fluoreszente Bioanalytik had already been normalized within one assay set, using the multi-chip analysis (RMA) algorithm [62, 63] and MAS5 . Generally the p-values for the individual hits were < 0.05, as long as the signal intensity of the MAS5 data was >15. All raw data sets are provided together with this manuscript and will be available to the public from www.richterlab.de.
To obtain relative expression differences between two samples, we used the MAS5 value for each gene, calculated the ratio (e.g. Q0_3d/Q56_3d) and sorted the full data set of 5900 probes regarding this value. The least-affected gene should then be at position 2950 with a quotient of 1. Thus the ratio column was divided by the median value obtained from position 2950. This approach was also used in the other data sets for normalization. All ratios were then converted to logarithmic values.
For the comparison of Q0-YFP and Q56-YFP samples, all possible combinations of the data sets were calculated and included. From the relative logarithmic expression differences we determined the average expression change and the standard deviation. The full set was then sorted according to this average relative expression change. Genes were included in the final hit list, if they showed expression differences of at least threefold (Q56 down versus Q0), 2.5 fold (Q56 up versus Q0) or 1.5 fold (Q30-YFP versus Q0-YFP, both directions). Expression changes below these values were not included in the hit list.
To define which signal intensities correspond to noise levels we used the S. pombe spots on the array as these probes will not detect specific genes and consequently define noise. This correlates well with the fluctuations in data with very low MAS5 values. Genes below this noise level (MAS5 values of 15) were generally not considered to be hits. During all these calculations no further genes were manually removed or edited.
All hits above noise level for either up- or down-regulation were included in the cluster analysis. To this end, we obtained the 20 highest ranking co-regulated genes for each of our hits from the SPELL database, which provides this information based on the analysis of more than 10,000 microarray data sets . Being listed together within these 21 genes (original gene plus 20 coregulators) was considered an incidence of co-regulation. A pairwise co-regulation matrix was built including the information from all the hits. Highly correlated genes were observed with up to 25 connections between each other. This matrix was used to draw a map with the open-source software CytoScape [65, 53]. The layout function “Edge-weighted Spring Embedded” was used for initial visualization of the clusters and the final maps were obtained by moving the nodes to prevent graphical overlap for better visualization.
In the likely case, that among the 20 co-regulated genes obtained from the SPELL database for a single hit genes showed up that were not in the hit list, a non-hit was included in the co-regulation matrix and if this non-hit was found co-regulated more often these co-regulated non-hits could obtain high levels of connectivity within our network. To test the significance of our network we picked those non-hit genes from our co-regulation matrix, which are best connected within the network and considered them to be part of the clusters in our network. If they are part of the clusters, they should also be regulated together with the cluster in our microarray data set, but possibly these genes were differentially expressed to an extent just below our inclusion threshold. To test, whether this is the case we obtained the experimental expression values for these “predicted further co-regulated genes” from our microarray data and tested, whether they are regulated in the appropriate direction in our experimental data set. If significantly more of these predicted genes are regulated together with the hit-based network than opposite to it, it would validate the predictive quality of the network and thus the validity of the clusters determined. For implementing this algorithm we used the integrated development environment Dev C++ 184.108.40.206 (Bloodshed Software, www.bloodshed.net). For all figures in this manuscript we used 20 co-regulated genes and included a “predicted co-regulator”, if the gene was found in correlation with another co-regulator above four times for the Q0/Q30-based network and above six times for the Q0/Q56-based network, which generally is better connected. Moreover we varied this parameter and the initial number of included hits to obtain information on the significance of this analysis.
Colony size analysis
Colony sizes of yeast strains after transformation with polyQ proteins were documented by taking photographs of whole petri dishes. A Canon EOS 60D digital camera mounted on a repro stand was used for this purpose.
In vivo NMR
31P-NMR was used to analyze phosphate levels in vivo. To this end yeast cells were washed off plates. We utilized cells directly derived from the transformation plates, where Q0-YFP and Q30-YFP could be obtained at very similar cell mass, while Q56-YFP expressing cells were harvested from several plates due to the toxicity of the transformation construct. Yeast cells were washed three times with 40 mM HEPES/KOH pH 7.0, 150 mM KCl to remove extracellular phosphate, which can originate directly from the agar plates. The cells were resuspended in 40 mM HEPES/KOH pH 7.0, 150 mM KCl to obtain OD595-levels of 300. The suspensions were then transferred to 5 mm NMR-tubes and pure D2O was added to 10 % final concentration in a final volume of 0.55 ml. 31P-NMR spectra were then measured at 25 °C using an AVANCE-III 500 NMR instrument (Bruker, Rheinstetten, Germany) equipped with a QNP cryo probe optimized for 31P-detection. Typically, 256 scans were assembled with a repetition time of 5 s. Spectra processing and analysis was performed with the software MestReNova 8.1. Prior to Fourier-transformation, the FID was multiplied with an exponential function (lb = 5). The spectra were calibrated to 0 ppm for external 85 % phosphoric acid.
Glucose metabolism was analyzed by 13C-NMR at 25 °C using a highly sensitive QNP cryo probe in the same NMR instrumentation. For this purpose, 20 mg/ml [U-13C6] glucose were added to the yeast cell suspensions and 1H-decoupled 13C NMR spectra were measured at intervals (0, 0.5, 1, 1.5, 2, 2.5 and 3 h after addition of the glucose tracer) using 256 scans and a repetition time of 5 s. Spectra were processed using the software MestReNova 8.1 using an exponential window function (lb = 1). Peak intensities were determined and the turnover rates were monitored on the basis of the glucose signal at 75.8 ppm and the ethanol signal intensity at 16.7 ppm. The kinetics of pyruvate accumulation were monitored based on the peak at 169.9 ppm and carbonate production based on the peak at 160.3 ppm.
Analysis of intracellular iron levels
Yeast cells were washed off plates and washed in HEPES buffer as described above. The cells were resuspended in pre-warmed 5 ml 40 mM HEPES/KOH pH 7.0, 150 mM KCl supplemented with 10 mM DETAPAC, pH 7.0, and 20 mM DFO, pH 8.0, and incubated at 30 °C for 15 min. DETAPAC blocks iron import, while DFO diffuses into cells and binds unincorporated free iron in an EPR-visible ferric form . Cells were then centrifuged at 4 °C, washed in ice-cold 40 mM HEPES/KOH pH 7.0, 150 mM KCl, and centrifuged again to pellet. Cell pellets were resuspended in a final volume of 300 μl of 40 mM HEPES/KOH pH 7.0, 150 mM KCl, 10 % glycerol, to gain final OD595-levels of 250–500. 200 μl were then transferred to a 4 mm quartz EPR tube, frozen on dry ice, and stored at −80 °C until assayed. EPR signals were measured with a Jeol JES FA 200 instrument. The spectrometer settings were as follows: temperature, −125 °C; microwave power, 10 mW; field center, 152 mT; field sweep, 100 mT; modulation amplitude, 0.25 mT; receiver gain, 1200; time constant, 0.3 s. A dewar in which samples were stored in liquid nitrogen before measurement was used to ensure constant temperature from sample to sample. Ferric chloride standards were prepared with DFO and EPR measurements were conducted as described above. The measured EPR signals were normalized according to cell density. Based on the spectra obtained from three independent samples we calculated the relative concentrations of free iron in the yeast cells.
Analysis of aconitase activity
Yeast cells were washed off plates as described in the EPR experiments. Aconitase activity was determined using a commercially available aconitase activity assay (Sigma-Aldrich, St. Louis, MO, USA). Cells were resuspended in aconitase assay buffer and disrupted. Mechanical disruption of 950 μl of cell suspension was achieved by adding 900 mg of 0.25–0.50 mm glass beads (Carl Roth, Karlsruhe, Germany) and subsequent shaking at 4 °C in 4 pulses of 2 min at 30 Hz in a MM400 bead mill (Retsch, Haan, Germany). Cell debris was removed by centrifugation (13300 rpm, 4 °C, 10 min) and yeast lysate was further treated according to the protocol. Isocitrate was processed in the assay to yield a colorimetrically detectable product at 450 nm and compared to a measured isocitrate standard curve. The aconitase activity was determined using the equation:
B [nmol], amount of isocitrate generated; SDF, sample dilution factor; T [min], time reaction incubated in minutes; V [ml], sample volume. The measured aconitase activities were normalized to the protein concentration determined by Bradford assay.
The ATP-levels were determined as described before [66, 67]. The protocol was slightly modified as harsher disruption methods had to be applied. In short: yeast cells were transformed, plated and incubated as described above. To enable analysis of the dominant colony population, large colonies were removed from pQ56 transformed plates. Plates were washed with 50 mM HEPES/KOH pH 7.8, 4 mM MgSO4 normalized to OD 5 and heated for 4 min at 95 °C. Mechanical disruption of 950 μl of cell suspension was achieved by adding 900 mg of 0.25–0.50 mm glass beads (Carl Roth, Karlsruhe, Germany) and subsequent shaking at 4 °C in 4 pulses of 2 min at 30 Hz in a MM400 bead mill (Retsch, Haan, Germany). After cell disruption the samples were transferred to ice and the total ATP content was analyzed by measuring the luciferase activity at 20 °C. 60 μl of the sample were added to 60 μl of a buffer containing 100 nM luciferase, 70 μM luciferin, 0.05 mg/ml BSA, 100 mM potassium phosphate buffer, pH 7.8, 25 mM glycylglycine, 0.02 mM EDTA, and bioluminescence was detected in a Tecan GENios™ microplate reader (Tecan Group, Männedorf, Switzerland).
Determination of total NAD+/NADH levels
Yeast cells were washed off plates with PBS as described above. NAD + /NADH levels were determined using a commercial available NAD/NADH Quantification Kit (Sigma-Aldrich, St. Louis, MO, USA). Briefly, cells were resuspended in NAD/NADH extraction buffer and the OD595-levels were determined before disrupting the cells using a bead mill. Cell debris was removed by centrifugation (17000 g, 4 °C) and the lysate further treated according to the protocol. The total NAD/NADH levels were quantified colorimetrically at 450 nm and calculated by using an NADH standard curve and corrected for the deployed amount of cells.
Detecting extracellular phosphatase activity
Extracellular phosphatase activity was analyzed as described before . Yeast cells were transformed, plated and incubated as described above. Plates were washed with no-phosphate buffer, spun down and washed twice with no-phosphate medium. The cells were incubated with an OD600 of 5 in 1 ml at 30 °C in no-phosphate medium. To determine the phosphatase activity at different time points 50 μl of the culture were added to 200 μl p-nitrophenylphosphate (5.62 mg/ml in 0.1 m sodium acetate, pH 4.2) and incubated at 20 °C for 15 min. By adding 200 μl of 10 % cold TCA the reaction was stopped and 200 μl were transferred to a new tube. Subsequently 200 μl of a saturated carbonate solution were added, mixed and spun down for 10 min at 3000 rpm.
The OD420 of 140 μl was recorded in quarz cuvettes in a Cary 50 UV/Vis spectrometer (Varian, Palo Alto, CA, USA). The extracellular phosphatase activity is given as the normalized OD420/OD600 x 1000.
For monitoring the fluorescence of GFP fused marker proteins an Axiovert 200 microscope (Zeiss Biomedical, Oberkochen, Germany) and the appropriate filter sets for GFP, mCherry fluorescence and DIC for transmitted light were deployed. The exposure time was not changed observing pQ0 and pQ56 transformed yeasts if not stated otherwise. For signal intensity quantification the mean intensity per pixel in the GFP channel was measured in areas densely populated with cells using ImageJ. The mean background intensity was similarly determined in empty image areas and subtracted from the mean intensity of fluorescent cells. The mean fluorescence was determined in triplicates.
Mitochondria were stained with MitoTracker rhodamine (Thermo Fisher, Waltham, MA, USA). Cells were washed off plates with PBS as described before and stained according to the manufacturer's protocol. In short, 250 nM MitoTracker was incubated for 15 min with pQ0 and pQ56 transformed yeast cells. Subsequently fluorescence was monitored without washing or further fixation of the cells.
Statistical analysis was performed using a double paired t-test with the OriginPro 8.6 software. Statistical analysis of microarray data was performed in Excel using the T.TEST function on the individual data sets.
Availability of supporting data
The data sets supporting the results of this article are included within the article and its additional files.
Stretch of polyglutamine residues
Nuclear magnetic resonance spectroscopy
Electron paramagnetic resonance spectroscopy
Diethylenetriamine pentaacetic acid
- Fe-S cluster:
Yeast extract peptone dextrose medium
- SD media:
Synthetic defined media
Yellow fluorescent protein
Polyglutamine-induced cellular arrest
Chiti F, Dobson CM. Protein misfolding, functional amyloid, and human disease. Annu Rev Biochem. 2006;75:333–66.
Taylor JP, Hardy J, Fischbeck KH. Toxic proteins in neurodegenerative disease. Science. 2002;296:1991–5.
Poirier M, Li H, Macosko J, Cai S, Amzel M, Ross C. Huntingtin spheroids and protofibrils as precursors in polyglutamine fibrilization. J Biol Chem. 2002;277:41032–7.
Ross C, Poirier M. Protein aggregation and neurodegenerative disease. Nat Med. 2004;10:S10–7.
Braun RJ, Büttner S, Ring J, Kroemer G, Madeo F. Nervous yeast: modeling neurotoxic cell death. Trends Biochem Sci. 2010;35:135–44.
Rubinsztein DC, Barton D, Davison B, Ferguson-Smith M. Analysis of the huntingtin gene reveals a trinucleotide-length polymorphism in the region of the gene that contains two CCG-rich stretches and a correlation between decreased age of onset of Huntington’s disease and CAG repeat number. Hum Mol Genet. 1993;2:1713–5.
Bates G. Huntingtin aggregation and toxicity in Huntington’s disease. The Lancet. 2003;361:1642–4.
Snell RG, MacMillan JC, Cheadle JP, Fenton I, Lazarou LP, Davies P, et al. Relationship between trinucleotide repeat expansions and phenotypic variation in Huntington’s disease. Nat Genet. 1993;4:393–7.
Zoghbi HY, Orr HT. Glutamine repeats and neurodegeneration. Annu Rev Neurosci. 2000;23:217–47.
Gatchel JR, Zoghbi HY. Diseases of unstable repeat expansion: mechanisms and common principles. Nat Rev Genet. 2005;6:743–55.
Williams AJ, Paulson HL. Polyglutamine neurodegeneration: protein misfolding revisited. Trends Neurosci. 2008;31:521–8.
Balch WE, Morimoto RI, Dillin A, Kelly JW. Adapting proteostasis for disease intervention. Science. 2008;319(February):916–20.
Bence NF, Sampat RM, Kopito RR. Impairment of the ubiquitin-proteasome system by protein aggregation. Science. 2001;292(May):1552–6.
Gidalevitz T, Ben-Zvi A, Ho KH, Brignull HR, Morimoto RI. Progressive disruption of cellular protein folding in models of polyglutamine diseases. Science. 2006;347(March):1471–5.
Olzscha H, Schermann SM, Woerner AC, Pinkert S, Hecht MH, Tartaglia GG, et al. Amyloid-like aggregates sequester numerous metastable proteins with essential cellular functions. Cell. 2011;144:67–78.
Lashuel H, Lansbury PT. Are amyloid diseases caused by protein aggregates that mimic bacterial pore-forming toxins? Q Rev Biophys. 2006;39:167–201.
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 U S A. 2001;98:13318–23.
Piccioni F, Pinton P, Simeoni S, Pozzi P, Fascio U, Vismara G, et al. Androgen receptor with elongated polyglutamine tract forms aggregates that alter axonal trafficking and mitochondrial distribution in motor neuronal processes. FASEB J. 2002;16:1418–20.
Büttner S, Habernig L, Broeskamp F, Ruli D, Vögtle FN, Vlachos M, et al. Endonuclease G mediates α-synuclein cytotoxicity during Parkinson’s disease. EMBO J. 2013;32:3041–54.
Kaiser CJO, Grötzinger SW, Eckl JM, Papsdorf K, Jordan S, Richter K. A network of genes connects polyglutamine toxicity to ploidy control in yeast. Nat Commun. 2013;4:1571.
Muchowski PJ, Schaffar G, Sittler A, Wanker EE, Hayer-hartl MK, Hartl FU. Hsp70 and Hsp40 chaperones can inhibit self-assembly of polyglutamine proteins into amyloid-like fibrils. Proc Natl Acad Sci U S A. 2000;97(July):7841–6.
Krobitsch S, Lindquist S. Aggregation of huntingtin in yeast varies with the length of the polyglutamine expansion and the expression of chaperone proteins. Proc Natl Acad Sci U S A. 2000;97:1589–94.
Treusch S, Lindquist S. An intrinsically disordered yeast prion arrests the cell cycle by sequestering a spindle pole body component. J Cell Biol. 2012;197:369–79.
Manogaran AL, Hong JY, Hufana J, Tyedmers J, Lindquist S, Liebman SW. Prion formation and polyglutamine aggregation are controlled by two classes of genes. PLoS Genet. 2011;7, e1001386.
Solans A, Zambrano A, Rodríguez M, Barrientos A. Cytotoxicity of a mutant huntingtin fragment in yeast involves early alterations in mitochondrial OXPHOS complexes II and III. Hum Mol Genet. 2006;15:3063–81.
Sokolov S, Pozniakovsky A, Bocharova N, Knorre D, Severin F. Expression of an expanded polyglutamine domain in yeast causes death with apoptotic markers. Biochim Biophys Acta. 2006;1757:660–6.
Hibbs M, Hess DC, Myers CL, Huttenhower C, Li K, Troyanskaya OG. Exploring the functional landscape of gene expression: directed search of large microarray compendia. Bioinformatics (Oxford, England). 2007;23:2692–9.
Hollander DJA, Ugurbil K, Brown TR, Bednar M, Redfield C, Shulman RG. Studies of anaerobic and aerobic glycolysis in sachharomyces cerevisiae. Biochemistry. 1986;25:203–11.
Hollander DJA, Ugurbil K, Brown TR, Shulman RG. Phosphorus-3 1 nuclear magnetic resonance studies of the effect of oxygen upon glycolysis in yeast. Biochemistry. 1981;20:5871–80.
Salhany JM, Yamane T, Shulman RG, Ogawa S. High resolution 31P nuclear magnetic resonance studies of intact yeast cells. Proc Natl Acad Sci U S A. 1975;72:4966–70.
Hughes RE, Lo RS, Davis C, Strand AD, Neal CL, Olson JM, et al. Altered transcription in yeast expressing expanded polyglutamine. Proc Natl Acad Sci U S A. 2001;98:13201–6.
Woodmansee AN, Imlay JA. Quantitation of intracellular free iron by electron paramagnetic resonance spectroscopy. Methods Enzymol. 2002;349:3–9.
Srinivasan C, Liba A, Imlay J, Valentine JS, Gralla EB. Yeast lacking superoxide dismutase(s) show elevated levels of “free iron” as measured by whole cell electron paramagnetic resonance. J Biol Chem. 2000;275:29187–92.
Lill R, Hoffmann B, Molik S, Pierik AJ, Rietzschel N, Stehling O, et al. The role of mitochondria in cellular iron-sulfur protein biogenesis and iron metabolism. Biochim Biophys Acta. 1823;2012:1491–508.
Castells-Roca L, Mühlenhoff U, Lill R, Herrero E, Bellí G. The oxidative stress response in yeast cells involves changes in the stability of Aft1 regulon mRNAs. Mol Microbiol. 2011;81:232–48.
Barakat R, Goubet I, Manon S, Berges T, Rosenfeld E. Unsuspected pyocyanin effect in yeast under anaerobiosis. MicrobiologyOpen. 2014;3:1–14.
Lai-zhang J, Mueller DM. Complementation of deletion mutants in the genes encoding the F 1 -ATPase by expression of the corresponding bovine subunits in yeast S. cerevisiae. Eur J Biochem. 2000;2418:2409–18.
Labbadia J, Morimoto RI. Huntington’s disease: underlying molecular mechanisms and emerging concepts. Trends Biochem Sci. 2013;38:378–85.
Orr HT, Zoghbi HY. Trinucleotide repeat disorders. Annu Rev Neurosci. 2007;30:575–621.
Contamine V, Picard M. Maintenance and integrity of the mitochondrial genome : a plethora of nuclear genes in the budding yeast. Microbiol Mol Biol Rev. 2000;64:281–315.
Netz DJ, Mascarenhas J, Stehling O, Pierik AJ, Lill R. Maturation of cytosolic and nuclear iron-sulfur proteins. Trends Cell Biol. 2014;24:303–12.
Outten CE, Albetel A-N. Iron sensing and regulation in Saccharomyces cerevisiae: Ironing out the mechanistic details. Curr Opin Microbiol. 2013;16:662–8.
Rutherford JC, Ojeda L, Balk J, Mühlenhoff U, Lill R, Winge DR. Activation of the iron regulon by the yeast Aft1/Aft2 transcription factors depends on mitochondrial but not cytosolic iron-sulfur protein biogenesis. J Biol Chem. 2005;280:10135–40.
Seguin A, Santos R, Pain D, Dancis A, Camadro J-M, Lesuisse E. Co-precipitation of phosphate and iron limits mitochondrial phosphate availability in Saccharomyces cerevisiae lacking the yeast frataxin homologue (YFH1). J Biol Chem. 2011;286:6071–9.
Maréchal A, Meunier B, Lee D, Orengo C, Rich PR. Yeast cytochrome c oxidase: a model system to study mitochondrial forms of the haem-copper oxidase superfamily. Biochim Biophys Acta. 1817;2012:620–8.
Lejeune F-X, Mesrob L, Parmentier F, Bicep C, Vazquez-Manrique RP, Parker JA, et al. Large-scale functional RNAi screen in C. elegans identifies genes that regulate the dysfunction of mutant polyglutamine neurons. BMC Genomics. 2012;13:91.
Becanovic K, Pouladi M, Lim RS, Kuhn A, Pavlidis P, Luthi-Carter R, et al. Transcriptional changes in Huntington disease identified using genome-wide expression profiling and cross-platform analysis. Hum Mol Genet. 2010;19:1438–52.
Kaltenbach LS, Romero E, Becklin RR, Chettier R, Bell R, Phansalkar A, et al. Huntingtin interacting proteins are genetic modifiers of neurodegeneration. PLoS Genet. 2007;3:e82.
Miller JP, Yates BE, Al-Ramahi I, Berman AE, Sanhueza M, Kim E, et al. A genome-scale RNA-interference screen identifies RRAS signaling as a pathologic feature of Huntington’s disease. PLoS Genet. 2012;8:e1003042.
Branco J, Al-Ramahi I, Ukani L, Pérez AM, Fernandez-Funez P, Rincón-Limas D, et al. Comparative analysis of genetic modifiers in Drosophila points to common and distinct mechanisms of pathogenesis among polyglutamine diseases. Hum Mol Genet. 2008;17:376–90.
Nollen E, Garcia SM, van Haaften G, Kim S, Chavez A, Morimoto RI, et al. Genome-wide RNA interference screen identifies previously undescribed regulators of polyglutamine aggregation. Proc Natl Acad Sci U S A. 2004;101:6403–8.
Zhu X, Beal MF, Wang X, Perry G, Smith MA. Mitochondria and neurodegenerative diseases. Journal of Alzheimers’s Disease. 2010;20:3233.
Tabrizi S, Workman J, Hart P, Mangiarini L, Mahal A, Bates G, et al. Mitochondrial dysfunction and free radical damage in the Huntington R6/2 transgenic mouse. Ann Neurol. 2000;47(January):80–6.
Oliveira JM. Nature and cause of mitochondrial dysfunction in Huntington’s disease: focusing on huntingtin and the striatum. J Neurochem. 2010;114:1–12.
Núñez MT, Urrutia P, Mena N, Aguirre P, Tapia V, Salazar J. Iron toxicity in neurodegeneration. Biometals. 2012;25:761–76.
Kim S-Y, Marekov L, Bubber P, Browne SE, Stavrovskaya I, Lee J, et al. Mitochondrial aconitase is a transglutaminase 2 substrate: transglutamination is a probable mechanism contributing to high-molecular-weight aggregates of aconitase and loss of aconitase activity in Huntington disease brain. Neurochem Res. 2005;30:1245–55.
Klockgether T, Evert B. Genes involved in hereditary ataxias. Trends Neurosci. 1998;21:413-18.
Jimenez A, Davies J. Expression of a transposable antibiotic resistance element in Sachharomyces. Nature. 1980;287:869–71.
Huh W, Falvo JV, Gerke LC, Carroll AS, Howson RW, Weissman JS, et al. Global analysis of protein localization in budding yeast. Nature. 2003;425(October):686–91.
Elble R. A simple and efficient procedure for transformation of yeasts. Biotechniques. 1992;13:18–20.
Wood V, Gwilliam R, Rajandream M, Lyne M, Lyne R, Stewart A, et al. The genome sequence of Schizosaccharomyces pombe. Nature. 2002;415(6874):871–80.
Speed TP, Bolstad BM, Irizarry RA, Astrand M. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003;19:185–93.
Irizarry R. Summaries of Affymetrix GeneChip probe level data. Nucleic Acids Res. 2003;31:15e–15.
Pepper SD, Saunders EK, Edwards LE, Wilson CL, Miller CJ. The utility of MAS5 expression summary and detection call algorithms. BMC bioinformatics. 2007;8:273.
Saito R, Smoot ME, Ono K, Ruscheinski J, Wang P-L, Lotia S, et al. A travel guide to Cytoscape plugins. Nat Methods. 2012;9:1069–76.
Yang N-C, Ho W-M, Chen Y-H, Hu M-L. A convenient one-step extraction of cellular ATP using boiling water for the luciferin–luciferase assay of ATP. Anal Biochem. 2002;306:323–7.
Gray MJ, Wholey W-Y, Wagner NO, Cremers CM, Mueller-Schickert A, Hock NT, et al. Polyphosphate is a primordial chaperone. Mol Cell. 2014;53:689–99.
Huang S, O’Shea EK. A systematic high-throughput screen of a yeast deletion collection for mutants defective in PHO5 regulation. Genetics. 2005;169:1859–71.
The authors are grateful to Prof. Dr. Klaus Köhler for access to the EPR device at his laboratory (Fachgebiet Anorganische Chemie, Technische Universität München). We further thank Christopher Stratil for critically reading the manuscript.
The authors declare that they have no competing interests.
Conceived and designed the experiments: KP, CJOK, AD, KR. Performed the experiments: KP, CJOK, AD, SWG, CH. Analyzed the data: KP, CJOK, AD, WE, KR. Wrote the paper: KP, CJOK, AD, KR. All authors read and approved the final manuscript.
Katharina Papsdorf and Christoph J. O. Kaiser contributed equally to this work.
Comparison between averaged Q 0 /Q 56 -YFP expression levels of microarrays. A) Relative expression differences with standard deviation. B) Relative expression differences with threshold levels for the inclusion in the initial hits list. C) Determination of noise levels for the exclusion of genes from the initial hits list. S. pombe MAS5 values from the microarray set Q0_2d. The blacsk line represents the applied noise threshold D) Correlation between Q0_3d and Q30_3d MAS5 signals. The lines represent the applied noise threshold. (PDF 144 kb)
Genes with reduced expression in Q 56 -YFP colonies versus Q 0 -YFP colonies. The table summarizes the averaged expression differences from the comparison of four data sets (Q0_2d, Q0_3d, Q56_3d and Q56_4d). Each combination was evaluated and the average expression difference was obtained. Standard deviation was calculated from these values. p-values were obtained by using the t-test. Hits with a p-value greater than 0.05 are indicated in grey. (DOCX 23 kb)
Predicted co-regulated genes. All predicted co-regulated genes integrated into the networks are listed together with their log2 regulation. (DOCX 19 kb)
Genes with reduced expression in Q 30 -YFP colonies versus Q 0 -YFP colonies. The table summarizes the expression differences of two data sets (Q0_3d, Q30_3d). Standard deviation and p-values were obtained as described above. Hits with a p-value greater than 0.05 are indicated in grey. (DOCX 16 kb)
Reduced expression between Q 30 -YFP and Q 0 -YFP yeast cells. A) Network of genes, which are lower expressed in Q30-YFP producing cells compared to Q0-YFP producing cells. Genes are colored in accordance to their log differences (100 % green: log2 > 1, 75 % green: 1 > log2 > 0.75, 50 % green: 0.75 > log2 > 0.5, 25 % green: 0.5 > log2 > 0.25, white: 0.25 > log2 > -0.25, 25 % red: -0.25 > log2 > -0.5, 50 % red: -0.5 > log2 > -0.75, 75 % red: -0.75 > log2 > -1, 100 % green: -1 > log2). Pink frames highlight the predicted co-regulated genes. B) Statistical analysis of predicted co-regulators within Figure S2A. The line shows the number of genes within the respective category, while the vertical bar chart indicates the percentage of these genes predicted by the SPELL-correlations. (PDF 180 kb)
Expression differences between Q 56 -YFP and Q 30 -YFP yeast cells. Cluster regulation in four possible combinations of the Q0/Q56 and Q0/Q30 experimental data sets. A) Consistency of the down-regulation of the genes in phosphate uptake and transport. Each cluster represents a comparison of different Q0 and Q30 (lower row) and Q56 experiments (upper row). B) Consistency of the up-regulation of the iron and the sulfur cluster. Each cluster represents a comparison of different Q0 and Q30 (lower row) and Q56 experiments (upper row). (PDF 4826 kb)
Peak assignments in the 31 P-NMR spectra. Peaks as observed in the 31P NMR spectra were assigned with the help of databases and literature. (DOCX 14 kb)
Genes with increased expression in Q 56 -YFP colonies versus Q 0 -YFP colonies. The table summarizes the averaged expression differences from the comparison of four data sets (Q0_2d, Q0_3d, Q56_3d and Q56_4d). Each combination was evaluated and the average expression difference was obtained. Standard deviation was calculated from these values. p-values were obtained by using the t-test. (DOCX 22 kb)
Genes with unaltered expression between Q 56 -YFP and Q 0 -YFP yeast cells. Network of genes regulated in Q56-YFP expressing cells compared to Q0-YFP expressing cells using the same color code as Fig. 2a and 4a. Genes are colored accordance to their log differences (100 % green: log2 > 2, 75 % green: 2 > log2 > 1.5, 50 % green: 1.5 > log2 > 1, 25 % green: 1 > log2 > 0.5, white: 0.5 > log2 > -0.5, 25 % red: -0.5 > log2 > -1, 50 % red: -1 > log2 > -1.5, 75 % red: -1.5 > log2 > -2, 100 % green: -2 > log2). Genes are clustered according to their GO-Terms. (PDF 307 kb)
Genes with increased expression in Q 30 -YFP colonies versus Q 0 -YFP colonies. The table summarizes the expression differences of two data sets (Q0_3d, Q30_3d). Standard deviation and p-values were obtained as described above. (DOCX 16 kb)
Enhanced expression between Q 30 -YFP and Q 0 -YFP yeast cells. A) Network of genes, which are higher expressed in Q30-YFP expressing cells compared to Q0-YFP expressing cells. Genes are colored in accordance to their log differences (100 % green: log2 > 1, 75 % green: 1 > log2 > 0.75, 50 % green: 0.75 > log2 > 0.5, 25 % green: 0.5 > log2 > 0.25, white: 0.25 > log2 > -0.25, 25 % red: -0.25 > log2 > -0.5, 50 % red: -0.5 > log2 > -0.75, 75 % red: -0.75 > log2 > -1, 100 % green: -1 > log2). Blue frames highlight the predicted co-regulated genes. B) Statistical analysis of predicted co-regulators within Additional file 6A. The line shows the number of genes within the respective category, while the vertical bar chart indicates the percentage of these genes predicted by the SPELL-correlations. (PDF 1709 kb)
EPR-measurement of yeast cells. A representative EPR measurement is shown. Detection of iron-content in Q0-YFP, Q30-YFP and Q56-YFP producing yeast cells. For detailed information see methods section. (PDF 363 kb)
Peak assignments in the 13 C-NMR spectra. Peaks as observed in the 13C NMR spectra were assigned with the help of databases and literature. (DOCX 14 kb)
Metabolism of [U- 13 C 6 ]-glucose in polyQ intoxicated yeast monitored by 13 C-NMR. A) NMR-analysis of [U-13C6]-glucose in the absence of yeast. B) NMR-analysis of [U-13C6]-glucose in the presence of Q0-YFP expressing yeasts or C) in the presence of Q56-YFP expressing yeasts. D) Kinetics of the metabolite production based on the peak at 124 ppm by pQ0 transformed yeasts (blue circles), pQ56 transformed yeasts (red square) or no yeasts present (black triangle). E) Kinetics of [U-13C6]-glucose consumption based on the peak at 60 ppm by pQ0 transformed yeasts (blue circles), pQ56 transformed yeasts (red square) or no yeasts present (black triangle). Without yeasts [U-13C6]-glucose is not metabolized. The chemical shift positions are assigned in Additional file 13. F) Assigned spectra in the range of -20 ppm and 230 ppm after 2 h of glucose metabolism for pQ0- and pQ56- transformed yeast cells. G) Detailed view on pyruvate (Pyr) and carbonate signals. Assigned spectra in the range of 156 ppm to 210 ppm expose the smaller peaks. Peak labeling corresponds to Additional file 13. (PDF 652 kb)
Deletion of mitochondrial fusion mediators fzo1 Δ and mgm1 Δ. The deletion of mitochondrial fusion mediators fzo1Δ and mgm1Δ results in loss of colony growth upon Q56-YFP expression. BY4741, fzo1Δ and mgm1Δ strains after transformation of either Q0-YFP (upper row) or Q56-YFP (lower row). The respective deletion is indicated on the top. Growth was documented after 4 days of incubation at 30 °C. Scale bar represents 10 mm. (PDF 14378 kb)
GFP-strains used in this manuscript. For the analysis of cellular structures the respective strains of the Yeast GFP fusion collection in the background of EY0986 (MAT a; his3Δ1; leu2Δ0; met15Δ0; ura3Δ0) were deployed and transformed with either p Q0 or p Q56 . Sul1 was not investigated further as it only shows very faint fluorescence. (DOCX 14 kb)
About this article
Cite this article
Papsdorf, K., Kaiser, C.J.O., Drazic, A. et al. Polyglutamine toxicity in yeast induces metabolic alterations and mitochondrial defects. BMC Genomics 16, 662 (2015). https://doi.org/10.1186/s12864-015-1831-7
- Iron-sulfur cluster
- Neurodegenerative disease