- 1Division of Pulmonary Medicine, Department of Medicine, Faculty of Medicine and Dentistry, University of Alberta, Edmonton, AB, Canada
- 2G.F. MacDonald Centre for Lung Health, Covenant Health, Edmonton, AB, Canada
Background: Up to 53% of individuals who had mild COVID-19 experience symptoms for >3-month following infection (Long-CoV). Dyspnea is reported in 60% of Long-CoV cases and may be secondary to impaired exercise capacity (VO2peak) as a result of pulmonary, pulmonary vascular, or cardiac insult. This study examined whether cardiopulmonary mechanisms could explain exertional dyspnea in Long-CoV.
Methods: A cross-sectional study of participants with Long-CoV (n = 28, age 40 ± 11 years, 214 ± 85 days post-infection) and age- sex- and body mass index-matched COVID-19 naïve controls (Con, n = 24, age 41 ± 12 years) and participants fully recovered from COVID-19 (ns-CoV, n = 14, age 37 ± 9 years, 198 ± 89 days post-infection) was conducted. Participants self-reported symptoms and baseline dyspnea (modified Medical Research Council, mMRC, dyspnea grade), then underwent a comprehensive pulmonary function test, cardiopulmonary exercise test, exercise pulmonary diffusing capacity measurement, and rest and exercise echocardiography.
Results: VO2peak, pulmonary function and cardiac/pulmonary vascular parameters were not impaired in Long- or ns-CoV compared to normative values (VO2peak: 106 ± 25 and 107 ± 25%predicted, respectively) and cardiopulmonary responses to exercise were otherwise normal. When Long-CoV were stratified by clinical dyspnea severity (mMRC = 0 vs mMRC≥1), there were no between-group differences in VO2peak. During submaximal exercise, dyspnea and ventilation were increased in the mMRC≥1 group, despite normal operating lung volumes, arterial saturation, diffusing capacity and indicators of pulmonary vascular pressures.
Interpretation: Persistent dyspnea after COVID-19 was not associated with overt cardiopulmonary impairment or exercise intolerance. Interventions focusing on dyspnea management may be appropriate for Long-CoV patients who report dyspnea without cardiopulmonary impairment.
Introduction
Up to 87% of individuals hospitalized as a result of COVID-19 and 53% of individuals with mild initial illness experience symptoms 3-month after infection; these individuals are colloquially known as COVID Long-Haulers (Long-CoV) (Carfì et al., 2020; Petersen et al., 2020; World Health Organization, 2021). Dyspnea is reported in approximately 60% of Long-CoV cases; little evidence is available directly linking persistent dyspnea after COVID-19 to known physiological mechanisms for increased dyspnea sensation (Davis et al., 2021).
Long-CoV is predominately thought to be due to lingering symptoms from impairments incurred during active infection. However, emerging evidence suggests a dissociation between pathophysiology and symptom burden, particularly in non-severe COVID where initial cardiopulmonary damage may be trivial. Lam et al. (2021) reported a Long-CoV phenotype characterized by persistent dyspnea, impaired 6-min walk distance and reduced health-related quality of life, but without accompanying pulmonary function abnormalities or increased neurologic, musculoskeletal or fatigue symptoms which could contribute to dyspnea and reduced cardiopulmonary fitness. Conversely, in a study of 103 hospitalized patients with overt cardiopulmonary impairment, Lerum et al. (2021) report that at 3-month follow-up, forced vital capacity (FVC), forced expiratory volume in 1 s (FEV1), diffusing capacity for carbon monoxide (DLCO), and ground glass opacities and parenchymal bands on chest CT improved, independent of improvements in 6-min walk distance or dyspnea score. These findings suggest that while cardiopulmonary function and symptoms may improve concomitantly, cardiopulmonary improvement may not relate to symptom improvement. In a large cohort (n = 26,823), Long-CoV symptom burden had a stronger association with belief in having had COVID-19, rather than serology confirmed history of infection, suggesting that symptoms may be falsely attributed to COVID-19 and may be a result of other or pre-existing conditions (Matta et al., 2022). Indeed, while several studies have reported cardiopulmonary fitness in Long-CoV, these patient samples contain a high proportion of participants with substantial smoking history, obesity, and cardiovascular comorbidities, all of which confound the ability to isolate and investigate COVID-19-related dyspnea (Clavario et al., 2020; Jahn et al., 2021; Motiejunaite et al., 2021). Though these study samples are consistent with Long-CoV risk factors, it remains unclear whether cardiopulmonary findings are incidental or due to COVID-19.
Here, we sought to examine the cardiorespiratory responses to exercise in isolated Long-CoV and determine whether persistent symptoms/dyspnea are associated with altered cardiorespiratory function. In other pathologies, such as chronic obstructive pulmonary disease (COPD), interstitial lung disease (ILD), and pulmonary hypertension (PH), exertional dyspnea can be due to the inability to maintain appropriate gas exchange, an exaggerated ventilatory response to exercise, elevated operating lung volumes (i.e., dynamic hyperinflation), and/or elevated pulmonary vascular pressures, all of which have detrimental impacts on exercise capacity (VO2peak) (Arena and Sietsema, 2011; Faisal et al., 2015; Obokata et al., 2018). We hypothesized that: VO2peak would be impaired in Long-CoV; and that reduced VO2peak would be associated with 1) elevated exertional dyspnea, 2) a greater ventilatory response to exercise, 3) altered breathing mechanics, or 4) reduced DLCO.
Materials and methods
Ethical approval
This study was approved by the University of Alberta Human Research Ethics Board (Pro00107436), registered as a clinical trial (NCT04732663) and conducted in accordance with the Declaration of Helsinki. All participants provided written, informed consent prior to enrollment.
Design
A cross-sectional study design was used to compare participants with self-reported Long-CoV to participants whom were no-longer symptomatic from COVID-19 (ns-CoV) and COVID-19 naïve controls (Con).
Participants
Long-CoV participants were at least 12-week from first molecular test positivity (mean time from test positivity to first research visit = 219 ± 82 days) and experiencing at least one persistent symptom (ageusia, anosmia, cough, diarrhea, dyspnea, fatigue, fever, headache, muscle pain) at the time of testing (World Health Organization, 2021). Participants were recruited and tested between March and August 2021 (latest possible infection, April 2021); no participants had received a vaccination against COVID-19 prior to their infection. Participants were excluded for: 1) diagnosis of PH predating COVID-19; 2) absolute contraindication to exercise testing or orthopedic limitation; 3) age <18 or >65 years; and 4) body mass index (BMI) > 30 kg/m2 to limit confounding obesity effects on cardiopulmonary function. Participants were recruited from the University of Alberta/Kaye Edmonton Post-COVID clinic. Con and ns-CoV participants were recruited by word of mouth from the greater Edmonton area and matched to Long-CoV participants for age, sex and BMI.
Participants were then stratified by modified Medical Research Council (mMRC) dyspnea grade to dyspneic-CoV (mMRC≥1) and retrospectively matched to non-dyspneic-CoV (mMRC = 0 from ns- or Long-CoV), and controls for sex, age and BMI to explore relationships between impaired cardiorespiratory physiology and persistent dyspnea.
Procedures
Participants reported to the laboratory at the University of Alberta—Edmonton, Canada—for three visits.
Visit 1
Participants first completed a medical history and the mMRC Dyspnea Scale, Post-COVID Functional Scale (PCFS; 0, no functional limitations; 1, negligible functional limitations; 2, slight functional limitations; 3, moderate functional limitations; 4, severe functional limitations) and EuroQoL-5D-5L Visual Analog Scale (EQ-5D VAS, scored 0–100; 0, the worst health you can imagine; 100, the best health you can imagine) questionnaires (Mahler and Wells, 1988; Herdman et al., 2011; Klok et al., 2020).
A pulmonary function test (PFT) was conducted in accordance with American Thoracic Society/European Respiratory Society guidelines, results are reported as percent of predicted (Graham et al., 2017; Graham et al., 2019; Hall et al., 2021). Participants then underwent an incremental cardiopulmonary exercise test as previously described (Phillips et al., 2021). Briefly, an incremental test was conducted on an electronically braked cycle ergometer (Ergoselect II 1200, Ergoline, Blitz, Germany) starting at 0 W (unloaded cycling) and progressing by 20 W every 2 min. Ratings of perceived exertion (RPE) for dyspnea and leg discomfort (modified Borg scale, 1–10) (Borg, 1982) and inspiratory capacity (IC) maneuvers were conducted every 2-min. A good quality test was based on attainment of three of the four following criteria: 1) volitional exhaustion; 2) respiratory exchange ratio (RER) greater than 1.10; 3) maximal heart rate within 10 beats per minute (bpm) of age predicted maximum heart rate (208—[0.7 x age]) (Tanaka et al., 2001); 4) increase in oxygen consumption <100 ml/min with an increase in power output (Stickland et al., 2012). Expired gas was analyzed using a metabolic measurement system (Encore229 Vmax, SensorMedics, Yorba Linda, United States); arterial oxygen saturation (SpO2) was estimated using finger pulse oximetry (N-595 Oximax, Nellcor, Boulder, United States), heart rate was measured using a 12-lead electrocardiogram (CardioSoft, GE Medical Systems, Milwaukee, United States) and blood pressure was taken by manual auscultation. Measures were continuously recorded and averaged in 30-s blocks; VO2 peak values were compared using prediction equations from population-based samples (Neder et al., 1999; Lewthwaite et al., 2020).
Visit 2
The multiple fraction of inspired oxygen (FIO2)-DLCO technique was used to measure DLCO, pulmonary capillary blood volume (VCap) and membrane diffusing capacity (DM) at rest, 40 W, and 80% of peak power, as previously described (Tedjasaputra et al., 2017; Tedjasaputra et al., 2018). Briefly, DLCO measurements were performed using a 6-s breath-hold; participants waited a minimum of 4-min between breath-holds at rest, and 2-min between breath-holds during exercise. Participants cycled continuously for the 40 W stage, and for 2-min at 80% of their maximum prior to conducting each breath-hold and returning to lower intensity cycling for active recovery.
Visit 3
A brief echocardiographic study was performed at rest and during cycling exercise at 40 W (Ergoselect 1200 Stress Echo Supine Ergometer, Blitz, Germany) in optimal echocardiography position (supine/tilted). Images were collected by a Canadian Registered Cardiac Sonographer using a commercially available ultrasound device (Vivid Q, GE Healthcare, Fairfield, United States). To reduce bias, blinded analysis was performed using the manual measurement tool of a commercially available analysis software (Us2. ai, Singapore, Singapore) by certified imaging specialists in accordance with American Society of Echocardiography guidelines and independently verified (Mitchell et al., 2019).
Protocol modification
Ongoing data analysis of the first n = 22 Long-CoV participants versus n = 9 ns-CoV and n = 16 Con yielded no significant differences, or trends for differences in exercising DLCO, VCap, DM, or echocardiography derived measures. Therefore, the investigators ceased Visit 2 and 3 testing in the remaining n = 6 Long-CoV, n = 5 ns-CoV and n = 8 Con. A sensitivity analysis of VO2peak (percent predicted) was conducted and no differences were found in each group between those completing all visits, versus those who completed only Visit 1.
Study size
A 3.5 ml/kg/min reduction in VO2peak is associated with a 12% and 17% decrease in survival in men and women respectively (Myers et al., 2002; Gulati et al., 2003). It was calculated that sixteen participants in each group would be sufficient to detect a 3.5 ± 3.5 ml/kg/min mean difference in VO2peak between groups (α = 0.05, β = 0.80).
Statistics
Data are presented as mean ± standard deviation unless otherwise stated. Statistical significance was set a priori at p < 0.05. Participant characteristics, pulmonary function, highest equivalent work rate and peak CPET data were compared between groups by one-way analysis of variance (ANOVA). Shapiro-Wilk tests were used to test normality in all outcomes prior to ANOVA. Bonferroni T-tests were used when ANOVA revealed a main effect and/or interaction. All statistical analyses were performed using IBM SPSS Statistics 24 (IBM Corporation, Armonk, United States).
Results
Participants
Recruitment, eligibility, enrollment and testing details are displayed in Figure 1. Participant characteristics are displayed in Tables 1, 2. Data are reported for twenty-eight Long-CoV, all of whom were seeking medical care for persistent symptoms. Participants were well matched between Long-CoV and Con for age, sex and BMI; fewer ns-CoV (n = 14) were recruited, however, the participant characteristics of the sample were not different from Long-CoV or Con (n = 24). Self-reported symptoms during acute infection are reported in Table 1; notably, all participants with a history of COVID-19 experienced symptoms during active infection. Fifty-eight percent of Long-CoV participants reported persistent dyspnea. Self-reported overall health (EQ-5D) and physical functioning (PCFS) were reduced in Long-CoV relative to ns-CoV and Con (Table 1). Participant characteristics after grouping by mMRC dyspnea grade are displayed in Table 2. Dyspneic-CoV participants were free of cardiovascular comorbidities and conditions that could provide an alternative explanation for dyspnea. Self-reported overall health (EQ-5D) and physical functioning (PCFS) were reduced in dyspneic-CoV relative to mMRC = 0 and Con groups.
TABLE 1. Participant characteristics, symptoms during acute infection, and time since molecular confirmed positive COVID-19 test for Long-CoV, no-longer symptomatic (ns)-CoV and Controls. BMI, body mass index; BSA, body surface area; PCFS, Post-COVID Functional Scale; EQ-5D VAS, EuroQol 5-Dimension visual analog scale.
TABLE 2. Participant characteristics, symptoms during acute infection, and time since molecular confirmed positive COVID-19 test for COVID-19 groups after stratification by mMRC dyspnea grade (0 = normal, ≥1 = dyspneic) and Controls. BMI, body mass index; BSA, body surface area; PCFS, Post-COVID Functional Scale; EQ-5D VAS, EuroQol 5-Dimension visual analog scale.
FIGURE 1. Participant recruitment, eligibility, enrollment and testing. aOngoing analysis showed no evidence of differences in exercising DLCO or echocardiography derived measurements. Therefore, investigators declined testing participants on Visits 2 and 3 as testing was not justifiable.
Pulmonary function testing
No differences were found between Long-CoV, ns-CoV, or Con for resting pulmonary function and means were all within normal limits (Table 3). When participants were split by mMRC grade, no differences in resting pulmonary function were found between groups, and means were within normal limits (Table 4). Three mMRC ≥1 participants, 4 mMRC = 0 participants, and 1 Con participant had an FEV1:FVC ratio below the LLN; 2 participants in the mMRC ≥1 group had a TLC below the LLN; and 1 participant in each of the mMRC = 0 and ≥1 groups had a DLCO below the LLN.
TABLE 3. Resting pulmonary function (pulmonary function test, PFT) and submaximal and peak exercise data for Long-CoV, no-longer symptomatic (ns)-CoV and Controls.
TABLE 4. Resting pulmonary function (pulmonary function test, PFT) and submaximal and peak exercise data for mMRC≥1 (dyspneic-CoV), mMRC = 0 and Controls.
Cardiopulmonary exercise testing
Long-CoV and ns-CoV had normal VO2peak (mean 106 ± 25 and 107 ± 25%predicted, respectively), but significantly lower VO2peak than Con (130%predicted) (Table 3; Figure 2A). Submaximal oxygen uptake and DLCO responses to exercise were similar for all groups (Figures 2A,C). Arterial saturation (Table 3) and the components of diffusing, DM and VCap (Table 3), were not different between groups, indicating normal gas exchange.
FIGURE 2. The (A) oxygen uptake (VO2)- and (B) ventilation (VE)-power relationships and the diffusing capacity for carbon monoxide (DLCO)-VO2 relationships in Long-CoV (circles, solid line), no-longer symptomatic (ns)-CoV (triangles, dashed line) and Con (squares, dotted line). † = p < 0.05 Long-CoV vs Con.
At the highest equivalent work rate (80 W), ventilation (VE, Long-CoV; 39 ± 8 vs. Con; 34 ± 5 L/min, p = 0.01) and the VE/VCO2 nadir were higher (i.e., less efficient ventilation) in Long-CoV relative to Con (Long-CoV; 29.8 ± 3.4 vs Con; 27.0 ± 3.6, p = 0.05), and respiratory rate (Long-CoV; 26 ± 6 vs. Con; 23 ± 5 breaths/minute, p = 0.07) trended to be higher, while end-tidal CO2 (PETCO2) was not different (Long-CoV; 36.3 ± 4.0 vs. Con; 38.0 ± 2.4 mmHg, p = 0.22) (Figure 2B, Figures 3A–C). At the highest equivalent work rate (80 W), no differences were found between Long-CoV, ns-CoV and Con for operating lung volumes and inspiratory reserve volume (Figures 4A,B). Dyspnea relative to VE and dyspnea relative to inspiratory reserve volume with incremental exercise were not different across groups (Figures 4B,C).
FIGURE 3. The (A) ventilatory efficiency (VE/VCO2)-, (B) end-tidal CO2−, and (C) respiratory rate (RR)-power relationships in Long-CoV (circles, solid line), no-longer symptomatic (ns)-CoV (triangles, dashed line) and Con (squares, dotted line). No significant differences.
FIGURE 4. (A) The operating lung volume-power (open, end-inspiratory lung volume, closed, end-expiratory lung volume); (B) dyspnea-inspiratory reserve volume (as percent of total lung capacity, TLC); (C) dyspnea-power (Borg rating of perceived exertion, RPE); and (D) dyspnea-ventilation relationships in Long-CoV (circles, solid line), no-longer symptomatic (ns)-CoV (triangles, dashed line) and Con (squares, dotted line). No significant differences.
After stratification by mMRC, no difference in VO2peak was found between COVID-19 groups (Table 4; Figure 5A). Oxygen uptake, and DLCO responses to exercise were similar for all groups (Figures 5A,C). Arterial saturation (Table 4) and the components of diffusing capacity, DM and VCap (Table 4), were not different between groups, indicating normal gas exchange. At the highest equivalent work rate (80 W), VE (dyspneic-CoV; 40 ± 9 vs. Con; 34 ± 5 L/min, p = 0.03) and dyspnea (dyspneic-CoV; 3 ± 2 vs. Con; 2 ± 1 Borg units, p = 0.04) were significantly higher for dyspneic-CoV as compared to Con, and respiratory rate (dyspneic-CoV; 27 ± 7 vs. Con; 23 ± 4 breaths/minute, p = 0.10) and the VE/VCO2 nadir (30.4 ± 4.0 vs 27.9 ± 2.3, p = 0.08) trended to be higher, while PETCO2 (mMRC ≥1; 35.9 ± 4.7 vs. Con; 38.2 ± 2.4 mmHg, p = 0.11) trended to be lower (Figures 6A–C). At the highest equivalent work rate, no differences were found between mMRC ≥1, mMRC = 0 and Con for operating lung volumes and inspiratory reserve volume (Figures 7A,B). Dyspnea relative to VE and dyspnea relative to inspiratory reserve volume with incremental exercise were not different across groups (Figures 7B,C).
FIGURE 5. The (A) oxygen uptake (VO2)- and (B) ventilation (VE)-power relationships and the diffusing capacity for carbon monoxide (DLCO)-VO2 relationships in mMRC > 1 (dyspneic-CoV, circles, solid line), mMRC = 0 (triangles, dashed line) and Con (squares, dotted line). † = p < 0.05 mMRC≥1 vs Con.
FIGURE 6. The ventilatory efficiency- (VE/VCO2, (A), end-tidal CO2- (B) and respiratory rate- (RR, (C) power relationships in mMRC≥1 (dyspneic-CoV, circles, solid line), mMRC = 0 (triangles, dashed line) and Con (squares, dotted line). No significant differences.
FIGURE 7. (A) The operating lung volume-power (open, end-inspiratory lung volume, closed, end-expiratory lung volume); (B) dyspnea-inspiratory reserve volume (as percent of total lung capacity, TLC); (C) dyspnea-power (Borg rating of perceived exertion, RPE); and (D) dyspnea-ventilation relationships in mMRC≥1 (dyspneic-CoV, circles, solid line), mMRC = 0 (triangles, dashed line) and Con (squares, dotted line). † = p < 0.05 mMRC≥1 vs Con.
Echocardiography
Results from rest and exercise echocardiography are displayed in Table 5, 6. No differences were found between Long-CoV and Con for right ventricle diastolic area, fractional area change, tricuspid annular plane systolic excursion (TAPSE) and the ratio of peak early mitral inflow velocity to early mitral annular diastolic velocity (E/e’). Pulmonary acceleration time was lower in Long-CoV relative to Con at rest, but similar during exercise. No differences were found between dyspneic-CoV and Con for right ventricle diastolic area or fractional area change, TAPSE, pulmonary acceleration time and E/e’ (left ventricle). Tricuspid regurgitant jet velocity could only be measured in 27% of participants; therefore, pulmonary artery systolic pressure was not calculated.
TABLE 5. Rest and exercise (40 W) echocardiography for Long-CoV, no-longer symptomatic (ns)-CoV and Controls. TAPSE; tricuspid annular plane systolic excursion.
TABLE 6. Rest and exercise (40 W) echocardiography data for mMRC≥1 (dyspneic-CoV), mMRC = 0 and Controls. TAPSE, tricuspid annular plane systolic excursion.
Discussion
Exercise testing is a valuable tool to evaluate functional capacity (i.e. VO2peak) and mechanisms of dyspnea and exercise intolerance (Radtke et al., 2019; Stickland et al., 2022). This study examined exercise responses in a sample of relatively young, non-obese and comorbidity-free COVID-19 participants to determine whether persistent symptoms/dyspnea were associated with altered cardiorespiratory function. In Long-CoV participants, ventilation for a given workload during submaximal exercise was elevated; however, resting cardiopulmonary function, VO2peak, and cardiopulmonary responses to exercise were otherwise normal. Similarly, when participants with a history of COVID-19 were grouped according to baseline dyspnea (mMRC grade 0 vs. ≥ 1), we found that dyspneic-CoV participants had greater dyspnea and ventilation for a given workload during submaximal exercise, while resting cardiopulmonary function, VO2peak and cardiopulmonary responses to exercise were normal. COVID-19 participants who report greater baseline dyspnea appear to have greater ventilation during submaximal exercise, which would contribute to an increased perception of exertional dyspnea. Importantly, during submaximal exercise, dyspnea appropriately matched ventilation, and the greater ventilation in dyspneic-CoV participants did not appear to be due to alterations in gas exchange, lung mechanics or any other impairment in cardiopulmonary physiology.
Pulmonary function
Previous studies have reported pulmonary function patterns consistent with restrictive defects in COVID-19, which normalize over time (Anastasio et al., 2021; Torres-Castro et al., 2021). In contrast to previous studies, the vast majority of participants in the present study were not hospitalized and Long- and dyspneic-CoV participants reported persistent dyspnea in the absence of restrictive defect or other pulmonary function impairment. The results of the present study reinforce findings by Lam et al., which suggest that a PFT may be insufficient to evaluate abnormal dyspnea after COVID-19 (Lam et al., 2021).
Mechanisms of exertional dyspnea
In obstructive and restrictive lung disease, heart failure (HF) and PH, dyspnea and exercise intolerance can be attributed to; impaired gas exchange, exaggerated ventilatory responses to exercise, dynamic hyperinflation (early critical tidal volume mechanical constraint), and/or elevated pulmonary vascular pressures (Arena and Sietsema, 2011; Guenette et al., 2014; Faisal et al., 2015; Obokata et al., 2018). Accordingly, each of these potential mechanisms of dyspnea were investigated.
Gas exchange
Given the extensive reporting of arterial hypoxemia and impaired DLCO in COVID-19, impaired gas exchange was a strong candidate mechanism of abnormal exertional dyspnea in Long- and dyspneic-CoV (Torres-Castro et al., 2021). Lerum et al. (2021) reported that 25% of patients hospitalized for COVID-19 had a DLCO below the LLN and ground-glass opacities (CT) 3-months after discharge; however, these abnormalities were not associated with dyspnea, arterial hypoxemia, or a reduction in 6-min walk distance. In the present study, only 2 (7%) COVID-19 participants had an SpO2 below 90% at peak exercise, and diffusing capacity (Figures 2C, 5C), DM and VCap were not different between Long- or dyspneic-CoV and controls (Tables 3, 4). These data indicate that gas exchange abnormalities are not requisite for persistent dyspnea in COVID-19 and do not explain abnormal exertional dyspnea during submaximal exercise in the COVID participants examined in the present study.
Ventilatory response
Ventilation, as measured at the highest equivalent submaximal workload, was elevated in Long- and dyspneic-CoV as compared to Con. Importantly, the dyspnea-ventilation relationship with incremental exercise was similar between Long- and dyspneic-CoV and controls, indicating that COVID-19 participants do not sense dyspnea differently during exercise relative to controls. Therefore, the greater dyspnea observed during submaximal exercise, but not at rest, in dyspneic-CoV could be explained by an exaggerated ventilatory response to exercise and the associated increased respiratory neural drive (Guenette et al., 2014). It is unclear why ventilation was increased in dyspneic-CoV, as the respiratory pattern does not completely align with other cardio-circulatory and respiratory diseases (Hansen and Wasserman, 1996; Arena et al., 2004; Neder et al., 2016). Dyspneic-CoV participants displayed a trend of mild tachypnea, but maintained normal tidal volume—in contrast to HF, PH and ILD where a rapid, shallow breathing pattern leads to increased anatomical deadspace ventilation and increased VE/VCO2 (Ponikowski et al., 2001; Wensel et al., 2013; Neder et al., 2016). Moreover, dyspneic-CoV participants displayed normal DLCO and VCap—contrasting with mild COPD where pulmonary vascular dysfunction marked by reduced DLCO and VCap leads to increased alveolar deadspace ventilation and increased VE/VCO2 (Tedjasaputra et al., 2018; Phillips et al., 2021). Finally, hyperventilation is observed across cardio-circulatory and respiratory diseases which can be secondary to muscle afferent feedback and deconditioning (Hansen and Wasserman, 1996; Arena et al., 2004; Neder et al., 2016). Peripheral deconditioning leading to an exaggerated ventilatory response to exercise is supported by a case series of eight dyspneic COVID-19 patients, 3-months after mild COVID-19 (Motiejunaite et al., 2021). In these patients, PFTs were normal, VE/VCO2 was elevated, and no patients reached their predicted VO2peak, which could be explained by greater peripheral fatigue (Motiejunaite et al., 2021). Our results extend and contrast these findings by reporting a similarly exaggerated ventilatory response to exercise, exertional dyspnea, and a trend of reduced ventilatory efficiency, but in the absence of VO2peak impairment, suggesting that peripheral deconditioning is an unlikely explanation for increased ventilation during submaximal exercise.
Operating lung volumes
Airway obstruction leading to dynamic hyperinflation and critical mechanical constraint is a well-known mechanism of dyspnea in COPD (O'Donnell et al., 2017). Long- and dyspneic-CoV participants did not demonstrate dynamic hyperinflation nor infringement on inspiratory reserve volume, as tidal volume was augmented appropriately, and end-expiratory and end-inspiratory lung volumes responded appropriately to exercise (Figures 4A, 7A). Thus, operating lung volume responses were normal and would not explain abnormal exertional dyspnea in the COVID participants examined in the current study.
Pulmonary vascular pressure
In contrast to previous reports, we did not identify structural or functional cardiac changes in Long- or dyspneic-CoV (Szekely et al., 2020). Aspects of the exercise response in dyspneic-CoV such as trends for tachypnea and hyperventilation are similar to what would be observed in conditions such as PH, but indicators of pulmonary vascular pressure (pulmonary acceleration time, E/e’) in the present sample were not elevated at rest or during exercise (Laveneziana et al., 2013). Moreover, metrics of right ventricular function (TAPSE, fractional area change) were not impaired in dyspneic or Long-CoV (Tables 5, 6). Taken together, these data imply that pulmonary vascular pressures and cardiac physiology are unlikely to contribute to persistent dyspnea in our Long- or dyspneic-CoV participants.
Caution is warranted when interpreting findings of increased ventilation and dyspnea during submaximal exercise. As increased ventilation in dyspneic-CoV does not appear to be a compensatory mechanism to maintain adequate gas exchange and the ventilation-dyspnea relationship was preserved, findings may partially be due to the fitness discrepancy, possibly predating COVID-19 infection. However, at lower work rates (40, 60 W) which more closely approximate the metabolic cost of walking up a slight hill (i.e. mMRC = 1), differences in submaximal exercise metrics were not detectable. This reinforces that factors beyond cardiopulmonary physiology appear to contribute to perceived symptom burden and exertional dyspnea after COVID-19.
Limitations
A limitation of the study is the cross-sectional design. VO2peak was lower in Long- and dyspneic-CoV relative to Con; however, we believe that this is due to selection bias, whereby exceptionally healthy participants volunteered for control groups, while symptomatic participants were recruited from Long-COVID clinics. However, we cannot rule out the possibility of a COVID-19-related drop in cardiopulmonary fitness. Given that pulmonary dysfunction and cardiopulmonary abnormalities were not observed in Long- or dyspneic-CoV, and that VO2peak in Long- and dyspneic-CoV was ∼106 and 98% of age-predicted normative data, respectively, a reduction in VO2peak as a result of COVID-19 infection is unlikely.
The results of this study should be extended with caution, particularly to those with severe initial infection and COVID-19 pneumonia requiring hospitalization. Only a small proportion of the sample included in the study were hospitalized during acute infection. Evidence suggests that after pneumonia and hospitalization, symptoms and functional impairment may persist for 6-months or longer, and that some symptoms may be attributable to underlying comorbidities (Metlay et al., 1997; Fumagalli et al., 2021).
Our sample consisted of a female to male ratio of ∼4:1 reporting Long-CoV and baseline dyspnea. The high female prevalence was consistent with our Long-CoV clinic, and others have reported female sex as a risk factor for Long-CoV (Sudre et al., 2021). More research is needed as it remains unclear whether potential sex differences in immune response, cardiopulmonary physiology, or other factors, such as care seeking behaviors, explain the female predominance of Long-CoV (Sudre et al., 2021).
Finally, emerging evidence suggests a potential role of autonomic dysfunction in cardiopulmonary impairment after COVID-19 (Shouman et al., 2021; Stute et al., 2021). Shouman et al. (2021) report impaired heart rate and blood pressure responses to physiologic stress, and orthostatic intolerance temporally related to COVID-19 infection that does not differ compared to orthostatic intolerance encountered following other viral disorders. Autonomic function was not measured in the present study, however, cardiopulmonary responses to exercise (including heart rate and blood pressure) and peak exercise capacity were normal in Long- and dyspneic-CoV.
Carotid chemoreceptor activity/sensitivity can play a role in cardiorespiratory regulation during exercise (Byers et al., 2019). We did not assess chemoreceptor activity/sensitivity in the present study; however, PETCO2 was not different across groups during exercise (Figures 3B, 6B), suggesting similar ventilatory responses and carotid chemoreceptor sensitivity across groups. Taken together, these data suggest that while autonomic function may be impaired in isolated Long-CoV cases, autonomic dysfunction is unlikely to explain elevated exertional dyspnea after COVID-19 in our sample.
Conclusion
The high incidence of Long-CoV, even after a non-severe COVID-19 infection, imposes a monumental healthcare burden as patients seek care for persistent dyspnea and other symptoms. Our data provide insight into persistent dyspnea after COVID-19 by demonstrating that in young, non-obese and comorbidity-free COVID-19 participants, dyspnea is not due to overt cardiopulmonary impairment or exercise intolerance, and that dyspnea sensation pathways are normal. Factors beyond cardiopulmonary physiology likely contribute to symptom burden, indicating that interventions focusing on dyspnea management may be appropriate for the phenotype of Long-CoV patients who report dyspnea.
Data availability statement
The raw data supporting the conclusion of this article will be made available by the authors, without undue reservation.
Ethics statement
The studies involving human participants were reviewed and approved by University of Alberta Human Research Ethics Board. The patients/participants provided their written informed consent to participate in this study.
Author contributions
RIB and MKS had full access to all of the data in the study and take responsibility for the integrity of the data and the accuracy of the data analysis. Concept and Design: All authors. Drafting of the manuscript: RIB, ARB, and MKS Critical revision of the manuscript for important intellectual content: All authors. Statistical analysis: RIB. Administrative and technical support: ARB, SdW, and DPF. Medical Oversight: RAV, RWD, GF, GYL, and MPS.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Nicole Plicka as well as Drs. Tracey Bryan, Sean Van Diepen, Cheryl Laratta, Angela Lau, and Eric Wong, and supporting staff at the University of Alberta/Kaye Edmonton post-COVID clinic.
Abbreviations
ANOVA, Analysis of variance; BMI, Body mass index; Con, Control group or participants; COPD, Chronic obstructive pulmonary disease; COVID-19, Corona Virus Disease 2019; Long-CoV, persistently symptomatic COVID-19; ns-CoV, no-longer symptomatic COVID-19; CT, Computed tomography; DLCO, Lung diffusing capacity for carbon monoxide; DM, Lung diffusing membrane capacity; E/e’, Ratio of early mitral filling velocity to early mitral annular tissue velocity; EQ-5D, EuroQoL 5-dimension; FEV1, Forced expiratory volume in 1-s; FiO2, Fraction of inspired oxygen; FVC, Forced vital capacity; HF, Heart failure; ILD, Interstitial lung disease; LLN, Lower limit of normal; mMRC, Modified Medical Research Council dyspnea grade; PCFS, Post-COVID functional scale; PETCO2, Pressure of end-tidal carbon dioxide; PFT, Pulmonary function test; PH, Pulmonary hypertension; TAPSE, Tricuspid annular plane systolic excursion; TLC, Total lung capacity; VCap, Pulmonary capillary blood volume; VE, Minute ventilation; VE/VCO2, Ventilation with respect to carbon dioxide production; VO2peak, peak oxygen uptake; W, Watts.
References
Anastasio F., Barbuto S., Scarnecchia E., Cosma P., Fugagnoli A., Rossi G., et al. (2021). Medium-term impact of COVID-19 on pulmonary function, functional capacity and quality of life. Eur. Respir. J. 58, 2004015. doi:10.1183/13993003.04015-2020
Arena R., Myers J., Aslam S. S., Varughese E. B., Peberdy M. A. (2004). Peak VO2 and VE/VCO2 slope in patients with heart failure: a prognostic comparison. Am. Heart J. 147 (2), 354–360. doi:10.1016/j.ahj.2003.07.014
Arena R., Sietsema K. E. (2011). Cardiopulmonary exercise testing in the clinical evaluation of patients with heart and lung disease. Circulation 123 (6), 668–680. doi:10.1161/circulationaha.109.914788
Borg G. A. V. (1982). Psychophysical bases of perceived exertion. Med. Sci. Sports Exerc. 14 (5), 377–381. doi:10.1249/00005768-198205000-00012
Byers B. W., Fuhr D. P., Moore L. E., Bhutani M., Wong E. Y. L., Stickland M. K. (2019). The effect of pulmonary rehabilitation on carotid chemoreceptor activity and sensitivity in chronic obstructive pulmonary disease. J. Appl. Physiol. 127 (5), 1278–1287. doi:10.1152/japplphysiol.00799.2018
Carfì A., Bernabei R., Landi F., and Gemelli Against COVID-19 Post-Acute Care Study Group (2020). Persistent symptoms in patients after acute COVID-19. JAMA 324 (6), 603–605. doi:10.1001/jama.2020.12603
Clavario P., De Marzo V., Lotti R., Barbara C., Porcile A., Russo C., et al. (2020). Assessment of functional capacity with cardiopulmonary exercise testing in non-severe COVID-19 patients at three months follow-up. medRxiv. doi:10.1101/2020.11.15.20231985
Davis H. E., Assaf G. S., McCorkell L., Wei H., Low R. J., Re'em Y., et al. (2021). Characterizing long COVID in an international cohort: 7 months of symptoms and their impact. EClinicalMedicine 38, 101019. doi:10.1016/j.eclinm.2021.101019
Faisal A., Alghamdi B. J., Ciavaglia C. E., Elbehairy A. F., Webb K. A., Ora J., et al. (2015). Common mechanisms of dyspnea in chronic interstitial and obstructive lung disorders. Am. J. Respir. Crit. Care Med. 193 (3), 299–309. doi:10.1164/rccm.201504-0841OC
Fumagalli A., Misuraca C., Bianchi A., Borsa N., Limonta S., Maggiolini S., et al. (2021). Long-term changes in pulmonary function among patients surviving to COVID-19 pneumonia. Infection [Epub ahead of print], 1–4. doi:10.1007/s15010-021-01718-2
Graham B. L., Brusasco V., Burgos F., Cooper B. G., Jensen R., Kendrick A., et al. (2017). 2017 ERS/ATS standards for single-breath carbon monoxide uptake in the lung. Eur. Respir. J. 49 (1), 1600016. doi:10.1183/13993003.00016-2016
Graham B. L., Steenbruggen I., Miller M. R., Barjaktarevic I. Z., Cooper B. G., Hall G. L., et al. (2019). Standardization of spirometry 2019 update. An official american thoracic society and european respiratory society technical statement. Am. J. Respir. Crit. Care Med. 200 (8), e70–e88. doi:10.1164/rccm.201908-1590st
Guenette J. A., Chin R. C., Cheng S., Dominelli P. B., Raghavan N., Webb K. A., et al. (2014). Mechanisms of exercise intolerance in global initiative for chronic obstructive lung disease grade 1 COPD. Eur. Respir. J. 44 (5), 1177–1187. doi:10.1183/09031936.00034714
Gulati M., Pandey D. K., Arnsdorf M. F., Lauderdale D. S., Thisted R. A., Wicklund R. H., et al. (2003). Exercise capacity and the risk of death in women: the St James women take heart project. Circulation 108 (13), 1554–1559. doi:10.1161/01.cir.0000091080.57509.e9
Hall G. L., Filipow N., Ruppel G., Okitika T., Thompson B., Kirkby J., et al. (2021). Official ERS technical standard: global lung function initiative reference values for static lung volumes in individuals of european ancestry. Eur. Respir. J. 57 (3), 2000289. doi:10.1183/13993003.00289-2020
Hansen J. E., Wasserman K. (1996). Pathophysiology of activity limitation in patients with interstitial lung disease. Chest 109 (6), 1566–1576. doi:10.1378/chest.109.6.1566
Herdman M., Gudex C., Lloyd A., Janssen M., Kind P., Parkin D., et al. (2011). Development and preliminary testing of the new five-level version of EQ-5D (EQ-5D-5L). Qual. Life Res. 20 (10), 1727–1736. doi:10.1007/s11136-011-9903-x
Jahn K., Sava M., Sommer G., Schumann D. M., Bassetti S., Siegemund M., et al. (2021). Exercise capacity impairment after COVID-19 pneumonia is mainly caused by deconditioning. Eur. Respir. J. 59, 2101136. doi:10.1183/13993003.01136-2021
Klok F. A., Boon G. J. A. M., Barco S., Endres M., Miranda Geelhoed J. J., Knauss S., et al. (2020). The post-COVID-19 functional status (PCFS) scale: a tool to measure functional status over time after COVID-19. Eur. Respir. J. 56, 2001494. doi:10.1183/13993003.01494-2020
Lam G. Y., Befus A. D., Damant R. W., Ferrara G., Fuhr D. P., Stickland M. K., et al. (2021). Exertional intolerance and dyspnea with preserved lung function: an emerging long COVID phenotype? Respir. Res. 22 (1), 222. doi:10.1186/s12931-021-01814-9
Laveneziana P., Garcia G., Joureau B., Nicolas-Jilwan F., Brahimi T., Laviolette L., et al. (2013). Dynamic respiratory mechanics and exertional dyspnoea in pulmonary arterial hypertension. Eur. Respir. J. 41 (3), 578–587. doi:10.1183/09031936.00223611
Lerum T. V., Aaløkken T. M., Brønstad E., Aarli B., Ikdahl E., Lund K. M. A., et al. (2021). Dyspnoea, lung function and CT findings 3 months after hospital admission for COVID-19. Eur. Respir. J. 57 (4), 2003448. doi:10.1183/13993003.03448-2020
Lewthwaite H., Benedetti A., Stickland M. K., Bourbeau J., Guenette J. A., Maltais F., et al. (2020). Normative peak cardiopulmonary exercise test responses in canadian adults aged ≥40 years. Chest 158 (6), 2532–2545. doi:10.1016/j.chest.2020.06.074
Mahler D. A., Wells C. K. (1988). Evaluation of clinical methods for rating dyspnea. Chest 93 (3), 580–586. doi:10.1378/chest.93.3.580
Matta J., Wiernik E., Robineau O., Carrat F., Touvier M., Severi G., et al. (2022). Association of self-reported COVID-19 infection and SARS-CoV-2 serology test results with persistent physical symptoms among french adults during the COVID-19 pandemic. JAMA Intern Med. 182 (1), 19–25. doi:10.1001/jamainternmed.2021.6454
Metlay J. P., Fine M. J., Schulz R., Marrie T. J., Coley C. M., Kapoor W. N., et al. (1997). Measuring symptomatic and functional recovery in patients with community-acquired pneumonia. J. Gen. Intern Med. 12 (7), 423–430. doi:10.1046/j.1525-1497.1997.00074.x
Mitchell C., Rahko P. S., Blauwet L. A., Canaday B., Finstuen J. A., Foster M. C., et al. (2019). Guidelines for performing a comprehensive transthoracic echocardiographic examination in adults: recommendations from the american society of echocardiography. J. Am. Soc. Echocardiogr. 32 (1), 1–64. doi:10.1016/j.echo.2018.06.004
Motiejunaite J., Balagny P., Arnoult F., Mangin L., Bancal C., d’Ortho M.-P., et al. (2021). Hyperventilation: a Possible explanation for long-lasting exercise intolerance in mild COVID-19 survivors? Front. Physiol. 11, 614590. doi:10.3389/fphys.2020.614590
Myers J., Prakash M., Froelicher V., Do D., Partington S., Atwood J. E. (2002). Exercise capacity and mortality among men referred for exercise testing. N. Engl. J. Med. 346 (11), 793–801. doi:10.1056/nejmoa011858
Neder J. A., Alharbi A., Berton D. C., Alencar M. C. N., Arbex F. F., Hirai D. M., et al. (2016). Exercise ventilatory inefficiency adds to lung function in predicting mortality in COPD. COPD J. Chronic Obstr. Pulm. Dis. 13 (4), 416–424. doi:10.3109/15412555.2016.1158801
Neder J. A., Nery L. E., Castelo A., Andreoni S., Lerario M. C., Sachs A., et al. (1999). Prediction of metabolic and cardiopulmonary responses to maximum cycle ergometry: a randomised study. Eur. Respir. J. 14 (6), 1304–1313. doi:10.1183/09031936.99.14613049
O'Donnell D. E., Elbehairy A. F., Webb K. A., Neder J. A. (2017). The link between reduced inspiratory capacity and exercise intolerance in chronic obstructive pulmonary disease. Ann. Am. Thorac. Soc. 14 (Suppl_1), S30–S39. doi:10.1513/AnnalsATS.201610-834FR
Obokata M., Olson T. P., Reddy Y. N. V., Melenovsky V., Kane G. C., Borlaug B. A. (2018). Haemodynamics, dyspnoea, and pulmonary reserve in heart failure with preserved ejection fraction. Eur. heart J. 39 (30), 2810–2821. doi:10.1093/eurheartj/ehy268
Petersen M. S., Kristiansen M. F., Hanusson K. D., Danielsen M. E., Á Steig B., Gaini S., et al. (2020). Long COVID in the Faroe Islands - a longitudinal study among non-hospitalized patients. Clin. Infect. Dis. 73 (11), e4058–e4063. doi:10.1093/cid/ciaa1792
Phillips D. B., Brotto A. R., Ross B. A., Bryan T. L., Wong E. Y. L., Meah V. L., et al. (2021). Inhaled nitric oxide improves ventilatory efficiency and exercise capacity in patients with mild COPD: a randomized‐control cross‐over trial. J. Physiol. 599 (5), 1665–1683. doi:10.1113/jp280913
Ponikowski P., Chua T. P., Anker S. D., Francis D. P., Doehner W., Banasiak W., et al. (2001). Peripheral chemoreceptor hypersensitivity: an ominous sign in patients with chronic heart failure. Circulation 104 (5), 544–549. doi:10.1161/hc3101.093699
Radtke T., Crook S., Kaltsakas G., Louvaris Z., Berton D., Urquhart D. S., et al. (2019). ERS statement on standardisation of cardiopulmonary exercise testing in chronic lung diseases. Eur. Respir. Rev. 28 (154), 180101. doi:10.1183/16000617.0101-2018
Shouman K., Vanichkachorn G., Cheshire W. P., Suarez M. D., Shelly S., Lamotte G. J., et al. (2021). Autonomic dysfunction following COVID-19 infection: an early experience. Clin. Auton. Res. 31 (3), 385–394. doi:10.1007/s10286-021-00803-8
Stickland M. K., Butcher S. J., Marciniuk D. D., Bhutani M. (2012). Assessing exercise limitation using cardiopulmonary exercise testing. Pulm. Med. 2012, 824091. doi:10.1155/2012/824091
Stickland M. K., Neder J. A., Guenette J. A., O'Donnell D. E., Jensen D. (2022). How we do it - using cardiopulmonary exercise testing to understand dyspnea and exercise intolerance in respiratory disease. Chest 161 (6), 1505–1516. doi:10.1016/j.chest.2022.01.021
Stute N. L., Stickford J. L., Province V. M., Augenreich M. A., Ratchford S. M., Stickford A. S. L. (2021). COVID‐19 is getting on our nerves: sympathetic neural activity and haemodynamics in young adults recovering from SARS‐CoV‐2. J. Physiol. 599 (18), 4269–4285. doi:10.1113/jp281888
Sudre C. H., Murray B., Varsavsky T., Graham M. S., Penfold R. S., Bowyer R. C., et al. (2021). Attributes and predictors of long COVID. Nat. Med. 27 (4), 626–631. doi:10.1038/s41591-021-01292-y
Szekely Y., Lichter Y., Taieb P., Banai A., Hochstadt A., Merdler I., et al. (2020). Spectrum of cardiac manifestations in COVID-19. Circulation 142 (4), 342–353. doi:10.1161/circulationaha.120.047971
Tanaka H., Monahan K. D., Seals D. R. (2001). Age-predicted maximal heart rate revisited. J. Am. Coll. Cardiol. 37 (1), 153–156. doi:10.1016/s0735-1097(00)01054-8
Tedjasaputra V., van Diepen S., Collins S. É., Michaelchuk W. M., Stickland M. K. (2017). Assessment of pulmonary capillary blood volume, membrane diffusing capacity, and intrapulmonary arteriovenous anastomoses during exercise. J. Vis. Exp. 120, 54949. doi:10.3791/54949
Tedjasaputra V., van Diepen S., Phillips D. B., Wong E. Y. L., Bhutani M., Michaelchuk W. W., et al. (2018). Pulmonary capillary blood volume response to exercise is diminished in mild chronic obstructive pulmonary disease. Respir. Med. 145, 57–65. doi:10.1016/j.rmed.2018.10.015
Torres-Castro R., Vasconcello-Castillo L., Alsina-Restoy X., Solis-Navarro L., Burgos F., Puppo H., et al. (2021). Respiratory function in patients post-infection by COVID-19: a systematic review and meta-analysis. Pulmonology 27 (4), 328–337. doi:10.1016/j.pulmoe.2020.10.013
Wensel R., Francis D. P., Meyer F. J., Opitz C. F., Bruch L., Halank M., et al. (2013). Incremental prognostic value of cardiopulmonary exercise testing and resting haemodynamics in pulmonary arterial hypertension. Int. J. Cardiol. 167 (4), 1193–1198. doi:10.1016/j.ijcard.2012.03.135
World Health Organization (2021). A clinical case definition of post COVID-19 condition by a delphi consensus. World Health Organization. Available at: https://www.who.int/publications/i/item/WHO-2019-nCoV-Post_COVID-19_condition-Clinical_case_definition-2021.1 (Accessed August 30, 2021).
Keywords: DLCO, long-covid, pulmonary function, pulmonary vascular, VO2, maximal oxygen consumption
Citation: Beaudry RI, Brotto AR, Varughese RA, de Waal S, Fuhr DP, Damant RW, Ferrara G, Lam GY, Smith MP and Stickland MK (2022) Persistent dyspnea after COVID-19 is not related to cardiopulmonary impairment; a cross-sectional study of persistently dyspneic COVID-19, non-dyspneic COVID-19 and controls. Front. Physiol. 13:917886. doi: 10.3389/fphys.2022.917886
Received: 11 April 2022; Accepted: 16 June 2022;
Published: 06 July 2022.
Edited by:
Yuanming Luo, First Affiliated Hospital of Guangzhou Medical University, ChinaReviewed by:
Matiram Pun, University of Calgary, CanadaDanilo C. Berton, Federal University of Rio Grande do Sul, Brazil
Copyright © 2022 Beaudry, Brotto, Varughese, de Waal, Fuhr, Damant, Ferrara, Lam, Smith and Stickland. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Michael K. Stickland, TWljaGFlbC5zdGlja2xhbmRAdWFsYmVydGEuY2E=