Impact Factor 3.367 | CiteScore 4.3
More on impact ›

Original Research ARTICLE

Front. Physiol., 22 August 2019 | https://doi.org/10.3389/fphys.2019.01075

The Magnitude of Diving Bradycardia During Apnea at Low-Altitude Reveals Tolerance to High Altitude Hypoxia

Pontus Holmström1*, Eric Mulder1, Angelica Lodin Sundström2, Prakash Limbu3 and Erika Schagatay1
  • 1Department of Health Sciences, Mid Sweden University, Sundsvall, Sweden
  • 2Department of Nursing Sciences, Mid Sweden University, Sundsvall, Sweden
  • 3Department of Clinical Physiology, Nepalese Army Institute of Health Sciences, Kathmandu, Nepal

Acute mountain sickness (AMS) is a potentially life-threatening illness that may develop during exposure to hypoxia at high altitude (HA). Susceptibility to AMS is highly individual, and the ability to predict it is limited. Apneic diving also induces hypoxia, and we aimed to investigate whether protective physiological responses, i.e., the cardiovascular diving response and spleen contraction, induced during apnea at low-altitude could predict individual susceptibility to AMS. Eighteen participants (eight females) performed three static apneas in air, the first at a fixed limit of 60 s (A1) and two of maximal duration (A2–A3), spaced by 2 min, while SaO2, heart rate (HR) and spleen volume were measured continuously. Tests were conducted in Kathmandu (1470 m) before a 14 day trek to mount Everest Base Camp (5360 m). During the trek, participants reported AMS symptoms daily using the Lake Louise Questionnaire (LLQ). The apnea-induced HR-reduction (diving bradycardia) was negatively correlated with the accumulated LLQ score in A1 (rs = −0.628, p = 0.005) and A3 (rs = −0.488, p = 0.040) and positively correlated with SaO2 at 4410 m (A1: r = 0.655, p = 0.003; A2: r = 0.471, p = 0.049; A3: r = 0.635, p = 0.005). Baseline spleen volume correlated negatively with LLQ score (rs = −0.479, p = 0.044), but no correlation was found between apnea-induced spleen volume reduction with LLQ score (rs = 0.350, p = 0.155). The association between the diving bradycardia and spleen size with AMS symptoms suggests links between physiological responses to HA and apnea. Measuring individual responses to apnea at sea-level could provide means to predict AMS susceptibility prior to ascent.

Introduction

Recreational trekking has become a popular sporting activity, and approximately 40 million people travel to high altitude (HA; ≥3000 m) each year (Neupane and Swenson, 2017). As both the individual rate of acclimatization (Schoene et al., 1984; Gallagher and Hackett, 2004) as well as the individual susceptibility to acute mountain sickness (AMS) varies greatly (Oliver et al., 2012), the prediction of individual susceptibility is imperative to reduce AMS incidence.

Several studies have previously investigated whether individual susceptibility to AMS can be predicted. These studies include a wide range of physiological markers, but with inconsistent findings. For example, Roach et al. (1998) found an association between AMS and SaO2, at 4200 m, wherein participants with SaO2 ≥ 84% were at less risk of developing AMS. In addition, Karinen et al. (2010) compared resting SaO2 and post-exercise SaO2 during ascent to 5300 m with the subsequent development of AMS. At 3500 and 4300 m both resting SaO2 and post-exercise SaO2 were lower in participants who developed AMS at higher altitudes. Conversely, some researchers have failed to find any association between AMS and SaO2 at altitude (Wagner et al., 2012).

Some investigators have found an association between both the hypoxic ventilatory response (HVR) and the hypercapnic ventilatory response (HCVR) with AMS (Nespoulet et al., 2012; Richalet et al., 2012). Nespoulet et al. (2012) measured carbon dioxide (CO2) and oxygen (O2) chemo-sensitivity during 20 min hypoxic eupnea in AMS susceptible individuals, during which both variables were found to be reduced in symptomatic compared to asymptomatic individuals. Others have not found any association between AMS and HVR (Milledge et al., 1991; Hohenhaus et al., 1995). For example, Hohenhaus et al. (1995) reported no difference in HVR and HCVR between AMS susceptible and symptom-free individuals, during isocapnic hypoxia.

Furthermore, short apneic duration at sea-level has been found to be an independent predictor of AMS (Austin and Sleigh, 1995). There was also an association between the gag reflex and AMS symptoms (Austin and Sleigh, 1995), and Hayat et al. (2006) found an association between hyperventilation capacity and development of AMS, indicating that general respiratory control may lead to lower risk of developing AMS. Moreover, Karinen et al. (2012) found that heart rate (HR) variability at altitude was associated with development of AMS, however, others have reported no association (Boos et al., 2017). Therefore, despite all past work, to date AMS susceptibility cannot be accurately predicted at sea-level prior to ascent.

Apneic diving naturally induces hypoxia, which initiates several protective responses. The cardiovascular diving response (CVD) is initiated by apnea which is enhanced by facial cooling and serves as the first line of defense against hypoxia during apnea, which occurs via sympathetic outflow to peripheral blood vessels leading to vasoconstriction and by parasympathetic outflow to the heart inducing bradycardia (Gooden, 1994; Andersson and Schagatay, 1998). The diving bradycardia slows the depletion of O2 during apnea, which extends the apnea duration. This response is more pronounced in professional apneic divers compared to non-divers (Schagatay and Andersson, 1998). The CVD is often quantified by HR-reduction, with reductions typically between 15 and 40% (Alboni et al., 2011). As a second line of defense against hypoxia during apnea, spleen contraction occurs which leads to a transient increase in circulating hematocrit (Hct) and hemoglobin concentration (Hb) and prolonged apneic duration (Schagatay et al., 2001; Bakovic et al., 2003). Expert apneic divers have been found to have exceptionally large spleens (Schagatay et al., 2012).

Accordingly, apneic diving and exposure to altitude are two naturally occurring activities in humans, possibly sharing protective mechanisms against an hypoxic insult. A recent study found that individuals with a pronounced diving response during maximal apnea also maintained higher SaO2 at simulated altitude (Starfelt et al., 2018). To evaluate the effects of low-altitude, relative to sea-level, on individual apnea-induced physiological responses, a pilot study was conducted in normobaric hypoxia, which showed that the physiological responses remain consistent between sea-level and simulated low-altitude conditions, without prior acclimatization (Holmström et al., 2018). Based on these promising results, our aim was to investigate if an apneic test that naturally induces hypoxia at sea-level could predict individual susceptibility to AMS during an expedition to HA. We hypothesized that the apnea-induced HR-reduction (diving bradycardia) and spleen contraction would be negatively associated with AMS symptom score in a Lake Louise Questionnaire (LLQ) used during a HA trek.

