BMC PhysiolBMC Physiology1472-6793BioMed CentralLondon1579042210798851472-6793-5-510.1186/1472-6793-5-5Research ArticleExercise responsive genes measured in peripheral blood of women with Chronic Fatigue Syndrome and matched control subjectsWhistlerToni1taw6@cdc.govJonesJames F1jaj9@cdc.govUngerElizabeth R1eru0@cdc.govVernonSuzanne D1sdv2@cdc.govViral Exanthems and Herpesvirus Branch, Centers for Disease Control and Prevention, Atlanta, GA 30333, USA2005243200555513920042432005Copyright © 2005 Whistler et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Background

Chronic fatigue syndrome (CFS) is defined by debilitating fatigue that is exacerbated by physical or mental exertion. To search for markers of CFS-associated post-exertional fatigue, we measured peripheral blood gene expression profiles of women with CFS and matched controls before and after exercise challenge.

Results

Women with CFS and healthy, age-matched, sedentary controls were exercised on a stationary bicycle at 70% of their predicted maximum workload. Blood was obtained before and after the challenge, total RNA was extracted from mononuclear cells, and signal intensity of the labeled cDNA hybridized to a 3800-gene oligonucleotide microarray was measured. We identified differences in gene expression among and between subject groups before and after exercise challenge and evaluated differences in terms of Gene Ontology categories.

Exercise-responsive genes differed between CFS patients and controls. These were in genes classified in chromatin and nucleosome assembly, cytoplasmic vesicles, membrane transport, and G protein-coupled receptor ontologies. Differences in ion transport and ion channel activity were evident at baseline and were exaggerated after exercise, as evidenced by greater numbers of differentially expressed genes in these molecular functions.

Conclusion

These results highlight the potential use of an exercise challenge combined with microarray gene expression analysis in identifying gene ontologies associated with CFS.

Background

In a state of health, physical exercise has a quantifiable effect on neuroendocrine, autonomic, and immune systems influencing metabolic and immune responses. However, in the initial phase of acute illness, there is an avoidance of physical stressors so energy can be dedicated to healing and a return to homeostasis. While physiologic disturbance in acute illness is transient, chronic illnesses, such as chronic fatigue syndrome (CFS), have prolonged disturbances that have a debilitating effect both physiologically and psychologically. Consequently, activities that are physiologic stressors, such as physical exercise, exacerbate the symptoms that define CFS.

CFS is a complex, multifactorial illness whose etiology and pathophysiology remain unclear [1]. CFS is defined by a characteristic symptom complex in the absence of other medical or psychiatric conditions with similar clinical characteristics [2,3]. Subtle differences in hypothalamic-pituitary-adrenal axis function [4], immune system function [5], and psychological profiles [6] between CFS patients and controls have been reported; however, no consistent distinguishing difference or frank abnormality has been confirmed [7,8], and it remains unclear whether CFS represents a unique disease or a common illness response to a variety of insults.

Perhaps the greatest methodological problem with studying CFS is that many individuals identified in population studies have been sick for at least 5 years [9]. During this time, the illness waxes and wanes, making it difficult to identify biomarkers or define pathogenesis. Physical, mental, and emotional stress exacerbate CFS and result in case-defining post-exertional fatigue [2] with measurable physiologic differences [10]. Therefore, exercise challenge of people with CFS is an effective method for calibrating CFS subjects and thus increasing the likelihood of uniformly identifying biomarkers and/or physiologic abnormalities.

We used gene expression profiling of peripheral blood to evaluate differences between CFS subjects and sedentary healthy controls both before and following an exercise challenge. Overall, we found the gene expression profiles to be quite similar, and of importance, most differences were present prior to exercise challenge. These differences were in G protein-coupled receptor and ion transport and ion channel activity ontologies. The latter was exaggerated after exercise as evidenced by differential expression of a greater number of genes involved in these molecular functions. Differences were also evident in exercise response, including chromatin and nucleosome assembly, cytoplasmic vesicles, membrane transport and G-protein coupled receptor ontologies. These differences may help explain the symptoms of CFS.

Results

