Figures
Abstract
Introduction
The components of minute ventilation, respiratory frequency and tidal volume, appear differentially regulated and thereby afford unique insight into the ventilatory response to exercise. However, respiratory frequency and tidal volume are infrequently reported, and have not previously been considered among military veterans with Gulf War Illness. Our purpose was to evaluate respiratory frequency and tidal volume in response to a maximal cardiopulmonary exercise test in individuals with and without Gulf War Illness.
Materials and methods
20 cases with Gulf War Illness and 14 controls participated in this study and performed maximal cardiopulmonary exercise test on a cycle ergometer. Ventilatory variables (minute ventilation, respiratory frequency and tidal volume) were obtained and normalized to peak exercise capacity. Using mixed-design analysis of variance models, with group and time as factors, we analyzed exercise ventilatory patterns for the entire sample and for 11 subjects from each group matched for race, age, sex, and height.
Results
Despite similar minute ventilation (p = 0.57, η2p = 0.01), tidal volume was greater (p = 0.02, η2p = 0.16) and respiratory frequency was lower (p = 0.004, η2p = 0.24) in Veterans with Gulf War Illness than controls. The findings for respiratory frequency remained significant in the matched subgroup (p = 0.004, η2p = 0.35).
Conclusion
In our sample, veterans with Gulf War Illness adopt a unique exercise ventilatory pattern characterized by reduced respiratory frequency, despite similar ventilation relative to controls. Although the mechanism(s) by which this pattern is achieved remains unresolved, our findings suggest that the components of ventilation should be considered when evaluating clinical conditions with unexplained exertional symptoms.
Citation: Lindheimer JB, Cook DB, Klein-Adams JC, Qian W, Hill HZ, Lange G, et al. (2019) Veterans with Gulf War Illness exhibit distinct respiratory patterns during maximal cardiopulmonary exercise. PLoS ONE 14(11): e0224833. https://doi.org/10.1371/journal.pone.0224833
Editor: Kathryn L. Weston, Teesside University/Qatar Metabolic Institute, UNITED KINGDOM
Received: May 7, 2019; Accepted: October 22, 2019; Published: November 12, 2019
This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.
Data Availability: All data were collected at the VA and the signed subject consent forms did not make provision for making individual data records publicly available, even in de-identified form. However the authors can provide the “metadata” – i.e. the numerical (aggregated data) results used to generate the figures. Requests for access can be sent to: Sharmaine Forde-Miller Privacy Officer East Orange Campus of the VA New Jersey Healthcare System 385 Tremont Ave East Orange, NJ 07018 PH: 973-676-1000 x1948 [email protected] or Joselyn McLaughlin, PhD Research Compliance Officer Directors's Office East Orange Campus of the VA New Jersey Healthcare System 385 Tremont Ave East Orange, NJ 07018 PH: 973-676-1000 x2035 [email protected].
Funding: This work was supported by Department of Veterans Affairs Office of Research & Development; Clinical Science Research & Development: 1I21CX000797 (MJF), 1K2CX001679 (JBL), www.research.va.gov. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Introduction
Increased minute ventilation (, l·min-1) secondary to physiological, pathophysiological or psychological stimuli is well-recognized. Less understood, however, are the independent changes in tidal volume (VT) and respiratory frequency (fR) that comprise
; variables that are often overlooked [1]. As reviewed by Tipton and colleagues [1], a small but growing number of studies suggest that that VT and fR are under differential regulation in response to stressors, such as exercise, with each contributing uniquely to exercise hyperpnoea. For example, Nicolò and colleagues [2] recently demonstrated that fR, unlike VT, responds rapidly to changes in workload during high-intensity cycling and recovery, is independent from metabolic factors (
and
), and is strongly associated with perceived exertion. These findings support earlier work [3–6] suggesting that VT is determined more by metabolic factors than fR. Understanding the differential control of VT and fR is of clinical relevance [1], and may provide unique insight for the interpretation of cardiopulmonary exercise testing (CPET).
A primary indication for CPET is the evaluation of undiagnosed exercise intolerance, particularly when symptoms are disproportionate with cardiopulmonary function [7]. Such a scenario has been described in military veterans with Gulf War Illness (GWI)–a chronic multisymptom illness affecting 25–32% of veterans deployed to Operations Desert Storm and Shield in 1990–1991 [8]. Although the etiology of GWI is incompletely understood, deployment related exposures (e.g., pesticides, nerve agents, oil well fires) are suspected to contribute [8]. As such, population-based studies indicate significantly greater respiratory symptoms between deployed and non-deployed Gulf War veterans [9] but similar pulmonary function [10]. Despite a clear indication for performing CPET (i.e., symptoms disproportionate to function), to our knowledge only two studies have utilized CPET for diagnostic purposes in GWI [11, 12], only one of which provided detailed CPET results [12, 13]. Numerous other studies have employed maximal or submaximal exercise testing in veterans with GWI as a stressor [13–21] to elucidate underlying mechanisms of GWI. These studies have provided considerable insight into the pathophysiology of GWI, but ventilatory patterns during CPET were beyond their purview and thus not reported. Therefore, our understanding of ventilatory patterns (, VT and fR) during exercise in GWI is nascent.
Cook and colleagues have previously found that veterans with GWI report enhanced pain sensitivity in response to exercise [16] and rate exercise as more effortful [13] relative to controls. Given that fR is strongly associated with perceived exertion and provides a valid index of physical effort in healthy adults [2–4], we hypothesized that veterans with GWI may have an exaggerated fR response to exercise that has not previously been recognized. An excessive increase in fR, and therefore , is an energetically inefficient strategy that may contribute to exercise-induced symptoms reported in this population. In this context, the present study was designed to compare maximal CPET performance between veterans with GWI and controls with respect to exercise ventilatory patterns (
, VT, and fR). Should between-group differences be detected as hypothesized, these findings may provide support for pursuing therapeutic strategies (i.e., breathing techniques) to ameliorate exertional symptoms and enhance exercise tolerance in veterans with GWI.
Materials and methods
Participants
Thirty-four individuals volunteered to participate in this study, including 20 cases of GWI (GWI+) and 14 controls (GWI-). GWI case status was determined using the Centers for Disease Control (CDC) and Kansas criteria [22]. In brief, cases must endorse moderate-to-severe symptoms in ≥ 3 domain areas (i.e., fatigue, pain, neurological/cognitive/mood, skin, gastrointestinal and respiratory) that began after 1990 and persisted for ≥ 1 year. Comorbid conditions (e.g., diabetes, heart disease, stroke, etc.) that may account for chronic symptoms were excluded per case definition [23], and were also excluded for control participants. In addition to meeting CDC and Kansas criteria, GWI+ cases were also required to be clinically fatigued as determined by the Fatigue Severity Scale [24]. Controls (GWI-) included both deployed veterans (n = 4) and civilians (n = 10). All participants completed a medical history that included self-reported physical activity (minutes per week) and the Veterans version of the Short Form 36 Health Survey [25] for descriptive purposes. Experimental procedures described below were reviewed and approved by the Department of Veterans Affairs New Jersey Health Care System’s Institutional Review Board, and all participants provided informed written consent.
Spirometry
Participants were asked to abstain from consuming meals or caffeine for at least 4 hours prior to arriving to our laboratory. Spirometry was performed in accordance with standard guidelines [26] using commercially available equipment (Cosmed Quark PFT; Rome, Italy) that was calibrated prior to each participant. Spirometric indices of forced vital capacity (FVC), forced expiratory volume in one second (FEV1), and FEV1/FVC were obtained and expressed as a percent of predicted [27].
Cardiopulmonary exercise testing
Resting 12-lead electrocardiogram (ECG) (Cosmed T12x; Rome, Italy) and blood pressure via manual auscultation were acquired during quiet rest. CPET was then performed on a cycle ergometer (Ergoline Ergoselect 100; Germany) using a ramp protocol (15 watts·min-1) and pedal cadence of 50–70 rpms until volitional exhaustion [7]. Before, during and after CPET, all participants were asked to rate their perceived exertion (6–20) and breathlessness (0–10) using the Borg scales. Maximal effort was defined as meeting two or more of the following criteria: 1) ≥ 1.1 respiratory exchange ratio, 2) ≥ 85% age-predicted heart rate, 3) ≥ 17 perceived exertion, or 4) plateau or decline despite increasing workload.
Pulmonary gas-exchange ( and
),
, VT, fR, and end-tidal CO2 (PETCO2) were obtained breath-by-breath using an oro-nasal mask (V2 Series; Hans Rudolph, Shawnee, KS) with bidirectional turbine and metabolic cart (Cosmed Quark CPET; Rome, Italy) that was calibrated prior to each participant. Heart rate was obtained from 12-lead ECG and integrated with breath-by-breath data. The ventilatory anaerobic threshold (VAT) was determined using a modified V-slope method [28], and common CPET variables were obtained at peak exercise (30-sec time average) and reported by group. Our primary variables of interest (
, VT, and fR) were calculated at distinct periods relative to each individual’s peak exercise capacity (
peak) to facilitate between-group comparisons across similar intensities (i.e., 0% (rest), 20%, 40%, 60%, 80%, and 100%
peak). A 15-breath average, aligned to the time of the central breath of peak
, was used to compute the times of relative intensities. Analyses were performed using a custom-built program (MATLAB; Mathworks; Natick, MA).
Statistical analysis
Participant characteristics, baseline pulmonary function and peak exercise variables were analyzed using independent t-tests for continuous variables and Fisher’s exact test for categorical variables. Normality assumptions were checked (Kolmogorov-Smirnov) for our primary dependent variables (, VT, fR) across groups (GWI+ vs. GWI-) and time-points (0%, 20%, 40%, 60%, 80%, and 100%
peak). We used separate mixed-design 2 (group) X 6 time ANOVA models to examine between-group differences for each dependent variable (α = 0.05). Degrees of freedom were adjusted (Greenhouse-Geisser) when the sphericity assumption was violated. Subsequent pairwise comparisons used a Bonferroni adjustment to control family-wise error. Effect sizes are expressed as partial eta squared (η2p) for ANOVA/ANCOVA models and Hedges’d for pairwise comparisons [29]. In addition, we also performed sub-group analysis using 11 subjects from each group to control for the main determinants of vital capacity—i.e., race, age, sex, and height. Demographics for matched pairs are provided in supporting information (S1 Table). Through matching, the average difference in vital capacity was within 200 mL. Analyses were performed using IBM SPSS Statistics (v. 25).
Results
Participant characteristics
Demographics, self-reported physical and mental health, spirometric indices, and resting heart rate, blood pressure, and arterial blood oxygen saturation are provided in Tables 1 and 2. Descriptive CPET information collected during exercise is provided in Tables 3 and 4. Data are presented as group means and standard deviations (mean ± SD).
https://doi.org/10.1371/journal.pone.0224833.t001
https://doi.org/10.1371/journal.pone.0224833.t002
https://doi.org/10.1371/journal.pone.0224833.t003
https://doi.org/10.1371/journal.pone.0224833.t004
Minute ventilation
Tidal volume
Full sample.
The VT ANOVA model yielded significant main effects for time (F2.11,67.38 = 221.02, p < 0.001, η2p = 0.87) and group (F1,32 = 6.01, p = 0.02, η2p = 0.16). The group-by-time interaction effect was non-significant (F2.11,67.38 = 2.44, p = 0.09, η2p = 0.07). Pairwise comparisons showed significantly higher VT values at 20%, 40%, 60%, 80%, and 100% of each individual’s in the GWI+ group compared to the GWI- group (all p < 0.05; Table 5; Fig 1). Hedges’ d effect sizes for these comparisons ranged from 0.68 to 0.92.
https://doi.org/10.1371/journal.pone.0224833.t005
Cases (GWI+) and controls (GWI-) are represented in filled and open circles, respectively. Values are mean ± SE.
https://doi.org/10.1371/journal.pone.0224833.g001
Matched subgroup.
The VT ANOVA model yielded a significant main effect for time (F1.82,36.34 = 145.84, p < 0.001, η2p = 0.88), but the group effect was no longer significant (F1,20 = 1.11, p = 0.31, η2p = 0.05). The group-by-time interaction effect remained non-significant (F1.82,36.34 = 0.86, p = 0.42, η2p = 0.04).
Respiratory frequency
Full sample.
The fR ANOVA model yielded a significant main effect for time (F2.46,78.78 = 70.63, p < 0.001, η2p = 0.69) and group (F1,32 = 9.87, p = 0.004, η2p = 0.24). The group-by-time interaction effect was non-significant (F2.46,78.78 = 1.3, p = 0.28, η2p = 0.04). Pairwise comparisons showed significantly lower fR values at rest and 20%, 40%, 60%, and 80% of each individual’s in the GWI+ group compared to the GWI- group (all p < 0.05; Table 6; Fig 1). Hedges’ d effect sizes for these comparisons ranged from -0.61 to -1.05.
https://doi.org/10.1371/journal.pone.0224833.t006
Matched subgroup.
The fR ANOVA model yielded a significant main effect for time (F3.26,65.2 = 54.01, p < 0.001, η2p = 0.73) and group (F1,20 = 10.88, p = 0.004, η2p = 0.35). The group-by-time interaction effect was non-significant (F3.26,65.2 = 1.54, p = 0.21, η2p = 0.07). Pairwise comparisons showed significantly lower fR values at rest, and 20%, 40%, 60%, and 80% of each individual’s in the GWI+ group compared to the GWI- group (all p < 0.05; Table 6; Fig 2). Hedges’ d effect sizes for these comparisons ranged from -0.91 to -1.48.
Cases (GWI+) and controls (GWI-) are represented in filled and open circles, respectively. Values are mean ± SE.
https://doi.org/10.1371/journal.pone.0224833.g002
Discussion
The present study investigated ventilatory patterns (, VT, and fR) during maximal cardiopulmonary exercise (CPET) among veterans with GWI and controls. As expected,
, VT, and fR increased with increasing exercise intensity across groups (Fig 1). Veterans with GWI adopted a unique exercise ventilatory pattern; however, this pattern was contrary to our hypothesis and characterized by greater depth (VT) and reduced frequency (fR) of breathing relative to controls. According to the ANOVA models, the effect of group significantly accounted for 15 and 25% of the variability in VT and fR, respectively. Moreover, post-hoc pairwise comparisons at each time-point during exercise (i.e. 20%, 40%, 60%, 80%, and 100%
peak) revealed moderate-to-large between-group effects for VT at all time-points during exercise (Hedges’ drange: 0.75–0.92). For fR, post-hoc pairwise comparisons revealed large and significant group effects for 20%, 40%, 60%, 80% of
peak (Hedges’ drange: 0.94–1.36), but not 100% of
peak. After performing a sub-group analysis of participants who were closely matched for major determinants of vital capacity (i.e., race, sex, age, height), the group effect remained significant for fR but not for VT (Fig 2). This suggests that fR is an especially robust aspect of
distinguishing Veterans with and without GWI. This disparate breathing strategy was observed despite non-significant between-group differences for
(full sample and matched subgroup), which is important to note given that the components of
(VT and fR) are frequently overlooked [1] despite their clinical relevance (e.g., tachypnea in congestive heart failure).
To our knowledge, of the 10 CPET studies involving Gulf War veterans [12–21], only one study has reported detailed cardiopulmonary function of Gulf War veterans as assessed during a maximal CPET [12, 13]. In this study, CPET performance was found to be similar between veterans with and without GWI at peak and submaximal intensities, apart from peak which was reduced in GWI. This effect appears to be driven by a reduced fR (g = -0.46) among ill Gulf War veterans in comparison to healthy Gulf War veterans [13], similar to that observed in the present study. Maximal exercise ventilatory responses during maximal testing from related populations (e.g., chronic fatigue syndrome [CFS] and fibromyalgia) are limited, but have generally found similar fR and VT between cases and controls in patients with CFS [30, 31], CFS and fibromyalgia [30], or fibromyalgia alone [32]. It is important to note that in these studies of patients with CFS and/or fibromyalgia, a majority of the sample was female which may not provide an appropriate comparison to the participants of the present study.
Published reference data for fR and VT are scant for cycling exercise but are available for maximal treadmill testing [33]. Differences in operating lung volumes between cycling and treadmill exercise, most notably larger VT during cycling [34], precludes the application of treadmill based prediction equations for interpreting cycling studies. Neder et al. reported average peak fR obtained during incremental cycle exercise stratified by age ranges (20–39, 40–59, and 60–80 years) and sex at four different peak levels (20, 40, 60, and 80 L∙min-1) from 120 sedentary adults (c.f. Table 2 [35]). Using these data, we estimated a predicted fR value for each participant based on age range, sex and approximate peak
. Predicted fR was 99.4±22.4% and 116.4±23.5% for cases (GWI+) and controls (GWI-), respectively. Although we would have expected predicted fR to be even lower among veterans with GWI, there are several important factors that may account for these observed differences including the lack of precision around age ranges and peak
levels, race and ethnicity of the reference sample, and instrumentation. Regarding the latter, Neder et al. [35] utilized a nose clip and mouthpiece to acquire ventilatory variables which is known to result in reduced fR and increased VT in comparison to an oronasal mask [36], which was used in the present study.
At exercise onset and during lower intensity workloads, increases in are primarily driven by VT whereby further increases in
at greater workloads and at peak exercise are achieved primarily through increased fR. This typical ventilatory response to exercise is well characterized and thought to reflect a strategy to maximize metabolic efficiency [37]. Excessive increases in either VT or fR deviate from this principle and result in increased work of breathing and increased physiologic dead space. Accurate measurement of work of breathing and dead space ventilation (VD/VT) require invasive procedures (i.e., esophageal pressures and arterial catheter) that were not employed in the present study; therefore, we cannot definitively rule out their contribution. However, abnormal findings such as increased VD/VT and work of breathing appear unlikely for several reasons. First, we computed non-invasive estimates of VD/VT [38] and work of breathing [39] and found similar results among our cases and controls that were within normal limits. Increased VD/VT and work of breathing are most common among individuals with chronic lung disease, but lung disease was exclusionary for this study. Further, we observed no evidence of ventilatory limitation at rest via spirometry (Table 2) nor during peak exercise (VE/MVV ≤ 80%; Table 4). That VT continued to increase throughout exercise also supports the lack of a ventilatory limitation, or more specifically the absence of a restrictive process that would limit VT expansion. Lastly, both groups demonstrated a similar ratio of VT to inspiratory capacity (VT:IC) throughout exercise which was within normal limits (VT:IC ≤ 0.8; Table 4). We conclude that there is little evidence to support an underlying subclinical pulmonary pathology among cases with GWI+ that may account for the reduced rate and greater depth of breathing observed in the present study. Moreover, our matched subgroup analysis further confirms the robustness of these findings. We posit that the observed ventilatory pattern might instead reflect a learned strategy to mitigate exercise-induced pain and/or unique metabolic demands associated with GWI.
As reviewed by Forster et al. [40], numerous experimental studies support a neural feed-forward (i.e., central command) mediation of exercise hyperpnoea, suggesting central nervous system control of breathing involves a behavioral or learned mechanism. Long-term modulation of the exercise ventilatory response has been demonstrated experimentally in both humans [41, 42] and goats [43] via a conditioning paradigm whereby exercise is paired with added external dead space in order to induce hypercapnia. Following 2–8 visits of conditioning and removal of added dead space, subsequent steady-state normocapnic exercise resulted in increased VT during early exercise. In a follow-up study to determine whether long-term modulation was driven primarily by hypercapnia or increased VT, Turner and Stewart [44] employed inspiratory resistive loading as a model of associative conditioning to avoid inducing hypercapnia, and observed long-term modulation of the early exercise ventilatory response of similar magnitude and duration to the aforementioned prior work. Data from these studies suggest that at least some component of the normal exercise ventilatory response may be learned, and increased VT, rather than hypercapnia, may provide the requisite stimulus. The persistence of these adaptations, however, appears temporary which likely reflects the limited exposure to conditioning. Wood et al. [42] support this conclusion and suggest that a permanent breathing pattern adaptation is unlikely to occur following only an acute exposure. Whether alterations to the exercise ventilatory response (i.e., increased VT) can be sustained in humans with chronic conditioning remains unclear.
Cook et al.[16] have previously found that Gulf War veterans with chronic pain perceive exercise as more painful and effortful than healthy veterans without pain, and also have greater thermal pain sensitivity post-exercise. Based on these results, we initially hypothesized that fR would be exaggerated in GWI given that non-metabolic stimuli, such as pain, are well known to elicit robust increases in fR in healthy adults [1]. However, perhaps the ventilatory response to pain is unique among those with chronic illness who have been managing chronic pain symptoms for many years? In support, fear of exacerbating pain symptoms secondary to exercise appears to be prevalent in civilians with medically unexplained widespread chronic musculoskeletal pain [45] and this perception may lead to a conscious or subconscious alteration in exercise ventilatory patterns (i.e., slower and deeper). There are now several studies that support the notion that pain is influenced by respiration whereby deep and slow breathing reduces pain intensity and perception [46–49]. Whether veterans with GWI adopt a deep and slow ventilatory pattern to mitigate exercise-induced pain cannot be determined from our data, but warrants future investigation.
Notwithstanding the limitations of a cross-sectional study, there are additional aspects of this study that should be considered when interpreting these findings. First, although the size of the participant sample is within the range of previous CPET studies involving Gulf War veterans, GWI is a heterogeneous chronic multisymptom illness and replication of this work with a larger sample size is necessary to enhance generalizability. As previously mentioned, the paucity of published cardiopulmonary data from exercise studies in GWI precludes a thorough comparison with work from other investigators. In light of the present findings, future studies that include a more sophisticated analysis of ventilatory patterns during exercise to determine operating lung volumes (e.g., serial inspiratory capacity maneuvers) as well as measurement of the work of breathing via esophageal pressures appear warranted to fully appreciate the mechanical and neural factors that contribute to this presentation of increased depth of breathing during exercise. Finally, we did not comprehensively assess symptomatic changes during exercise. Considering the frequent use of exercise to elucidate illness pathophysiology and that the gold standard for GWI case assignment relies on self-report methods, an intuitive step when measuring physiological responses to exercise is to also measure perceptual responses. Thus, incorporating subjective measures for symptoms that are characteristic of the illness (e.g., pain, fatigue) and are sensitive to acute exercise-induced changes may aid future attempts to link changes in physiological parameters to pathophysiological indices of chronic multisymptom illnesses such as GWI [50].
Conclusions
Excessive increase in either VT or fR in response to exercise is an inefficient ventilatory strategy that contributes to poor performance and exertional symptoms; however, VT and fR are infrequently reported and often overlooked. In the present study, we observed a unique exercise ventilatory pattern among veterans with GWI characterized predominantly by an attenuated fR response to maximal exercise. The present study was not designed to assess the mechanism(s) of this response, but warrants further investigation. Confirmation of these findings may also provide a rationale on which to guide exercise prescription for veterans with GWI (e.g., wearable devices for fR [51]). In conclusion, researchers and clinicians should consider looking beyond and separately consider VT and fR responses to exercise, which appear to offer unique insight with respect to exercise hyperpnoea.
Supporting information
S1 Table. Matched-pair demographics.
For each case (GWI+) and control (GWI-) pair, we report i) race and ethnicity, ii) birth sex, iii) age in years, iv) height in centimeters, and v) forced vital capacity (FVC). Absolute differences for each variable for each pair (Δ = case—control) are reported for age, height and FVC. Overall, we were able to match evenly by sex and race/ethnicity. Other factors were very similar as indicated by the low average differences observed for age (3.7 years), height (6.5 cm), and FVC (0.2 L). Note that acceptable and repeatable spirometry were not available for one Veteran (pair 8).
https://doi.org/10.1371/journal.pone.0224833.s001
(DOCX)
Acknowledgments
The authors gratefully acknowledge the participants who volunteered for this study as well as the research personnel who helped acquire and analyze data, especially Bishoy Samy. The contents do not represent the views of the U.S. Department of Veterans Affairs or the United States Government.
References
- 1. Tipton MJ, Harper A, Paton JF, Costello JT. The human ventilatory response to stress: rate or depth? The Journal of Physiology. 2017.
- 2. Nicolo A, Marcora SM, Bazzucchi I, Sacchetti M. Differential control of respiratory frequency and tidal volume during high-intensity interval training. Exp Physiol. 2017;102(8):934–49. pmid:28560751.
- 3. Nicolo A, Bazzucchi I, Haxhi J, Felici F, Sacchetti M. Comparing continuous and intermittent exercise: an "isoeffort" and "isotime" approach. PLoS One. 2014;9(4):e94990. pmid:24736313.
- 4. Nicolo A, Marcora SM, Sacchetti M. Respiratory frequency is strongly associated with perceived exertion during time trials of different duration. J Sports Sci. 2016;34(13):1199–206. pmid:26503587.
- 5. Bell HJ, Duffin J. Rapid increases in ventilation accompany the transition from passive to active movement. Respir Physiol Neurobiol. 2006;152(2):128–42. Epub 2005/09/13. pmid:16153897.
- 6. Thornton JM, Guz A, Murphy K, Griffith AR, Pedersen DL, Kardos A, et al. Identification of higher brain centres that may encode the cardiorespiratory response to exercise in humans. J Physiol. 2001;533(Pt 3):823–36. pmid:11410638.
- 7. ATS/ACCP Statement on Cardiopulmonary Exercise Testing. 2003;167(2):211–77. pmid:12524257.
- 8. White RF, Steele L, O'Callaghan JP, Sullivan K, Binns JH, Golomb BA, et al. Recent research on Gulf War illness and other health problems in veterans of the 1991 Gulf War: Effects of toxicant exposures during deployment. Cortex. 2016;74:449–75. pmid:26493934
- 9. Kelsall HL, Sim MR, Forbes AB, McKenzie DP, Glass DC, Ikin JF, et al. Respiratory health status of Australian veterans of the 1991 Gulf War and the effects of exposure to oil fire smoke and dust storms. Thorax. 2004;59(10):897–903. pmid:15454658
- 10. Karlinsky JB, Blanchard M, Alpern R, Eisen SA, Kang H, Murphy FM, et al. Late prevalence of respiratory symptoms and pulmonary function abnormalities in Gulf War I Veterans. Arch Intern Med. 2004;164(22):2488–91. pmid:15596641
- 11. Medinger AE, Chan TW, Arabian A, Rohatgi PK. Interpretive algorithms for the symptom-limited exercise test: assessing dyspnea in Persian Gulf war veterans. Chest. 1998;113(3):612–8. pmid:9515833
- 12. Nagelkirk P, Cook D, Peckerman A, Kesil W, Sakowski T, Natelson B, et al. Aerobic capacity of Gulf War veterans with chronic fatigue syndrome. Military medicine. 2003;168(9):750. pmid:14529252
- 13. Cook DB, Nagelkirk PR, Peckerman A, Poluri A, LaManca JJ, Natelson BH. Perceived exertion in fatiguing illness: Gulf War veterans with chronic fatigue syndrome. Medicine and science in sports and exercise. 2003;35(4):569–74. pmid:12673138
- 14. Broderick G, Ben-Hamo R, Vashishtha S, Efroni S, Nathanson L, Barnes Z, et al. Altered immune pathway activity under exercise challenge in Gulf War Illness: An exploratory analysis. Brain, Behavior, and Immunity. 2013;28:159–69. http://dx.doi.org/10.1016/j.bbi.2012.11.007. pmid:23201588
- 15. Broderick G, Kreitz A, Fuite J, Fletcher MA, Vernon SD, Klimas N. A pilot study of immune network remodeling under challenge in Gulf War Illness. Brain, Behavior, and Immunity. 2011;25(2):302–13. pmid:20955779
- 16. Cook DB, Stegner AJ, Ellingson LD. Exercise Alters Pain Sensitivity in Gulf War Veterans With Chronic Musculoskeletal Pain. The Journal of Pain. 2010;11(8):764–72. pmid:20338824
- 17. Donta ST C DJ EJ CC, et al. Cognitive behavioral therapy and aerobic exercise for gulf war veterans' illnesses: A randomized controlled trial. JAMA. 2003;289(11):1396–404. pmid:12636462
- 18. Rayhan RU, Raksit MP, Timbol CR, Adewuyi O, VanMeter JW, Baraniuk JN. Prefrontal lactate predicts exercise-induced cognitive dysfunction in Gulf War Illness. American Journal of Translational Research. 2013;5(2):212–23. pmid:23573365
- 19. Rayhan RU, Stevens BW, Raksit MP, Ripple JA, Timbol CR, Adewuyi O, et al. Exercise Challenge in Gulf War Illness Reveals Two Subgroups with Altered Brain Structure and Function. PLoS ONE. 2013;8(6):e63903. pmid:23798990
- 20. Smylie AL, Broderick G, Fernandes H, Razdan S, Barnes Z, Collado F, et al. A comparison of sex-specific immune signatures in Gulf War illness and chronic fatigue syndrome. BMC Immunology. 2013;14(1):29. pmid:23800166
- 21. Whistler T, Fletcher M, Lonergan W, Zeng XR, Lin JM, LaPerriere A, et al. Impaired immune function in Gulf War Illness. BMC Medical Genomics. 2009;2(1):12.
- 22.
IOM. Chronic Multisymptom Illness in Gulf War Veterans: Case Definitions Reexamined. Washington, DC: The National Academies Press; 2014. 130 p.
- 23. Steele L. Prevalence and patterns of Gulf War illness in Kansas veterans: association of symptoms with characteristics of person, place, and time of military service. Am J Epidemiol. 2000;152(10):992–1002. Epub 2000/11/25. pmid:11092441.
- 24. Krupp LB, LaRocca NG, Muir-Nash J, Steinberg AD. The fatigue severity scale. Application to patients with multiple sclerosis and systemic lupus erythematosus. Arch Neurol. 1989;46(10):1121–3. pmid:2803071.
- 25. Falvo MJ, Serrador JM, McAndrew LM, Chandler HK, Lu SE, Quigley KS. A retrospective cohort study of U.S. service members returning from Afghanistan and Iraq: is physical health worsening over time? BMC Public Health. 2012;12:1124. Epub 2013/01/01. pmid:23272950.
- 26. Miller MR, Hankinson J, Brusasco V, Burgos F, Casaburi R, Coates A, et al. Standardisation of spirometry. EurRespirJ. 2005;26(2):319–38.
- 27. Hankinson J, Odencrantz J, Fedan K. Spirometric Reference Values from a Sample of the General U.S. Population. American Journal of Respiratory and Critical Care Medicine. 1999;159(1):179–87. pmid:9872837
- 28. Beaver WL, Wasserman K, Whipp BJ. A new method for detecting anaerobic threshold by gas exchange. J Appl Physiol (1985). 1986;60(6):2020–7. pmid:3087938.
- 29. Fritz CO, Morris PE, Richler JJ. Effect size estimates: current use, calculations, and interpretation. J Exp Psychol Gen. 2012;141(1):2–18. Epub 2011/08/10. pmid:21823805.
- 30. Cook DB, Nagelkirk PR, Poluri A, Mores J, Natelson BH. The influence of aerobic fitness and fibromyalgia on cardiorespiratory and perceptual responses to exercise in patients with chronic fatigue syndrome. Arthritis Rheum. 2006;54(10):3351–62. pmid:17009309
- 31. Inbar O, Dlin R, Rotstein A, Whipp BJ. Physiological responses to incremental exercise in patients with chronic fatigue syndrome. Med Sci Sports Exerc. 2001;33(9):1463–70. Epub 2001/08/31. pmid:11528333.
- 32. Sanudo B, Galiano D. Using cardiovascular parameters and symptom severity to prescribe physical activity in women with fibromyalgia. Clin Exp Rheumatol. 2009;27(5 Suppl 56):S62–6. Epub 2010/03/12. pmid:20074442.
- 33. Loe H, Steinshamn S, Wisloff U. Cardio-respiratory reference data in 4631 healthy men and women 20–90 years: the HUNT 3 fitness study. PLoS One. 2014;9(11):e113884. Epub 2014/11/27. pmid:25426954.
- 34. Tanner DA, Duke JW, Stager JM. Ventilatory patterns differ between maximal running and cycling. Respir Physiol Neurobiol. 2014;191:9–16. Epub 2013/11/12. pmid:24211317.
- 35. Neder JA, Dal Corso S, Malaguti C, Reis S, De Fuccio MB, Schmidt H, et al. The pattern and timing of breathing during incremental exercise: a normative study. Eur Respir J. 2003;21(3):530–8. Epub 2003/03/29. pmid:12662013.
- 36. Hirsch JA, Bishop B. Human breathing patterns on mouthpiece or face mask during air, CO2, or low O2. J Appl Physiol Respir Environ Exerc Physiol. 1982;53(5):1281–90. Epub 1982/11/01. pmid:6816769.
- 37. Sheel AW, Romer LM. Ventilation and respiratory mechanics. Compr Physiol. 2012;2(2):1093–142. pmid:23798297.
- 38. Pianosi P, D'Souza SJ, Esseltine DW, Charge TD, Coates AL. Ventilation and gas exchange during exercise in sickle cell anemia. Am Rev Respir Dis. 1991;143(2):226–30. Epub 1991/02/01. pmid:1990932.
- 39. Aaron EA, Johnson BD, Seow CK, Dempsey JA. Oxygen cost of exercise hyperpnea: measurement. J Appl Physiol (1985). 1992;72(5):1810–7. Epub 1992/05/01. pmid:1601790.
- 40. Forster HV, Haouzi P, Dempsey JA. Control of Breathing During Exercise. Comprehensive Physiology: John Wiley & Sons, Inc.; 2011.
- 41. Turner DL, Sumners DP. Associative conditioning of the exercise ventilatory response in humans. Respir Physiol Neurobiol. 2002;132(2):159–68. Epub 2002/08/06. pmid:12161329.
- 42. Wood HE, Fatemian M, Robbins PA. A learned component of the ventilatory response to exercise in man. J Physiol. 2003;553(Pt 3):967–74. pmid:14514870.
- 43. Martin PA, Mitchell GS. Long-term modulation of the exercise ventilatory response in goats. J Physiol. 1993;470:601–17. pmid:8308746.
- 44. Turner D, Stewart JD. Associative conditioning with leg cycling and inspiratory resistance enhances the early exercise ventilatory response in humans. Eur J Appl Physiol. 2004;93(3):333–9. pmid:15375661.
- 45. Turk DC, Robinson JP, Burwinkle T. Prevalence of fear of pain and activity in patients with fibromyalgia syndrome. The Journal of Pain. 2004;5(9):483–90. pmid:15556826
- 46. Zautra AJ, Fasman R, Davis MC, Craig AD. The effects of slow breathing on affective responses to pain stimuli: an experimental study. Pain. 2010;149(1):12–8. Epub 2010/01/19. pmid:20079569.
- 47. Chalaye P, Goffaux P, Lafrenaye S, Marchand S. Respiratory effects on experimental heat pain and cardiac activity. Pain Med. 2009;10(8):1334–40. Epub 2009/08/13. pmid:19671085.
- 48. Busch V, Magerl W, Kern U, Haas J, Hajak G, Eichhammer P. The effect of deep and slow breathing on pain perception, autonomic activity, and mood processing—an experimental study. Pain Med. 2012;13(2):215–28. Epub 2011/09/24. pmid:21939499.
- 49. Grant JA, Rainville P. Pain sensitivity and analgesic effects of mindful states in Zen meditators: a cross-sectional study. Psychosom Med. 2009;71(1):106–14. Epub 2008/12/17. pmid:19073756.
- 50. Lindheimer JB, Meyer JD, Stegner AJ, Dougherty RJ, Van Riper SM, Shields M, et al. Symptom variability following acute exercise in myalgic encephalomyelitis/chronic fatigue syndrome: a perspective on measuring post-exertion malaise. Fatigue: Biomedicine, Health & Behavior. 2017;5(2):69–88.
- 51. Nicolo A, Massaroni C, Passfield L. Respiratory Frequency during Exercise: The Neglected Physiological Measure. Front Physiol. 2017;8:922. Epub 2018/01/13. pmid:29321742.