Skip to main content

Protein signatures of oxidative stress response in a patient specific cell line model for autism

Abstract

Background

Known genetic variants can account for 10% to 20% of all cases with autism spectrum disorders (ASD). Overlapping cellular pathomechanisms common to neurons of the central nervous system (CNS) and in tissues of peripheral organs, such as immune dysregulation, oxidative stress and dysfunctions in mitochondrial and protein synthesis metabolism, were suggested to support the wide spectrum of ASD on unifying disease phenotype. Here, we studied in patient-derived lymphoblastoid cell lines (LCLs) how an ASD-specific mutation in ribosomal protein RPL10 (RPL10[H213Q]) generates a distinct protein signature. We compared the RPL10[H213Q] expression pattern to expression patterns derived from unrelated ASD patients without RPL10[H213Q] mutation. In addition, a yeast rpl10 deficiency model served in a proof-of-principle study to test for alterations in protein patterns in response to oxidative stress.

Methods

Protein extracts of LCLs from patients, relatives and controls, as well as diploid yeast cells hemizygous for rpl10, were subjected to two-dimensional gel electrophoresis and differentially regulated spots were identified by mass spectrometry. Subsequently, Gene Ontology database (GO)-term enrichment and network analysis was performed to map the identified proteins into cellular pathways.

Results

The protein signature generated by RPL10[H213Q] is a functionally related subset of the ASD-specific protein signature, sharing redox-sensitive elements in energy-, protein- and redox-metabolism. In yeast, rpl10 deficiency generates a specific protein signature, harboring components of pathways identified in both the RPL10[H213Q] subjects’ and the ASD patients’ set. Importantly, the rpl10 deficiency signature is a subset of the signature resulting from response of wild-type yeast to oxidative stress.

Conclusions

Redox-sensitive protein signatures mapping into cellular pathways with pathophysiology in ASD have been identified in both LCLs carrying the ASD-specific mutation RPL10[H213Q] and LCLs from ASD patients without this mutation. At pathway levels, this redox-sensitive protein signature has also been identified in a yeast rpl10 deficiency and an oxidative stress model. These observations point to a common molecular pathomechanism in ASD, characterized in our study by dysregulation of redox balance. Importantly, this can be triggered by the known ASD-RPL10[H213Q] mutation or by yet unknown mutations of the ASD cohort that act upstream of RPL10 in differential expression of redox-sensitive proteins.

Background

More than 1% of the general population display ASDs (autism spectrum disorders) which encompass early childhood autism, Asperger-Syndrome and pervasive-developmental-disorders not-otherwise-specified [1]. ASD patients are characterized by a clinically assessed triad of symptoms of impairments in reciprocal social interaction and communication, and of restricted interests and/or repetitive behaviors [2]. The heritability of ASD is estimated to range between 37% to 80% for the disorder [3–5] depending on sample size and structure of the studies. Genetic studies indicate that a combination of common but lowly penetrant and rare but highly penetrant variants can explain a large proportion of the susceptibility for ASD [6, 7]. Although single gene disorders with increased risk for ASD, ASD-specific mutations and chromosomal copy number variations are estimated to explain up to 20% of the liability for the disorder [8], the majority of cases, as well as the underlying biological pathomechanisms, remain elusive. This is in part due to the very wide difference of penetrance of each one of the three autistic core impairments recognizable in autistic patients [9]. Therefore, autism research now has a focus on the investigation of common cellular pathophysiologies, within which individual, patient specific malfunctions would generate a personalized autistic phenotype. So far, results obtained from genetic, functional and biomarker screens have been gathered to assemble the heterogeneous risk factors into pathways affected in ASD [10, 11]. Among these are pathways of energy metabolism, neurotransmitter signaling and/or oxidative stress response and protein synthesis in the context of synaptic plasticity [12–15]. These studies confirmed that the functional impact of genetic variants is not only determined by the known functional role of the gene product, but also by its network context [16]. Thus, studies investigating putatively causative genetic variants in cases of otherwise undefined genetic status may prove beneficial for understanding malfunctions in cellular pathways altered in ASD [17].

We have previously characterized the genetics of two rare X-linked, recessive amino acid substitutions in the RPL10 gene (Ribosomal Protein of the Large Ribosomal Subunit 10) co-segregating with ASD in three unrelated German families [18, 19]. At present, only a few other reports on the genetics of human RPL10 are available. Gong et al. did not identify any additional RPL10 mutations when screening 141 patients, which could probably be attributed to the low frequency of RPL10 mutations [20]. At the level of mRNA expression, Gong et al. reported no differential regulation of RPL10 in ASD [20]. This is supported by our showing that the RPL10 mutations neither affect X-chromosomal inactivation nor alter the overall mRNA expression levels of RPL10[19]. However, a suggestive etiological role of RPL10 is emerging from genome-wide studies: for example, Piton et al. identified a novel exonic silent mutation C > T [A211A] [21] of yet unknown pathological impact. Silent mutations have been reported to influence translation kinetics and protein folding [22, 23]. While the pathophysiological role of the silent RPL10 C > T [A211A] mutation is not known, we note that the corresponding amino-acid position falls within a RPL10 C-terminal sequence tract delineated by the RPL10[L206M] and the RPL10[H213Q] mutations identified in ASD patients described by us [18, 19]. Piton et al. did not identify non-synonymous mutations in RPL10 [21]. These results are underlined by the finding that the functionality of RPL10 is readily compromised by single amino-acid mutations affecting one of the many functions of this translational regulator [24–26] rather than by changes in copy number [27]. Again, this is in agreement with both the rare occurrence of RPL10 single amino-acid mutations in ASD specific screens and the results of two genome-wide copy number variation (CNV) studies, which reported de-novo copy number duplication spanning RPL10 in two independent male patients but not in controls [14, 28].

Collectively, our and other studies suggest that RPL10 dysfunction in ASD pathophysiology may be triggered by either altered gene dosage or loss of function within the cellular compartment of RPL10 activity. Indeed, we and others have shown in the eukaryotic model system yeast that selected mutations in rpl10 or altered gene-dosage of rpl10 and other ribosomal proteins selectively shift the mRNA expression profile to generate altered expression levels of distinct protein sets [24, 27]. At the functional level, this is thought to be mediated by altering the role of Rpl10 during initiation of translation [29, 30] and mRNA recruitment and/or in the kinetics of accommodating the incoming tRNAs during elongation of protein synthesis [26]. Interestingly, we could show that the rare RPL10 variants, L206M and H2013Q, reduce the cellular translational capacity in a yeast model, that is, decrease of translating polysomes and increase in translationally less active monosomes and free ribosomal subunits [18]. It is known that wild-type yeast cells employ this strategy in response to a change in metabolic conditions, either in response to a change in nutrient availability [31] or in response to oxidative stress [32]. While reduction of translational capacity is associated with reduction of global protein synthesis, it is also accompanied by increased expression of distinct stress response proteins [32].

Studies in peripheral cells, such as lymphoblastoid cell lines (LCL), have been successfully employed to study dysregulation of protein expression in ASD [33]. In a systematic review Rossignol and Frye reported studies in LCLs that monitored abnormalities in ASD within distinct pathways, with the strongest evidence for immune dysregulation/inflammation and oxidative stress [34]. In this study here, making use of available LCLs we investigated the effect of the human ASD-specific RPL10[H213Q] mutation [Genebank NM_006013, c.639C > G] on cellular pathophysiology, in particular on protein expression of LCLs, and performed a comparative study employing an unrelated ASD patient set.

The effect of the RPL10[H213Q] mutation on protein expression was analyzed by a comparative two-dimensional-DIGE (difference gel electrophoresis) approach studying RPL10[H213Q] mutation carriers, RPL10 wild-type relatives and unrelated controls. Comparing RPL10[H213Q] induced proteomic changes to protein patterns present in unrelated ASD patients without any mutations in RPL10, we did not find individually overlapping proteins, but rather identified a common pattern among the differentially expressed proteins related to the response to oxidative stress. Finally, in a proof of principle study we aimed at testing the hypothesis that rpl10 deficiency induces a protein expression pattern reminiscent of an oxidative stress response. When using a yeast model system, we obtained evidence that indeed a gene-dosage reduction of RPL10 generated a protein signature, which was overlapping with the pattern observed in wild-type cells subjected to oxidative stress.

This approach has allowed us to assess if and how the proteomic network is altered by the autism-associated RPL10[H213Q] mutation [Genebank: NM_006013, c.639C > G] and to which degree the resulting proteomic signature resembles alterations in protein expression present in the general ASD population without RPL10 mutations. Finally we are able to discuss our data in the context of substantial literature reports on peripheral markers of ASD pathophysiology [34] and references therein.

Methods

Selection of patients and controls

