Carotid body dysregulation contributes to the enigma of long COVID

Participants

Supplementary Figure 1 shows recruitment, screening, cases excluded and the final sample size. Sixty-four individuals contacted the group about the study, and all were sent a participant information sheet. Of these, 32 replied and completed phone screening. Six were excluded due to screen failure (met exclusion criteria). Twenty-six participants were recruited (long COVID n=16 and controls n=10). Two participants (1 long COVID and 1 control) did not complete the study since they could not tolerate the mask for spirometry, 1 long COVID participant was excluded during the study visit due to identification of cardiac disease (valve disease or non-benign arrythmia) and 1 participant was excluded between visits due to diagnosis of gout (control). The final sample size was 14 participants with long COVID and 8 controls (recovered from initial viral infection within 4 weeks without ongoing symptoms). Data sets from six healthy participants from a previous study (NHS REC numbers: 17/SW/0171 and 18/SW/0241), completed before the SARS-CoV-2 pandemic using the same methods, equipment, and location (clinical room) were added to the control group. These participants were identified based on age, sex, and body mass index (BMI) so that they matched the participants in the long COVID study.

Participant demographics are shown in Table 1. Age, BMI, height, and body mass were similar between the control and long COVID groups (P>0.05). Clinic systolic blood pressure (SBP), diastolic BP (DBP) and heart rate (HR), were not different between groups (P>0.05).

All participants who had COVID-19 were first infected from April 2020 to December 2021, when the dominant variants of SARS-CoV-2 were the original virus, Alpha, Delta, and Omicron. The average number of days since first infection to the study visit, in the long COVID and control groups, were similar (Table 1, P=0.1006). As per the exclusion criteria, none of the participants were hospitalised (admitted) for COVID-19 and had a previous positive PCR test at the time of the initial infection or were seropositive for SARS-CoV-2 antibodies (see inclusion/exclusion criteria in methods). The initial infection caused mild to moderate symptoms in all participants. Table 2 outlines the most common ongoing symptoms previously reported^9,31,32 and the percentage of the long COVID participants in this study reporting these symptoms. Long COVID participants reported at least 3 symptoms, with the most common symptoms being dyspnoea at rest and on exertion, extreme fatigue, “brain fog” and chest pain.

Ten out of the 14 long COVID participants had been seen in a cardiology outpatient long COVID clinic at a UK based University Hospital NHS Trust, where structural cardiac disease (e.g., myocarditis, heart failure, ischaemia) had been excluded. All the participants reporting chest pain were reviewed in a cardiology clinic, which was diagnosed as non-cardiac in origin. None of the long COVID participants (or controls) had been diagnosed with pre-existing cardiovascular or pulmonary diseases. Four of the long COVID participants had been prescribed ivabradine in the cardiology outpatient long COVID clinic. These participants stopped taking ivabradine 48 hours before their study visits (plasma half-life: 2 hours and effective half-life: 11 hours).

Severe asthma was an exclusion criterion, however 2 controls and 6 long COVID participants had been previously diagnosed with mild asthma (before COVID). All participants had a normal resting 12-lead ECG, which was completed as part of the study screening process. Twelve of the 14 long COVID participants self-reported having a chest x-ray, all of which were normal. Finally, Supplementary Table 1 provides information of drugs prescribed in both groups.

Resting spirometry

Spirometry data are summarized in Table 1 and are presented as z-scores with a normal z-score being between ± 1.64³³. There were no differences between the two groups for FEV₁ (forced expiratory volume in 1 second), FVC (forced vital capacity), FEV₁/FVC and FEF_75% (forced expiratory flow after 75% of the FVC has been exhaled; P>0.05). No participants had an FEF_75% smaller than −1.64 z-score, therefore excluding small airway disease. One control subject and 3 long COVID participants had a reduced FEV₁ and FVC with a normal FEV₁/FVC ratio, suggestive of a restrictive ventilatory defect. One long COVID participant had a reduced FEV₁ and FVC, with an FEV₁/FVC ratio consistent with an obstructive ventilatory defect (z-score = −3.25) indicating severe airflow obstruction (z-score = −3.85)³³.

Sit to standing blood pressure test

To check for orthostatic intolerance in both groups, a sit-to-standing test (like that used in the HYVET study^34,35) was completed (3 mins of standing). Sit-to-stand was completed in all long COVID participants and 8 of the controls. BP was measured at seated rest, immediately after standing, 1 min, 2 mins and 3 mins of standing.