Materials and Methods

Participants

Twenty-two participants, who were already participating in a trekking expedition to mount Everest base camp (EBC; 5360 m), volunteered to participate in the Study. Included in the final analysis were 18 participants; 8 females and 10 males (mean ± SD age was 44 ± 14 years; height: 174 ± 8 cm; weight: 72 ± 12 kg; and vital capacity [VC]: 4.8 ± 0.8 L). Four participants were excluded for taking Acetazolamide which is used to prevent and reduce symptoms of AMS.

Twelve participants had visited HA (≥3000 m) previously, but none in the previous 2 months. Thus, the participants were recreational un-acclimatized trekkers with no prior experience in apneic diving. Human ethical review boards in Sweden and Nepal approved test protocols and all participants gave written informed consent to participate in accordance with the Declaration of Helsinki.

Procedures

The experimental procedure consisted of two parts: (1) an apnea test conducted in Kathmandu, Nepal, at an altitude of 1400 m and (2) a high altitude trek to EBC.

The participants were asked to report to the lab after at least 1 h without eating or drinking caloric beverages. They entered the lab directly after climbing stairs up to the fourth floor, their height and weight was measured, and they were seated to rest on the chair where the apnea test was conducted and they signed a consent form. Before the apnea test started, after 2 min of seated rest, baseline VC (L) was recorded in duplicate (Vitalograph, Ltd., Compact II, Buckingham England) and the larger volume used for analysis. During the following rest, they received written and oral instructions about the test and after a minimum of 10 min rest, a 2-min countdown for the apnea test started.

Apnea Test

The test consisted of three consecutive static apneas in air, with the first at a fixed time limit of 1 min: apnea one (A1), followed by two maximal voluntary duration apneas: apnea two (A2) and apnea three (A3), all of which were separated by 2 min rest (Figure 1; Engan et al., 2014). While normal breathing preceded the first two apneas, participants were instructed to hyperventilate for 15 s before A3, combined with a maximal voluntary ventilation (MVV) test. Participants were handed a nose-clip at 30 s before apneas and notified when 10 s remained and a second by second countdown was initiated.

FIGURE 1
www.frontiersin.org

Figure 1. Timing of apneas (A1: apnea 1; A2: apnea 2; A3: apnea 3) with spleen measure (doted line), continuous measure of SaO2 (arterial oxygen saturation, black line) and HR (heart rate, gray line) and lung volume measures. 15 s of hyperventilation preceded apnea 3 (A3) combined with a MVV (maximal voluntary ventilation) test.

Participants were instructed to exhale fully and take a deep but not maximal breath preceding each apnea, while inspiratory volume was measured. Participants were told to keep the chest relaxed and refrain from swallowing or exhaling during apneas. If SaO2 levels decreased below 60% during apneas, the participants would be told to resume breathing in order not to risk syncope, which may occur below 50% SaO2. Upon termination of apnea, the nose-clip was removed, and the participants were instructed to resume normal breathing. Apneic duration was monitored with a stopwatch.

Measurements

After 5 min of rest, blood pressure, and forehead temperatures were measured in duplicate to exclude hypertension (200/100) or fever (>38°C). The spirometer were also used for inspiratory and expiratory volumes as well as MVV measurements in the apnea tests. SaO2, HR, partial pressure of end-tidal CO2 and respiratory rate were continuously logged throughout the apnea test, using a combined pulse oximeter-capnograph (Medair Lifesense LS1-9R, Medair AB, Delsbo, Sweden) attached to the middle finger, from two min prior to A1 until 10 min after A3 (Figure 1). Data were stored via memory unit (Trendsense, Nonin Medical, Inc., Medair AB, Hudiksvall, Sweden). Spleen size was measured via ultrasonic imaging (M-Turbo Ultrasound system, FUJIFILM SonoSite, Inc., Bothell, WA, United States) each minute for 10 min during the rest period and immediately after the termination of each apnea, during which measurements were taken each minute from the end of each apnea until 10 min after A3 (Figure 1). During each minute, three maximal spleen measures were recorded (length, thickness, and width).

High Altitude Trek

The morning after the apnea test, the participants flew to Lukla (2840 m) and started the EBC trek. They hiked in three groups, organized by a local trekking agency and followed the same ascent profile to Dingboche (4410 m; Figure 2) reaching Gorak Shep (5160 m) on the 9th day. The trek included two rest days, and arrived back in Kathmandu on the morning of the 14th day. At least one experienced Sherpa guide led each group. During the ascent, the participants documented their AMS symptoms in a booklet, each morning using the LLQ. The participants had been thoroughly informed on how to use the LLQ booklet in conjunction with the apnea test. SaO2 and HR was also measured each morning with a portable pulse oximeter (Nonin Medical, Inc., Medair AB, Delsbo, Sweden), which was supervised by the guides. After the hike, participants returned the LLQ booklet to the research staff.

FIGURE 2
www.frontiersin.org

Figure 2. Sleeping altitudes for data collection of Lake Louise Questionnaire (LLQ) score (days 1–7) included in the analysis during ascent to EBC (n = 18).

Lake Louise Questionnaire

The LLQ is a self-administered questionnaire consisting of five items quantifying the presence of the most frequent AMS symptoms, including: (1) headache, (2) gastrointestinal symptoms, (3) fatigue, (4) dizziness, and (5) difficulty sleeping. Each item is graded on a scale from 0 (no symptom) to 3 (severe symptom) as described by Roach et al. (1993). In addition, medication taken during the trek and morning values of SaO2 and HR were noted in the LLQ booklet.

Analysis