LCLs from RPL10[H213Q] mutation carrier families 277 and 440 [18, 19] were available for the two male index patients (hemizygous) and the mothers (heterozygous), an unaffected carrier sister of family 440 (heterozygous) and an unaffected wild-type half-brother of family 277. Two randomly selected unrelated controls (male and female) were chosen as wild-type controls for this RPL10 mutation study. For the study of ASD-specific protein patterns LCLs from 10 male RPL10 wild-type ASD patients (ASD setup) were matched for ethnicity (Caucasian) and diagnosis (International Classification of Diseases 10th revision (ICD-10): F84.0 with IQ >65) to the index patients of families 277 and 440 [see Additional file 1: Table S1]. Ten cell lines from male uncharacterized Caucasian control individuals with RPL10 wild-type were used as controls (CTRL setup). RPL10 mutational status was confirmed by direct sequencing. LCLs had been established as published in Neitzel [35].

The investigation conformed to the principles outlined in the Declaration of Helsinki. All subjects or caregivers provided written informed consent for the participation in this study and publication of this report which had been approved by the German local institutional ethics committees in Frankfurt (Fachbereich Medizin der Johann Wolfgang Goethe-Universität Frankfurt am Main, Ethik-Kommission) and Heidelberg (Ethikkommission I der Medizinischen Fakultät Heidelberg).

Growth conditions for lymphoblastoid cell lines

Viable cells (1 × 106) from a stationary phase culture were inoculated into 5 ml RPMI medium supplied with 10% fetal bovine serum (FBS), 100 U/ml penicillin, 100 μg/ml streptomycin, 2 mM L-glutamine (all from GIBCO®, Life Technologies, Paisley, UK). Volume was doubled each time cultures reached a density of 1 × 106 cells/ml. Cells were harvested at a final density of 1 × 106 cells/ml in 40 ml by centrifugation at 700 × g for five minutes and washed twice with PBS.

Yeast growth conditions

Yeast strains BY4743 (wild type (WT); genotype: MATa/MAT α his3 Δ0/his3 Δ0; leu2 Δ/leu2 Δ0; met15 Δ0/MET15; LYS2/lys2 Δ0; ura3 Δ0/ura3 Δ0) and rpl10 deficient strain (BY4743; genotype: MATa/MAT α his3 Δ 0/his3 Δ0; leu2 Δ /leu2 Δ0; met15 Δ0/MET15; LYS2/lys2 Δ0; ura3 Δ0/ura3 Δ0; rpl10::kanMX4/RPL10) (all EUROSCARF, Frankfurt, Germany) were grown in 100 ml synthetic complete media with a starting concentration of OD600 = 0.1. Oxidative stress was induced by applying 0.8 mM H2O2 after the cells had reached an OD600 of 0.4. Cells were harvested at an OD600 = 0.7 by centrifugation at 1,500 × g for five minutes and washed twice with H2O.

Protein preparation

Protein preparation was performed at 4°C. LCL pellets were resuspended in 750 μl DIGE lysis buffer (7 M urea, 2 M thiourea, 30 mM Tris pH 8.5 and 4% 3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate (CHAPS), all from Roth, Karlsruhe, Germany) and subjected to five successive freeze and thaw cycles. Yeast cell pellets were resuspended in 1 ml DIGE lysis buffer for yeast (9 M urea, 4% CHAPS, 30 mM Tris pH 8.8, all from Roth) and mechanically lysed using 0.45 mm acid-washed glass beads. The cytosolic protein fractions from LCLs or yeast cells were purified by centrifugation at 13,000 × g for 20 minutes. The protein concentration of the supernatant was measured using the Protein Bradford Assay (Biorad, Munich, Germany).

Two-dimensional difference gel electrophoresis

Two-dimensional-DIGE experiments were performed as previously published [36]. LCL and yeast protein extracts were processed the same way using the DIGE minimal labeling kit (GE-Healthcare, Uppsala, Sweden) according to the manufacturer’s recommendation. All individual samples were cyanine 5 (Cy5)-labelled and compared to a Cy3-labelled equimolar mixture of all samples (internal standard). Gels were scanned at a resolution of 100 μm on the Typhoon 9200 Variable Mode Imager (GE-Healthcare). DeCyder Software version 5.0 (GE-Healthcare) was used for spot detection and primary analyses; only spots with an area larger than 400 pixels and an intensity slope smaller than 2 were included. At least 20 spots per gel were matched manually and an additional 20 spots per gel were verified by manual inspection after automated spot matching. On average, 1,200 spots per gel were successfully matched and included in further calculations. The relative volume of each individual spot was normalized against the respective internal standard. T-test on the log2 ratios (sample/internal standard) comparing the two groups was performed. For further statistical analyses ‘R’ software version 2.15 was used.

Selection of spots of interest for identification

Spots of interest were defined based on two criteria: first, candidate spots that were differentially regulated and second, spots that had a variance which was different between the two analyzed groups (this analysis was only performed in the ASD versus CTRL setup).

To identify candidates relevant to the ASD etiology including also the ones with small effects we decided to focus on T-test P-values only including spots regulated below the common threshold of 1.96 to 2 fold standard deviations.

After excluding the spots where the variance was similar for both (that is, the patient and control set), we focused on those spots where the variance was either larger or smaller between the two biological groups. We were also interested in minor changes and, thus, applied F-test statistics at low threshold (F-test P-value 0.1).

In summary, protein spots were selected for tandem mass spectroscopy (MS/MS) identification if: 1) the P-value of the students T-test was below 0.05; or 2) if an altered variance was detected (F-test P-value <0.10); and 3) if the protein spot was visible on a preparative two-dimensional gel.

Protein identification and validation

For two-dimensional-preparative gels 500 μg of total protein were used and Coomassie stained following published protocols [37]. Picked protein spots were sent to the Max Planck Institute, Department of Biochemistry, Martinsried, Germany and identified by MS/MS.

Computational analysis

All statistical calculations were performed using DeCyder 5.0 or ‘R’ 2.15. Hierarchical cluster analysis was performed by applying the heatmap.2() function from the R-package ‘limma’. Spots to be included in the cluster analysis were selected based on their P-values in the Bayes linear model predicting the biological groups using the R function lmfit(). Three models were visualized including candidates with P-values ≤ 0.05, ≤0.03 or ≤0.01. Best clustering was observed when selecting spots with a Bayesian model P-values ≤ 0.05 for ASD versus controls and ≤ 0.03 for RPL10[H213Q] mutation carriers versus non-carriers.

Network analysis was performed using Cytoscape version 2.8.2 [38] with GeneMania plug-in version 3.2 [39]. Gene Ontology (GO)-term enrichment analysis was performed using the BiNGO plug-in version 2.44 [40]. GO-terms for biological processes (GOTERM_BP_FAT) for each identified candidate (human or yeast) were retrieved from the Database for Annotation, Visualization and Integrated Discovery (DAVID) [41].

The selection of the functional sub-networks was based on key terms searched within the list of retrieved GO-terms. Candidates were included in the network named ‘Energy metabolism’ if they were associated with a GO-term including the words ‘energy’ or ‘glycolysis’. For ‘redox metabolism’ we applied the search terms ‘oxidative stress’, ‘reduction’, ‘oxidation reduction’ or ‘peroxide’. The ‘Protein and mRNA metabolism’ network was defined by the search terms ‘mRNA binding’, ‘mRNA transport’, ‘mRNA degradation’, ‘protein metabolism/catabolism/anabolism’, ‘protein folding’, ‘translation’, ‘protein localization’, ‘peptidase activity’, ‘ER’, ‘endoplasmatic reticulum’ or ‘protein catabolism’. A list of the retrieved candidates and selected GO-terms is available in Additional file 2: Table S7. We manually added RPL10 to each list of genes if it was not a member by itself.

Selected candidates of the respective sub-networks were analyzed and visualized based on co-expression using the GeneMania plug-in mining expression data from the default 20 out of 161 databases for gene-expression available within GeneMania. For details see [42].

Results

The impact of an ASD-associated RPL10 mutation on differential protein expression was investigated comparing LCLs of carriers of a RPL10[H213Q] mutation to wild-type RPL10 carriers including two unrelated healthy controls [see Additional file 1: Table S1 and Figure S1]. RPL10 mutant carriers were compared to RPL10 wild-type carriers disregarding ASD diagnosis to identify protein expression patterns related to the mutation and not to the ASD itself.

Two-dimensional difference gel electrophoresis experiments

LCLs generally showed a large variance in protein expression, which can be attributed to the clonality effects during transformation as described in [43]. However, within the set of all RPL10[H213Q] mutation carriers 62 significantly regulated candidate spots were observed [see Additional file 3: Table S2] of which 7 were identified using MS/MS (Table 1). Among the up-regulated candidates we identified the lipid β-oxidation related enzyme ECH1 (FC (fold-change) = 1.76, P-value = 0.0003) and the translation elongation factor EEF1D (FC = 1.29, P-value = 0.0085). Furthermore, we observed differential expression of two central glycolytic enzymes. We report a significant down-regulation of three spots of different charge identified as TPI1 (triosephosphate isomerase, all spots with a FC < −1.30 and P < .01) and a trend for GAPDH (glyceraldehyde 3-phosphate dehydrogenase, FC = −1.45, P-value = 0.0528). For TPI1 this suggests an alteration also on the posttranslational level. Additionally, a mRNA regulatory enzyme, HNRNPK (FC = −2.79, P-value = 0.0017), was identified (Table 1, Additional file 3: Table S2).