Exercise response genes were evaluated using a random variance t test in a paired, class comparison analysis of control subjects before and after exercise, and 21 genes were identified as being differentially expressed (Table 2). The probability of identifying these 21 genes by chance if there are no real differences between the classes was 0.056 as determined by the multivariate permutation test. Among the 21 genes, 16 could be categorized in the Gene Ontology (GO) of biological process and 15 in molecular function (results not shown). The most significant categories or "themes" of these exercise-responsive genes as assessed by an EASE score of <0.10, pertained to the biological process of transport (both vesicle-mediated and protein transport). 5 of the 21 genes were involved in this process.

Since these 21 genes reflect a healthy subject's peripheral blood gene expression response to exercise challenge, we reasoned that the expression of these would be altered in CFS subjects. To have a visual representation of these differences, the gene list from Table 2 (differentially expressed genes in control subjects, compared before and after exercise challenge) was used in a two-way hierarchical cluster analysis (Figure 1). The response of 10 of the 21 genes was quite similar in terms of magnitude and direction for both CFS and control subjects (Figure 1, marked in blue). For the other 11 genes, the magnitude of the exercise change was considerably smaller in CFS subjects (Figure 1, subject clusters 2 and 4) than in control subjects (clusters 1 and 3). With regard to the GO categories of these 21 genes, 10 genes were associated with binding and 8 with metabolism, all of which were equally distributed between the two response types (denoted by # and in Figure 1). However, 5 genes classified in vesicle-mediated and protein-transport ontologies differed between CFS and control subjects (Figure 1, highlighted in yellow).

No differentially expressed genes were identified by class comparison analysis (at a significance level of p > 0.005) for CFS subjects before and after exercise, for CFS subjects before exercise compared with controls before exercise, or for CFS subjects after exercise compared with controls after exercise.

Because differentially expressed genes were identified by class comparison prior to the exercise challenge, we reasoned that a comparison of gene expression in CFS and control subjects, using genes categorized by ontology, would more efficiently reveal perturbed physiological pathways. Figures 2 and 3 show the results of these analyses. GO terms with defined parent-child relationships are grouped together in the figures and are color-matched. Only two exercise-responsive GO categories [phospholipid binding (orange) and chromatin architecture (pale green)] were common to controls and CFS subjects (Figure 2a and 2b). However, for the chromatin architecture category, the CFS comparison highlighted 7 overlapping ontologies (containing 59 unique genes), compared with 1 ontology of 33 genes in the control comparison. The 33 genes overlap with the 59 identified in the CFS comparison. The phospholipids-binding ontologies were identical in gene composition. Exercise-related changes that were identified as significant only in control subjects were associated with genes involved with vesicle (yellow, Figure 2a), dehydrogenase (grey, Figure 2a), ATPase (pink, Figure 2a), and transporter (blue, Figure 2a), activities. Exercise-related changes that were seen only in CFS subjects were related to G-protein-coupled receptor signaling (purple, Figure 2b).

Gene ontology comparison was also used to evaluate differences between control and CFS subjects before (i.e. baseline, Figure 3a) and after (Figure 3b) exercise. Baseline differences between CFS subjects and controls that continued after exercise involved GO terms relating to ion transport (blue, Figure 3a). After exercise, these differences appear to be amplified, as evidenced by increased numbers of genes present in these GO categories and also by inclusion of more GO terms pertaining to ATPase transmembrane movement of ions (pink, Figure 3b). G-protein-coupled receptor binding (purple, Figure 3a), part of the broad functional category of signal transduction, differed between CFS subjects and controls prior to exercise. This baseline difference between controls and CFS subjects was not significant after exercise. Interestingly, complement activation (dark green, Figure 3b) was one of the exercise-induced differences between subjects and controls that was present only after challenge. Genes in most of the ontologies identified as different between CFS and control subjects had lower expression levels in CFS subjects.

Discussion