According to HYVET study guidelines, a fall in SBP>15 and DBP>7 mmHg was used as a diagnosis of orthostatic intolerance. During standing, the mean SBP and DBP did not fall by more than 15 and 7 mmHg in either group at any time point ³⁴ (Supplementary Figure 2; Supplementary Table 2 for mean change data). Two control participants had a fall in SBP >15 or DBP >7 mmHg, and 1 participant from the long COVID group met this criterion (supplementary figure 2). Mean data at each timepoint were analysed using a mixed model ANOVA. Interestingly, we found that the increase in HR was larger in the control group versus the long COVID group (main Time*Group effect; P=0.0258, see supplementary results for ANOVA details and figure 2).

Resting cardiopulmonary data

The resting HR, BP, minute ventilation, tidal volume and breathing frequency were similar between groups (Table 3, P>0.05). Resting partial pressure of the end tidal CO₂ (P_ETCO₂) was lower in the long COVID group (range: 22 to 34 mmHg versus controls (range: 30-36 mmHg, Table 3, P=0.0378). The minute ventilation was similar between groups despite a lower P_ETCO₂ in the long COVID group. This was coupled with a higher ventilation for the volume of CO₂ expired (V_E/VCO₂ ratio, P=0.0246); suggesting that there is some hyperventilation occurring at rest in the long COVID group. There were no differences in the expiratory and inspiratory times, or the ratio of inspiration to expiratory times between groups (P>0.05; Table 3 for exact P-values). Interestingly, the tidal volume/inspiratory time ratio (V_T/Ti; an index of inspiratory flow³⁶) was 16% higher in the long COVID group versus the controls (P=0.0483) indicative of increased inspiratory drive.

Cardiopulmonary exercise testing

Cardiopulmonary exercise testing (CPET) was completed on a cycle ergometer to peak oxygen consumption (VO₂ peak) to assess exercise tolerance and identify any breathing or cardiovascular abnormalities that might not be evident at rest. We also wanted to understand whether any dysfunctional breathing (defined here as inappropriate hyperventilation or poor breathing efficiency characterised by elevated minute ventilation/volume of CO₂ expired (V_E/VCO₂) slopes¹⁷) during exercise could be linked to carotid chemoreflex hyperactivity in the long COVID group versus control.

To assess whether maximal effort was achieved and as an objective assessment of the quality of the test, a respiratory exchange ratio>1.15, maximum predicted HR>85%, rating of perceived exertion (RPE; 6-20 Borg Scale) of 17-20 and a plateau in VO₂ were used. Overall, 86% of the control group and 86% of the long COVID group achieved these criteria. Thus, differences between groups are likely not due to differences in effort or quality of the CPET. In fact, the control group reported a lower RPE at peak exercise (17; 17-18 (median, IQR)) versus the long COVID group (19; 18-19, P=0.0081; Mann Whitney U). Table 4 shows the mean ± SD (range) for all CPET variables. None of the participants desaturated during the test (defined as SpO₂%<95%, where peak SpO₂% was similar between groups), and there were no exercise induced cardiac ischaemic changes observed on the 12-lead ECG.

The VO₂ peak measured in L/min, mL/kg/min and as a % of peak predicted were lower in the long COVID group versus controls (Table 4; P<0.05). Absolute VO₂ at anaerobic threshold (AT) was also lower (P=0.0025), but the % of peak VO₂ at which the AT occurred was similar between groups (P=0.6490), indicating normal metabolic function. Along these lines there were no other indications of defects in muscle metabolic pathways indicating no muscle metabolic limitations to VO₂ peak in either group (i.e. normal peak respiratory exchange ratio (RER) and a normal VO₂/work rate relationship^37,38 (Table 4). Maximum HR was lower in the long COVID group versus controls, but there was no difference between groups when HR was measured as the % of predicted maximum (P=0.1326). The slope of the HR plotted against the %VO₂ peak was lower in the long COVID group (Table 4; P=0.0114) versus controls indicating a blunted HR response to exercise and potentially some degree of chronotropic incompetence. Heart rate at AT was also lower in the long COVID group versus controls (P=0.0453). Finally, the peak oxygen pulse (VO₂/HR) was lower in the long COVID group versus controls, indicating a cardiovascular limitation to exercise. Despite this, the HR recovery at 1 min was similar between groups (P=0.4653); indicating that parasympathetic engagement post-exercise in the long COVID group was similar to control.

The VCO₂ at peak was lower in the long COVID group (Table 4; P=0.0477), whereas minute ventilation, tidal volume and breathing frequency at peak were similar between groups (Table 4, P>0.05). Supplementary figure 3 shows the minute ventilation versus the VCO₂ and VO₂ at rest, AT, and peak exercise in both groups. The mean breathing reserve (max minute ventilation/predicted maximal voluntary ventilation (calculated as FEV1 * 40) was similar between groups (P=0.8708). Thus it is likely that the lower peak VO₂ in the long COVID group is not a result of pulmonary mechanical limitations^39,40.

The V_E/VCO₂ ratio at peak (ANOVA; P=0.0051) and at AT (P=0.0477) was higher in the long COVID group versus control (see Supplementary figure 3 for ANOVA details and V_E/VCO₂ ratio plotted against the VO₂ at rest, AT, and peak exercise). The V_E/VCO₂ slope (Table 4, P=0.0008) and the V_E/VCO₂ nadir (P=0.0020) during exercise were also higher in the long COVID group versus controls indicating lower breathing efficiency (a higher minute ventilation to remove a given volume of CO₂) indicating hyperventilation at any point during exercise (Figure 1 for example slopes (raw data) in a control and long COVID participant). Thirteen percent of the control group had a V_E/VCO₂ slope higher than the normal range (20-30^41,42) versus 88% in the long COVID group. Finally, the P_ETCO₂ was lower at the AT and at peak exercise in the long COVID group versus the control group (time*group effect; P=0.0119, mixed model ANOVA), but the magnitude of rise of the P_ETCO₂ from rest to AT and from rest to peak exercise was the same in both groups (P=0.0904, supplementary figure 4, mixed model ANOVA).

• Download figure
• Open in new tab

Figure 1.

A) Example carotid chemoreflex test in a control (solid line) and long COVID (dotted line) participants. Data shows the raw SpO₂ response to nitrogen exposure (bottom panel) and the resultant tidal volume, breathing frequency, and minute ventilation responses. The minute ventilation (average of the largest 2 consecutive breaths) is plotted against the nadir SpO₂ for each nitrogen exposure²¹ (panel B). The resultant slope of the linear regression is the hypoxic ventilatory response. Panel C shows the minute ventilation (V_E) plotted against the volume of expired carbon dioxide (VCO₂) for each breath during exercise in the same 2 participants. The slope of the regression is the V_E/VCO₂ slope and is used as a measure of breathing efficiency.

Hypoxic ventilatory response

Figure 1A shows examples carotid chemoreflex tests in a control and a long COVID participants and the resultant hypoxic ventilatory responses (calculated from the minute ventilation) in these participants (Figure 1B). The hypoxic ventilatory response (minute ventilation response to reductions in SpO₂%) was elevated in the long COVID (−0.44 ± 0.23 l/min/ SpO₂%, R²=0.77±0.20) group compared to controls (−0.17 ± 0.13 l/min/SpO2%, R²=0.54±0.38, Figure 2A, P=0.0007); thus, for a given decrease in SpO₂, the participants with long COVID had a greater increase in minute ventilation. This was driven by a greater increase in tidal volume (Figure 2B) in the participants with long COVID rather than a greater increase in breathing frequency compared to the controls (Figure 2C). Taken together this indicates that the participants with long COVD have an elevated carotid chemoreflex sensitivity to hypoxia.

• Download figure
• Open in new tab

Figure 2.

The hypoxic ventilatory response based on A) minute ventilation, B) tidal volume (V_T) and C) breathing frequency (f_B). There was a higher hypoxic ventilatory response based on minute ventilation and tidal volume in the participants with long COVID (n=14) versus control (n=14). Panels D, E and F show the heart rate, systolic and diastolic blood pressure hypoxic responses. Data are the slope of the linear regression where SpO₂% was plotted against the heart rate, and BP during each hypoxic exposure. There was no difference in the HR or BP response to hypoxia between groups.