Table 1 Proteins identified by tandem mass spectrometry following two-dimensional DIGE analysis

Hierarchical cluster analysis using all detected protein spots in RPL10[H213Q] carriers did not distinguish carriers from non-carriers. However, selecting spots that were contributing to a Bayesian linear model (BLM) predicting mutation status with a significance of P-value ≤0.03 allowed correct clustering of individuals according to mutation status. This criterion was fulfilled by 33 spots, among them TPI1, ECH1 and HNRNPK. GAPDH was slightly above the threshold (BLM P-value = 0.0340) and, thus, was excluded but should still be considered as a biologically relevant target of differential protein expression in a RPL10[H213Q] background. In addition, using cluster analysis we were able to demonstrate that the protein profile of RPL10[H213Q] carriers is distinct from wild-type carriers, as the unaffected brother of family 277 clearly maps to the cluster representing the controls. Moreover, the protein signatures of the mutation carriers themselves fall into two distinct groups, representing the individual families’ genetic background (Figure 1A). We conclude that BLM cluster analysis is able to both separate carriers from non-carriers as well as to report similarities within the carriers of the two families.

Figure 1
figure 1

Two-dimensional-DIGE results. (A) Hierarchical clustering of RPL10[H213Q] mutation carriers (MUT) and wild-type allele carriers (WT) including spots with P-value ≤0.03 only. These 33 spots separated MUT from WT and, in addition, separately clustered families. Two distinct protein clusters (A, B) are visible with cluster B being more affected in family 440. Asterisks mark samples from ASD patients. (B) Best hierarchical clustering between ASD and controls (CTRL) was observed when selecting spots with P-value ≤0.05 (N = 52). Two distinct clusters are down- (A) and up-regulated (B), respectively, in the ASD cohort. For better visualization missing values in A and B omitted during calculations are shown in black. (C) In yeast, proteins differentially regulated under hydrogen peroxide (H2O2) induced oxidative stress overlap with a set of proteins (marked in grey) differentially regulated in the rpl10 deficient yeast strain L10/ΔL10 (all P-values ≤0.05, comparison for both conditions are to wild-type yeast). (D) Co-expression networks of identified candidates based on biological functions related to GO-terms. RPL10 was manually spiked into each network and is co-regulated in every module. Gene-expression data are retrieved from the GeneMania database. For details on methods and network construction see the Materials and Methods section. Symbols used in protein names in B and C: ‘&’ = a mixture of both proteins was identified; ‘|’ = MS data were not able to clearly distinguish between named protein isoforms. ASD, autism spectrum disorders; DIGE, difference gel electrophoresis; GO, Gene Ontology database; MS, mass spectroscopy.

To distinguish between RPL10[H213Q]-derived changes and ASD-specific variations independent of the RPL10 mutation, we advanced to proteome analysis of LCLs derived from 10 male patients, all with a comparable ASD diagnosis and similar IQ and 10 unrelated randomly selected male controls [see Additional file 1: Table S1 and Figure S1]. Of the 1,215 detected protein spots 1,195 passed the quality criteria. We investigated protein patterns with respect to change in protein expression level (T-test P-value ≤0.05) and variance of protein expression (F-test ≤0.1). Sixty-six candidates showed a significantly altered expression level of which 17 spots were identified using tandem MS/MS (Table 1), representing 15 different proteins (PCBP1, ALDOC, ERP29, ECHS1, ATP5B, PSME2, ACTB, TAGLN2, HNRNPA2B1, two isoforms, ATP5A1, GAPDH, PGK1, two isoforms, HSPD1, PSMA1, GLUD1; Table 1, Additional file 1: Figure S2, for full names see Additional file 4: Table S4). Of the more than 200 protein spots that had a differential expression variance we identified 10 proteins. Of these, four proteins had also been identified as being differentially expressed (ACTB, HNRNPA2B1, PSME2, PSMA1) and six additional proteins exclusively showed an altered variance in protein expression (ATP5H, EEF1D, HNRNPK, HNRPDL, HSPA5, PRDX2, see also Table 1 and Additional file 5: Table S3).

Clustering of individuals according to affected status was not possible using all spots. The best hierarchical clustering between ASD samples and controls was observed when selecting spots that contributed to the underlying BLM with a significance of P-value ≤0.05 (Figure 1B), suggesting a specific set of proteins to be etiologically relevant. Twelve of the 52 spots that fulfilled these criteria had already been identified in the differentially regulated set mentioned above. Again, a distinct protein signature emerged which is differentiating cases and controls. However, we also note that one of the control individuals has a protein profile reminiscent of the ASD patient set (Figure 1B).

Candidate proteins (N = 21) collectively identified in the ASD versus CTRL set mapped to etiopathologically relevant pathways previously described in neurodevelopmental disorders: energy metabolism (ALDOC, ATP5A1, ATP5B, ATP5H, ECHS1, GAPDH, PGK1, GLUD1), mRNA and protein metabolism (EEF1D, ERP29, HSPA5, HNRNPA2B1, HNRPDL, HNRNPK, HSPD1, PCBP1, PSMA1, PSME2), cytoskeleton (ACTB, TAGLN2) and redox scavenger (PRDX2). Isoforms of two proteins, HNRNPA2B1 and PGK1, were detected. For example, the glycolytic enzyme PGK1 showed two spots down-regulated (FC < −1.11, and P-values <0.05), suggesting again an alteration at the level of post-translational modification.

Interestingly, none of the differentially expressed proteins identified in the RPL10 family cohort was overlapping with the differentially expressed proteins in the general ASD cohort. However, an overlap within the ASD cohort showing altered variance was observed. This overlapping set included HNRNPK isoforms, EEF1D and ECH1, which is the mitochondrial homolog to the peroxisomal ECHS1.

Validation of protein signatures

In an attempt to validate our proteome results with other techniques, we used the investigated LCLs for RT-PCR of candidate mRNAs [see Additional file 6: Table S8] and Western blot analysis [see Additional file 7: Table S9]. Of the five significantly regulated candidates identified in RPL10[H213Q] versus RPL10-WT setup (ECH1, EEF1D, GAPDH, HNRNPK and TPI1), we were able to confirm increased protein expression for ECH1 (FC = 1.32; P-value = 0.0250 U-test), but no changes in mRNA expression. Interestingly, we identified alterations at the protein level in the RPL10[H213Q] mutation carriers of candidate proteins identified in the ASD versus CTRL setup, namely overall expression of ACTB, HNRPDL, HSPD1 and PGK1 and specific isoforms of ALDOC and ECHS1. These candidates were not significantly changed in the two-dimensional DIGE experiment but the direction of FC was comparable between the two setups.

For the ASD versus CTRL setup we were able to confirm reduction of ACTB (FC = −1.19, P-value = 0.041 U-test) and observed a suggestive up-regulation for ERP29 (FC = 1.26, P-value = 0.096) at the protein level. At the mRNA level we observed changes in the EEF1D and the HNRNPK gene-expression. At the level of protein expression in the two-dimensional DIGE setup both candidates showed an altered variance but were not differentially regulated. Our results suggest a predominant regulation of the differentially expressed proteins at the translational level. This is supported by a study of Schwanhäusser et al. who quantified global mammalian gene and protein expression and found that cellular abundance of proteins is predominantly controlled at the level of translation [45].

Validation of protein signatures at functional levels

In search for validation of the candidate proteins on the functional level we started to look for shared molecular properties. The identification of a redox scavenger differing in variance (PRDX2; F-test P-value = 0.0095) and several differentially expressed mitochondrial enzymes (for example, ATP5A1, ATP5B, ECHS1; all T-test P-values <0.03) led us to the hypothesis of a redox-sensitive protein network. Indeed, a literature search revealed that 22 out of the total 25 candidates identified in this study (combining RPL10[H213Q]-regulated and ASD-specific proteins) are described as oxidative stress response proteins, or as regulators of the levels of reactive oxygen species (ROS) (for references see Additional file 1: Table S5).

