Abstract
Chronic fatigue syndrome (CFS) is a debilitating and complex disorder characterized by unexplained fatigue not improved by rest. An area of investigation is the likely connection of CFS with defective mitochondrial function. In a previous work, we investigated the proteomic salivary profile in a couple of monozygotic twins discordant for CFS. Following this work, we analyzed mitochondrial proteins in the same couple of twins. Nano-liquid chromatography electrospray ionization mass spectrometry (nano-LC-MS) was used to study the mitochondria extracted from platelets of the twins. Subsequently, we selected three proteins that were validated using western blot analysis in a big cohort of subjects (n=45 CFS; n=45 healthy), using whole saliva (WS). The selected proteins were as follows: aconitate hydratase (ACON), ATP synthase subunit beta (ATPB) and malate dehydrogenase (MDHM). Results for ATPB and ACON confirmed their upregulation in CFS. However, the MDHM alteration was not confirmed. Thereafter, seeing the great variability of clinical features of CFS patients, we decided to analyze the expression of our proteins after splitting patients according to clinical parameters. For each marker, the values were actually higher in the group of patients who had clinical features similar to the ill twin. In conclusion, these results suggest that our potential markers could be one of the criteria to be taken into account for helping in diagnosis. Furthermore, the identification of biomarkers present in particular subgroups of CFS patients may help in shedding light upon the complex entity of CFS. Moreover, it could help in developing tailored treatments.
Similar content being viewed by others
Introduction
Chronic fatigue syndrome (CFS) is a debilitating and complex disorder characterized by unexplained fatigue not improved by rest, that lasts longer than 6 months, and that may be aggravated by physical or mental activity. In CFS, fatigue is also characterized by rheumatological, and neuropsychiatric symptoms such as joint pain, aching or stiff muscles, concentration and memory problems, sleep disruption, sore throat, headaches and so on.1, 2 In 2011, the International Consensus Criteria have been defined to identify CFS by revising the Canadian definition released in 2003;3 however, at present, no laboratory tests can be used to objectively diagnose CFS. Therefore, the diagnosis of CFS remains a clinical endeavor and can be made after the exclusion of medical conditions that may explain the prolonged fatigue as well as a number of psychiatric diagnoses (for example, bipolar disorder, eating disorders, psychotic disorders and melancholic depression).
This difficulty in defining CFS, due to the lack of diagnostic markers, complicates epidemiological studies. Prevalence varies from as low as 0.2% to as high as 6.41%,4, 5 with an incidence rate that differed strongly with age, with a first peak in the age group 10â19 years and a second peak in the age group 30â39 years.6, 7
Up till now, the etiology of CFS remains unclear; it is likely that CFS is a heterogeneous disorder that could be explained by a multifactorial model where genetic, infectious, neuroendocrine and psychological factors converge.
One of the most common proposed causative agent is a viral infection, and vaccines have been also indicated to contribute to the development of CFS; however, data are controversial.8, 9, 10, 11, 12 In addition, several immune abnormalities have been reported in CFS patients; however, many of the observed abnormalities were not confirmed by others, resulting in inconsistent findings across studies.13, 14 Hence, there is a clear need to enhance CFS understanding.
Another area of investigation is the likely connection of CFS with defective mitochondrial function. Recently, studies on metabolism and CFS suggest irregularities in energy metabolism and oxidative stress metabolism.15 In particular, it is increasingly recognized that CFS can be underpinned by mitochondrial dysfunctions;16, 17, 18 moreover, fatigue is a prominent symptom in patients with mitochondrial disease.19 This is an attractive hypothesis, as a metabolic alteration at the mitochondrial level could entail an energy deficiency and explain the peculiar fatigue.
Hitherto, it has often been suggested that CFS may have a genetic predisposition; therefore, twin studies supply the best matched controls in a caseâcontrol research. Thereby, in a previous work, we investigated the proteomic salivary profile in a couple of monozygotic twins discordant for CFS.20 In particular, our results highlighted the involvement of inflammatory response in CFS pathogenesis. Following the encouraging results obtained with this preliminary study, we sought to improve our knowledge about mechanisms of CFS by analyzing mitochondrial proteins in the same couple of monozygotic twins discordant for CFS. The expectation was that our findings could be able to provide evidence of a mitochondrial role in CFS pathogenesis.
With this purpose, we isolated mitochondria from platelets. Platelets are circulating blood cells that are essential not only in hemostasis and coagulation but also in innate and adaptive immunity.21 Platelets have functionally active mitochondria, whose involvement in cell physiology and in pathological conditions has been increasingly demonstrated.21, 22 The recent assessments of mitochondrial function in platelets offer several advantages for diagnostic purposes, such as their abundance in peripheral blood, the relative ease of obtaining platelets, their limited complexity,23 the considerable number of mitochondria contained in platelets and the ability to make repeated measurements for each patient.21
Platelets have already been used in investigations of mitochondrial involvement in some disorders, for example, in sepsis,24 neurodegenerative diseases25, 26, 27 or diabetes.28, 29 However, these studies were limited to mitochondria DNA analysis and electron transfer activities, whereas no investigation has been performed on the total mitochondrial protein content. Therefore, we analyzed with nano-liquid chromatography electrospray ionization mass spectrometry (nano-LC-MS) the mitochondria extracted from platelets of the same couple of monozygotic twins discordant for CFS who have been the subjects of our previous work. Subsequently, the most promising biomarkers were validated by western blot analysis in a big cohort of patients, using whole saliva (WS), whose sampling is a noninvasive, simple, safe and stress-free procedure that can be easily applied to large groups of subjects.
Seeing the vast heterogeneity of the disease, the aim was to suggest biomarkers to stratify CFS cases, rather than selecting some biomarkers for identifying CFS. This can help to customize intervention.
Materials and methods
Study design
This translational caseâcontrol study was subdivided into three different steps. The first âdiscoveryâ phase was aimed at characterizing the mitochondrial proteomic profile of a patient suffering from CFS in comparison with his healthy monozygotic twin. This purpose was achieved by using a proteomics approach: nano-LC-MS. The second âexploratoryâ phase was aimed at investigating the changing proteins; we applied the Ingenuity Pathway Analysis (IPA) Knowledge base. This platform enabled to visualize the potential interactions between the identified biomarker signatures in mitochondria helping us to select some candidates for the âvalidationâ phase. Selected proteins were first searched in WS of the twin brothers and, once the differences have been confirmed, the validation was widened to a larger number of patients and controls. Finally, we have amplified our analysis looking closely at clinical data, with the specific aim to investigate how the biomarkers correlate with patientsâ clinical parameters. This leads us to identify mitochondrial proteins useful in subtyping patients with CFS.
Patients
The recruitment of the subjects was managed by the Division of Rheumatology, Department of Clinical and Experimental Medicine, University of Pisa, Pisa, Italy. This study has been approved by the local Ethics Committee, and an informed consensus was obtained for diagnostic or clinical purposes.
The twins of the previous study20 were enrolled to collect peripheral venous blood samples for the purpose of âdiscoveryâ phase, and saliva for the validation. Samples, both blood and WS, were collected three times from the twins. In addition, a total of 45 patients affected by CFS were consecutively recruited for the âvalidationâ phase; mean age 38.6±11.9 years (mean±s.d.; M±s.d.; Table 1). Forty-five healthy subjects, with similar mean age (36.4±10.9, M±s.d.; 25 females, 20 males) and demographic characteristics, were included as controls (see Supplementary Methods 1). The patients were classified based on the classification criteria of Fukuda et al.30
The patients performed a rheumatologic visit with routine clinical evaluation of medical history. The subjects completed the Fibromyalgia Impact Questionnaire (FIQ), in which pain and fatigue severity were assessed with 10-cm Visual Analog Score Scale (VAS) and short form-36. Fatigue was assessed also by means of the Functional Assessment of Chronic Illness Therapy-Fatigue Scale (FACIT). We considered the VAS from the FIQ asking whether the patient felt rested upon awaking during the last week. In order to assess the presence of psychiatric concomitant disorders, patients had a psychiatric evaluation based on the administration of the Structured Clinical Interview for Diagnostic and Statistical Manual of Mental Disorders, 4th Edition, axis-I disorders.31
To better classify the subjects the following hematochemical parameters were evaluated:
-
Immunological parameters: antinuclear antibodies, extractable nuclear antigen antibodies, rheumathoid factor (Ra test) and anticardiolipin antibodies.
-
Hormone profile: cortisole and adrenocorticotropic hormone, serotonin, growth hormone, insulin-like growth factor-1 and vitamin D3.
-
Cytokine profile: tumor necrosis factor (TNF)-α, interferon (IFN)-γ, interleukin (IL)-6, IL-10, IL-2, IL-8 and IL-1b.
Table 1 summarizes clinical data.
WS collection
Salivary samples were collected from patients (n=45) and controls (n=45) with a saliva collector sponge (Surescreen Diagnostics, Derby, UK). WS samples were collected early in the morning (between 0800 and 1100 hours) according to a standard protocol.32 No evidence of oral pathologies or inflammatory processes was observed. The collected saliva was immediately centrifuged at 17â000âg for 20âmin at 4â°C to yield clear samples. The resulting supernatants were stored at â80â°C. Protein amounts of samples were determined using the Bio-Rad DC-protein assay (Bio-Rad, Hercules, CA, USA). Bovine serum albumin was used as a standard.
Platelet collection and mitochondrial enrichment
Blood samples were drawn from an antecubital vein using a 19-gauge needle and the first 2.5âml were discarded. All samples were collected in Vacutainer tubes (Becton Dickinson, San Jose, CA, USA), containing 3.8% sodium citrate as anticoagulant.
Peripheral venous blood sample (25âml) was drawn from the subjects and placed in proper collecting tubes with heparin. The samples were immediately processed to obtain platelets. Blood was centrifuged for 15âmin at 300âg to obtain a platelet-rich plasma and a pellet-containing red cells, granulocytes and lymphocytes. Platelets were precipitated from platelet-rich plasma by centrifugation at 3000âg for 10âmin and resuspended in 0.9% NaCl; an aliquot of 300âμl was stored for normalization, whereas the remaining was used for mitochondrial enrichment.
Isolation of mitochondrial fractions from human platelets was performed by differential centrifugation. The platelets were centrifuged for 10âmin at 3000âg. A volume of cold lysis buffer (10âmM Tris-HCl pH 7.6), depending on the size of the pellet, was added and kept 7âmin on ice stirring occasionally. After centrifugation for 10âmin at 1000âg we obtained unbroken cells (P1). The supernatant, containing mitochondria, was immediately centrifuged for 20âmin at 9000âg obtaining mitochondrial fraction M1. P1 was further resuspended in cold lysis buffer, and kept for 7âmin on ice stirring occasionally. After a centrifugation for 10âmin at 1000âg, the resulting supernatant was immediately centrifuged for 20âmin at 9000âg obtaining mitochondrial fraction M2. M1 and M2 were immediately resuspended in cold buffer (10âmM Tris-HCl pH 7.4, 2âmM EDTA, 250âmM Sucrose) and combined. Protein amounts were determined using the Bio-Rad DC-protein assay (Bio-Rad). Bovine serum albumin was used as a standard.
The integrity of mitochondria was assayed by measuring the activity of cytochrome oxidase. Western blot (WB) analysis was performed to evaluate the purity of the mitochondria as described in Supplementary Method 2. Blood was collected three times, and mitochondria obtained from each experiments were pooled and analyzed in triplicate by LC-MS-based label-free quantification.
LC-MS-based label-free quantification
Sample preparation for LC-MS
Platelets, kept for normalization, and mitochondria were analyzed by label-free LC-MS analysis. Samples were resuspended in the lysis buffer (50âmM Tris pH 7.8, 150âmM NaCl, 1% sodium dodecylsulphate (SDS)) supplemented with protease inhibitors and phosphatase inhibitors (Roche, Mannheim, Germany; 1 tablet per 10âml buffer) and sonicated. After centrifugation (18â000âg for 15âmin at 4â°C), the protein amount of the lysate was determined using the bicinchoninic acid protein assay according to the manufacturerâs instructions (Pierce BCA Protein Assay Kit, Thermo Scientific, Bremen, Germany). Cysteines were reduced with 10âmM dithiothreitol at 56â°C for 30âmin and the free thiols were carbamidomethylated with 30âmM iodoacetamide at room temperature for 30âmin in the dark.
Sample cleaning and proteolysis were based on the filter-aided sample preparation protocol.33, 34 Briefly, cell lysates corresponding to 135âμg of protein were diluted up to 500âμl with freshly prepared 8âM urea/100âmM Tris-HCl (pH 8.5) buffer.35 Diluted samples were placed on the centrifugal devices (PALL microsep, 30âKDa cutoff) and were centrifuged at 13â800âg at room temperature for 20âmin. All the following centrifugation steps were performed under similar conditions. To eliminate residual SDS, three washing steps were carried out using 100âμl of 8.0âM urea/100âmM Tris-HCl buffer, pH 8.5 and finally for the buffer exchange, the devices were washed three times with 100âμl of 50âmM NH4HCO3 buffer, pH 7.8. To the concentrated proteins, 100âμl of proteolysis buffer comprising trypsin (Promega, Madison, WI, USA; 1:25 w/w, enzyme to protein), 0.2âM guanidine hydrochloride and 2âmM CaCl2 in 50âmM NH4HCO3, pH 7.8, were added and incubated at 37â°C for 14âh. The generated tryptic peptides were recovered by centrifugation with 50âμl of 50âmM NH4HCO3, followed by 50âμl of ultrapure water. Finally, the peptides were acidified with 10% trifluoroacetic acid to pH<3, and the digests were quality-controlled as described previously.36
LC-MS analysis
Each condition was measured in triplicate (500âng each) using an Ultimate 3000 nanoRSLC system coupled to an Orbitrap Elite mass spectrometer (both Thermo Scientific). Samples were analyzed in a randomized order to minimize systematic errors. Briefly, peptides were preconcentrated on a 100âμm à 2âcm C18 trapping column for 10âmin using 0.1% trifluoroacetic acid at a flow rate of 20âμlâminâ1, followed by separation on a 75âμm à 50âcm C18 analytical column (both Pepmap, Thermo Scientific) with a 90-min LC gradient ranging from 3 to 42% of buffer B (84% acetonitrile, 0.1% formic acid) at a flow rate of 250ânlâminâ1. MS survey scans were acquired in the Orbitrap from m/z 300 to 1500 at a resolution of 60â000 using the polysiloxane ion at m/z 371.101236 as lock mass.37 The 15 most intense signals were subjected to collision-induced dissociation in the ion trap, taking into account a dynamic exclusion of 30âs. Collision-induced dissociation spectra were acquired with a normalized collision energy of 35% and an activation time of 10âms. Automatic gain control target values were set to 1 à 106 for Orbitrap MS and 1x104 for ion trap MS/MS scans, and maximum injection times were set to 100âms for both full MS and MS/MS scans.
Label-free data analysis
Label-free quantification of the acquired MS data was performed using the Progenesis LC-MS software from Nonlinear Dynamics (Newcastle upon Tyne, UK) version 4.1. MS data-processing, including alignment of raw data, selection of the reference LC-MS run and peak picking, was done automatically by Progenesis. The features within retention time and m/z windows from 0 to 90âmin and 300â1500âm/z with charge states +2, +3 and +4 were considered for peptide statistics, analysis of variance and principal component analysis. Spectra were exported as peak lists, searched against a concatenated target/decoy version of the human Uniprot database, (downloaded on 11 December 2013, containing 20â273 target sequences) using Mascot 2.4 (Matrix Science, Boston, MA, USA). Trypsin was selected as enzyme with a maximum of two missed cleavage sites, carbamidomethylation of Cys was set as fixed and oxidation of Met was selected as variable modification. MS and MS/MS tolerances were set to 10âp.p.m. and 0.5âDa, respectively. Peptide Shaker version 0.28.0 (ref. 38) was used for processing peptide-spectrum matches obtained from Mascot and filtering the results at a false discovery rate of 1% on the peptide-spectrum match, peptide and protein level. Results were exported using the advanced Peptide Shaker features, allowing direct re-import of the quality-controlled data into Progenesis. Only proteins that were quantified with ⩾2 unique peptides were considered for further analysis. The signal intensities of each protein present in the mitochondrial fraction were normalized on the total content of the same proteins obtained from the analysis of total platelets.
Signaling pathway analysis
Proteins differentially expressed were functionally analyzed through the use of QIAGENâs IPA (Ingenuity System, QIAGEN, Redwood City, CA, USA) to select candidates for validation.
All differentially expressed proteins with a fold variation in CFS greater than 2.0 with respect to control were included in bioinformatic analysis to identify molecular functions most strongly associated with the protein list. Supplementary Method 3 gives more details about this analysis.
WB for validation
To validate the expression changes in WS, we used WB analysis. The selected proteins were as follows: aconitate hydratase (ACON), ATP synthase subunit beta (ATPB) and malate dehydrogenase (MDHM). For ACON and MDHM, we loaded 30âμg of proteins and 10âμg for ATPB.
Proteins were resolved by 12% SDS-PAGE gels and transferred on nitrocellulose membranes as previously described.39 ACON and MDHM antibodies have been used at 1:1000 dilution (Cell Signaling Technology, Danvers, MA, USA), whereas the ATPB antibody (Santa Cruz Biotechnology, Dallas, TX, USA) has been used at 1:200 dilution. A horseradish peroxidase-conjugated goat anti-rabbit (Stressgen, Farmingdale, NY, USA) was used as a secondary antibody at 1:10â000 dilution. Immunoblots were developed using the enhanced chemiluminescence detection system. The chemiluminescent images were acquired using LAS4010 (GE Health Care, Uppsala, Sweden). The antigen-specific bands were quantified using the Image Quant-L (GE Health Care) software. In order to normalize the optical density of immune-reactive bands, the optical density of total proteins was calculated. Therefore, immediately after WB, the membranes were stained with 1âmM RuBP as previously described.39
Statistical analysis
All data are presented as mean±s.e.m. Comparisons between groups were performed using the MannâWhitney U-test for non-normal data. A P-value<0.05 was considered significant. Linear regression analysis was used to determine the correlation among levels of different biomarkers. Logistic regression was used to determine the weight given to each marker and then to calculate a specific formula to provide a combined risk index. To estimate whether this marker combination might increase the marker's performance in CFS diagnosis, receiver operating characteristic (ROC) curves were plotted, and the areas under curves (AUCs) were calculated with their 95% confidence intervals using standard techniques to evaluate sensitivity and specificity of each marker and their combination. Statistical analyses were performed with GraphPad Prism (San Diego, CA, USA) and SPSS (SPSS/PC Statistical Package for the Social Science, update for 10.1., SPSS, Chicago, IL, USA, 2000).
Clinical correlations
To determine the statistical correlations among putative biomarkers and serological measures, the Spearmanâs rank correlation coefficient was calculated, a non-parametric measure of correlation based on data ranks. A P-value<0.05 was considered significant. Clinical correlations were performed with SPSS.
Results
Proteomic analysis
In order to profile the mitochondrial proteins in CFS, we purified mitochondria from platelets of the same couple of monozygotic twins, discordant for CFS, which was investigated in our previous work.20 To evaluate the quality of the isolated mitochondria, enzymatic assay and WB analysis were employed (see Supplementary Results 1), confirming that we reached satisfactory mitochondrial enrichment (Supplementary Figure). We obtained a mitochondrial fraction of 450âμg of proteins starting roughly from 5 à 109 platelets.
Moreover, with the aim to overcome the limit of reproducibility in the mitochondrial preparation, experiments were performed three times. Mitochondria obtained from each of the three experiments were pooled together and analyzed in triplicate.
The protein composition of mitochondria was analyzed by nano-LC-MS and a differential analysis between CFS patient and his healthy brother was performed. Overall, we identified 1007 proteins with a number of peptides used for quantitation ⩾2. Out of 1007 proteins detected using the label-free approach, 194 mitochondrial proteins were significantly modified in CFS. Among these 194 proteins, 41 had a fold variation in CFS with respect to control greater than 2.0, 34 were upregulated in CFS and 7 were downregulated. Table 2 contains data of these 41 proteins. Table containing the remaining 153 proteins is supplied as Supplementary Table S1.
These 41 proteins were included in the analysis with IPA. Each identified protein was converted to its gene and mapped to its corresponding gene object in the IPA knowledge base.
We utilized IPA to retrieve the known functions of each protein. As one protein may have multiple functions, we selected the functions with the highest statistical significance. The first three top-ranked biological functions were as follows: metabolism of isocitric acid (P-value=1.83eâ08), metabolism of NADH (P-value=4.00eâ07) and metabolism of nucleic-acid component or derivative (P-value=8.15eâ07). These significantly over-represented categories, together with the related molecules found in our study, are depicted in Figures 1a. For the validation we selected a protein for each category by considering both fold variation and the P-value of proteins from the proteomics analysis (see Supplementary Results 2 and Supplementary Tables S2).
Expression of biomarkers in WS
We selected three mitochondrial proteins for the validation on WSâACON, ATPB and MDHM. First of all, we assessed whether their variation in CFS with respect to control was confirmed in WS of the twin brothers. As depicted in Figure 1b, the upregulation in CFS was corroborated by WB for all the proteins. ACON, ATPB and MDHM showed an increase of 3.1, 1.5 and 4.7, with a significant P-value of 0.006, 0.04 and 0.003, respectively (Table 3).
Thereafter, we moved to the analysis of these proteins in a vast cohort of CFS patients (n= 45) and healthy subjects (n= 45). We started with a global analysis of the expression in WS of these proteins before moving to the examination of different CFS subtypes.
Table 3 and Figure 1c summarize these results. The changes for ACON and ATPB were consistent with the results from nano-LC-MS, but not for MDHM.
Then, we classified patients according to their clinical features from questionnaires, that is, we examined the following: FACIT, FIQ, VAS_pain, VAS_fatigue and VAS_sleep. CFS patients were split into two groups considering as reference the twin suffering from CFS. One group with the clinical features whose values were in common with (or greater than) those of the ill twin (group A), and the second with the patients with lower values (group B). Intensities were higher for each of the three markers in group A with respect to group B (Figure 2). Details are given in SR3.
ROC curves
ROC curves were calculated to assess the clinical potential of our selected and validated proteins to distinguish CFS from healthy subjects in WS. The AUCs were calculated for each protein individually, showing whether each marker alone can discriminate CFS from healthy subjects. ATPB and ACON were better differentiating biomarkers than MDHM. ATPB had an AUC of 0.700, with sensitivity=54% and specificity=78%. ACON had an AUC of 0.738, with sensitivity=61% and specificity=78%. Moreover, using a logistic regression analysis, we investigated whether the discriminative power of each marker could be potentially increased with the combination of various markers. With this purpose, we tested all the different combinations in order to select the best association of biomarkers useful to discriminate CFS from healthy subjects. Thus, we found that the discriminative power slightly increased if ATPB was combined with ACON, as shown in Figure 1d. Indeed, the AUC increased to 0.793 and the sensitivity to 85%, but the specificity was 72%. Finally, we also calculated ROC curves for testing the ability of biomarkers to differentiate patients from controls according to the clinical parametersâonly FACIT gave good results. When we split patients according to FACIT, for ATPB and ACON we obtained a ROC curve (group A versus controls) with an AUC of 0.780 and 0.722, respectively, whereas the combined ROC curve had an AUC of 0.844 (sensitivity=93% and specificity=70%).
Clinical correlations
Finally, statistical correlations were carried out calculating the Spearmanâs rank correlation coefficient. The considered clinical parameters were as follows: age, age of onset, duration of disease and the hematochemical parameters (for example, TNF-α, IL-6 and so on; see the âPatientsâ section). ACON correlated negatively with IL-2, with a Spearmanâs rank correlation coefficient of â0.458 (P-value=0.010). The levels of ATPB correlated positively with TNF-α and IFN-γ. The Spearmanâs rank correlation coefficient was 0.352 (P-value=0.048) and 0.516 (P-value=0.004) for TNF-α and IFN-γ, respectively. MDHM also correlates with TNF-α; coefficient 0.447, P-value=0.013.
Discussion
The discovery and validation of more objective markers could provide an important support for CFS diagnosis and therapy. As evidence suggests that CFS can be underpinned by mitochondrial dysfunctions, we analyzed mitochondrial proteins obtained from platelets. Platelets were from a couple of monozygotic twins discordant for CFS, which were the subjects of our previous study.20 Nano-LC-MS highlighted the changes in expression of 41 proteins with a P-value<0.05 and fold of difference greater than 2. Among these proteins, we selected for the validation three proteins from each of the three most significant âbiological functionsâ associated with our proteins pointed out after IPA analysis.
These 3 proteins were: ACON for âmetabolism of isocitric acidâ, MDHM for âmetabolism of NADHâ and ATPB for the âmetabolism of nucleic-acid component or derivativeâ. Moreover, validation of our candidate protein biomarkers was restricted by the availability of antibodies that can allow detection of proteins in saliva samples.
Despite this study analyzed only two subjects, the differences of protein expression could be linked to the disease itself, as the twins differed only regarding the presence of CFS.20 By this way, our study on a couple of twins could permit to focus on some biomarkers that can be searched in a wide number of patients in order to confirm their usefulness in CFS diagnosis. We decided to explore the WS of CFS patients preferring to adopt a less invasive procedure in collecting the samples. Although WS is a less complex fluid, it has been proven to be a promising diagnostic tool, which reflects the physiological function of the body.40, 41First of all, it was mandatory to check whether the significant increase of these proteins in CFS, with respect to control, could be verified in WS of the twins brothers. Actually, we confirmed the results from nano-LC-MS.
The pathway analysis highlighted the âmetabolism of NADHâ as one of the most important biological functions involved in CFS. The NAD+/NADH ratio has an omnipresent role in regulating the intracellular redox status and, therefore, represents a function of the metabolic state. Some studies suggest that NADH concentrations are significantly lower in CFS patients compared with healthy controls.42, 43, 44 It is likely that high levels of MDHM and isocitrate dehydrogenase observed in CFS may be an adaptive change to deficiencies of NADH which, in turn, has a critical role in mitochondrial ATP production.
ATPB is the subunit beta of the ATP synthase, the key enzyme for ATP production. In addition to the increase of ATPB in CFS, we also detected significant increases of both subunits alpha and gamma using nano-LC-MS (P-value of 0.005 and 0.001 respectively; Supplementary Table S1). Tiredness and energy are tightly connected; therefore, an attractive hypothesis is that there is a metabolic dysfunction with the result that enough energy is not being produced. In this perspective the increase can be explained as an attempt to increase the ATP production. Indeed, there are studies suggesting that the decrease in mitochondrial ATP synthesis in CFS patients is not caused by a defect in the enzyme complexes catalyzing oxidative phosphorylation, but in another factor.45 In order for ADP to be recycled, and for ATP to pass into the cytosol from the matrix, the exchange of cytoplasmic ADP with mitochondrial ATP across the mitochondrial inner membrane is needed. This function is provided by a translocator protein. Consistent with other lines of work,17, 46 mitochondrial ADP/ATP translocase expression is modified in our study: translocase 2 and 3 decreased in CFS (P-value 0.009 and 0.002, respectively; Supplementary Table S1).
Together with ATPB, the most promising biomarker seems to be ACON. This is the mitochondrial form of aconitase, an enzyme that catalyses the stereo-specific isomerization of citrate to isocitrate via cis-aconitate in the tricarboxylic acid cycle. ACON is commonly used as a biomarker for oxidative stress and has been suggested to serve as a sensor of redox status.47 Indeed, the reaction between ACON and superoxide generates free hydroxyl radical. It is likely that the reaction between superoxide and ACON may enhance mitochondrial oxidative damage associated with the pathophysiology of CFS.48 As a matter of fact, the excessive production of ROS and nitrogen species is usually related to CFS,16 and we found a significant upregulation of superoxide dismutase (P-value=0.0007; Table 2).
After the validation on the WS of twins, we have broadened the validation to a big cohort of CFS patients. Results for ATPB and ACON corroborated the upregulation in CFS, and their clinical potential to discriminate patients from controls was evaluated. ROC curves highlighted good levels of specificity but low sensitivity. If the biomarkers were combined, we obtained a better performance in diagnosis, and the sensitivity increased to 85%, despite a little decrease in specificity.
On the other hand, the MDHM alteration found in twins was not confirmed in the global cohort of patients, but we could observe its increase when we considered only the patients who had similar or higher clinical parameter values than the twin affected by CFS. This finding brought us to inspect whether other putative biomarkers could also demarcate distinct CFS subgroups, seeing the great variability of clinical features of our patients as demonstrated by Table 1. Therefore, we have amplified our analysis looking closely at clinical data, with the specific aim to investigate how the biomarkers were linked with patientsâ clinical characteristics. We analyzed the expression of our proteins after splitting patients according to clinical parameters. For each marker, the values were actually higher in the group of patients who had clinical features similar to the ill twin. It seemed that the increase in expression of these proteins was corresponding to a worse patient outcome.
In addition, the combination of ATPB with ACON gave a ROC curve with sensitivity of 93% and specificity of 70%, when we checked patients with values of FACIT greater than equal to those of CFS twin.
It is noteworthy also to point out that CFS patients with comorbid bipolar disorder, and not with other psychiatric comorbidities, showed a significant upregulation of ATPB expression (P-value=0.03, data not shown). In our cohort of patients only eight CFS patients had comorbid bipolar disorder (Table 1); anyway, this aspect deserves future investigation, considering that hampered mitochondrial function in bipolar patients has been reviewed.49, 50
Finally, we investigated the correlation of our proteins with hematochemical parameters. Cytokine abnormalities are common in CFS, but the findings are mixed and no individual marker is used as an objective biomarker for the diagnosis or management of CFS. It is noteworthy that the increase of ATPB in CFS correlated positively with the increase in IFN-γ and TNF-α. Immunological dysregulation has been always proposed as a significant component of the CFS, and CFS patients had usually high levels of pro-inflammatory cytokines, such as IFN-γ and TNF-α.14, 51, 52 In addition, MDHM correlated positively with TNF-α. These correlations are instructive, as all these variables can be thought as markers. It is interesting to point out that ACON levels correlated negatively with IL-2, in line with the fact that the twin with CFS had lower levels of IL-2 with respect to his healthy twin brother.20
In conclusion, these results suggest that our potential markers could be one of the criteria to be taken into account for helping in diagnosis. Furthermore, even if these proteins may not necessarily be useful for a universal diagnosis of CFS, given the huge heterogeneity of the illness, they offer a strong potential for characterizing subsets. The identification of biomarkers present in particular subgroups of CFS patients may help in shedding light upon the complex entity of CFS. Moreover, it could help in developing tailored treatments.
By this way, this work may contribute to the continuous development of a more clear framework for a complex and nuanced disease such as CFS.
References
Werker CL, Nijhof SL, van de Putte EM . Clinical Practice: chronic fatigue syndrome. Eur J Pediatr 2013; 172: 1293â1298.
Yancey JR, Thomas SM . Chronic fatigue syndrome: diagnosis and treatment. Am Fam Physician. 2012; 86: 741â746.
Carruthers BM, van de Sande MI, De Meirleir KL, Klimas NG, Broderick G, Mitchell T et al. Myalgic encephalomyelitis: international consensus criteria. J Intern Med 2011; 270: 327â338.
Johnston S, Brenu EW, Staines D, Marshall-Gradisnik S . The prevalence of chronic fatigue syndrome/myalgic encephalomyelitis: a meta-analysis. Clin Epidemiol 2013; 5: 105â110.
Brurberg KG, Fønhus MS, Larun L, Flottorp S, Malterud K . Case definitions for chronic fatigue syndrome/myalgic encephalomyelitis (CFS/ME): a systematic review. BMJ Open 2014; 4: e003973.
Bakken I, Tveito K, Gunnes N, Ghaderi S, Stoltenberg C, Trogstad L et al. Two age peaks in the incidence of chronic fatigue syndrome/myalgic encephalomyelitis: a population-based registry study from Norway 2008â2012. BMC Med 2014; 12: 167.
Pheby D, Sneddon P, Heinrich I . Severe ME/CFS in adults - a report from the CHROME database. Bull IACFS/ME 2009; 17: 143â167.
Donegan K, Beau-Lejdstrom R, King B, Seabroke S, Thomson A, Bryan P . Bivalent human papillomavirus vaccine and the risk of fatigue syndromes in girls in the UK. Vaccine 2013; 31: 4961â4967.
Rosenblum H, Shoenfeld Y, Amital H . The common immunogenic etiology of chronic fatigue syndrome: from infections to vaccines via adjuvants to the ASIA syndrome. Infect Dis Clin North Am 2011; 25: 851â863.
Magnus P, Brubakk O, Nyland H, Wold BH, Gjessing HK, Brandt I et al. Vaccination as teenagers against meningococcal disease and the risk of the chronic fatigue syndrome. Vaccine 2009; 27: 23â27.
Shapiro JS . Does varicella-zoster virus infection of the peripheral ganglia cause chronic fatigue syndrome? Med Hypotheses 2009; 73: 728â734.
Exley C, Swarbrick L, Gherardi RK, Authier FJ . A role for the body burden of aluminium in vaccine-associated macrophagic myofasciitis and chronic fatigue syndrome. Med Hypotheses 2009; 72: 135â139.
Nijs J, Nees A, Paul L, De Kooning M, Ickmans K, Meeus M et al. Altered immune response to exercise in patients with chronic fatigue syndrome/myalgic encephalomyelitis: a systematic literature review. Exerc Immunol Rev 2014; 20: 94â116.
Brenu EW, Huth TK, Hardcastle SL, Fuller K, Kaur M, Johnston S et al. Role of adaptive and innate immune cells in chronic fatigue syndrome/myalgic encephalomyelitis. Int Immunol 2014; 26: 233â242.
Armstrong CW, McGregor NR, Butt HL, Gooley PR . Metabolism in chronic fatigue syndrome. Adv Clin Chem 2014; 66: 121â172.
Morris G, Maes M . Mitochondrial dysfunctions in myalgic encephalomyelitis/chronic fatigue syndrome explained by activated immuno-inflammatory, oxidative and nitrosative stress pathways. Metab Brain Dis 2014; 29: 19â36.
Booth NE, Myhill S, McLaren-Howard J . Mitochondrial dysfunction and the pathophysiology of Myalgic Encephalomyelitis/Chronic Fatigue Syndrome (ME/CFS). Int J Clin Exp Med 2012; 5: 208â220.
Myhill S, Booth NE, McLaren-Howard J . Chronic fatigue syndrome and mitochondrial dysfunction. Int J Clin Exp Med 2009; 2: 1â16.
Gorman GS, Elson JL, Newman J, Payne B, McFarland R, Newton JL et al. Perceived fatigue is highly prevalent and debilitating in patients with mitochondrial disease. Neuromuscul Disord 2015; 25: 563â566.
Ciregia F, Giusti L, Da Valle Y, Donadio E, Consensi A, Giacomelli C et al. A multidisciplinary approach to study a couple of monozygotic twins discordant for the chronic fatigue syndrome: a focus on potential salivary biomarkers. J Transl Med 2013; 11: 243.
Garcia-Souza LF, Oliveira MF . Mitochondria: biological roles in platelet physiology and pathology. Int J Biochem Cell Biol 2014; 50: 156â160.
Zharikov S, Shiva S . Platelet mitochondrial function: from regulation of thrombosis to biomarker of disease. Biochem Soc Trans 2013; 41: 118â123.
Burkhart JM, Vaudel M, Gambaryan S, Radau S, Walter U, Martens L et al. The first comprehensive and quantitative analysis of human platelet protein composition allows the comparative analysis of structural and functional pathways. Blood 2012; 120: e73âe82.
Gründler K, Angstwurm M, Hilge R, Baumann P, Annecke T, Crispin A et al. Platelet mitochondrial membrane depolarization reflects disease severity in patients with sepsis and correlates with clinical outcome. Crit Care 2014; 18: R31.
Silva AC, Almeida S, Laço M, Duarte AI, Domingues J, Oliveira CR et al. Mitochondrial respiratory chain complex activity and bioenergetic alterations in human platelets derived from pre-symptomatic and symptomatic Huntington's disease carriers. Mitochondrion 2013; 13: 801â809.
Carelli V, Ghelli A, Ratta M, Bacchilega E, Sangiorgi S, Mancini R et al. Leber's hereditary optic neuropathy: biochemical effect of 11778/ND4 and 3460/ND1 mutations and correlation with the mitochondrial genotype. Neurology 1997; 48: 1623â1632.
Mann VM, Cooper JM, Krige D, Daniel SE, Schapira AH, Marsden CD . Brain skeletal muscle and platelet homogenate mitochondrial function in Parkinson's disease. Brain 1992; 115: 333â342.
Tang WH, Stitham J, Jin Y, Liu R, Lee SH, Du J et al. Aldose reductase-mediated phosphorylation of p53 leads to mitochondrial dysfunction and damage in diabetic platelets. Circulation 2014; 129: 1598â1609.
Avila C, Huang RJ, Stevens MV, Aponte AM, Tripodi D, Kim KY et al. Platelet mitochondrial dysfunction is evident in type 2 diabetes in association with modifications of mitochondrial anti-oxidant stress proteins. Exp Clin Endocrinol Diabetes 2012; 120: 248â251.
Fukuda K, Straus SE, Hickie I, Sharpe MC, Dobbins JG, Komaroff A . The chronic fatigue syndrome: a comprehensive approach to its definition and study. International Chronic Fatigue Syndrome Study Group. Ann Intern Med 1994; 121: 953â959.
First MB, Spitzer RL, Gibbon M, JBW Williams . Structured Clinical Interview for DSM-IV Axis I Disorders, Research Version, Patient Edition. (SCID-I/P). Biometrics Research, New York State Psychiatric Institute: New York, NY, USA, 2002.
Giusti L, Baldini C, Bazzichi L, Ciregia F, Tonazzini I, Mascia G et al. Proteome analysis of whole saliva: a new tool for rheumatic diseases-the example of Sjögrenâs syndrome. Proteomics 2007; 7: 1634â1643.
Manza LL, Stamer SL, Ham AJ, Codreanu SG, Liebler DC . Sample preparation and digestion for proteomic analyses using spin filters. Proteomics 2005; 5: 1742â1745.
WiÅniewski JR, Zougman A, Nagaraj N, Mann M . Universal sample preparation method for proteome analysis. Nat Method 2009; 6: 359â362.
Kollipara L, Zahedi RP . Protein carbamylation: in vivo modification or in vitro artefact? Proteomics 2013; 13: 941â944.
Burkhart JM, Schumbrutzki C, Wortelkamp S, Sickmann A, Zahedi RP . Systematic and quantitative comparison of digest efficiency and specificity reveals the impact of trypsin quality on MS-based proteomics. J Proteomics 2012; 75: 1454â1462.
Olsen JV, de Godoy LM, Li G, Macek B, Mortensen P, Pesch R et al. Parts per Million Mass Accuracy on an Orbitrap Mass Spectrometer via Lock Mass Injection into a C-trap. Mol Cell Proteomics 2005; 4: 2010â2021.
Vaudel M, Burkhart JM, Zahedi RP, Oveland E, Berven FS, Sickmann A et al. PeptideShaker enables reanalysis of MS-derived proteomics data sets. Nat Biotechnol 2015; 33: 22â24.
Giusti L, Mantua V, Da Valle Y, Ciregia F, Ventroni T, Orsolini G et al. Search for peripheral biomarkers in patients affected by acutely psychotic bipolar disorder: a proteomic approach. Mol Biosyst 2014; 10: 1246â1254.
Lima DP, Diniz DG, Moimaz SAS, Sumida DH, Okamoto AC . Saliva: reflection of the body. Intl J Infect Dis 2010; 14: e184âe188.
Javaid MA, Ahmed AS, Durand R, Tran SD . Saliva as a diagnostic tool for oral and systemic diseases. J Oral Biol Craniofac Res 2016; 6: 66â75.
Mikirova N, Casciari J, Hunninghake R . The assessment of the energy metabolism in patients with chronic fatigue syndrome by serum fluorescence emission. Altern Ther Health Med 2012; 18: 36â40.
Castro-Marrero J, Cordero MD, Segundo MJ, Sáez-Francà s N, Calvo N, Román-Malo L et al. Does oral coenzyme Q10 plus NADH supplementation improve fatigue and biochemical parameters in chronic fatigue syndrome? Antioxid Redox Signal 2015; 22: 679â685.
Castro-Marrero J, Sáez-Francà s N, Segundo MJ, Calvo N, Faro M, Aliste L et al. Effect of coenzyme Q10 plus nicotinamide adenine dinucleotide supplementation on maximum heart rate after exercise testing in chronic fatigue syndrome - a randomized, controlled, double-blind trial. Clin Nutr 2015; 35: 826â834.
Vermeulen RC, Kurk RM, Visser FC, Sluiter W, Scholte HR . Patients with chronic fatigue syndrome performed worse than controls in a controlled repeated exercise study despite a normal oxidative phosphorylation capacity. J Transl Med 2010; 8: 93.
Hsiao CP, Wang D, Kaushal A, Saligan L . Mitochondria-related gene expression changes are associated with fatigue in patients with nonmetastatic prostate cancer receiving external beam radiation therapy. Cancer Nurs 2013; 36: 189â197.
Bubber P, Hartounian V, Gibson GE, Blass JP . Abnormalities in the tricarboxylic acid (tca) cycle in brain of schizophrenia patients. Eur Neuropsychopharmacol 2011; 21: 254â260.
Vasquez-Vivar J, Kalyanaraman B, Kennedy MC . Mitochondrial aconitase is a source of hydroxyl radical. An electron spin resonance investigation. J Biol Chem 2000; 275: 14064â14069.
Shao L, Martin MV, Watson SJ, Schatzberg A, Akil H, Myers RM et al. Mitochondrial involvement in psychiatric disorders. Ann Med 2008; 40: 281â295.
Callaly E, Walder K, Morris G, Maes M, Debnath M, Berk M . Mitochondrial dysfunction in the pathophysiology of bipolar disorder: effects of pharmacotherapy. Mini Rev Med Chem 2015; 15: 355â365.
Brenu EW, van Driel ML, Staines DR, Ashton KJ, Ramos SB, Keane J et al. Immunological abnormalities as potential biomarkers in chronic fatigue syndrome/myalgic encephalomyelitis. J Transl Med 2011; 9: 81.
Klimas NG, Broderick G, Fletcher MA . Biomarkers for chronic fatigue. Brain Behav Immun 2012; 26: 1202â1210.
Acknowledgements
This work has been supported by the Italian âFondazione ItPAâ onlus with the âPremio di Mobilità ItPA, 2014â given to FC. LK, RPZ and AS gratefully acknowledge the financial support by the Ministerium für Innovation, Wissenschaft und Forschung des Landes Nordrhein-Westfalen, the Senatsverwaltung für Wirtschaft, Technologie und Forschung des Landes Berlin and the Bundesministerium für Bildung und Forschung.
Author information
Authors and Affiliations
Corresponding author
Ethics declarations
Competing interests
The authors declare no conflict of interest.
Additional information
Supplementary Information accompanies the paper on the Translational Psychiatry website
Rights and permissions
This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articleâs Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/
About this article
Cite this article
Ciregia, F., Kollipara, L., Giusti, L. et al. Bottom-up proteomics suggests an association between differential expression of mitochondrial proteins and chronic fatigue syndrome. Transl Psychiatry 6, e904 (2016). https://doi.org/10.1038/tp.2016.184
Received:
Revised:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/tp.2016.184