Gene expression profiling affords a unique opportunity to characterize CFS at a systems biology level. Changes in gene expression underlie many biologic processes and may provide insight into disease-specific gene expression and the response of genes to environmental stimuli. In a proof-of-concept study, we found that CFS patients had different blood mononuclear cell gene expression patterns than non-fatigued controls [11] and that CFS is a heterogeneous illness as evidenced by different gene expression profiles for patients reporting gradual onset of their illness compared with those reporting sudden onset of illness [12]. In addition, differential display polymerase chain reaction on a small number of CFS and control subjects identified candidate biomarkers in the peripheral blood [13,14].

CFS is defined by a post-exertional fatigue that does not subside 24 hours following physical stress. In contrast, exercise in healthy, untrained people induces changes in cellular homeostasis in 1 to 4 hours and a return to basal levels within 24 hours, as measured in muscle [15]. Analysis of peripheral blood gene expression in the healthy control subjects confirmed this observation since the majority of gene expression levels were the same before and 24 hours following exercise challenge. This implied that expression either returned to basal levels or was unchanged as a result of the exercise challenge. And indeed, many of the 21 exercise-induced, differentially expressed genes in control subjects were characterized by GOs that reflect a diverse set of molecular functions necessary for cell function and viability. (These ontologies overlapped with those identified in the GO comparison analysis given in Figure 2a). Figure 1 clearly illustrates the reciprocal pattern of gene expression in the 21 genes for most of the control subjects. In contrast, 11 of the genes were unchanged in CFS subjects before and after exercise; with 5 being classified in a transport-related ontology. Because this difference in gene expression is so dramatic, it implicates a fundamental perturbation in the biochemical activity of lymphocyte and monocyte peripheral blood fractions from CFS subjects compared with control subjects that does not affect classical immunologic markers (i.e, CD45) that have been shown to be unaffected in CFS patients [16,17]. Rather, low expression of these genes may have subtle effects on immune function. Immune dysfunction has been inconsistently implicated in CFS pathogenesis [18].

Class comparison was used to identify these 21 differentially expressed genes, which indicated the possible disturbance of biologic pathways (Figure 1). To explore this possibility, we used the GO comparison that is based on the knowledge that gene expression levels are dependent variables in biological processes, cellular components, and molecular functions. In this way, multiple genes in the same category reinforce each other and enhance the power for identifying the significance of the category. The GO categories considered significantly different (p < 0.005) when comparing CFS subjects with controls after exercise challenge were those pertaining to ion transporter activity (a total of 87 genes applied to this category in the comparison of CFS and controls after exercise) and ATPase activity coupled to transmembrane movement (42 genes). When the CFS and control classes are compared prior to exercise, ion transport activity and voltage-gated, ion channel activity are identified (38 and 44 genes within the GO categories, respectively).

It is evident that ion transport and ion channel activity segregate cases from controls and that exercise seems to intensify these differences. Several other conditions have been reported in which fluctuating fatigue occurs that are known to be caused by abnormal ion channels. These conditions include genetically determined channelopathies and acquired conditions such as neuromyotonia, myasthenic syndromes, multiple sclerosis, and polyneuropathies [19,20]. There are other transmembrane functions associated with differences between controls and CFS patients, including signal transducer activity through receptor binding/activity (Figure 3a). Signal transduction of transmembrane receptors occurs by a number of mechanisms, including structural changes, ion channels, and changes of transmembrane potentials. The G-protein-coupled receptors play an important role in the membrane trafficking machinery [21]. The most obvious exercise-induced changes in CFS cases pertain to gene regulation at the point of chromatin structure; whether these changes reflect the differences seen in the mRNA transcripts relating to membrane trafficking differences between cases and controls has not yet been determined.

One interesting correlate of this study was the finding that the complement pathway showed significant differences between CFS and control subjects after exercise. This has been reported previously in the analysis of these same exercise challenge-derived specimens. Sorensen et al. [22] measured levels of complement split products in the sera of these subjects and found differences between CFS and control subjects in C4a after exercise challenge. Complement activation was identified as an ontology that was significantly different between CFS and control subjects after exercise. The correlates on the data are interesting as their study measured protein levels (i.e., gene product levels) and this study measured the transcript levels.