There was no difference in the HR response to hypoxia (Figure 2D) between groups (control; −0.66 ± 0.35 beats/min/ SpO₂% versus long COVID; −0.58 ± 0.60 beats/min/SpO2%, P=0.5978). Additionally, there was no difference in the SBP or DBP response to hypoxia between groups (Figure 2E and 2F). The BP response to hypoxia was variable between individuals with some participants showing a depressor response to hypoxia and others showing a pressor (most common, where the slope of the regression is negative in figure 2E and 2F) or biphasic response to hypoxia.

The hypoxic ventilatory response was correlated with the V_E/VCO₂ slope during exercise (r=0.54, P=0.0037); indicating that the high carotid chemoreflex sensitivity in the long COVID participants may partially explain the higher V_E/VCO₂ slope (poorer breathing efficiency) during exercise. Therefore, the carotid chemoreflex may play an important role in driving hyperventilation or poor breathing efficiency during exercise. Targeting the carotid body or the carotid chemoreflex could be helpful in improving ongoing symptoms in participants with long COVID and improving exercise tolerance.

The hyperventilatory state of long COVID

In this population of participants with long COVID, we show that despite similar lung function (resting spirometry and ventilatory reserve during exercise) to the control group, these participants hyperventilate at rest and during exercise. At rest, this is evidenced by a similar minute ventilation compared to control despite a lower resting P_ETCO₂. Resting levels of arterial CO₂ provides the key stimulus for respiratory drive, mainly via the central chemoreceptors.⁴³ Reductions in P_ETCO₂ indicate decreases in PaCO₂ ^44,45 which is a stimulus to lower ventilation; however, this has not occurred in the long COVID participants since their level of ventilation was the same as the controls suggesting a resetting of chemoreceptor set-point for breathing. If resting P_ETCO₂ levels were restored to normal control levels, it is possible that minute ventilation would be greater in the long COVID group. Additionally, the V_E/VCO₂ ratio at rest is higher in the long COVID participants, showing that they are breathing more to remove the same volume of CO₂ as controls.