The magnitude of the diving response was interpreted through the diving bradycardia which was quantified by the apnea-induced HR-reduction from baseline. Baseline HR was computed by calculating the mean HR (bpm) from 90 to 30 s prior to A1 and used as a reference value. To calculate the mean apnea-induced HR-reduction, the mean HR across the apnea duration, minus the first 30 s (A-30) was computed, thus eliminating the initial tachycardia and decline phase (Andersson and Schagatay, 1998). Maximal HR-reduction was also calculated by identifying the HRnadir during each apnea. Quantification of the HR-reduction was done by calculating the percentage change between the baseline HR and the mean HR during the apnea (mean HR-reduction) and the HRnadir (maximal HR-reduction; Figure 3).

FIGURE 3
www.frontiersin.org

Figure 3. Recording of an individual participants diving bradycardia during 1 min apnea. Note the initial tachycardia, followed by a decline phase. The diving response was defined as the percentage reduction in HR (heart rate) during the period from 30 s into the apnea until the end of the apnea, marked as Mean HR, which was compared to the period 90–30 s before the apnea (baseline).

Baseline SaO2 was computed by calculating the mean SaO2 from 90 to 30 s prior to A1. To account for the circulatory delay from the lung to the finger, the arterial O2 desaturation was calculated as the percentage change between baseline SaO2 and the SaO2nadir from initiation of apnea until 30 s after apnea.

Measurements of the maximal splenic length (L), width (W), and thickness (T) were used to calculate spleen volume according to the Pilström equation: Vspleen = Lπ(WT-T2)/3.

The formula describes the difference between two ellipsoids divided by two, based on the observed average shape of the spleen (Schagatay et al., 2005; Lodin-Sundström and Schagatay, 2010). The baseline (resting) spleen volume was obtained by averaging the two consecutive maximal values from the 10-min period before the first apnea and changes that occurred during apneas and recovery were calculated by comparing resting volumes with those immediately after apnea and throughout the recovery.

As all participants followed the same ascent profile to 4410 m, LLQ scores could be summed up to this altitude and used in the analysis. Quantification of AMS symptoms throughout the ascent via the accumulated LLQ score minimizes the influence of the many different conditions that mimic AMS, e.g., exhaustion, dehydration, infection, hangover, migraine and hypoglycemia, which all may lead to some degree of headache (Luks et al., 2017). The primary symptom included in the diagnostic criteria of AMS is headache (Roach et al., 1993), which can potentially lead to an incorrect diagnosis. Furthermore, AMS susceptibility was also evaluated by dividing the participants into subgroups, using median split based on individual accumulated LLQ scores: low LLQ group (Lo; LLQ ≤ 6p); medium LLQ group (Me; LLQ 7–12p) and high LLQ group (Hi; LLQ ≥ 13p).

Statistical Analysis

Data are reported as mean ± SD, except for LLQ data, which are reported as median (range). Normality distribution of the data was assessed using Kolmogorov–Smirnov and Shapiro–Wilk tests (p > 0.05). Differences in HR-reduction, resting spleen volume and contraction, O2 desaturation during each apnea, VC, resting HR and resting SaO2 (independent variables) between subgroups were evaluated by independent samples t-tests. Within-subject’s comparisons on all variables were performed using paired samples t-tests with Bonferroni corrections for multiple comparisons. Magnitude of observed effects were estimated using the standardized mean difference (Cohen’s d, effect size [ES]), computed as the mean difference divided by the pooled SD. Effect sizes was presented along with 95% confidence intervals (CI). An ES of 0.0–0.3 was considered a small effect, 0.4–0.7 medium and >0.8 a large effect (Lee, 2016). Bivariate multiple correlation tests with Spearman’s rank correlation coefficients (rs) were conducted to analyze correlation between dependent (LLQ score) and independent variables. Pearson correlation coefficient (r) were used to analyze correlation between parametric variables. A binomial logistic regression using a conditional forward stepwise selection method was performed to analyze the prediction of LLQ score (dependent variable) by HR-reduction, resting spleen volume, VC, arterial O2 desaturation during apnea and gender (independent variables) and odds ratio (OR) estimated, with 95% CI. Dichotomization of AMS symptoms in the logistic regression, were established based on median split of LLQ score, into high and low LLQ score (low ≤ 12p; high ≥ 13p). Only one variable: apnea-induced bradycardia during A1, was retained in the regression model. p < 0.05 was considered significant and tendencies were denoted p < 0.1.

Results

LLQ Score and Apnea Duration

Apneic duration for A1 was 56 ± 7 s with seven participants not able complete 60 s. The maximal voluntary apnea durations were 70 ± 16 s (range: 35–97 s) for A2 and 92 ± 25 s (range: 55–151 s) for A3.

While most participants reported only mild AMS symptoms, one subject was evacuated by helicopter from 4410 m due to severe AMS. There was an individual variation in accumulated LLQ (0–22) LLQ points with a median score of 9 at 4410 m. Individual apneic durations did not correlate with LLQ score (NS for all apneas).

Apnea-Induced Bradycardia

The mean apnea-induced bradycardia was, for A1: 11 ± 12% (HRnadir: 70 ± 9 bpm), for A2: 13 ± 14% (HRnadir: 65 ± 9 bpm), and for A3: 10 ± 13% (HRnadir: 68 ± 9 bpm) (p < 0.001). The individual diving bradycardia varied between −34 and +20% and reached baseline resting values during each of the recovery periods (NS).

The mean apnea-induced bradycardia was negatively correlated with the LLQ score (Figures 4A–C). There was also a correlation between LLQ score and the maximal diving bradycardia for A1 (rs = −0.542, p = 0.02) which became insignificant during A2 (rs = −0.328, p = 0.184) and A3 (rs = −0.354, p = 0.149). The diving bradycardia during apnea was also correlated with SaO2 at 4410 m; for A1: r = 0.655 (p = 0.003; Figure 4D), for A2: r = 0.471 (p = 0.049) and for A3: r = 0.635 (p = 0.005). Logistic regression revealed an OR of 1.146 ([CI 1.01–1.30]; p = 0.034) for A1.

FIGURE 4
www.frontiersin.org

Figure 4. The figures depict correlation plots of accumulated LLQ score at 4410 m and magnitude of the HR (heart rate) reduction (diving bradycardia) during Apnea 1 (A), Apnea 2 (B) and Apnea 3 (C). (D) Shows HR-reduction (diving bradycardia) during Apnea 1 and SaO2 at 4410 m.

Spleen Volume and Contraction During Apnea