The class comparison analysis performed in this study accounted for multiple testing and the over fitting problems of microarray data analysis. The lack of statistical significance in the 3 other class comparison analyses performed (CFS cases compared before and after exercise, comparison of cases to controls at baseline, and the comparison of cases to controls 24 hours after exercise) reflects low experimental sensitivity, most likely due to a small number of subjects, rather than an absence of biological effect. This is accounted for in the gene ontology comparison tool where classes are compared by GO category rather than with regard to individual genes.

The next line of research will detail larger numbers of subjects in the expression arrays. The emphasis in such studies will be on developing a gene expression-based multivariate function, or predictor, that accurately predicts the class membership of a new sample on the basis of the expression levels of key genes. Class discovery tools will also be applied to CFS subjects' expression profiles in an attempt to further describe discrete subsets of this disease on the basis of gene expression as we have done for gradual and sudden onset of illness [12]. However, the methods used in this study will be applied to these data sets too, as these analytical tools will prove to be very helpful in defining the pathophysiology of CFS. It is hoped that this broader, more fully encompassing approach to CFS research will open many doors to the understanding of this syndrome and perhaps of fatigue and un-wellness in general.

MethodsStudy subjects

This study adhered to human experimentation guidelines of the U.S. Department of Health and Human Services and the Helsinki Declaration. All participants were volunteers who gave informed consent. Study parameters have been described [22]. In brief, women who attended a CFS outpatient fatigue clinic volunteered for the study. All of these patients met the current working definition for CFS [2] that includes the following criteria: fatigue lasting longer than 6 months or more, no other illness that could explain the fatigue, 4–8 concurrent symptoms, fatigue that is not relieved by rest, and fatigue that interferes with occupational, educational and social or personal activities.

Healthy controls were recruited through advertisements. They were similar in age, sex, and activity level (sedentary to moderately active) to the CFS patients. All women, were scheduled for exercise challenge 5 to 10 days after the first day of their menstrual cycle. All subjects were asked not to use inhaled or systemic corticosteroids, anti-histamines, or anti-inflammatory medication for 7 days prior to the exercise challenge. Subjects performed a submaximal (70% predicted maximum work load), steady-state exercise for 20 minutes on a stationary bicycle ergometer. All subjects met these challenge criteria. Blood samples were obtained immediately before and 24 hours following exercise. The selection of subjects for inclusion into the gene expression pilot study, focused on those without allergies, to whom an exercise challenge was given: 5 women with CFS and 5 female controls.

RNA isolation

Immediately following blood collection, peripheral blood mononuclear cells were isolated, using Ficoll gradients and stored in RNAqueous™ lysis buffer (Ambion Inc., Austin, TX) at -70°C. Total RNA was extracted using the RNAqueous™ kit (Ambion), with quality and quantity assessed by agarose gel electrophoresis as described previously [23].

Preparation and hybridization of labeled cDNA

Biotinylated cDNA synthesis from 1 μg of total RNA, microarray hybridization, and detection were performed as previously described [23]. Atlas™ Human 3.8I oligonucleotide glass microarrays (Clontech Laboratories, Palo Alto, CA) were used. The slides were archived and images captured using the GSD-501™ scanner (Invitrogen Life Technologies, Carlsbad, CA).

Data analysisPreprocessing data

The scanned TIFF images were processed using ArrayVision™ (Imaging Research Inc., Ontario, Canada) as previously described [23]. A median background value was calculated around each of the 3757 features and subtracted from the mean feature signal to give the net signal for the respective gene. Values were log2-transformed, and a global normalization was used to median center the log intensity values. Genes whose expression differed by at least 1.5 fold from the median in at least 20% of arrays were retained (3699 genes), thus excluding genes showing minimal variation across the set of arrays. If an expression value was missing or filtered out in more than 50% of the arrays, that feature was not included, leaving a total of 3682 genes for analysis.

Analysis

We used several analytical approaches to look for differential expression among four predefined classes: control subjects before exercise, control subjects after exercise, CFS subjects before exercise, and CFS subjects after exercise. Comparisons included paired analyses of controls before versus after exercise and CFS cases before versus after exercise. In addition, we examined pre-exercise (baseline) controls versus pre-exercise CFS subjects and post-exercise controls versus post-exercise CFS subjects.

Analytical approach 1