To investigate a possible link between alterations in redox-sensitive proteins, energy metabolism, altered RPL10 function and change in protein expression pattern, we performed a proof-of-principle study in yeast [see Additional file 1: Figure S1]. First, we confirmed previous findings that in response to oxidative stress, exemplified by H2O2 stress, yeast wild-type cells differentially express redox scavengers, redox-sensitive enzymes as well as glycolytic enzymes [46, 47] (Figure 1C and Additional file 8: Table S6). We observed regulation of yeast Tdh3 (one of the three yeast GAPDH homologues), where two isoforms were differentially expressed (FC = 1.34, P-value = 0.0027 and FC = −1.07, P-value = 0.0118). Interestingly, proteome analysis of the yeast rpl10 deficient strain under standard conditions and not subjected to oxidative stress showed that 13 of 45 significantly (P-value ≤0.05) regulated protein spots were up- or down-regulated in an identical fashion as a subset of the overall set of H2O2 stress-regulated candidates in wild-type yeast cells [see Additional file 8: Table S6, Figure 1C]. We were able to identify in 9 of these 13 spots a total of 11 proteins: Asc1, Cys4, Frd1, Rim1, a mixture of Rps5 and Cdc19, Ssb1 or Ssb2, Tma19, Vma1, and Tma29 [see Additional file 4: Table S4 and Additional file 8: Table S6, Figure 1C]. In addition to the glycolytic enzyme Cdc19 (pyruvate kinase; FC = −1.04, P-value = 0.0189), we also noted differential regulation of a Tdh3 isoform, however with a suggestive significance only (FC = −1.03, P-value = 0.0777, Additional file 8: Table S6).

The overlap in glycolytic proteins in the RPL10[H213Q] mutant carriers, the ASD patient set, the oxidatively stressed yeast cells and the yeast cells with rpl10 deficiency suggested a gene ontology association study. Indeed, a first investigation on biological processes revealed enrichment for energy metabolism, protein and mRNA metabolism, and redox metabolism [see Additional file 2: Table S7 and Additional file 1: Figure S2]. Subsequently, based on their GO-terms, we selected those candidates of related function – protein and mRNA metabolism, energy metabolism and redox metabolism – and performed a gene co-expression network analysis based on available transcriptome-data from GeneMania [48] (for details see Methods and Additional file 2: Table S7). For both yeast and human mRNA expression data, we observed that the individual networks of the selected candidates built ensembles that are co-regulated at mRNA levels with RPL10 (Figure 1D). In addition, the tight co-regulation at mRNA levels of the RPL10[H213Q] and the ASD candidate proteins strengthens the reported functional overlap between the molecular phenotypes generated by the RPL10 mutation and the ASD cohort, respectively.

Discussion

An estimated 1 out of 110 individuals in the US is currently affected by ASD [1]. For decades up until now ASD has been clinically defined by behavioral observations following standard manuals. In spite of tremendous effort, the etiology of ASD is still unclear and large consortia have obtained genetic evidence that the majority of ASD cases do not result from simple gene or chromosomal disorder. A meta-analysis of research trends in ASD conducted by Rossignol and Frye indicates that in addition to dysfunctions of the CNS thought to be causative for the cognitive and behavioral features of ASD, pathophysiologies that transcend specific CNS function accompany most ASD cases [49]. These authors combined reports obtained in neurocognitive areas with focus on CNS malfunction (neuroimaging, neuropathology, theory of mind and genetics) with reports on the study of pathophysiological processes that could more directly result in cellular dysfunctions and the subsequent development of psychiatric disorders, such as ASD. In particular, there are well documented differences in ASD patients compared to controls in immune response, especially that of the gastrointestinal tract [50], in oxidative stress response/redox regulation and general detoxification [51], as well as in energy metabolism and associated mitochondrial function [52]. One common theme of all pathophysiological parameters reported in this meta-study is their relatively small, albeit significant, deviation from the respective level in typically developed controls. While different systemic abnormalities may accompany most ASD cases, they must be small to allow for organismal viability. In sum, this study suggests that ASD, rather than being a purely CNS disorder, displays a wide range of systemic abnormalities making the study of peripheral cells and systemic fluids an attractive target to investigate the etiology of ASD.

The results of our study mirror in several ways the documentations of Rossignol and Frye [34]. First, we report that an ASD-associated RPL10[H213Q] mutation alters protein expression patterns similar but not identical to the one observed in an unrelated ASD cohort without this mutation. This suggests a common downstream pathophysiology resulting from different unidentified upstream genetic variants. Second, the RPL10[H213Q] expression signature differentiating carriers from non-carriers, is modified by the genetic background of the carriers within families, as observed in cluster analysis. We hypothesize that distinct different protein expression levels in individuals may define distinct pathophysiological phenotypes. Furthermore, while RPL10[H213Q] heterozygous carriers, all women, and RPL10[H213Q] hemizygous carriers, all ASD cases, share common motifs in altered protein expression, the contribution of the remaining wild-type RPL10 [19] to the establishment of a healthy or subclinical phenotype of the female carriers has not been studied here. Thus, it remains to be tested if the protein signatures identified here ultimately correlate with severity of ASD as reported for mitochondrial biomarkers (for an extensive meta-review compare Rossignol and Frye [49]).

Third, the protein signatures observed in the RPL10[H213Q] mutant cohort and in the unrelated ASD cohort are characterized by both having different expression levels and containing different redox-sensitive members of identical pathways. This is in agreement with observations made by Rossignol and Frye [34], namely that depending on the study cohort, overlap is observed in same pathophysiologies.

To assess the phenomenon of oxidative stress and redox control in the context of ASD pathophysiology it is important to recall that endogenous oxidative stress results from an excess of ROS which no longer can be reduced to non-pathological levels by anti-oxidant scavenger and repair mechanisms. Ultimately, all excess ROS are generated by malfunctions in energy metabolism, primarily fueled by glycolysis in the cytoplasm and the tricarboxylic acid (TCA) cycle and oxidative phosphorylation (OXPHOS) in mitochondria. The successive breakdown of carbon fuels during glycolysis and the TCA cycle supplies electrons to the electron transport chain for mitochondrial OXPHOS to produce ATP and H2O. During this process leakage of electrons results in the formation of ROS, the superoxide radical (O2–), hydrogen peroxide (H2O2) and the hydroxyl radical (OH). If increasing amounts of ROS escape detoxification, excess oxidative power attacks lipids, proteins and DNA (reviewed in [53]). Lipid peroxidation is detrimental to the function of receptors and ion channels in all cellular membranes and, thus, also to the balance of membrane potential in the neuron [54]. Loss of control over membrane potential, in turn, has been reported to affect receptors’ function [55] and may thus lead to altered signal transduction between pre- and post-synaptic membranes. Protein oxidation leads to loss of enzyme function and/or loss of structural integrity [56]. Oxidation of DNA is a major source of somatic mutations [57]. Each of these pathologies is harmful to all cells and has been shown to have a severe negative impact on neural cells, in particular in the context of inflammation and synaptic function [54, 58]. On the level of protein expression, oxidative stress is known to quantitatively and qualitatively alter mRNA translation in model systems [32, 59].

Altered translation has been until now one of the best documented cellular processes executing correct synaptic plasticity by modulating density and composition of synaptic proteins, the loss of which is one of the hallmark pathologies in ASD brain tissue [60]. Our observation that altered availability of the translational regulator RPL10 generates a cellular redox-sensitive stress response signature in the yeast model, suggests that a mutation in human RPL10 could trigger a similar phenomenon in neural cells. This is supported by our observation of a redox-sensitive protein signature in LCLs of RPL10[H213Q] patients.

How then might an ASD-mutation modulate oxidative stress in the periphery and the CNS of ASD patients? First, a mutation in a translational regulator, such as the ribosomal protein RPL10, might result in unscheduled or increased recruitment of redox-sensitive proteins that are a mimic of the translational response driven by other, more upstream ASD genetic variations (for example, provided by different genetic backgrounds). Basically, in all these cases this scenario would supply the cell with a similar compromised redox microenvironment as observed in our study. If there is a related translational response for other or even most ASD cases remains to be established. Second, several affected genes in our study encode enzymes that were reported to produce or eliminate ROS and thus, when mutated lead to loss of control of ROS production associated with impaired gene expression [51]. Third, any enzyme directly active in energy metabolism, when mutated, would fulfill the criteria for an autism candidate gene, as it ultimately would favor imbalance in mitochondrial ROS production. This is underscored by a series of enzymes of glycolysis, the TCA cycle and OXPHOS found to be affected in neuropsychiatric diseases [61, 62]. GAPDH and TPI, identified in this study to be differentially regulated in ASD, have been reported to be differentially expressed in brains of schizophrenia patients [63]. The close neighborhood of the RPL10 protein to both proteins in the energy metabolism co-expression network may thus underscore the etiological relevance of RPL10 in neurodevelopmental disorders. Finally, a protein identified as dysfunctional in ASD may be directly sensitive to cellular redox status, either by oxidative stress induced alterations in posttranslational modifications [64, 65] or by shifting the equilibrium of redox-sensitive disulphide bonds in regulatory proteins [66]. The former proposal is supported by the observation of differential expression of posttranslationally altered proteins in our study, and the latter one by redox-sensitive disulphide bonds reported for ASD candidate proteins [66, 67].

For a discussion of putative functional roles of individual proteins identified with the redox-sensitive proteins signature in the etiology of ASD we will focus on the ASD candidate protein RPL10, which is co-expressed with candidates in all three GO-term enriched networks identified in both human and yeast studies (Figure 1D). A second such protein is GAPDH/TDH3 which was highly connected within the co-expression networks of energy metabolism (human and yeast) and redox response (human and yeast) as well as mRNA and protein expression (yeast only, as human GAPDH is so far not GO-term annotated with translation or mRNA regulation).