Resting spleen volume was 191 ± 50 mL (range: 86–261 mL), which was reduced after all apneas, to a mean volume of 160 ± 50 mL (−16%) after A1, 147 ± 36 mL (−23%) after A2 and 149 ± 39 mL (−22%) after A3 (Figure 5). Five minutes after A3, the volume had increased to 180 ± 51 mL (−5%; p = 0.299).

FIGURE 5
www.frontiersin.org

Figure 5. The spleen volume change after three apneic episodes (apnea 1–3) spaced by 2 min of recovery and the volume after 5 min recovery. ∗∗ indicates p < 0.01 from baseline.

Baseline (resting) spleen volume correlated negatively with LLQ score (rs = −0.479, p = 0.044; Figure 6A). Maximal volume reduction across all apneas correlated with resting volume (rs = −0.625, p = 0.006), but the maximal volume reduction did not correlate with LLQ score (rs = −0.350, p = 0.155; Figure 6B). Neither was there any correlation between LLQ score and relative apnea-induced spleen volume reduction (%).

FIGURE 6
www.frontiersin.org

Figure 6. Shows a correlation plot between accumulated LLQ score at 4400 m and baseline (resting) spleen volume (A) and maximal volume reduction induced by apnea (B) in Kathmandu at 1400 m.

Comparing the differences in diving bradycardia and spleen volume between sub-groups according to their LLQ score, the diving bradycardia during A1 was more pronounced in the Lo group compared to the Hi group (p = 0.014, ES: 1.73 [0.3–2.9]; Figure 7A) and resting spleen volume was greater in the Lo LLQ group compared to Hi LLQ group (p = 0.027, ES: 1.59 [0.19–2.75]; Figure 7B).

FIGURE 7
www.frontiersin.org

Figure 7. Shows the magnitude of the heart rate (HR) reduction (diving bradycardia) (A) and resting spleen volume (B) for each of the groups based on accumulated acute mountain sickness (AMS) symptoms score: Low LLQ score, Medium LLQ score, and High LLQ score (p < 0.05).

Arterial Oxygen Saturation During Apneas

Resting SaO2 at low-altitude was 96.3%, which was reduced by 4 ± 1 to 92 ± 2% during A1 (p < 0.0001), by 6 ± 4 to 90 ± 4% during A2 (p = 0.02) and by 9 ± 6 to 87 ± 6% during A3 (p = 0.0014). Arterial O2 desaturation during A1 tended to correlate with LLQ score (rs = −0.443, p = 0.065; Figure 8A). During A2, the correlation became insignificant (rs = −0.352, p = 0.152), with no correlation between O2 desaturation and LLQ score during A3 (rs = −0.121, p = 0.633).

FIGURE 8
www.frontiersin.org

Figure 8. Shows correlation plots of accumulated LLQ score at 4400 m and arterial O2 desaturation during apnea 1 (A), vital capacity (VC; B), and arterial O2 saturation (SaO2) at 4400 m (C).

Baseline Cardiovascular Variables at 4410 m

Resting SaO2 at 4410 m was 79 ± 9% and correlated negatively with LLQ score (rs = −0.757, p < 0.001; Figure 8C). Resting HR at 1400 m was 77 ± 10 and 80 ± 12 bpm at 4410 m (Table 1). The Lo group had higher SaO2 at altitude compared to Hi group (Table 1).

TABLE 1
www.frontiersin.org

Table 1. Baseline and functional characteristics of the AMS subgroups based on LLQ score.

Lung Volume

VC was 4.8 ± 1 L and varied between 3 and 6 L. Individual VC correlated negatively with LLQ score (Figure 8B). VC was also associated with apnea duration (A1: r = 0.478, p = 0.045; A2: r = 0.566, p = 0.014; A3: r = 0.679, p = 0.002). The Lo LLQ group had larger lungs compared to Me and Hi groups (Table 1). The results were similar when lung volumes were expressed related to participants height.

Discussion

In the present study we investigated the association between apnea-induced physiological responses with the development of AMS symptoms during a trek to mount Everest base camp (EBC). The principal findings were that: (I) the apnea-induced bradycardia for A1 and A3 was negatively associated with accumulated LLQ score; (II) the apnea-induced bradycardia was also associated with SaO2 at high altitude (HA); (III) resting spleen volume was negatively associated with accumulated LLQ score with ascent.

Cardiovascular Diving Response

The finding that the apnea-induced bradycardia is associated with symptoms of AMS suggests that a powerful diving response, which is known to conserve O2 and characterizes successful free divers (Schagatay and Andersson, 1998), may also be characteristic of individuals who are tolerant to HA hypoxia. This novel finding provides a promising approach for further studies aiming to develop a simple sea-level test to predict AMS. The association between the diving response and AMS likely result from higher SaO2, leading to maintained O2 availability for the brain and the heart. This was clearly shown by the association between the apnea-induced bradycardia during both 1 min apnea (A1) and maximal apneas (A3) with accumulated LLQ score, showing that an individual with a pronounced diving response develop less AMS symptoms. However, we did not find any association between the second apnea (A2) and accumulated LLQ score, which was unexpected as we observed a similar magnitude in the HR-reduction for A1 and A3 and a similar O2 desaturation during A1. Likely, this may reflect variations in individual HR responses across the apnea series wherein four participants responded with a small HR increase during A2. A noteworthy observation was that resting HR did not correlate with LLQ score; hence, it is the actual apnea-induced bradycardia magnitude that is relevant, with a powerful response reducing the individual risk of developing AMS symptoms.

The metabolic rate decreases during apnea through the CVD, which is strongly associated with apneic diving performance (Schagatay and Andersson, 1998) most likely due to a decrease in O2 consumption (Andersson et al., 2002). We suggest that during work at HA, the diving response may have a similar effect. During apnea, the brain and the heart are typically prioritized organs for oxygenation and this situation would also be beneficial at HA. Our investigation also show a less pronounced arterial O2 desaturation at HA with a concomitant powerful diving response, which is in agreement with a previous study using simulated altitude, reporting that individuals with a pronounced HR-reduction and maintained levels of SaO2 during apnea also maintained high SaO2 during simulated altitude (Starfelt et al., 2018). As such, individuals with a powerful diving response during apnea do not only develop less AMS symptoms, but also show less arterial O2 desaturation at HA, and thus are less hypoxic and more tolerant to HA hypoxia.