The more commonly used approach to microarray gene expression analysis is to establish if gene expression profiles differ between subjects in predetermined classes and to identify the genes responsible for the differences. The Class Comparison analysis method in BRB ArrayTools [24] uses this approach, applying a random variance t-test to the data. This is an improvement over the standard t-test as it permits sharing information among genes about within-class variation without assuming that all genes have the same variance [25]. It is an advisable modification when the numbers of samples per class are small. Genes were considered statistically significant if the parametric p-value was less than 0.005. Next a global, multivariate permutation test was performed to determine whether the expression profiles differed between the classes by permuting the array classification. We used this test to provide a median false discovery rate of 10%.

The EASE software package [26] was then used to evaluate the biologic significance of the ontology of genes identified as differentially expressed by class comparison. EASE performs a statistical analysis of gene categories in a gene list relative to all genes on the array and calculates a conservative variant of the standard Fisher exact probability called the EASE score. The most significant categories as assessed by EASE score are deemed "themes" of the gene list. These themes correspond to the systematic and standardized nomenclature developed by the GO (Gene Ontology) Consortium [27]. The three organizing principles of GO are molecular function, biological process and cellular component, and presently 22,665 human gene products have been annotated. Associating genes with related GO terms assists in the interpretation of expression patterns. A GO category includes genes described by that term and those included in any subset (or children) of that GO term. A gene may be categorized in more than one ontology category.

To demonstrate the distinct gene clusters, we performed a two-way hierarchical cluster analysis on the genes identified by this analytical approach, using the algorithm described by Eisen et al. [28].

Analytical approach 2

We used the Gene Ontology comparison module of BRB Array tools [24] as a second analytical approach. In this method, the genes are assigned GO terms prior to analysis. For each GO term the total number of genes on the array belonging to that category was determined. A random variance t-test or F-test [25] was used to determine the p-value for differences between the predefined classes for each gene in the GO term. Two statistics are computed to summarize the p-values for all genes in the GO term: the Fisher statistic and the Kolmogorov-Smirnov statistic [24]. These p-values provide a list of GO categories that have more genes differentially expressed among the classes than expected by chance. We considered a GO category differentially regulated if either significance level was less than 0.005. All GO categories with between 5 and 100 genes represented on the array were considered.

Competing interests

The author(s) declare that they have no competing interests.

Authors' contributions

TW contributed to the design of the experimental approach, the analysis of the data and drafted the manuscript. JFJ contributed to the conception and design of this study, conducted the exercise challenges and assisted in drafting the manuscript. ERU contributed to the design of the experimental approach, and to the manuscript preparation. SDV was instrumental in the design of the experimental approach, analysis, presentation, discussion of these data and manuscript preparation. All authors read and approved the final manuscript.

Acknowledgements

We thank Kelly Bowen for her technical assistance. Bristol Sorensen for co-coordinating the exercise challenge data included in this manuscript, Eric Aslakson for help with the diagrams and Dr. William C. Reeves for help preparing the manuscript. Analyses were performed using BRB ArrayTools developed by Dr. Richard Simon and Amy Peng Lam.