Altered V_E/VCO₂ ratio (and P_ETCO₂) could result from changes to gas exchange between the alveoli and pulmonary circulation. However, we found that during exercise the magnitude of the rise in P_ETCO₂ was the same in the long COVID participants versus the controls. During exercise (with a normal gas exchange) P_ETCO₂ rises and peaks around anaerobic threshold. When gas exchange is impaired by abnormal ventilation or perfusion of the lung, the P_ETCO₂ will show smaller rises, no rise, or even decrease as seen in COPD, heart failure and pulmonary hypertension⁴⁶. As such, the similar rise in P_ETCO₂ with exercise in long COVID patients and control participants suggests no differences in gas exchange between these groups. However, we acknowledge that this needs to be confirmed with gas diffusion studies.

We also found evidence of elevated resting V_T/Ti, an index of inspiratory flow, in the long COVID group versus the control group. The V_T/Ti is used as surrogate of primary neural drive for inspiration³⁶ and thus could indicate a higher ‘inspiratory drive’ at rest in the long COVID participants which may contribute to their feelings of dyspnoea. High inspiratory flow can also indicate changes in airway mechanics⁴⁷, however, this is unlikely in the long COVID group, since resting lung function (spirometry) was similar between groups. During exercise, hyperventilation or poor breathing efficiency was also evidenced by elevated V_E/VCO₂ slopes and an elevated V_E/VCO₂ ratio at anaerobic threshold and peak exercise. Since the carotid chemoreflex plays an important part in the control of breathing both at rest and exercise (including inspiratory drive⁴⁸); it is possible that it is contributing to hyperventilation in long COVID.

Carotid chemoreflex

The mammalian carotid chemoreflex is a protective reflex that contributes to ventilatory and cardiovascular control⁴⁹. In addition to signalling the need for increased ventilation during hypoxia (and other stimuli¹⁹), the carotid body chemoreceptors contribute to resting ventilatory drive, as demonstrated by their tonic activity in humans and animal models^48,50. It is well established that certain disease states exhibit exaggerated carotid chemoreflex sensitivity²⁰ and tonicity¹⁸. Given that elevated carotid chemoreflex sensitivity measured by the hypoxic ventilatory response predicts symptom burden, exercise intolerance, dyspnoea, and hyperventilation in heart failure^21,51, we aimed to assess carotid chemoreflex sensitivity in long COVID. The hypoxic ventilatory response was 145% higher in the long COVID group versus age-, BMI- and sex-matched controls indicating amplified carotid chemoreceptor sensitivity in the long COVID group. The higher hypoxic ventilatory response was observed despite a lower resting P_ETCO₂ in the long COVID group even though hypocapnia normally depresses carotid chemoreflex sensitivity to hypoxia⁵². The hypoxic ventilatory response was in fact similar to that measured in a group of patients with heart failure reduced ejection fraction in our lab (−0.44 ±0.23 L/min/SpO2% versus −0.48± 0.30 L/min/SpO2%; mean ± SD) using the same equipment before March 2020 (see supplementary Figure 5). The elevated hypoxic ventilatory response in the long COVID participants was driven by amplified responses in tidal volume rather than breathing frequency, which supports previous reports in humans where the increase in minute ventilation during mild hypoxia is driven mainly by increases in tidal volume⁵³.