Several investigations have found an association between AMS and SaO2 measured at HA (Roach et al., 1998; Karinen et al., 2010; Oliver et al., 2012; Richalet et al., 2012) which is in line with our findings of lower SaO2 in participants with more AMS symptoms. For instance, in a large epidemiological study including 1300 participants, Richalet et al. (2012) reported that both resting SaO2 and arterial O2 desaturation during exercise is greater in participants who develop AMS compared to symptom-free individuals. These findings support the view that SaO2 distinguishes those sensitive to altitude hypoxia and those more tolerant.

When participants were divided by AMS symptoms, we found a similar pattern that the group who developed most symptoms also desaturated most, whereas the group who developed least AMS symptoms desaturated less and hence were less hypoxic. Interestingly, we found the same pattern in the three LLQ groups for the diving bradycardia. The group who developed less AMS symptoms had the most pronounced apnea-induced bradycardia, followed by the intermediate LLQ group and high LLQ group. This further shows the contribution of the diving response to maintain high SaO2 levels, even in conditions of severe poikilocapnic hypoxia. The low SaO2 in AMS susceptible individuals gives credence to an inefficient O2-conserving effect through a relative marginal diving bradycardia in the etiology of AMS.

The association between the diving bradycardia and AMS symptoms suggests that the initiating mechanism of the diving response shares a common pathway between HA and apneic diving. At HA, a high HVR is associated with superior performance (Schoene et al., 1984) and with protective properties against AMS (Nespoulet et al., 2012; Richalet et al., 2012). HVR has also been found to have a protective function during apnea, wherein HVR predicts apneic duration and SaO2nadir (Feiner et al., 1995), and apneic divers have higher HVR during poikilocapnic hypoxia compared to non-apneic divers (Costalat et al., 2014), resulting from higher chemoreceptor responsiveness in the former group. However, during apnea the diving bradycardia is not initiated by hypoxia, but by cessation of respiration and augmented by face immersion via stimulation of cutaneous facial cold receptors (Schutema and Holm, 1988). This was also apparent in the present study, although face immersions were not used, thus the magnitude of the apnea-induced bradycardia was almost similar during A1 compared to A2 and A3, which did not cause a severe drop in SaO2 during A1. However, PaO2 regulates the bradycardia, where decreased PaO2 below normal values will cause a more powerful response (Gooden, 1994), which might imply that individuals with a higher chemoreceptor responsiveness during apnea, eliciting a more powerful diving response, may also have a higher HVR at HA.

A powerful diving response may also reflect a dynamic vascular system, involving a general ability of the body to transiently prioritize blood flow to the most important regions. Endothelial function [flow-mediated dilatation (FMD)] is essential for vasoregulation. At HA, preservation of the vascular system through maintained FMD is crucial for O2 delivery to vital organs. Redistribution of blood flow through peripheral vasoconstriction occurs at HA during exercise, whereby the body down-prioritizes muscle oxygenation (Calbet et al., 2003) most likely in favor of the most hypoxia-sensitive organs, as occurs during diving. In addition, cerebral blood flow increases at HA in response to low PO2, as a means to swiftly shuttle O2 to the brain (Ainslie and Subudhi, 2014). Upon acute exposure to HA, FMD is reduced in those who develop AMS compared to asymptomatic individuals, where reduction in FMD occurs well before the onset of AMS with increased arterial O2 desaturation (Bruno et al., 2016). We speculate that a pronounced diving response contributes to/reflects a dynamic cardiovascular system at HA, which would reduce O2-consumption by prioritizing and redistributing blood flow. This mechanism could have contributed to the higher SaO2 observed at HA in symptom-free individuals in the present study, and thereby be protective against AMS.

Spleen Volume and Contraction

An additional new finding of the present study is that resting spleen volume is negatively associated with AMS symptoms development, i.e., individuals with larger spleens develop less AMS symptoms. Previous investigations have shown that different hypoxic situations stimulate the spleen to contract and release red blood cells (RBCs), both during apnea (Schagatay et al., 2001) and during normobaric poikilocapnic hypoxia (Lodin-Sundström and Schagatay, 2010), resulting in increased circulating RBC which increases the blood gas transportation. The primary driver to induce splenic contraction is hypoxia (Lodin-Sundström and Schagatay, 2010), however, other mechanisms have been proposed, such as hypercapnia (Richardson et al., 2012) and sympathetic activation (Frances et al., 2008), which may be a possible cause for the response which occurred during A1 in the present study, during which SaO2 was still broadly maintained above 92%.

We hypothesized that the function to increase blood gas transportation capacity would also be efficient at HA, in which a greater spleen volume and ability to contract could be beneficial in coping with HA hypoxia, by enhancing O2-carrying capacity. The findings of an association between individual spleen volume and AMS symptoms supports our hypothesis. The observation is in line with earlier findings in free divers, which suggests that spleen volume was largest in the athletes with the best results in a free diving world championship (Schagatay et al., 2012). Also, native Bajau divers, who are continuously exposed to hypoxia, were recently reported to have larger spleens compared to a native population within close proximity but not exposed to hypoxia associated with apneic diving (Ilardo et al., 2018). We also reported larger spleens in the most hypoxic COPD patients, compared to less hypoxic patients (Schagatay et al., 2015). Furthermore, we found that spleen contraction was enhanced after a HA expedition to the summit of mount Everest (Engan et al., 2014). These studies suggest that spleen volume and function could be affected by training/exposure. However, it cannot be determined in the present study if participants tolerant to AMS had larger spleens due to genetics, or as a result of training. None of the participants in our sample had visited HA (≥3000 m) in the 2 months prior to the trek to EBC.

The mechanism by which a large spleen could be protective against AMS is that a higher circulating Hb would enhance blood O2 carrying-capacity, and the larger the spleen, the more RBC can be expelled, which is supported by a strong association between spleen volume and volume contraction in our study. Therefore, an association could also be expected with the maximal volume reduction, however, only a weak tendency was observed in our study. We expected that the participants most sensitive to HA, with the greatest SaO2 reductions at altitude, would have a greater volume contraction. However, there were no association between SaO2 at HA and the magnitude of spleen volume reduction. On the other hand, we found that the group with least AMS symptoms, had both larger resting spleen volume and maximal spleen volume reduction compared to the group who developed most AMS symptoms, indicating the superior contractile function of the bigger spleens, allowing a greater increase in RBC, which is associated with tolerance to HA hypoxia.