The glycolytic enzyme GAPDH is a redox-sensitive key enzyme of glycolysis and together with its immediate glycolytic neighbor TPI1 (Figure 1D) has previously been reported to be differentially regulated in schizophrenia [63]. This functional interaction of human RPL10 mutations with glycolytic enzymes is confirming previous findings from yeast studies [24, 25]. GAPDH and other glycolytic enzymes are receiving increasing interest as regulatory proteins in differentiation and disease [68]. It has been shown that isoforms of GAPDH and downstream glycolytic enzymes are expressed in different metabolic stages in yeast [69] as well as in cellular differentiation in mammalian cells [70]. These studies propose that posttranslational modifications of GAPDH and other glycolytic enzymes, some targeted by oxidative modification, make these glycolytic enzymes to platforms within the successive stages of energy metabolism, which can be used to switch between high and low levels of energy metabolism [69]. For example, several studies have shown, that reduction of GAPDH function reroutes the glycolytic flow to the pentose phosphate pathway, used to generate NADPH, which serves as the reducing power to restore glutathione (GSH) levels in the GSH/GSSH redox couple, the major intracellular antioxidant/detoxification species [71]. Alterations in GAPDH function, induced by mutation, altered expression level or change in posttranslational modification, may thus affect energy metabolism and redox status, which is of particular importance in high energy consuming cells like neurons.

Ribosomal protein RPL10 is also of specific interest to further studies on the etiology of ASD. So far, translational regulators described in more detail as mediating dysfunctional translation in ASD are FMRP [72] and eukaryotic initiation factor 4E, eIF4E [73]. FMRP is a translational repressor binding to and, thereby, preventing protein synthesis of synaptic mRNAs. Mutated FMRP leads to unscheduled protein synthesis at the synapse, resulting in compromised synaptic plasticity [74]. Also, mutationally altered eIF4E leads to hyperconnectivity of neural circuits due to increased synaptic protein synthesis. EIF4E is a key component of the eukaryotic translation initiation machinery, which when over expressed leads to increased expression of neuroligins, which are postsynaptic proteins that are causally linked to ASD [73]. RPL10 is a component of the translating ribosome and different modifications of RPL10 may alter different ribosome functions, for example mutations in RPL10 have been described that alter translation initiation [30] or protein expression levels [24]. At present, we do not know if there are modifications of RPL10 resulting from upstream regulators of translation, for example upstream regulators of synaptic protein expression and/or oxidative stress response. If so, it would be interesting to learn whether they are known ASD candidate proteins.

In sum, our results on differential redox-sensitive protein expression patterns in patients with RPL10[H213Q] mutations and unrelated RPL10 wild-type ASD patients, in particular in the context of alterations in energy metabolism and protein translation, have led us to formulate the following hypothesis (Figure 2). The altered redox-sensitive protein signature, as observed in peripheral cells derived from our patient sets, mirrors a cellular response to oxidative stress. We suppose that this signature is largely compatible with normal cellular function in peripheral tissues under conditions of moderate ROS levels. However, the brain is the organ with the highest oxygen consumption and metabolism [53] and, thus, we propose that in the presence of high oxygen turnover in neural cells, the altered redox-sensitive protein signatures may fail to buffer oxidative stress in energy metabolism, redox metabolism and protein metabolism. To arrive at a molecular understanding of the role of an individual, RPL10-regulated and differentially expressed redox-sensitive protein, GAPDH may serve as an example (Figure 2). GAPDH is a key glycolytic enzyme and major player in rerouting reductive power to balance increased oxidative stress. In the process of adapting to different levels of cellular oxidative stress, GAPDH may tolerate only a distinct range of change in expression level. While the altered expression level of GAPDH observed in lymphoblastoid cells may be compatible with effective redox control in this cell type, it may fail to provide the required response to the high oxygen metabolism present in a neuron. If this insufficiency in the presence of high oxygen stress either results from a loss of regulatory flexibility and/or inadequate network function as a consequence of the differentially expressed state, remains to be established.

Figure 2
figure 2

ASD mutation-induced redox-sensitive protein signatures compromise core synaptic functions in autism. (A) The ribosome recruits mRNAs, encoding proteins for specific pathways (black, green, blue mRNAs). Subsets of mRNAs are under the control of the translational regulator RPL10 (orange), which itself is controlled by upstream regulators (grey). This forms a protein signature characteristic for example for mRNA and protein metabolism (black), redox-metabolism (green) and energy metabolism (blue). (B) In the presence of (mis)functional RPL10 variants (for example, RPL10[H213Q]) or other ASD variants (red crosses) regulating RPL10 the protein signatures are selectively shifted (red rectangles) resulting in alterations (change in number of the framed protein spots) of pathway-specific redox-sensitive proteins as observed here. The altered protein signature thus indicates a response to oxidative stress elicited by the ASD mutation(s). (C) In peripheral tissues, for example, a lymphoblast, deficiencies in RPL10 or upstream regulators thereof alter the redox-sensitive signature in a way which in the presence of moderate ROS (reactive oxygen species) production (small red bolt) is still able to balance cellular stress response. In particular, altered expression of the glycolytic enzyme GAPDH (blue) may be employed for redox buffering by rerouting glycolysis into the PPP (pentose phosphate pathway). This will increase reductive power to balance oxidative power caused by ROS. (D) Under conditions of high oxygen metabolism (large red bolt) as present in neural cells, redox buffering will still support basal synaptic functions, but may no longer be able to prevent oxidative damage of lipids and proteins required for cellular fine-tuning (undulations) of synaptic functions and plasticity (purple, double headed arrow). We hypothesize that the effect of different ASD-mutation induced protein signatures would drive individual autistic phenotypes with differential failure to secure correct executing of synaptic plasticity. Question marks indicate unknown individual steps. ASD, autism spectrum disorders.

Consequently, we propose that ASD-mutation induced altered redox-sensitive protein signatures fail to balance the high oxidative stress in the neuron. The resulting excessive oxidative damage may compromise correct functioning of synaptic proteins and membrane potential to execute synaptic plasticity.

Conclusions

Here, we show for the first time in a patient-specific cell line model for ASD that ASD patients, either with defined or with undefined mutations at the level of protein expression, display modulated expression of redox-sensitive components in those cellular pathways that execute synaptic plasticity in the brain: energy metabolism, redox metabolism and protein metabolism. Individual differences in these expression patterns may contribute to the phenotypic continuum observed in the wide spectrum of ASD patients.

Availability of supporting data

The data sets supporting the results of this article are included within the article and its additional files [see Additional files 1, 2, 3, 4, 5, 6, 7 and 8].

Abbreviations

ASD:

autism spectrum disorder

BLM:

Bayesian linear model

CHAPS:

3-[(3-cholamidopropyl)dimethylammonio]-1-propanesulfonate

CNS:

central nervous system

CTRL:

control

Cy3:

cyanine 3

Cy5:

cyanine 5

DAVID:

Database for Annotation Visualization and Integrated Discovery

DIGE:

difference gel electrophoresis

FBS:

fetal bovine serum

FC:

fold change

GO:

Gene Ontology Database

ICD-10:

International Classification of Diseases Version 10

IQ:

intelligence quotient

LCL:

lymphoblastoid cell line

MS:

mass spectrometry

MUT:

mutation

OXPHOS:

oxidative phosphorylation

PBS:

phosphate-buffered saline

RPMI:

Roswell Park Memorial Institute cell culture medium

RT-PCR:

reverse transcriptase polymerase chain reaction

ROS:

reactive oxygen species

TCA:

tricarboxylic acid

WT:

wild-type.