The hypoxic ventilatory response was correlated with the V_E/VCO₂ slope, so that a person with high chemoreflex sensitivity had a higher V_E/VCO₂ slope (or vice versa). This suggests that elevated carotid chemoreflex sensitivity may partially explain reduced breathing efficiency and hyperventilation during exercise in participants with long COVID^14,15,17. In fact, breathing efficiency was poor in some of the participants in this study; 64% of the long COVID group had a V_E/VCO₂ slope >34 (versus only 7% (one person) in the control group), which is a powerful prognostic indicator of future outcomes for people with heart failure⁵⁴. It is possible that in long COVID the carotid chemoreflex has a similar function as that described heart failure, partially explaining the dysregulated breathing and feelings of breathlessness in patients with ongoing symptoms after their initial mild infection. Elevated V_E/VCO₂ slopes are also caused by increased dead space ventilation due to 1) inadequate aeration of the alveoli and 2) poor perfusion of the aerated lung spaces, affecting gas exchange. However, in this study there was no evidence of mechanical lung or small airway issues (due to a similar prevalence of normal breathing reserve in the two groups and similar resting spirometry between the groups). No obvious evidence of gas exchange issues were observed either, because the P_ETCO₂ during exercise (to AT) increased by a similar magnitude in the participants with long COVID versus the controls.

Finally, the resting P_ETCO₂ was <30 mmHg in 4 participants with long COVID. Low P_ETCO₂ causes similar symptoms as those often experienced by patients with long COVID including feelings of brain fog (due to poor cerebral perfusion) and paraesthesia ⁵⁵. The level of P_ETCO₂ however did not reach apnoeic threshold levels ⁵⁶ in any patients because there were no apnoeas at rest and no evidence of periodic breathing patterns⁵⁷. Low resting P_ETCO₂ and hyperventilation is a hallmark of hyperventilatory syndrome which has been cited as a cause of ongoing in symptoms in some cohorts of patients with long COVID⁵⁸. It is possible that problems with breathing at rest and during exercise, leading to lower P_ETCO₂ levels and manifesting as hyperventilation syndrome could be partially driven by the carotid chemoreflex in patients with long COVID.

Mechanisms

Possible mechanisms of increased carotid chemoreceptor sensitivity after SARS-CoV-2 infection include local changes within the parenchyma of the carotid body and/or dysfunction occurring in medullary regions that process afferent sensory information and control efferent ventilatory responses.

In the acute phase of the infection, local viral and immune cell invasion of the carotid body²⁵ could disrupt normal functioning, potentially via immune cells destroying infected glomus cells explaining silent hypoxia in some patients²³. In the long-term these cells are likely replaced but there is the possibility of long-term local inflammation and disruption of ACE1 and ACE2 balance leading to elevated carotid chemoreceptor drive. Additionally, since the carotid bodies are highly sensitive to disturbances in perfusion, any blood flow disruption caused by microthrombi and endothelial dysfunction²⁵ could elevate carotid body activity⁵⁹. It is also possible that dysfunction of the petrosal ganglion, which carries afferent sensory input into the brainstem, causes hyper-reactivity in the carotid chemoreflex. Of note, the petrosal ganglion participates in mediating taste signals in the brain. However, only 7% of our long COVID participants reported ongoing impaired sense of taste and smell ¹¹.

Additionally, inflammation or neuronal damage in the medullary regions mediating afferent signals from the carotid body could contribute to increased chemoreceptor drive. Along these lines there is evidence from animal models that the original SARS-CoV infiltrates medullary brain regions⁶⁰, but there is no strong evidence of this occurring with existing variants of SARS-CoV-2.

Finally, some evidence suggests the SARS-CoV-2 may affect mitochondrial function⁶¹. The glomus cells in the carotid body have unique mitochondria, which are important for O₂ sensing⁶². Any mitochondrial dysfunction in glomus cells could lead to augmented chemoreflex activity⁶². Since mitochondrial dysfunction could occur in any organ, the exercise intolerance and hyperventilation observed in long COVID could also be driven, for example, by the metaboreceptors in skeletal muscle due to poor mitochondrial function and this needs to be examined.