We believe this lack of association between spleen volume reduction and AMS symptoms is partly due to the discrepancy between AMS subgroups, i.e., the spleen volume reduction was smaller in the intermediate group compared to the high LLQ group, which was unexpected. Looking at the group with low AMS symptoms, we observe a resting volume of 226 mL, a maximal volume reduction of 77 mL and a pronounced diving response (18% HR-reduction), a clear protective characteristic of HA tolerant individuals resulting in higher SaO2 and subsequently less AMS symptoms. Conversely, the two other groups exhibits smaller and almost similar resting spleen volumes, of 169 and 176 mL, but different volume reduction and levels of SaO2. Presumably, the smaller relative contraction of the intermediate group compared to the Hi group is due to other superior means to cope with hypoxia, e.g., a powerful diving response of 13% HR-reduction, leading to less severe arterial O2 desaturation at HA. Meanwhile, participants who develop more AMS symptoms, exhibit less functional defense systems against hypoxia, e.g., have a smaller magnitude in apnea-induced bradycardia of <3%, leading to more severe arterial O2 desaturation during apnea and at HA. This would thus lead to a stronger/maximal spleen volume contraction, which would thus serve as a second line of defense against hypoxia, and only be triggered maximally if the hypoxia is severe enough. If the spleen were small, however, it would not have as large effect in protecting against hypoxic insult, as the RBC boost into circulation would still be relatively small.

Lung Volume

We also found that VC was associated with AMS symptoms development, which has been observed previously (Hackett et al., 1982), wherein VC measured in Kathmandu was inversely associated with SaO2 at HA and AMS symptoms. This is also in line with previous findings in free divers who reported greater lung volume in elite apneic divers compared to non-apneic divers (Schagatay et al., 2012). A greater lung is a vital part of apneic diving performance, as it leads to greater O2 storage capacity, hence increased O2 availability during apnea. In regard to altitude sensitivity, individuals with relatively smaller lungs are more likely to desaturate more (Hackett et al., 1982), thus be at increased risk of developing AMS. It has also been reported that lung volume characterizes individuals that tolerate HA hypoxia well. Tibetans, for instance, exhibit both greater VC and forced VC compared to lowland Han Chinese (Gilbert-Kawai et al., 2014), a physiological adaptation brought about generations of continued exposure to HA hypoxia.

Study Limitations

The current study design is associated with a number of limitations, which should be acknowledged. Although the strength of the study relates to the experimental procedure and establishing a well-functioning laboratory in the field for advanced physiological testing during apnea. A study limitation is that AMS symptoms are self-reported. The LLQ depends solely on individual subjective perception of the occurrence and severity of symptoms (Roach et al., 1993).

Although all participants were in the same trekking group, hiking together along the same ascent profile, it may still be difficult to standardize an exact ascent rate, water and food intake and additional drug use, e.g., ibuprofen to relieve headache, factors which may have effects on AMS symptoms development.

SaO2 measurements via pulse oximetry in the field may also have its limitations. This may lead to misinterpretations in certain situations with peripheral vasoconstriction due to cold extremities or hypoxia at HA or during apnea, which may influence the data output (Luks and Swenson, 2011). However, efforts were made to minimize these factors causing measurements error. All measurements were conducted indoors; maintaining peripheral skin temperatures, and by using a pulse oximeter with a light signal indicating when blood flow was sufficient to obtain good data.

Furthermore, the low sample size of the current study may also have affected the outcome of the logistic regression model. Therefore, to continue to develop the prediction of AMS susceptibility through this easy-to use and non-invasive apnea test we propose, these results should be further studied in a larger population.

Conclusion

We conclude that the apnea-induced bradycardia and resting spleen volume were inversely associated with AMS symptoms development at HA. These novel findings highlights important links between physiological responses to HA and apneic diving which could contribute to furthering the understanding of the great inter-individual differences in susceptibility to AMS. Measuring individual cardiovascular, splenic and hematological responses induced during apnea at sea-level could possibly be an alternative approach to predict AMS susceptibility in native lowlanders prior to ascent.

Data Availability

The raw data supporting the conclusions of this manuscript will be made available by the authors, without undue reservation, to any qualified researcher.

Author Contributions

PH planned and organized the laboratory and field study tests and procedures, collected and analyzed the data, and wrote the manuscript. EM planned and organized the laboratory and field study tests and procedures, collected the data, and proofread the manuscript. AS organized the field laboratory, collected the data, and proofread the manuscript. PL collected the data and proofread the manuscript. ES conceived the idea, planned and organized the laboratory and field study tests and procedures, collected the data, wrote the manuscript, and proofread the manuscript.

Funding

This project was financially supported by the Swedish National Centre for Research in Sports (CIF) and Mid Sweden University.

Conflict of Interest Statement

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.

Acknowledgments

We thank the participants for their efforts, and the people helping us in Nepal, especially Mr. Ram Sapkota (Mountain Delights Treks & Expeditions), Mr. Nicke Sundström (Pathfinder Travel), Mr. Pradeep Pakhrin (Managing Director of Fuji Hotel), and Dr. Arjun Karki (Director of Patan Academy), for their contributions. We acknowledge the International Society for Mountain Medicine Conference in Kathmandu, Nepal in November 2018, where the part of the data was presented and published as an abstract.

References

Ainslie, P. N., and Subudhi, A. W. (2014). Cerebral blood flow at high altitude. High Alt. Med. Biol. 15, 133–140. doi: 10.1089/ham.2013.1138

PubMed Abstract | CrossRef Full Text | Google Scholar

Alboni, P., Alboni, M., and Gianfranchi, L. (2011). Diving bradycardia: a mechanism of defence against hypoxic damage. J. Cardiovasc. Med. 12, 422–427. doi: 10.2459/JCM.0b013e328344bcdc

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, J., Linér, M., Runow, E., and Schagatay, E. (2002). Diving response and arterial oxygen saturation during apnea and exercise in breath-hold divers. J. Appl. Physiol. 93, 882–886. doi: 10.1152/japplphysiol.00863.2001

PubMed Abstract | CrossRef Full Text | Google Scholar

Andersson, J., and Schagatay, E. (1998). Arterial oxygen desaturation during apnea in humans. Undersea Hyperb. Med. 25, 21–25.

PubMed Abstract | Google Scholar

