Clinical history of two subjects with ME/CFS displays significant muscle and cognitive impairments. In our current study, we included a 71 years old white male ME/CFS patient. Sometimes in 1988, when the subject was a 40 years old, he was initially diagnosed with acute flu-like symptom. The disease eventually progressed with severe symptomology of sore throat, acute joint pain, swollen cervical lymph nodes, mental fog, memory impairment, extreme dizziness, and severe fatigue leaving the patient bedridden for 90% of the time. Although, the underlying cause behind these symptoms was not determined, but subsequent antigen testing detected positive titers for EBV, HSV1, and CMV infection. In 1991, the patient was admitted to Dr. Peterson’s clinic in a cluster of patients with similar clinical symptoms. A single-photon emission computerized tomography (SPECT) scan imaging analyses revealed that there was chronic Myalgic encephalomyelitis with bilateral hypoperfusion.
Another patient in our current study is currently a 67 years old female who had similar initial flu-like symptoms in 1988, which rapidly progressed to severe sinus infection with methicillin-resistant staphylococcus aureus (MRSA) infection. Eventually the patient received 6 turbinate surgeries. Although, the subject did not have any sinus problems, but has a history of depression, severe difficulty in sleeping (staying asleep specifically) and experiences fatigue, pain (general), IBS symptomology, severe memory impairment, and major brain fog. Recently the patient experienced extremely high blood pressure, with no noted cardiac dysfunction.
Both patients are classic cases of ME/CFS with extreme muscle fatigue and severe cognitive impairment. They routinely visit Dr. Peterson’s clinic for their treatment regimen. Blood collection was performed with their consent and based on the guideline of approved IRB protocol.
Autophagy-driven degradation of metabolically inactive protein is an important quality control process of cellular metabolism that directly regulates the function of muscle tissue [33] and also coordinates cognitive function[34]. Interestingly, both of these biological processes are significantly damaged in chronic fatigue syndrome, which intrigued us to study the impairment of autophagy in ME/CFS. Since autophagy impairment can be evaluated with the upscaled expressions of autophagy proteins in serum, next, we wanted to study the expressions of autophagy proteins in the serum of ME/CFS patients.
Estimation of protein aggregation propensities in serum samples of ME/CFS patients. Severe impairment in autophagy might result the increased tendency of protein aggregation in serum[35–37]. Therefore, first, we performed a thioflavin T (ThT)-based protein aggregation study in patient serum. The analysis is a fluorimetric tracking analysis that generates a sigmoidal evolution of protein from its liner structure to aggregated morphology. This type of kinetic study was previously done[38] in a cell-free system to analyze the kinetics of amyloid-β protein aggregation[39]. In our present study, we adopted similar fluorimetric tracking study to explore if the serum proteins of ME/CFS patients have higher propensities to form aggregates. The accelerated rate of aggregation does not necessarily reflect the impairment of physiological mechanism of protein degradation; however, it might indicate that there is a higher tendency to form aberrant protein aggregates in serum samples of ME/CFS patients. Accordingly, serum samples of both 74 years old healthy control (Fig. 1A) and a 71 years old male ME/CFS (Fig. 1B) patient displayed increasing but non-liner pattern of aggregation as indicated with nonlinear fitted curves of Florescence intensities (486 nm/450 nm) with the function of increasing time from 0 to 90 minutes of ThT addition. The plotting equation is \(\text{y}={a}_{1}{x}^{2}+{a}_{2}x+{a}_{3}\) with goodness of fit was determined as R2 = 0.9780 (~ 1 for liner) for healthy control and 0.9570 for patient. Interestingly, we observed that serum proteins of patient displayed significantly faster rate of protein aggregation compared to healthy control as indicated with steeper slope of 70.70 ± 1.956 in patient compared to relatively flat slope of 17.72 ± 0.3806 in control (Fig. 1E). Moreover, the significance analysis revealed that there is a strong difference of protein aggregation rate in ME/CFS patient [F29,29 =15.99 (***p < 0.0001)] compared to its gender-matched healthy control (Supplementary Table 1). Next, we performed a similar comparison study between a 68 years old healthy (Fig. 1C) female and a 67 years old female ME/CFS patient (Fig. 1D) that resulted a similar outcome of protein aggregation upon ThT addition. In that case-control study, we also observed a significantly faster aggregation of serum proteins in female patient compared to healthy control (Supplementary Table 1) [F29,29 =15.99 (***p < 0.0001)]. The non-liner fitted curve generated a steeper slope of 59.68 ± 3.005 in patient versus 19.36 ± 1.599 in healthy control (Fig. 1F). Taken together, our ThT assay revealed that serum samples of two ME/CFS patients have higher and faster propensity of protein aggregation compared to their age-matched and gender-balanced control.
Upregulations of autophagy markers in the serum of ME/CFS patients. Upregulations of autophagy markers in serum have been documented in many neurodegenerative and metabolic disorders. To evaluate these markers in ME/CFS, we performed a human autophagy antibody array analysis in the serum of two case-control subjects. Total protein concentration was assessed by Bradford procedure followed by loading of equal amount of proteins (1.4 mg per sample) in array blots encoded with 20 autophagy-related proteins (Supplementary Fig. 1C). After that, the signal was detected with IRDye800-conjugated streptavidin and the image was captured in Li-Cor Odyssey Sa imager as discussed in method section. Surprisingly, we observed that expressions of ATG5, ATG7, ATG13, p62, Rheb, and α-syn were strongly upregulated in serum sample of 71-year-old ME/CFS patient (Fig. 2B; Supplementary Fig. 1B) compared to age-matched healthy subject (Fig. 2A; Supplementary Fig. 1A). The result was further confirmed with a densitometric analyses (Fig. 2C, 2D, 2E, 2F, 2G, and 2H are ATG 5, 7, 13, p62, Rheb, and α-syn respectively) with the significance of means calculated based on three independent experiments at ***p < 0.0001. To further confirm the result, we performed ELISA analyses of ATG5, ATG13, p62 and α-syn in serum sample of male ME/CFS patient with commercially available kits as described in method section. The detection sensitivity of these ELISA kits varies from lot to lot and also based on serum concentrations of proteins. Therefore, to demonstrate the accurate detection of these proteins, we adopted a dilution series of serum samples for ELISA-based detection of ATG5, ATG13, p62, and α-syn. Accordingly, a dose responsive non-liner fitting curve analysis followed by measuring X intercept demonstrate that 1:4 dilution of serum provides most efficient detection of these serum-derived factors (Supplementary Fig. 2A-C). Subsequent ELISA analyses revealed that ATG5 (Fig. 2I), ATG13 (Fig. 2J), p62 (Fig. 2K), and α-syn (Fig. 2L) were strongly upregulated in serum samples of 71-year-old ME/CFS patients with most significant elevation in ATG13 level as evaluated with an unpaired t-test (***p < 0.0005 = 0.0003 vs. healthy after 5 independent experiments). Although, ELISA is most reliable technology to quantify serum proteins ,sometimes low signal-to-noise ratio of HRP detection method limits the accurate quantification[40]. In our manuscript, we introduced a novel near infrared-based ELISA method, which produces extremely sensitive detection of low abundance protein in serum mainly due to the high signal-to-noise ratio of near-infrared probe[41]. Upon binding to biotinylated secondary antibody, IRDye-800 conjugated streptavidin produces extremely sharp signal at the end point that was detected in Li-Cor Odyssey Sa plate reader. Evidently, strong green signals representing upscaled abundance of ATG5 (Fig. 2M), ATG13 (Fig. 2N), p62 (Fig. 2O), and α-syn (Fig. 2P) were observed in serum sample of 71-year-old ME/CFS male compared to low green signals indicative of low levels of these proteins in age-matched control. The specificity of the signal was validated with no signal in 2◦antibody control whereas equal loading was verified with total albumin in respective serum samples. The result was further quantified with a densitometric analysis relative to total albumin concentrations (Supplementary Fig. 3A). Collectively, these results suggest that serum sample of 71-year-old male ME/CFS patient displayed upscaled expressions of ATG5, ATG13, p62, and α-syn, with maximum difference in ATG13.
Next, we performed similar antibody array analysis in the serum sample of 67-years-old female patient, which indicated significant elevations of LC3a, LC3b, p62, α-syn, and ATG-13 proteins (Fig. 3B: Supplementary Fig. 1D) compared to age-matched healthy control (Fig. 3A; Supplementary Fig. 1C). The result was further quantified with densitometric analyses (Fig. 3C-3G). To nullify the possibility of unequal loading of protein in array, we estimated the total protein concentration with Bradford method[42] and 1.5 mg total protein was loaded on array membrane for each sample. Since antibody array is a semi-quantitative method to estimate serum protein, next we performed a quantitative ELISA analysis of LC3a, ATG5, ATG13, p62 and α-syn. Although, LC3a ELISA analyses (Fig. 3H) did not display any significance, we observed strong elevations of ATG13 (Fig. 3I), p62 (Fig. 3J), and α-syn (Fig. 3K) in 67-years-old female ME/CFS patient compared to healthy control with maximum difference in ATG13 level (***p < 0.001). Results were confirmed after five independent experiments.
To further confirm, we performed near-infrared ELISA assay for LC3a (Fig. 3L), ATG13 (Fig. 3M), p62 (Fig. 3N), and α-syn (Fig. 3O) with albumin as loading control. Interestingly, similar to our quantitative ELISA method, no difference in LC3a expression was observed (Fig. 3L; Supplementary Fig. 3B). However, expressions of ATG13, p62, and α-syn were consistently higher in patient compared to healthy control (Fig. 3M-3O: Supplementary Fig. 3B). In both case-control studies, we observed that serum levels of p62, α-syn and ATG13 were consistently elevated suggesting that these patients might have significant impairment in the formation of autophagosomes.
Evaluation of lysosomal function in ME/CFS patients. Upon enclosure in autophagosomes, cellular components are directed and fused to lysosomes for hydrolytic degradation[43]. Therefore, next, we wanted to evaluate if there is any lysosomal impairment in these patients. Lysosome is the home of cathepsin, a class of acid proteases that plays an essential role in the degradation of autophagic material and maintaining the cellular homeostasis of metabolism[44]. Therefore, lysosomal impairment can be evaluated by measuring the level of different cathepsins and other proteases. We adopted a human protease proteome profiler array, which can quantitatively detect 40 different blood-borne proteases (layout was shown in supplementary Fig. 4) including Cathepsins, Kallikreins, and matrix metalloproteinases. Interestingly, we did not observe any difference in the expression of lysosomal cathepsins between healthy control (Fig. 4A) and CFS male patient (Fig. 4B) as well as between female control (Fig. 4C) and female ME/CFS patient (Fig. 4D). Although, strong signals were observed for lysosomal proteases such as Cathepsin A (Fig. 4E), D (Fig. 4F), and X (Fig. 4G) in both case-control group, no significant difference was observed. We observed upregulation of cathepsins in male patient, however that difference did not pass t-test for the statistical significance. In contrast, we observed significant upregulation of extra-lysosomal proteases such as MMP-9 (Fig. 4H) and proteinase-3 (Fig. 4K) in male ME/CFS, but not female ME/CFS patient. Collectively, these results suggest that although there was significant impairment in the formation of autophagosomes, there was no abnormality in the expression of lysosomal proteases and a majority of other extra-lysosomal proteases.
Exploring the effect of autophagy marker ATG13 in the microglial stress response. Since we observed strong upregulations of ATG13 in serum samples of both patients, next we wanted to see if ATG13 plays any role in the pathogenesis of ME/CFS. Increase ROS production has been implicated in the pathogenesis of ME/CFS[45] and several other neuroinflammatory diseases such as MS[46], PD[47], and AD[48]. In order to explore the effect of ATG13 on microglial ROS production, we created an in vitro cell culture model (Fig. 5A), in which HMC3 human microglial cells were treated with serum-constituted DMEM media (1:1 v/v) followed by measuring ROS load in different time periods starting from 30 mins to 2 hrs with DCFDA staining procedure[46]. As a control, healthy serum-constituted media was added on HMC3 cells and evaluated for ROS production. Interestingly, our immunofluorescence analyses (Fig. 5B) revealed that patient, but not healthy serum, time-dependently increased the ROS production. The DCFDA-labelled ROS signal appeared as early as 30 minutes and reached maximum at 120 min. This semi-quantitative analyses was further confirmed with a real-time kinetic study (Fig. 5C), in which ROS production was measured at a ratio of 485/535 nm as a function of increasing time. Interestingly, we observed that serum sample of male patient increased microglial ROS production with increasing time, whereas healthy serum did not evoke any ROS production. Next, we performed similar fluorescence imaging (Fig. 5D) and kinetic (Fig. 5E) analyses with serum samples of female ME/CFS patient and age-matched healthy control. Accordingly, we observed that patient serum elevated ROS production with increasing time. Taken together, these results established a new cell culture model to evaluate the ROS-inducing capacity in serum samples of ME/CFS patients. With the aid of this model, we were able to quantify ROS-production in HMC3 human microglial cells upon stimulation with ME/CFS serum. Next, we wanted to study if ATG13 in ME/CFS serum enhances ROS production in HMC3 microglial cells. To study that, we adopted a neutralization technique (Fig. 5A; lower panel), in which patient sera were added on a 96 well plate, which is coated with ATG13 neutralization antibody. After 30 minutes, ATG13-depleted sera were harvested from the top and then proceeded with a similar sets of experiments to evaluate microglial ROS production. Interestingly, ATG13 depletion from the sera of both male (Fig. 5F) and female patients (Fig. 5G), significantly attenuated the production of ROS in HMC3 cells suggesting a direct role of serum-derived ATG13 in stimulating oxidative stress in microglial cells. Microglial production of nitric oxide (NO) also plays an important role in the pathogenesis of neuroinflammatory diseases[46, 49]. Next, we wanted to study if ATG13 could directly stimulate nitrite production in microglia. Surprisingly, Stimulation of HMC3 microglial cells with increasing doses (1, 2, 5, and 10 ng/mL) of ATG13 protein alone induced nitrite production (Fig. 6A). The nitrite production was compared with positive control LPS. Although, ATG13-stimulated nitrite production was almost 2-times lower than that of LPS stimulation, combination of 1 ng/mL ATG13 with 0.5 µg/mL LPS stimulated the production of nitrite on a par with 1 µg/mL LPS. Accordingly, 1:4 (25 µL serum + 75 µL serum-free DMEM/F12 media) and 1:2 (50 µL serum + 50 µL serum-free DMEM/F12 media) dilution (v/v) of serum from male ME/VFS patient dose-dependently stimulated nitrite production in HMC3 microglial cells. On the other hand, we observed that depletion of ATG13 significantly ablated the ability of the patient’s serum to stimulate nitrite production in HMC3 cells. Taken together, these results suggest that elevated ATG13 in the serum of ME/CFS patients directly stimulates oxidative stress and nitrite production in microglial cells.
Regulation of ATG13 in the serum samples of ME/CFS patients. Next, we wanted to measure ATG13 level and levels of other autophagy-related proteins such as ATG5 and p62 in the sera of 12 ME/CFS patients and 12 healthy controls. The cohort of 24 samples from Dr. Peterson’s clinic was randomly selected and displayed with their age, ethnicity, and gender in a table (Table 1). Briefly, serum samples were diluted (1:4 dilution = 25 µL serum + 75 µL sample diluent) and then assayed with competitive ELISA technique for ATG13 (Fig. 7A), ATG5 (Fig. 7B), LC3a (Fig. 7C), and p62 (Fig. 7D). Quantification was performed based on the standard curve equation as shown in supplementary Fig. 5. Interestingly, we observed that ATG13 levels are higher in 10 patients (out of 12) compared to 10 healthy controls (total 12). No difference was observed while comparing the levels of ATG5, LC3a, and p62. Taken together, our results suggest that dysregulation of autophagy resulting the release of ATG13 in the serum of ME/CFS patients stimulates the oxidative stress response and nitric oxide production in microglial cells.
Table 1
Gender, ethnicity, age, disease status, and serum ATG13 (± SD) levels in a cohort of 24 subjects (12 healthy and 12 ME/CFS). High or low was determined with reference to mean (= 4.8030). Expression above mean is considered to be high and below mean is considered to be low. Red = disease and Blue = control. OI = orthostatic intolerance; SFN = small fiber polyneuropathy; POTs = Postural Orthostatic Tachycardia.
S/N | Status | Age | Gender | Ethnicity | Disease (Y/N) | ATG13 ng/mL | ATG13 (High/low) |
1 | Control | 58 | F | W | N | 5.2438 | high |
2 | Control | 60 | F | W | N | 3.8755 | low |
3 | Control | 44 | F | W | N | 1.3909 | low |
4 | CFS+, OI+, POTS+, SFN+ | 72 | F | W | Y | 11.299 | high |
5 | SFN + only, no CFS | 83 | M | W | Y | 1.6975 | low |
6 | Control | 66 | F | W | N | -4.298 | low |
7 | CFS+, OI+, POTS+ | 40 | F | W | Y | 45.380 | Very high |
8 | Control | 25 | M | W | N | 3.933 | low |
9 | CFS+, OI+, SFN+ | 66 | F | W | Y | 7.559 | high |
10 | CFS+, OI+, SFN+ | 65 | M | W | Y | 4.999 | high |
11 | Control | 59 | F | W | N | 2.061 | low |
12 | Control | 59 | F | W | N | 22.547 | high |
13 | Control | 67 | M | W | N | -7.765 | low |
14 | Control | 55 | M | W | N | -2.682 | Low |
15 | CFS + only | 73 | M | W | Y | 15.436 | high |
16 | Control | 68 | F | W | N | -1.029 | low |
17 | CFS + only | 33 | M | W | Y | 13.490 | high |
18 | Control | 32 | M | W | N | -6.861 | low |
19 | CFS+, OI+, SFN+ | 28 | F | W | Y | 7.8578 | high |
20 | CFS+, OI+, SFN+ | 65 | F | W | Y | 16.087 | high |
21 | Control | 65 | M | W | N | -0.6706 | low |
22 | CFS+, OI+ | 56 | F | W | Y | 6.51124 | high |
23 | CFS+, OI+ | 61 | F | W | Y | -9.1545 | low |
24 | CFS + only | 72 | F | W | Y | 7.377 | high |