Limitations

Firstly, the control group are a combination of prospectively recruited participants who recovered from a SARS-CoV-2 infection in <4 weeks, and participants who were retrospectively identified from our previous studies (all measures were taken using the same equipment, procedures, and site). Thus 6 of the controls did not have a previous SARS-CoV-2 infection. We cannot rule out that SARS-CoV-2 infection has long lasting effects of the cardiopulmonary system even if individuals recover with no ongoing symptoms. Secondly, we measured the carotid chemoreflex sensitivity using poikilocapnic hypoxia, thus P_ETCO₂ decreased after each hypoxic exposure (due to ventilatory adjustments). However, it is unlikely that this affected the ventilatory response to hypoxia because 1) the P_ETCO₂ decreased only after the breaths that were used as measurements for ventilation and thus could not affect the data and 2) we waited for the P_ETCO₂ to return to baseline after each hypoxic exposure. Thirdly, this was a single site study with a small sample size and thus needs to be evaluated in a larger population of patients with long COVID. Finally, although we see evidence of elevated carotid chemoreceptor sensitivity in patients with long COVID, we need to evaluate whether dampening down the hyperreflexia helps to improve symptoms.

Implications

Here we show for the first time that patients with long COVID have elevated carotid chemoreceptor sensitivity and that this is correlated with poor breathing efficiency or hyperventilation during exercise. Interventions that temper carotid body excitability could be explored as a treatment option for long COVID. Previously our group had shown that P2X3 receptors in the carotid body can be targeted to reduce carotid chemoreflex hyperreflexia in an animal model of hypertension¹⁸ and heart failure⁶³, and could be a viable target in humans with long COVID. Gefapixant, an oral P2X3 receptor antagonist, has recently demonstrated efficacy and an acceptable safety profile in chronic cough in phase 3 clinical trials⁶⁴. P2X3 receptors could therefore be a viable target in humans with long COVID.

Design

This was a single-site case-control study.

Participants

Ethical approval for the study was granted by South Central Hampshire NHS Research Ethics Committee (21/SC/0260) and the Health Research Authority. Participants gave their written informed consent. All participants were asked to abstain from intense exercise and alcohol consumption 24 hours before the study. All experimental protocols conformed to the Declaration of Helsinki. Inclusion criteria for all participants were; aged 18-80 years and a positive SARS-CoV-2 antibody test before vaccination, or a positive COVID-19 PCR antigen swab test. Long COVID participants had received a diagnosis of long COVID, where symptoms developed during or after an infection consistent with COVID-19 and continued for more than 12 weeks (and could not be explained by an alternative diagnosis) as per NHS (UK) National Institute for Health and Care Excellence guidelines⁶⁵.

Participants without long COVID had symptoms lasting less than 4 weeks after their initial infection. See the online-only Supplement for exclusion criteria.

Experimental protocol

Participants attended the NIHR Bristol Clinical Research facility for 2 studies, completed at the same time of day, and the laboratory conditions were at a set temperature (22°C). In visit one, informed consent, in-depth medical history, clinic blood pressure assessment, sit-to-stand test for orthostatic intolerance, lung function tests (spirometry), 12-lead ECG, pregnancy tests and a symptom limited incremental cardiopulmonary exercise test were completed. The second visit involved resting ventilation and cardiovascular measurements followed by carotid chemoreflex assessment via the hypoxic ventilatory response.

Procedures

Clinic blood pressure: Participants rested in a chair for 10 minutes prior to clinic BP being assessed (Omron, 705IT, Omron Healthcare, Kyoto, Japan). Clinic BP was assessed in-line with European Hypertension Society Guidelines⁶⁶

Orthostatic intolerance was assessed using a sit-to-stand test using the HYVET protocol³⁴ where BP was measured when sitting, immediately upon standing, and after 1, 2 and 3 minutes of standing.

Resting spirometry: Resting spirometry was used to assess lung function, to ensure no mechanical lung function abnormalities were present. Spirometry (Ergostick, CPET system, LoveMedical, UK) was completed in line with the joint American Thoracic and European Respiratory Society guidelines ⁶⁷. The Global Lung Function Initiative (GLI) network reference values were used to calculate the percentage of predicted values and z-scores⁶⁸.