EvengardBSchacterleRSKomaroffALChronic fatigue syndrome: new insights and old ignoranceJ Intern Med19992464554691058371510.1046/j.1365-2796.1999.00513.xFukudaKStrausSEHickieISharpeMCDobbinsJGKomaroffAThe chronic fatigue syndrome: a comprehensive approach to its definition and study. International Chronic Fatigue Syndrome Study GroupAnn Intern Med19941219539597978722ReevesWCLloydAVernonSDKlimasNJasonLABleijenbergGIdentification of ambiguities in the 1994 chronic fatigue syndrome research case definition and recommendations for resolutionBMC Health Serv Res20033251470220210.1186/1472-6963-3-25PapanicolaouDAAmsterdamJDLevineSMcCannSMMooreRCNewbrandCHNeuroendocrine aspects of chronic fatigue syndromeNeuroimmunomodulation20041165741475805210.1159/000075315LyallMPeakmanMWesselySA systematic review and critical evaluation of the immunology of chronic fatigue syndromeJ Psychosom Res20035579901293250510.1016/S0022-3999(02)00515-9Moss-MorrisRPetrieKJDiscriminating between chronic fatigue syndrome and depression: a cognitive analysisPsychol Med20013146947911305855CleareAJThe neuroendocrinology of chronic fatigue syndromeEndocr Rev20032423625212700181AfariNBuchwaldDChronic fatigue syndrome: a reviewAm J Psychiatry20031602212361256256510.1176/appi.ajp.160.2.221NisenbaumRJonesAJonesJReevesWLongitudinal analysis of symptoms reported by patients with chronic fatigue syndromeAnn Epidemiol2000104581101836810.1016/S1047-2797(00)00119-8GaabJHusterDPeisenREngertVHeitzVSchadTHypothalamic-pituitary-adrenal axis reactivity in chronic fatigue syndrome and health under psychological, physiological, and pharmacological stimulationPsychosom Med2002649519621246120010.1097/01.PSY.0000038937.67401.61VernonSDUngerERDimulescuIMRajeevanMReevesWCUtility of the blood for gene expression profiling and biomarker discovery in chronic fatigue syndromeDis Markers20021819319912590173WhistlerTUngerERNisenbaumRVernonSDIntegration of gene expression, clinical, and epidemiologic data to characterize Chronic Fatigue SyndromeJ Transl Med20031101464193910.1186/1479-5876-1-10PowellRRenJLewithGBarclayWHolgateSAlmondJIdentification of novel expressed sequences, up-regulated in the leucocytes of chronic fatigue syndrome patientsClin Exp Allergy200333145014561451915410.1046/j.1365-2222.2003.01745.xSteinauMUngerERJonesJFRajeevanMSDifferential dispaly PCR of peripheral blood for biomarker discovery in chronic fatigue syndromeMolecular Medicine20048275075510.1007/s00109-004-0586-4PilegaardHOrdwayGASaltinBNeuferPDTranscriptional regulation of gene expression in human skeletal muscle during recovery from exerciseAm J Physiol Endocrinol Metab2000279E806E81411001762StrausSEFritzSDaleJKGouldBStroberWLymphocyte phenotype and function in the chronic fatigue syndromeJ Clin Immunol1993133040809527010.1007/BF00920633PeakmanMDealeAFieldRMahalingamMWesselySClinical improvement in chronic fatigue syndrome is not associated with lymphocyte subsets of function or activationClin Immunol Immunopathol1997828391900004610.1006/clin.1996.4284LyallMPeakmanMWesselySA systematic review and critical evaluation of the immunology of chronic fatigue syndromeJ Psychosom Res20035579901293250510.1016/S0022-3999(02)00515-9ChaudhuriABehanPOFatigue in neurological disordersLancet20042097898810.1016/S0140-6736(04)15794-2ChaudhuriAWatsonWSPearnJBehanPOThe symptoms of chronic fatigue syndrome are related to abnormal ion channel functionMed Hypotheses20005459631079072510.1054/mehy.1998.0822MunroSOrganelle identity and the organization of membrane trafficNat Cell Biol200464694721517045310.1038/ncb0604-469SorensenBStreibJEStrandMMakeBGiclasPCFleshnerMComplement activation in a model of chronic fatigue syndromeJ Allergy Clin Immunol20031123974031289774810.1067/mai.2003.1615OjaniemiHEvengardBLeeDUngerEVernonSImpact of RNA extraction from limited samples on microarray resultsBiotechniques20033596897314628670BRB ArraytoolsWrightGWSimonRA random variance model for detection of differential gene expression in small microarray experimentsBioinformatics200319244824551466823010.1093/bioinformatics/btg345EASE software packageAshburnerMBallCABlakeJABotsteinDButlerHCherryJMGene ontology: tool for the unification of biology. The Gene Ontology ConsortiumNat Genet20002525291080265110.1038/75556EisenMBSpellmanPTBrownPOBotsteinDCluster analysis and display of genome-wide expression patternsProc Natl Acad Sci U S A1998951486314868984398110.1073/pnas.95.25.14863Figures and Tables