Austin, D., and Sleigh, J. (1995). Prediction of acute mountain sickness. BMJ 311, 989–990. doi: 10.1136/bmj.311.7011.989

PubMed Abstract | CrossRef Full Text | Google Scholar

Baković, D., Valic, Z., Eterović, D., Vukovic, I., Obad, A., Marinović-Terzić, I., and Dujić, Z. (2003). Spleen volume and blood flow response to repeated breath-hold apneas. J. Appl. Physiol. 95, 1460–1466. doi: 10.1152/japplphysiol.00221.2003

PubMed Abstract | CrossRef Full Text | Google Scholar

Boos, C. J., Vincent, E., Mellor, A., O’hara, J., Newman, C., Cruttenden, R., et al. (2017). The effect of sex on heart rate variability at high altitude. Med. Sci. Sports Exerc. 49, 2562–2569. doi: 10.1249/MSS.0000000000001384

PubMed Abstract | CrossRef Full Text | Google Scholar

Bruno, R. M., Ghiadoni, L., and Pratali, L. (2016). Vascular adaptation to extreme conditions: the role of hypoxia. Artery Res. 14, 15–21. doi: 10.1016/j.artres.2016.02.003

CrossRef Full Text | Google Scholar

Calbet, J. A., Boushel, R., Rådegran, G., Søndergaard, H., Wagner, P. D., and Saltin, B. (2003). Why is Vo2 max after altitude acclimatization still reduced despite normalization of arterial O2 content? Am. J. Physiol. Regul. Integr. Comp. Phys. 284, R304–R316.

PubMed Abstract | Google Scholar

Costalat, G., Pichon, A., Coquart, J., Bauer, F., and Lemaître, F. (2014). Cardio-ventilatory responses to poikilocapnic hypoxia and hypercapnia in trained breath-hold divers. Respir. Physiol. Neurobiol. 192, 48–54. doi: 10.1016/j.resp.2013.12.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Engan, H., Lodin-Sundström, A., Schagatay, F., and Schagatay, E. (2014). The effect of climbing Mount Everest on spleen contraction and increase in hemoglobin concentration during breath-holding and exercise. High Alt. Med. Biol. 15, 52–57. doi: 10.1089/ham.2013.1061

PubMed Abstract | CrossRef Full Text | Google Scholar

Feiner, J. R., Bickler, P. E., and Severinghaus, J. W. (1995). Hypoxic ventilatory response predicts the extent of maximal breath-holds in man. Respir. Physiol. 100, 213–222. doi: 10.1016/0034-5687(94)00132-j

PubMed Abstract | CrossRef Full Text | Google Scholar

Frances, M. F., Dujic, Z., and Shoemaker, J. K. (2008). Splenic constriction during isometric handgrip exercise in humans. Appl. Physiol. Nutr. Metab. 33, 990–996. doi: 10.1139/H08-087

PubMed Abstract | CrossRef Full Text | Google Scholar

Gallagher, A., and Hackett, P. (2004). High-altitude illness. Emerg. Med. Clin. 22, 329–355.

Google Scholar

Gilbert-Kawai, E. T., Milledge, J. S., Grocott, M. P., and Martin, D. S. (2014). King of the mountains: Tibetan and Sherpa physiological adaptations for life at high altitude. Physiology 29, 388–402. doi: 10.1152/physiol.00018.2014

PubMed Abstract | CrossRef Full Text | Google Scholar

Gooden, B. A. (1994). Mechanism of the human diving response. Integr. Physiol. Behav. Sci. 29, 6–16. doi: 10.1007/bf02691277

CrossRef Full Text | Google Scholar

Hackett, P., Rennie, D., Hofmeister, S., Grover, R., Grover, E., and Reeves, J. (1982). Fluid retention and relative hypoventilation in acute mountain sickness. Respiration 43, 321–329. doi: 10.1159/000194501

PubMed Abstract | CrossRef Full Text | Google Scholar

Hayat, A., Hussain, M. M., Aziz, S., Siddiqui, A. H., and Hussain, T. (2006). Hyperventilatory capacity–a predictor of altitude sickness. J. Ayub Med. Coll. Abbottabad 18, 17–20.

PubMed Abstract | Google Scholar

Hohenhaus, E., Paul, A., McCullough, R., Kucherer, H., and Bartsch, P. (1995). Ventilatory and pulmonary vascular response to hypoxia and susceptibility to high altitude pulmonary oedema. Eur. Respir. J. 8, 1825–1833. doi: 10.1183/09031936.95.08111825

CrossRef Full Text | Google Scholar

Holmström, P., Mulder, E., and Schagatay, E. (2018). “Effects of 300 m and 1500 m altitude on cardiovascular responses during apnea,” in Proceedings of the 12th World Congress on Mountain Medicine, Kathmandu.

Google Scholar

Ilardo, M. A., Moltke, I., Korneliussen, T. S., Cheng, J., Stern, A. J., Racimo, F., et al. (2018). Physiological and genetic adaptations to diving in sea nomads. Cell 173, 569.e15–580.e15. doi: 10.1016/j.cell.2018.03.054

PubMed Abstract | CrossRef Full Text | Google Scholar

Karinen, H. M., Peltonen, J. E., Kähönen, M., and Tikkanen, H. O. (2010). Prediction of acute mountain sickness by monitoring arterial oxygen saturation during ascent. High Alt. Med. Biol. 11, 325–332. doi: 10.1089/ham.2009.1060

PubMed Abstract | CrossRef Full Text | Google Scholar

Karinen, H. M., Uusitalo, A., Vähä-Ypyä, H., Kähönen, M., Peltonen, J. E., Stein, P. K., et al. (2012). Heart rate variability changes at 2400 m altitude predicts acute mountain sickness on further ascent at 3000–4300 m altitudes. Front. Physiol. 3:336. doi: 10.3389/fphys.2012.00336

CrossRef Full Text | Google Scholar

Lee, D. K. (2016). Alternatives to P value: confidence interval and effect size. Korean J. Anesthesiol. 69, 555–562.

PubMed Abstract | Google Scholar

Lodin-Sundström, A., and Schagatay, E. (2010). Spleen contraction during 20 min normobaric hypoxia and 2 min apnea in humans. Aviat. Space Environ. Med. 81, 545–549. doi: 10.3357/asem.2682.2010

PubMed Abstract | CrossRef Full Text | Google Scholar