References

  1. Weintraub K: The prevalence puzzle: autism counts. Nature. 2011, 479: 22-24. 10.1038/479022a.

    Article  CAS  PubMed  Google Scholar 

  2. World Health Organization: International Classification of Mental and Behavioral Disorders. Clinical Descriptions and Diagnostic Guidelines. 1992, Geneva: World Health Organization, 10

    Google Scholar 

  3. Constantino JN, Todorov A, Hilton C, Law P, Zhang Y, Molloy E, Fitzgerald R, Geschwind D: Autism recurrence in half siblings: strong support for genetic mechanisms of transmission in ASD. Mol Psychiatry. 2013, 18: 137-138. 10.1038/mp.2012.9.

    Article  CAS  PubMed  Google Scholar 

  4. Hallmayer J, Cleveland S, Torres A, Phillips J, Cohen B, Torigoe T, Miller J, Fedele A, Collins J, Smith K, Lotspeich L, Croen LA, Ozonoff S, LaJonchere C, Grether JK, Risch N: Genetic heritability and shared environmental factors among twin pairs with autism. Arch Gen Psychiatry. 2011, 68: 1095-1102. 10.1001/archgenpsychiatry.2011.76.

    Article  PubMed Central  PubMed  Google Scholar 

  5. Lichtenstein P, Carlström E, Råstam M, Gillberg C, Anckarsäter H: The genetics of autism spectrum disorders and related neuropsychiatric disorders in childhood. Am J Psychiatry. 2010, 167: 1357-1363. 10.1176/appi.ajp.2010.10020223.

    Article  PubMed  Google Scholar 

  6. Anney R, Klei L, Pinto D, Almeida J, Bacchelli E, Baird G, Bolshakova N, Bolte S, Bolton PF, Bourgeron T, Brennan S, Brian J, Casey J, Conroy J, Correia C, Corsello C, Crawford EL, Jonge M, Delorme R, Duketis E, Duque F, Estes A, Farrar P, Fernandez BA, Folstein SE, Fombonne E, Gilbert J, Gillberg C, Glessner JT, Green A: Individual common variants exert weak effects on the risk for autism spectrum disorderspi. Hum Mol Genet. 2012, 21: 4781-4792. 10.1093/hmg/dds301.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  7. Devlin B, Scherer SW: Genetic architecture in autism spectrum disorder. Curr Opin Genet Dev. 2012, 22: 229-237. 10.1016/j.gde.2012.03.002.

    Article  CAS  PubMed  Google Scholar 

  8. Abrahams BS, Geschwind DH: Advances in autism genetics: on the threshold of a new neurobiology. Nat Rev Genet. 2008, 9: 341-355. 10.1038/nrg2346.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  9. Freitag CM, Staal W, Klauck SM, Duketis E, Waltes R: Genetics of autistic disorders: review and clinical implications. Eur Child Adolesc Psychiatry. 2010, 19: 169-178. 10.1007/s00787-009-0076-x.

    Article  PubMed Central  PubMed  Google Scholar 

  10. Xu L, Li J, Huang Y, Zhao M, Tang X, Wei L: AutismKB: an evidence-based knowledgebase of autism genetics. Nucleic Acids Res. 2011, 40: D1016-D1022.

    Article  PubMed Central  PubMed  Google Scholar 

  11. Banerjee-Basu S, Packer A: SFARI Gene: an evolving database for the autism research community. Dis Model Mech. 2010, 3: 133-135. 10.1242/dmm.005439.

    Article  PubMed  Google Scholar 

  12. Anney RJ, Kenny EM, O’Dushlaine C, Yaspan BL, Parkhomenka E, Buxbaum JD, Sutcliffe J, Gill M, Gallagher L: Gene-ontology enrichment analysis in two independent family-based samples highlights biologically plausible processes for autism spectrum disorders. Eur J Hum Genet. 2011, 19: 1082-1089. 10.1038/ejhg.2011.75.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  13. Van der Zwaag B, Franke L, Poot M, Hochstenbach R, Spierenburg HA, Vorstman JA, Van Daalen E, Jonge MV, De Verbeek NE, Brilstra EH, Van’t Slot R, Ophoff RA, Van Es MA, Blauw HM, Veldink JH, Buizer-Voskamp JE, Beemer FA, van den Berg LH, Wijmenga C, Van Amstel HKP, Van Engeland H, Burbach JPH, Staal WG: Gene-network analysis identifies susceptibility genes related to glycobiology in autism. PLoS One. 2009, 4: e5324-10.1371/journal.pone.0005324.

    Article  PubMed Central  PubMed  Google Scholar 

  14. Pinto D, Pagnamenta AT, Klei L, Anney R, Merico D, Regan R, Conroy J, Magalhaes TR, Correia C, Abrahams BS, Almeida J, Bacchelli E, Bader GD, Bailey AJ, Baird G, Battaglia A, Berney T, Bolshakova N, Bölte S, Bolton PF, Bourgeron T, Brennan S, Brian J, Bryson SE, Carson AR, Casallo G, Casey J, Chung BH, Cochrane L, Corsello C: Functional impact of global rare copy number variation in autism spectrum disorders. Nature. 2010, 466: 368-372. 10.1038/nature09146.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  15. Frustaci A, Neri M, Cesario A, Adams JB, Domenici E, Dalla Bernardina B, Bonassi S: Oxidative stress-related biomarkers in autism: systematic review and meta-analyses. Free Radic Biol Med. 2012, 52: 2128-2141. 10.1016/j.freeradbiomed.2012.03.011.

    Article  CAS  PubMed  Google Scholar 

  16. Barabási A, Gulbahce N, Loscalzo J: Network medicine: a network-based approach to human disease. Nat Rev Genet. 2011, 12: 56-68. 10.1038/nrg2918.

    Article  PubMed Central  PubMed  Google Scholar 

  17. Kou Y, Betancur C, Xu H, Buxbaum JD, Ma’ayan A: Network- and attribute-based classifiers can prioritize genes and pathways for autism spectrum disorders and intellectual disability. Am J Med Genet C: Semin Med Genet. 2012, 160: 130-142.

    Article  CAS  Google Scholar 

  18. Klauck SM, Felder B, Kolb-Kokocinski A, Schuster C, Chiocchetti A, Schupp I, Wellenreuther R, Schmötzer G, Poustka F, Breitenbach-Koller L, Poustka A: Mutations in the ribosomal protein gene RPL10 suggest a novel modulating disease mechanism for autism. Mol Psychiatry. 2006, 11: 1073-1084. 10.1038/sj.mp.4001883.

    Article  CAS  PubMed  Google Scholar 

  19. Chiocchetti A, Pakalapati G, Duketis E, Wiemann S, Poustka A, Poustka F, Klauck SM: Mutation and expression analyses of the ribosomal protein gene RPL10 in an extended German sample of patients with autism spectrum disorder. Am J Med Genet A. 2011, 155: 1472-1475. 10.1002/ajmg.a.33977.

    Article  CAS  Google Scholar 

  20. Gong X, Delorme R, Fauchereau F, Durand CM, Chaste P, Betancur C, Goubran-Botros H, Nygren G, Anckarsäter H, Rastam M, Gillberg IC, Kopp S, Mouren-Simeoni MC, Gillberg C, Leboyer M, Bourgeron T: An investigation of ribosomal protein L10 gene in autism spectrum disorders. BMC Med Genet. 2009, 10: 7-

    Article  PubMed Central  PubMed  Google Scholar 

  21. Piton A, Gauthier J, Hamdan FF, Lafrenière RG, Yang Y, Henrion E, Laurent S, Noreau A, Thibodeau P, Karemera L, Spiegelman D, Kuku F, Duguay J, Destroismaisons L, Jolivet P, Côté M, Lachapelle K, Diallo O, Raymond A, Marineau C, Champagne N, Xiong L, Gaspar C, Rivière J, Tarabeux J, Cossette P, Krebs M, Rapoport JL, Addington A, Delisi LE: Systematic resequencing of X-chromosome synaptic genes in autism spectrum disorder and schizophrenia. Mol Psychiatry. 2011, 16: 867-880. 10.1038/mp.2010.54.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  22. Cortazzo P, Cerveñansky C, Marín M, Reiss C, Ehrlich R, Deana A: Silent mutations affect in vivo protein folding in Escherichia coli. Biochem Biophys Res Commun. 2002, 293: 537-541. 10.1016/S0006-291X(02)00226-7.

    Article  CAS  PubMed  Google Scholar 

  23. Lazrak A, Fu L, Bali V, Bartoszewski R, Rab A, Havasi V, Keiles S, Kappes J, Kumar R, Lefkowitz E, Sorscher EJ, Matalon S, Collawn JF, Bebok Z: The silent codon change I507-ATC-ATT contributes to the severity of the ΔF508 CFTR channel dysfunction. FASEB J. 2013, 27: 4630-4645. 10.1096/fj.13-227330.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  24. Pachler K, Karl T, Kolmann K, Mehlmer N, Eder M, Loeffler M, Oender K, Hochleitner EO, Lottspeich F, Bresgen N, Richter K, Breitenbach M, Koller L: Functional interaction in establishment of ribosomal integrity between small subunit protein rpS6 and translational regulator rpL10/Grc5p. FEMS Yeast Res. 2004, 5: 271-280. 10.1016/j.femsyr.2004.07.009.

    Article  CAS  PubMed  Google Scholar 

  25. Oender K, Loeffler M, Doppler E, Eder M, Lach S, Heinrich F, Karl T, Moesl R, Hundsberger H, Klade T, Eckl P, Dickinson JR, Breitenbach M, Koller L: Translational regulator RpL10p/Grc5p interacts physically and functionally with Sed1p, a dynamic component of the yeast cell surface. Yeast. 2003, 20: 281-294. 10.1002/yea.963.

    Article  CAS  PubMed  Google Scholar 

  26. Petrov AN, Meskauskas A, Roshwalb SC, Dinman JD: Yeast ribosomal protein L10 helps coordinate tRNA movement through the large subunit. Nucleic Acids Res. 2008, 36: 6187-6198. 10.1093/nar/gkn643.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  27. Bauer JW, Brandl C, Haubenreisser O, Wimmer B, Weber M, Karl T, Klausegger A, Breitenbach M, Hintner H, von der Haar T, Tuite MF, Breitenbach-Koller L: Specialized yeast ribosomes: a customized tool for selective mRNA translation. PLoS One. 2013, 8: e67609-10.1371/journal.pone.0067609.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  28. Levy D, Ronemus M, Yamrom B, Lee Y, Leotta A, Kendall J, Marks S, Lakshmi B, Pai D, Ye K, Buja A, Krieger A, Yoon S, Troge J, Rodgers L, Iossifov I, Wigler M: Rare de novo and transmitted copy-number variation in autistic spectrum disorders. Neuron. 2011, 70: 886-897. 10.1016/j.neuron.2011.05.015.

    Article  CAS  PubMed  Google Scholar 

  29. De Keersmaecker K, Atak ZK, Li N, Vicente C, Patchett S, Girardi T, Gianfelici V, Geerdens E, Clappier E, Porcu M, Lahortiga I, Lucà R, Yan J, Hulselmans G, Vranckx H, Vandepoel R, Sweron B, Jacobs K, Mentens N, Wlodarska I, Cauwelier B, Cloos J, Soulier J, Uyttebroeck A, Bagni C, Hassan BA, Vandenberghe P, Johnson AW, Aerts S, Cools J: Exome sequencing identifies mutation in CNOT3 and ribosomal genes RPL5 and RPL10 in T-cell acute lymphoblastic leukemia. Nat Genet. 2013, 45: 186-190.

    Article  CAS  PubMed  Google Scholar 

  30. Baierlein C, Hackmann A, Gross T, Henker L, Hinz F, Krebber H: Monosome formation during translation initiation requires the SR protein Npl3. Mol Cell Biol. 2013, 5: 271-280.

    Google Scholar 

  31. Fuge EK, Braun EL, Werner-Washburne M: Protein synthesis in long-term stationary-phase cultures of Saccharomyces cerevisiae. J Bacteriol. 1994, 176: 5802-5813.

    PubMed Central  CAS  PubMed  Google Scholar 

  32. Grant CM: Regulation of translation by hydrogen peroxide. Antioxid Redox Signal. 2011, 15: 191-203. 10.1089/ars.2010.3699.

    Article  CAS  PubMed  Google Scholar 

  33. Gross C, Bassell GJ: Excess protein synthesis in FXS patient lymphoblastoid cells can be rescued with a p110β-selective inhibitor. Mol Med. 2012, 18: 336-345.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  34. Rossignol DA, Frye RE: A review of research trends in physiological abnormalities in autism spectrum disorders: immune dysregulation, inflammation, oxidative stress, mitochondrial dysfunction and environmental toxicant exposures. Mol Psychiatry. 2012, 17: 389-401. 10.1038/mp.2011.165.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  35. Neitzel H: A routine method for the establishment of permanent growing lymphoblastoid cell lines. Hum Genet. 1986, 73: 320-326. 10.1007/BF00279094.

    Article  CAS  PubMed  Google Scholar 

  36. Laimer M, Kocher T, Chiocchetti A, Trost A, Lottspeich F, Richter K, Hintner H, Bauer JW, Onder K: Proteomic profiling reveals a catalogue of new candidate proteins for human skin aging. Exp Dermatol. 2010, 19: 912-918. 10.1111/j.1600-0625.2010.01144.x.

    Article  CAS  PubMed  Google Scholar 

  37. Candiano G, Bruschi M, Musante L, Santucci L, Ghiggeri GM, Carnemolla B, Orecchia P, Zardi L, Righetti PG: Blue silver: a very sensitive colloidal Coomassie G-250 staining for proteome analysis. Electrophoresis. 2004, 25: 1327-1333. 10.1002/elps.200305844.

    Article  CAS  PubMed  Google Scholar 

  38. Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, Christmas R, Avila-Campilo I, Creech M, Gross B, Hanspers K, Isserlin R, Kelley R, Killcoyne S, Lotia S, Maere S, Morris J, Ono K, Pavlovic V, Pico AR, Vailaya A, Wang PL, Adler A, Conklin BR, Hood L, Kuiper M, Sander C, Schmulevich I, Schwikowski B, Warner GJ: Integration of biological networks and gene expression data using Cytoscape. Nat Protoc. 2007, 2: 2366-2382. 10.1038/nprot.2007.324.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  39. Montojo J, Zuberi K, Rodriguez H, Kazi F, Wright G, Donaldson SL, Morris Q, Bader GD: GeneMANIA Cytoscape plugin: fast gene function predictions on the desktop. Bioinformatics. 2010, 26: 2927-2928. 10.1093/bioinformatics/btq562.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  40. Maere S, Heymans K, Kuiper M: BiNGO. A Cytoscape plugin to assess overrepresentation of gene ontology categories in biological networks. Bioinformatics. 2005, 21: 3448-3449. 10.1093/bioinformatics/bti551.

    Article  CAS  PubMed  Google Scholar 

  41. DAVID: Database for Annotation, Visualization and Integrated Discovery (DAVID ) v6.7. [http://david.abcc.ncifcrf.gov]

  42. Genemania. [http://genemania.org]

  43. Plagnol V, Uz E, Wallace C, Stevens H, Clayton D, Ozcelik T, Todd JA: Extreme clonality in lymphoblastoid cell lines with implications for allele specific expression analyses. PLoS One. 2008, 3: e2966-10.1371/journal.pone.0002966.

    Article  PubMed Central  PubMed  Google Scholar 

  44. SWISS-2DPAGE: Two-dimensional polyacrylamide gel electrophoresis database. [http://world-2dpage.expasy.org/swiss-2dpage]

  45. Schwanhäusser B, Busse D, Li N, Dittmar G, Schuchhardt J, Wolf J, Chen W, Selbach M: Global quantification of mammalian gene expression control. Nature. 2011, 473: 337-342. 10.1038/nature10098.

    Article  PubMed  Google Scholar 

  46. Le Moan N, Clement G, Le Maout S, Tacnet F, Toledano MB: The Saccharomyces cerevisiae proteome of oxidized protein thiols: contrasted functions for the thioredoxin and glutathione pathways. J Biol Chem. 2006, 281: 10420-10430. 10.1074/jbc.M513346200.

    Article  CAS  PubMed  Google Scholar 

  47. Magherini F, Tani C, Gamberi T, Caselli A, Bianchi L, Bini L, Modesti A: Protein expression profiles in Saccharomyces cerevisiae during apoptosis induced by H2O2. Proteomics. 2007, 7: 1434-1445. 10.1002/pmic.200600796.

    Article  CAS  PubMed  Google Scholar 

  48. Warde-Farley D, Donaldson SL, Comes O, Zuberi K, Badrawi R, Chao P, Franz M, Grouios C, Kazi F, Lopes CT, Maitland A, Mostafavi S, Montojo J, Shao Q, Wright G, Bader GD, Morris Q: The GeneMANIA prediction server: biological network integration for gene prioritization and predicting gene function. Nucleic Acids Res. 2010, 38: W214-W220. 10.1093/nar/gkq537.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  49. Rossignol DA, Frye RE: Mitochondrial dysfunction in autism spectrum disorders: a systematic review and meta-analysis. Mol Psychiatry. 2012, 17: 290-314. 10.1038/mp.2010.136.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  50. Ashwood P, Wakefield AJ: Immune activation of peripheral blood and mucosal CD3+ lymphocyte cytokine profiles in children with autism and gastrointestinal symptoms. J Neuroimmunol. 2006, 173: 126-134. 10.1016/j.jneuroim.2005.12.007.

    Article  CAS  PubMed  Google Scholar 

  51. Chauhan A, Chauhan V: Oxidative stress in autism. Pathophysiology. 2006, 13: 171-181. 10.1016/j.pathophys.2006.05.007.

    Article  CAS  PubMed  Google Scholar 

  52. Rossignol DA, Bradstreet JJ: Evidence of mitochondrial dysfunction in autism and implications for treatment. Am J Biochem Biotechnol. 2008, 4: 208-217.

    Article  CAS  Google Scholar 

  53. Halliwell B, Gutteridge , John MC: Free Radicals in Biology and Medicine. 2010, Oxford: Oxford Univ. Press, 4

    Google Scholar 

  54. Volinsky R, Kinnunen , Paavo KJ: Oxidized phosphatidylcholines in membrane-level cellular signaling: from biophysics to physiology and molecular pathology. FEBS J. 2013, 280: 2806-2816. 10.1111/febs.12247.

    Article  CAS  PubMed  Google Scholar 

  55. Jana A, Hogan EL, Pahan K: Ceramide and neurodegeneration: susceptibility of neurons and oligodendrocytes to cell damage and death. J Neurol Sci. 2009, 278: 5-15. 10.1016/j.jns.2008.12.010.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  56. Butterfield DA, Perluigi M, Reed T, Muharib T, Hughes CP, Robinson , Renã AS, Sultana R: Redox proteomics in selected neurodegenerative disorders: from its infancy to future applications. Antioxid Redox Signal. 2012, 17: 1610-1655. 10.1089/ars.2011.4109.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  57. Jena NR: DNA damage by reactive species: mechanisms, mutation and repair. J Biosci. 2012, 37: 503-517. 10.1007/s12038-012-9218-2.

    Article  CAS  PubMed  Google Scholar 

  58. Nakamura T, Cho D, Lipton SA: Redox regulation of protein misfolding, mitochondrial dysfunction, synaptic damage, and cell death in neurodegenerative diseases. Exp Neurol. 2012, 238: 12-21. 10.1016/j.expneurol.2012.06.032.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  59. Morano KA, Grant CM, Moye-Rowley WS: The response to heat shock and oxidative stress in Saccharomyces cerevisiae. Genetics. 2012, 190: 1157-1195. 10.1534/genetics.111.128033.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  60. Kelleher RJ, Bear MF: The autistic neuron: troubled translation?. Cell. 2008, 135: 401-406. 10.1016/j.cell.2008.10.017.

    Article  CAS  PubMed  Google Scholar 

  61. Marazziti D, Baroni S, Picchetti M, Landi P, Silvestri S, Vatteroni E, Catena Dell’Osso M: Psychiatric disorders and mitochondrial dysfunctions. Eur Rev Med Pharmacol Sci. 2012, 16: 270-275.

    CAS  PubMed  Google Scholar 

  62. Martins-de-Souza D, Gattaz WF, Schmitt A, Rewerts C, Maccarrone G, Dias-Neto E, Turck CW: Prefrontal cortex shotgun proteome analysis reveals altered calcium homeostasis and immune system imbalance in schizophrenia. Eur Arch Psychiatry Clin Neurosci. 2009, 259: 151-163. 10.1007/s00406-008-0847-2.

    Article  PubMed  Google Scholar 

  63. Martins-de-Souza D, Maccarrone G, Wobrock T, Zerr I, Gormanns P, Reckow S, Falkai P, Schmitt A, Turck CW: Proteome analysis of the thalamus and cerebrospinal fluid reveals glycolysis dysfunction and potential biomarkers candidates for schizophrenia. J Psychiatr Res. 2010, 44: 1176-1189. 10.1016/j.jpsychires.2010.04.014.

    Article  PubMed  Google Scholar 

  64. Ehlers MD: Activity level controls postsynaptic composition and signaling via the ubiquitin-proteasome system. Nat Neurosci. 2003, 6: 231-242. 10.1038/nn1013.

    Article  CAS  PubMed  Google Scholar 

  65. Figueiredo-Pereira ME, Yakushin S, Cohen G: Disruption of the intracellular sulfhydryl homeostasis by cadmium-induced oxidative stress leads to protein thiolation and ubiquitination in neuronal cells. J Biol Chem. 1998, 273: 12703-12709. 10.1074/jbc.273.21.12703.

    Article  CAS  PubMed  Google Scholar 

  66. Yap L, Garcia JV, Han DS, Cadenas E: Role of nitric oxide-mediated glutathionylation in neuronal function: potential regulation of energy utilization. Biochem J. 2010, 428: 85-93. 10.1042/BJ20100164.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  67. Gu F, Chauhan V, Chauhan A: Impaired synthesis and antioxidant defense of glutathione in the cerebellum of autistic subjects: alterations in the activities and protein expression of glutathione-related enzymes. Free Radic Biol Med. 2013, 65C: 488-496.

    Article  Google Scholar 

  68. Breitenbach M, Ralser M, Perrone GG, Iglseder B, Rinnerthaler M, Dawes IW: Oxidative stress and neurodegeneration: the yeast model system. Front Biosci (Landmark Ed). 2013, 18: 1174-1193. 10.2741/4171.

    Article  CAS  Google Scholar 

  69. Ralser M, Wamelink MM, Kowald A, Gerisch B, Heeren G, Struys EA, Klipp E, Jakobs C, Breitenbach M, Lehrach H, Krobitsch S: Dynamic rerouting of the carbohydrate flux is key to counteracting oxidative stress. J Biol. 2007, 6: 10-10.1186/jbiol61.

    Article  PubMed Central  PubMed  Google Scholar 

  70. Tosato V, Grüning N, Breitenbach M, Arnak R, Ralser M, Bruschi CV: Warburg effect and translocation-induced genomic instability: two yeast models for cancer cells. Front Oncol. 2012, 2: 212-

    PubMed Central  PubMed  Google Scholar 

  71. Aung-Htut MT, Ayer A, Breitenbach M, Dawes IW: Oxidative stresses and ageing. Subcell Biochem. 2012, 57: 13-54.

    Article  CAS  PubMed  Google Scholar 

  72. Darnell JC, Van Driesche SJ, Zhang C, Hung KY, Mele A, Fraser CE, Stone EF, Chen C, Fak JJ, Chi SW, Licatalosi DD, Richter JD, Darnell RB: FMRP stalls ribosomal translocation on mRNAs linked to synaptic function and autism. Cell. 2011, 146: 247-261. 10.1016/j.cell.2011.06.013.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  73. Gkogkas CG, Khoutorsky A, Ran I, Rampakakis E, Nevarko T, Weatherill DB, Vasuta C, Yee S, Truitt M, Dallaire P, Major F, Lasko P, Ruggero D, Nader K, Lacaille J, Sonenberg N: Autism-related deficits via dysregulated eIF4E-dependent translational control. Nature. 2013, 493: 371-377.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

  74. Greenough WT, Klintsova AY, Irwin SA, Galvez R, Bates KE, Weiler IJ: Synaptic regulation of protein synthesis and the fragile X protein. Proc Natl Acad Sci USA. 2001, 98: 7101-7106. 10.1073/pnas.141145998.

    Article  PubMed Central  CAS  PubMed  Google Scholar 

Download references

Acknowledgements

We thank all the individuals and their families who have participated in research for their cooperation and the professionals for collecting data. We thank C. Wirth and S. Karolus for excellent technical assistance. This study was partially funded by grant Po 255/17-4 of the Deutsche Forschungsgemeinschaft to Fritz Poustka and grant DEBRA Austria Projekt P_147200_09 to Johann W. Bauer and Lore Breitenbach-Koller.

Dedication

We dedicate this paper to our late friend and mentor Annemarie Poustka.

Author information

Authors and Affiliations

Authors

Corresponding authors

Correspondence to Sabine M Klauck or Lore Breitenbach-Koller.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors' contributions

AGC co-designed and co-conceived most experiments and together with DH performed most experiments. AGC performed statistical analyses and wrote the manuscript. AGC and DH are joint first authors. MBo performed preliminary experiments for yeast analyses. TK performed the literature search. CMF and FP recruited and diagnosed the patient sample. BS provided expertise on software for analyses of two-dimensional-DIGE data. JWB and HH provided expertise for use and the antibodies for validation experiments. JK and FL performed mass spectrometric experiments and analyzed mass spectrometry data. SW and MBr contributed to study design and provided experimental expertise. LBK and SMK conceived and designed experiments, oversaw the study, interpreted results and wrote the manuscript. SMK and LBK are joint senior authors and project co-leaders. All authors read and approved the final manuscript.

Andreas G Chiocchetti, Denise Haslinger contributed equally to this work.

Electronic supplementary material

13229_2013_106_MOESM1_ESM.doc

Additional file 1: Additional Methods. mRNA expression and Western blot analysis. Figure S1. Experimental setup. Figure S2. Hierarchical structure of gene enrichment analysis results. Table S1. List of samples used in this study. Table S5. Literature search for redox-sensitive candidates. (DOC 649 KB)

Additional file 2: Table S7: Gene enrichment analysis of identified candidate proteins. (XLS 118 KB)

13229_2013_106_MOESM3_ESM.xls

Additional file 3: Table S2: Two-dimensional-DIGE results of RPL10[H213Q] carriers versus wild-type allele carriers. (XLS 694 KB)

Additional file 4: Table S4: List of candidates discussed in this publication. (XLS 69 KB)

Additional file 5: Table S3: Two-dimensional-DIGE results of ASD versus controls. (XLS 723 KB)

Additional file 6: Table S8: mRNA analysis of identified candidate proteins. (XLS 97 KB)

Additional file 7: Table S9: Western blot analysis of identified candidate proteins. (XLS 1 MB)

Additional file 8: Table S6: Two-dimensional-DIGE results of yeast proteomic analysis. (XLS 704 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Open Access This article is published under license to BioMed Central Ltd. This is an Open Access article is distributed under the terms of the Creative Commons Attribution License ( https://creativecommons.org/licenses/by/2.0 ), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Chiocchetti, A.G., Haslinger, D., Boesch, M. et al. Protein signatures of oxidative stress response in a patient specific cell line model for autism. Molecular Autism 5, 10 (2014). https://doi.org/10.1186/2040-2392-5-10

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/2040-2392-5-10

Keywords