Hierarchical clustering of exercise responsive genes in control subjects. The 21 differentially expressed genes identified by a class comparison test of control subjects (before and after exercise challenge (Table 2)) were clustered using a two-way hierarchical algorithm. In the matrix each row represents the hybridization results for a single gene, and each column represents a subject. Transcript levels are depicted as above (red) or below (green) the mean. The dendograms illustrates average-linkage hierarchical clustering of subjects (top) and genes (left). Refseq IDs for each gene is given on the right of the matrix. Those with similar exercise responses in both CFS and control subjects are at the top of the matrix, and the remainder of genes (highlighted in blue) show a diminished exercise response in CFS cases. Refseq IDs highlighted in yellow classify to the GO categories of protein or vesicle-mediated transport (Biological Process). Refseq IDs followed by: are classified to the GO category of binding (Molecular Function); and/or # are classified in the GO category of metabolism (Biological Process). The subjects group into 4 clusters which approximate to: 1) Control subjects before exercise (Con0); 2) CFS cases before exercise (CFS0); 3) Control subjects after exercise (Con24); and 4) CFS cases after exercise (CFS24).

Significant gene ontology categories defining exercise-related changes in control (a) and CFS subjects (b). The three organizing principles of GO (represented as grey shaded boxes) are molecular function, biological process, and cellular component. Related ontologies and/or subgroups of the ontologies are denoted by similarly colored squares in all tables. Ontologies presented in these figures were significant at a p-value <0.005 by one or both of the LS and KS permutation tests.

Gene Ontology categories identified as significantly different between controls and CFS cases at baseline (pre-exercise) (a) and post-exercise (b). The three organizing principles of GO (represented by grey shaded squares) are molecular function, biological process, and cellular component. Related ontologies and/or subgroups of the ontologies are denoted by similarly colored squares in all tables. Ontologies presented in these figures were significant at a p-value <0.005 by one or both of the LS and KS permutation tests.

Demographic and challenge related information of subjects included in this study.

Exercise Challenge InformationSubjects

CFSControls
Mean Age41.836.1
Age Range35–4827–45
Average workload performed during exercise challenge in watts.13601664
Mean change in score of daily symptom diary (± SE) calculated for the exercise challenge.#1.61 ± 0.65*0.54 ± 0.31

#A 10 symptom daily diary with a 0 to 4 scale was filled out each day for 2 weeks before and 1 week after challenge [22].

*Significant relative to control patients at P < 0.05

List of genes differentially expressed in exercised control subjects. The parametric p-value is a measure of the significance of the random variance t-statistic test used to identify differentially expressed genes in the class comparisons.

GenBank IDDescriptionParametric p-value
NM_003431Zinc finger protein 124 (HZF-16)0.0001
NM_003466Paired box gene 80.0001
NM_003477Pyruvate dehydrogenase complex, lipoyl-containing component X; E3-binding protein0.0002
NM_003591Cullin 20.0003
NM_003985Tyrosine kinase, non-receptor, 10.0004
NM_003490Synapsin III0.0004
NM_003764Syntaxin 110.0004
NM_003494Dysferlin, limb girdle muscular dystrophy 2B (autosomal recessive)0.0006
NM_003715Vesicle docking protein p1150.0006
NM_003692Transmembrane protein with EGF-like and two follistatin-like domains 10.0006
NM_003853Interleukin 18 receptor accessory protein0.0006
NM_003716Ca2+-dependent activator protein for secretion0.0006
NM_003558Phosphatidylinositol-4-phosphate 5-kinase, I beta0.0008
NM_003854Interleukin 1 receptor-like 20.0009
NM_003543H4 histone family, member H0.0010
NM_003488A kinase (PRKA) anchor protein 10.0013
NM_003487TATA box binding protein (TBP)-associated factor, RNA polymerase II, N, 68 kD0.0014
NM_003528H2B histone family, member Q0.0019
NM_003693Acetyl LDL receptor; scavenger receptor expressed by endothelial cells0.0021
NM_003473Signal transducing adaptor molecule (SH3 domain and ITAM motif) 10.0029
NM_004653SMC (mouse) homolog, Y chromosome0.0035