Luks, A. M., and Swenson, E. R. (2011). Pulse oximetry at high altitude. High Alt. Med. Biol. 12, 109–119. doi: 10.1089/ham.2011.0013

PubMed Abstract | CrossRef Full Text | Google Scholar

Luks, A. M., Swenson, E. R., and Bärtsch, P. (2017). Acute high-altitude sickness. Eur. Respir. Rev. 26:160096.

Google Scholar

Milledge, J., Beeley, J., Broome, J., Luff, N., Pelling, M., and Smith, D. (1991). Acute mountain sickness susceptibility, fitness and hypoxic ventilatory response. Eur. Respir. J. 4, 1000–1003.

PubMed Abstract | Google Scholar

Nespoulet, H., Wuyam, B., Tamisier, R., Saunier, C., Monneret, D., Remy, J., et al. (2012). Altitude illness is related to low hypoxic chemoresponse and low oxygenation during sleep. Eur. Respir. J. 40, 673–680. doi: 10.1183/09031936.00073111

PubMed Abstract | CrossRef Full Text | Google Scholar

Neupane, M., and Swenson, E. R. (2017). High-altitude pulmonary vascular diseases. Adv. Pulm. Hypertens. 15, 149–157. doi: 10.21693/1933-088x-15.3.149

CrossRef Full Text | Google Scholar

Oliver, S. J., Sanders, S. J., Williams, C. J., Smith, Z. A., Lloyd-Davies, E., Roberts, R., et al. (2012). Physiological and psychological illness symptoms at high altitude and their relationship with acute mountain sickness: a prospective cohort study. J. Travel Med. 19, 210–219. doi: 10.1111/j.1708-8305.2012.00609.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Richalet, J.-P., Larmignat, P., Poitrine, E., Letournel, M., and Canouï-Poitrine, F. (2012). Physiological risk factors for severe high-altitude illness: a prospective cohort study. Am. J. Respir. Crit. Care Med. 185, 192–198. doi: 10.1164/rccm.201108-1396oc

PubMed Abstract | CrossRef Full Text | Google Scholar

Richardson, M. X., Engan, H. K., Lodin-Sundström, A., and Schagatay, E. (2012). Effect of hypercapnia on spleen-related haemoglobin increase during apnea. Diving Hyperb. Med. 42, 4–9.

PubMed Abstract | Google Scholar

Roach, C. R., Bartsch, P., Hackett, P., and Oelz, O. (1993). The Lake Louise acute mountain sickness scoring system. Hypoxia Mol. Med. 4, 272–274.

Google Scholar

Roach, R. C., Greene, E. R., Schoene, R. B., and Hackett, P. H. (1998). Arterial oxygen saturation for prediction of acute mountain sickness. Aviat. Space Environ. Med. 69, 1182–1185.

PubMed Abstract | Google Scholar

Schagatay, E., and Andersson, J. (1998). Diving response and apneic time in humans. Undersea Hyperb. Med. 25, 13–19.

PubMed Abstract | Google Scholar

Schagatay, E., Andersson, J., Hallén, M., and Pålsson, B. (2001). Selected contribution: role of spleen emptying in prolonging apneas in humans. J. Appl. Physiol. 90, 1623–1629. doi: 10.1152/jappl.2001.90.4.1623

PubMed Abstract | CrossRef Full Text | Google Scholar

Schagatay, E., Haughey, H., and Reimers, J. (2005). Speed of spleen volume changes evoked by serial apneas. Eur. J. Appl. Physiol. 93, 447–452. doi: 10.1007/s00421-004-1224-0

PubMed Abstract | CrossRef Full Text | Google Scholar

Schagatay, E., Hubinette, A., Lodin-Sundström, A., Engan, H., and Stenfors, N. (2015). Exercise induced hemoconcentration following spleen contraction in subjects with COPD. COPD Res. Pract. 1:13.

Google Scholar

Schagatay, E., Richardson, M., and Lodin-Sundström, A. (2012). Size matters: spleen and lung volumes predict performance in human apneic divers. Front. Physiol. 3:173. doi: 10.3389/fphys.2012.00173

CrossRef Full Text | Google Scholar

Schoene, R., Lahiri, S., Hackett, P., Peters, R., Milledge, J., Pizzo, C., et al. (1984). Relationship of hypoxic ventilatory response to exercise performance on Mount Everest. J. Appl. Physiol. 56, 1478–1483. doi: 10.1152/jappl.1984.56.6.1478

PubMed Abstract | CrossRef Full Text | Google Scholar

Schutema, K., and Holm, B. (1988). The role of different facial areas in eliciting human diving bradycardia. Acta Physiol. Scand. 132, 119–120. doi: 10.1111/j.1748-1716.1988.tb08306.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Starfelt, V., Rodriguez, L., Patrician, A., Olander, C., Lodin-Sundström, A., and Schagatay, E. (2018). “Association between the diving response during apnea at sea-level and the arterial oxygen saturation at simulated high altitude,” in Proceedings of the XII World Congress on Mountain Medicine, Kathmandu.

Google Scholar

Wagner, D. R., Knott, J. R., and Fry, J. P. (2012). Oximetry fails to predict acute mountain sickness or summit success during a rapid ascent to 5640 meters. Wilderness Environ. Med. 23, 114–121. doi: 10.1016/j.wem.2012.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: acute mountain sickness, breath-hold diving, hypoxia, prediction, cardiovascular diving response, spleen

Citation: Holmström P, Mulder E, Sundström AL, Limbu P and Schagatay E (2019) The Magnitude of Diving Bradycardia During Apnea at Low-Altitude Reveals Tolerance to High Altitude Hypoxia. Front. Physiol. 10:1075. doi: 10.3389/fphys.2019.01075

Received: 04 June 2019; Accepted: 05 August 2019;
Published: 22 August 2019.

Edited by:

Costantino Balestra, Haute École Bruxelles-Brabant (HE2B), Belgium

Reviewed by:

Jochen D. Schipke, University Hospital of Düsseldorf, Germany
Danilo Cialoni, Dan Europe Foundation, Italy
Kate Lambrechts, Haute École Bruxelles-Brabant (HE2B), Belgium

Copyright © 2019 Holmström, Mulder, Sundström, Limbu and Schagatay. 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: Pontus Holmström, Pontus.holmstrom@miun.se