12-lead ECG: Resting 12-lead ECG was performed and checked by a Cardiologist at the Bristol Heart Institute for any ECG abnormalities, and to clear participants to exercise.

Cardiopulmonary exercise testing

After acclimatisation to the facemask and sitting on the cycle ergometer, participants completed a 5-minute steady state resting period followed by 3 minutes of unloaded cycling. Participants then completed a continuous ramp incremental exercise test to volitional exhaustion where work rate increased by 15-30 W depending on their physical ability.

Exercise tests were completed on an electronically braked cycle ergometer (Ergoselect 5; Ergoline, Germany). Cardiorespiratory data were recorded using a metabolic measurement system (Ergostik; LoveMedical, UK) with integrated 12-lead ECG and finger pulse oximetry for heart rate and SpO₂% monitoring. Brachial arterial blood pressure was measured via an integrated automated auscultatory blood pressure cuff (LoveMedical, UK. Ratings of perceived exertion (6-20 Borg scale) and dyspnoea scores (modified Borg scale) were obtained at rest, and every minute during exercise, and at the end of exercise. Peak cardiopulmonary data were averaged over the last 30 seconds of exercise. The anaerobic threshold was measured via the V-slope method and Dual Criterion methods as recommended by the American Thoracic Society guidelines ⁶⁹. The minute ventilation and VCO₂ values from initiation to peak exercise were used to measure V_E/VCO₂ slope via least squares linear regression⁷⁰.

Carotid chemoreflex assessment

Resting carotid chemoreceptor sensitivity tests were completed with participants in a semi-supine position with continuous monitoring of beat-to-beat blood pressure (Finapres), SpO₂(ear-lobe pulse oximeter; Radical-7; Masimo Corp, USA), heart rate (lead II of 3-lead ECG, AD Instruments, Australia). Simultaneously, ventilation was measured via a facemask attached to a one-way non-rebreathing circuit (Hans Rudolph, Inc., USA). The inhalation part of the circuit delivered room air or 100% nitrogen gas (transiently) for carotid chemosensitivity testing. The exhalation arm of the circuit was connected to a gas analyser (Ad Instruments) and flow head (MLT3000L; AD Instruments) fitted with a differential pressure transducer (FE141 Spirometer; AD Instruments) for the measurement of inspired and expired fractions of O₂ and CO₂, tidal volume, breathing frequency and minute ventilation. All data were continuously monitored and recorded with a data acquisition system (Powerlab 16/30; AD Instruments) and stored for subsequent analyses using associated software (LabChart 8.0 Pro; AD Instruments).

The transient hypoxic ventilatory response test was used to measure the sensitivity of the carotid chemoreflex ^21,71. After a period of quiet rest breathing room air (10 mins baseline), the researcher added extra-nitrogen to the room air being delivered to the face-mask.

Nitrogen gas administration was controlled silently using a high-pressure electric valve. The nitrogen blended into the room air was delivered for 2-8 breaths, followed by a 3-minute recovery period or until ventilation and haemodynamic variables returned to baseline levels. This was repeated 6-8 times to obtain a range of oxygen saturations (SpO₂: ∼70–100%).

The average of the two largest consecutive breaths in the 1 minute proceeding the nitrogen exposure was used to calculate the ventilatory response to reductions in SpO2%²⁰. The hypoxic ventilatory response was evaluated as the slope of the linear regression relating the minute ventilation to the nadir of SpO₂ for each nitrogen exposure ^20,71 (Figure 1B). The response of tidal volume and breathing frequency to reductions in SpO₂ were also evaluated in the same way as that for minute ventilation.

The peak heart rate and blood pressure was determined following each hypoxic challenge using a 3-beat rolling average and plotted against the nadir oxygen saturation. The hypoxic heart rate (beats/min/%) and blood pressure (mmHg/%) were calculated as the slope of the simple linear regression obtained from baseline and the hypoxic challenges.

Statistical analysis

Statistical analyses were completed in GraphPad Prism (V9.5.1). Participant demographics, resting spirometry, cardiopulmonary exercise and hypoxic ventilatory response data were analysed using an independent samples t-test or Mann Whitney U test if data were not normally distributed or were non-parametric. Where data are compared across multiple time points between the groups, a mixed model ANOVA was used with a Bonferroni correction for pairwise comparisons. Data are reported as mean ± standard deviation or median (interquartile range) unless otherwise stated. A significance level of α<0.05 was used for all analyses.