Abstract
Background Obstructive sleep apnoea (OSA) causes, among other things, intermittent blood oxygen desaturations, increasing the sympathetic tone. Yet the effect of desaturations on heart rate variability (HRV), a simple and noninvasive method for assessing sympathovagal balance, has not been comprehensively studied. We aimed to study whether desaturation severity affects the immediate HRV.
Methods We retrospectively analysed the electrocardiography signals in 5-min segments (n=39 132) recorded during clinical polysomnographies of 642 patients with suspected OSA. HRV parameters were calculated for each segment. The segments were pooled into severity groups based on the desaturation severity (i.e. the integrated area under the blood oxygen saturation curve) and the respiratory event rate within the segment. Covariate-adjusted regression analyses were performed to investigate possible confounding effects.
Results With increasing respiratory event rate, the normalised high-frequency band power (HFNU) decreased from 0.517 to 0.364 (p<0.01), the normalised low-frequency band power (LFNU) increased from 0.483 to 0.636 (p<0.01) and the mean RR interval decreased from 915 to 869 ms (p<0.01). Similarly, with increasing desaturation severity, the HFNU decreased from 0.499 to 0.364 (p<0.01), the LFNU increased from 0.501 to 0.636 (p<0.01) and the mean RR interval decreased from 952 to 854 ms (p<0.01). Desaturation severity-related findings were confirmed by considering the confounding factors in the regression analyses.
Conclusion The short-term HRV response differs based on the desaturation severity and the respiratory event rate in patients with suspected OSA. Therefore, a more detailed analysis of HRV and desaturation characteristics could enhance OSA severity estimation.
Abstract
Higher short-term HRV is related to more severe oxygen desaturations and a higher rate of respiratory events. Considering HRV and desaturation characteristics in the diagnosis of OSA could be useful when assessing the cardiac consequences of OSA. https://bit.ly/3yZSYTR
Introduction
Obstructive sleep apnoea (OSA) is one of the most prevalent sleep disorders increasing the risk for cardiovascular diseases and other severe sequelae [1–3]. In OSA, nocturnal respiratory events, i.e. apnoeas and hypopnoeas, cause physiological consequences such as oxygen desaturations and oscillations of the heart rate [4]; therefore, they also affect heart rate variability (HRV) [5, 6]. HRV is an effective, noninvasive measure to evaluate the state of the autonomic nervous system, especially the activities of its sympathetic and parasympathetic branches regulating the vital functions of the body [7]. Yet, electrocardiogram (ECG) usage is very limited in the diagnostics of OSA, although it is routinely recorded during polysomnography (PSG) [8].
Together with neuro-cardiac interactions, the sympathetic nervous system (SNS) and parasympathetic nervous system (PNS) regulate the functioning of the heart and, thus, HRV [7]. Normally, stress increases SNS activity, which is associated with increased oxidative stress, systemic inflammation, decreased long-term (≥24 h) HRV and a higher risk for coronary artery disease [9, 10]. Conversely, the PNS dominates the sympathovagal balance during rest, elevating long-term HRV [10]. In OSA, the frequent respiratory events shift the sympathovagal balance towards SNS dominance by causing intermittent hypoxaemia, hypercapnia, intrathoracic pressure changes and recurrent arousals [3]. Furthermore, OSA patients have reduced long-term HRV [5, 11] but the events also lead to increased ultra-short-term time-domain HRV [12]. Long-term HRV cannot be used to assess the immediate physiological responses to respiratory events and hypoxaemia, although decreased long-term HRV is generally associated with poor health [10, 13].
Currently, the severity assessment of OSA is predominantly based on the average number of respiratory events per hour of sleep, i.e. the apnoea-hypopnoea index (AHI), determined from PSG [8]. The AHI does not consider the severity of individual respiratory events or their physiological consequences. Therefore, novel computational parameters using integrated area under the blood oxygen saturation curve have been developed [14, 15] and have a stronger association with the consequences of OSA than the AHI [15–18]. Longer respiratory events cause more severe desaturations and increased heart rate after the event [4, 19] and are associated with higher ultra-short-term (<5 min) HRV [12]. Nonetheless, the effect of desaturations on HRV has not been comprehensively studied in patients with OSA. Because HRV parameters can be easily derived from ECG measurements, we hypothesised that they could be used more extensively to assess the state of the autonomic nervous system in OSA patients. Therefore, in the future they could potentially evaluate the risk for cardiovascular diseases alongside conventional diagnostic parameters.
In this study, our main hypothesis was that, in the short term (5 min), more severe desaturations are associated with higher time-domain HRV, increased low-frequency band power and blunted high-frequency band power in frequency-domain HRV. In addition, we hypothesised that the severity of desaturations affects the short-term HRV more strongly than the rate of respiratory events. Using a large clinical population of patients with a suspicion of OSA, we aimed to investigate whether the severity of desaturations and the rate of respiratory events affect time- and frequency-domain HRV.
Methods
PSG data
We retrospectively studied the type I PSGs of 901 consecutive patients referred to a sleep study owing to clinical suspicion of OSA. The PSGs were recorded at the Princess Alexandra Hospital (Brisbane, Australia) between 2015 and 2017 using the Grael Acquisition System (Compumedics, Abbotsford, Australia). The recordings were manually scored by experienced sleep technicians in compliance with the American Academy of Sleep Medicine 2012 scoring criteria [8]. A respiratory event was scored as an apnoea if the respiratory airflow decreased by ≥90% from the baseline for ≥10 s. An event was scored as a hypopnoea if the respiratory airflow decreased ≥30% from the baseline for ≥10 s and there was arousal or a ≥3% drop in blood oxygen saturation related to the respiratory event. The Institutional Human Research Ethics Committee of Princess Alexandra Hospital approved the data collection and re-use (HREC/16/QPAH/021 and LNR/2019/QMS/54313). After patient exclusion (criteria described in figure 1), we included a total of 642 patients in this study (table 1).
Inclusion and exclusion criteria for the patient and 5-min segment selection. PSG: polysomnography; OSA: obstructive sleep apnoea; ECG: electrocardiogram; Ev: event rate; Ev0=0 events per 5 min, EvQ1=1 event per 5 min<EvQ2≤3 events per 5 min<EvQ3≤5 events per 5 min<EvQ4; DesSev: desaturation severity; DesSev0=0%, DesSevQ1≤0.305%<DesSevQ2≤0.906%<DesSevQ3≤2.463%<DesSevQ4.
Demographic characteristics of the study population
HRV analysis
ECGs were recorded during the PSGs using lead II [8] with a sampling frequency of 256 Hz. We detected the R peaks from the ECGs using the Kubios HRV Premium 3.4.1 software (Kubios Oy, Kuopio, Finland) with the default settings [20]. The Kubios detects R peaks with the algorithm based on the Pan–Tompkins method [21] using the amplitude threshold and expected time between adjacent R peaks as the decision rules. They both are adaptively adjusted after each detected R peak. To improve the time resolution, the R wave is interpolated to 2000 Hz before extracting the time of the R peak. The software also uses automatic correction for artefacts due to ectopic peaks and missed peak detections; detected artefact peaks are replaced using cubic spline interpolation. We divided the resulting R peak time series into non-overlapping 5-min segments during sleep (n=54 453). The exclusion criteria for the segments are described in figure 1.
Further, we divided the remaining 5-min segments (n=39 132) into severity groups based on the respiratory event rate (Ev), including any event starting during the segment, i.e. apnoeas and hypopnoeas. We considered the segments having zero apnoeas and hypopnoeas as the reference group (Ev0, n=20 873). The remaining segments (n=18 259) were divided into quartiles EvQ1 to EvQ4 with thresholds of 1, 3 and 5 events in the 5-min segment (EvQ1=1<EvQ2≤3<EvQ3≤5<EvQ4, table 1). A similar division was separately performed based on the severity of the desaturations starting during the segment. The severity of desaturations was assessed using the desaturation severity (DesSev) parameter that describes the severity of hypoxic load by considering the depth and duration of the individual desaturation events as the integrated area under the blood oxygen saturation curve. The DesSev was calculated by dividing the sum of individual desaturation areas by the duration of the segment (5 min) [14]. The segments having DesSev=0 (DesSev0, n=23 191) were considered as the reference group. The remaining segments (n=15 941) were pooled into quartiles DesSevQ1 to DesSevQ4 with thresholds of 0.305%, 0.906% and 2.463% (0% <DesSevQ1≤0.305%<DesSevQ2≤0.906%<DesSevQ3≤2.463%<DesSevQ4, table 1).
Next, we computed the values of short-term (5 min) time- and frequency-domain HRV parameters for each 5-min segment. The time-domain HRV parameters consisted of the mean RR interval, the standard deviation of the RR intervals from which artefacts had been corrected (SDNN), the root mean square of the successive differences (RMSSD) and the proportion of adjacent RR intervals differing by more than 50 ms (pRR50) [7]. Before the frequency-domain analysis, the 5-min RR interval segments were resampled by cubic spline interpolation and detrended by the smoothness priors method with λ=500 [22]. Then, we estimated the power spectral densities (PSDs) with Welch's method [20]. From PSDs, we calculated the frequency-domain HRV parameters: power in the high-frequency band (HF) (0.15–0.40 Hz), the low-frequency band (LF) (0.04–0.15 Hz), and their ratio (LF/HF ratio) [7]. We additionally calculated the normalised HF and LF band powers, HFNU=HF/(HF+LF) and LFNU=LF/(HF+LF), respectively [7]. The HF and HFNU band powers are thought to represent PNS activity, whereas both SNS and PNS activities contribute to the LF and LFNU band powers [7]. Finally, we calculated the medians of HRV parameter values separately for the Ev and DesSev groups and compared them between 1) Ev groups, 2) DesSev groups and 3) the corresponding Ev and DesSev groups.
Because the severity groups could contain numerous 5-min segments from the same patient, we assumed the groups not to be independent. Therefore, we used a Wilcoxon signed-rank test to evaluate the statistical significance of differences in cases 1), 2) and 3). Because the Wilcoxon signed-rank test conducts a pairwise comparison, we computed a total of 5000 randomly chosen permutations of the HRV parameter pairs in each of cases 1), 2) and 3). Only one HRV parameter was compared between two severity groups at a time. We defined the statistical significance as the median of the 5000 p-values and a significance level of p<0.01 was used for these analyses owing to large sample sizes and multiple testing.
To consider the effect of potential confounding factors on HRV parameters, we also performed covariate-adjusted regression analysis for the 5-min segments. DesSev was the continuous variable and sex, age, body mass index (BMI), daytime sleepiness based on the Epworth Sleepiness Scale (ESS) questionnaire, history of arrhythmias, hypertension, hypothyroidism and chronic obstructive pulmonary disease (COPD) were adjusting covariates. The HRV, statistical and regression analyses were performed with MATLAB R2018b (MathWorks Inc., Natick, MA, USA).
Results
Frequency-domain HRV
The HF, LF and LFNU band powers, and the LF/HF ratio increased while the HFNU decreased with the increasing number of respiratory events within the 5-min segments (table 2, figure 2). All frequency-domain HRV parameter values in different Ev groups were different from the reference group Ev0 (p<0.001). The LF band power increased between each Ev group towards the EvQ4 group (p<0.001).
Frequency-domain HRV parameter values in the severity groups based on the Ev and the DesSev in the 5-min segments
Median power spectral densities (PSDs) in the severity groups pooled based on the respiratory event rate (Ev) within the 5-min segments. Vertical dashed lines separate the frequency bands used in heart rate variability analyses: the low-frequency band (0.04–0.15 Hz) and the high-frequency band (0.15–0.40 Hz). Ev0=0 events per 5 min, EvQ1=1 event per 5 min<EvQ2≤3 events per 5 min<EvQ3≤5 events per 5 min<EvQ4.
The LF and LFNU band powers increased and the HFNU power decreased towards the more severe DesSev groups (table 2, figure 3). Similarly, the HF band power and the LF/HF ratio showed a trend for an increase towards the most severe DesSev group (DesSevQ4). LF and LFNU band powers and LF/HF ratio values in all DesSev groups were higher than in the DesSev0 group (p<0.001). The frequency-domain HRV parameter values in the DesSevQ4 group were higher than all other DesSev groups (p<0.001).
Median power spectral densities (PSDs) in the groups pooled based on the desaturation severity (DesSev) within the 5-min segments. Vertical dashed lines separate the frequency bands used in heart rate variability analyses: the low-frequency band (0.04–0.15 Hz) and the high-frequency band (0.15–0.40 Hz). DesSev0=0%, DesSevQ1≤0.305%<DesSevQ2≤0.906%<DesSevQ3≤2.463%<DesSevQ4.
All HRV parameter values in the DesSev0 group, except HFNU, were higher than in the Ev0 group (p<0.001, table 2). In the DesSevQ2 and DesSevQ3 groups, the LF band power was lower than in the EvQ2 and EvQ3 groups (p<0.001), respectively.
Time-domain HRV
A higher Ev was associated with the shorter mean RR interval and higher other short-term time-domain HRV parameter values (table 3). Within the most severe event group, EvQ4, the mean RR intervals were shorter (p<0.001) and the SDNN, RMSSD and pRR50 parameter values were greater (p<0.001) than in the reference group Ev0 and all other Ev groups. The increase in the SDNN values was significant (p≤0.002) between all Ev groups.
Time-domain HRV parameter values in the severity groups based on the Ev and the DesSev in the 5-min segments
Similarly, shorter mean RR intervals and higher other short-term time-domain HRV parameter values were observed with more severe desaturations (table 3). Mean RR intervals in all DesSev groups were shorter and SDNN values greater than in the reference group DesSev0 (p<0.001). The mean RR intervals were significantly shorter and SDNN and RMSSD values higher in the DesSevQ3 and DesSevQ4 groups than in all other DesSev groups (p<0.001).
In the reference group DesSev0, all time-domain HRV parameter values were greater (p<0.001) than in the Ev0 group (table 3). The mean RR interval and SDNN values in DesSev groups differed from the corresponding Ev groups in the majority of comparisons (p<0.001). The decrease in mean RR intervals from the DesSev0 to DesSevQ4 group was greater than that of the corresponding Ev groups.
Desaturation characteristics in the severity groups
Within the Ev groups, the median number and depth of the desaturations within the 5-min segments increased towards the more severe Ev groups (p<0.001, table 4). However, the desaturations were the longest in the EvQ3 group (p<0.001). Similarly, the median number, duration and depth of the desaturations consistently increased towards the more severe DesSev groups (p<0.001). In general, the desaturations were longer and deeper (p<0.001) within the DesSev groups compared to the corresponding Ev groups.
Desaturation characteristics in the 5-min segments in the severity groups
Covariate-adjusted regression analysis
Covariate-adjusted regression analysis showed that even after adjustments, increasing DesSev decreased the HFNU power and mean RR interval and increased the LFNU power and other time-domain HRV parameters (table 5). Overall, male sex was associated with lower HFNU power, higher LFNU power and higher time-domain HRV parameter values. Age blunted the total PSD power and increased the mean RR intervals, and higher BMI was associated with lower LFNU power, mean RR intervals and the SDNN. The ESS score was only marginally associated with HRV. Of the comorbidities, a history of arrhythmias had the greatest effect on HRV parameters, decreasing the LFNU power and increasing the HFNU power and time-domain HRV parameter values.
Covariate-adjusted regression analysis to investigate the effect of covariates on selected heart rate variability parameters
Discussion
In this study, we investigated whether the severity of the desaturation events (DesSev) and the rate of the respiratory events (Ev) affect the short-term time- and frequency-domain HRV parameters. To our knowledge, the connection between desaturation severity and HRV has not previously been studied in OSA patients. Supporting our hypothesis, we observed that the more severe desaturations were associated with decreased normalised HF band power and higher absolute and normalised LF band powers, LF/HF ratio and short-term time-domain HRV parameter values. In addition, a higher Ev was associated with higher absolute LF band power, LF/HF ratio and mean RR intervals compared to corresponding desaturation severity groups.
Previous studies show that, compared with the AHI, the severity of the nocturnal hypoxic load has a stronger association with OSA symptoms and comorbidities such as daytime sleepiness, heart failure and cardiovascular disease-related mortality [15–18]. Against our hypothesis, the higher rate of respiratory events (apnoeas and hypopnoeas) within 5-min segments was related to generally higher absolute time-domain HRV parameter values and SNS activation based on frequency-domain HRV compared to the corresponding DesSev groups (tables 2 and 3). Greater changes were also seen in HRV between reference and other groups with Ev quartiles (tables 2 and 3). This finding seems fairly counterintuitive because, generally, DesSev quartiles include more desaturations and they are significantly longer and deeper compared to the corresponding Ev quartiles (table 4). Nevertheless, a higher Ev seems to lead to stronger SNS activation (table 2). By contrast, intermittent respiratory events lead to cyclical heart rate variation [23], potentially increasing the short-term HRV, whereas longer desaturations may result in the heart rate decreasing more steadily or even stabilising at a certain level. However, other OSA-related factors could partly explain these differences. Shorter respiratory events lead to less severe desaturations [19] but increase the arousability [24, 25]. Arousals increase the heart rate and the SNS activity [26] and, therefore, they most probably contribute to our findings. In our present study, increasing both DesSev and Ev in the 5-min segments led to higher short-term HRV, which is logical because desaturations are a consequence of the respiratory events. The driving force behind these differences in HRV between Ev and DesSev groups remains unclear and more comprehensive research, including other potential factors such as sleep stages and arousals, is required.
The frequency-domain HRV analysis showed that all parameter values, except HFNU, increased with the increasing Ev and increasing DesSev. There was a higher relative increase in absolute LF compared with absolute HF band power towards the more severe DesSev and Ev groups, which elevated the LF/HF ratio (table 2, figures 2 and 3). In addition, the LFNU band power increased and HFNU power decreased towards more severe DesSev and Ev groups (table 2). These findings indicate a shift towards sympathetic dominance. Previous studies show that patients with OSA have higher SNS activity compared to healthy controls [6, 27, 28]. In addition, they demonstrate an overactive SNS compared to PNS in terms of higher LF band power and LF/HF ratio with more severe OSA, and blunted HF band power [5, 6, 29, 30]. Our present results suggest that increasing hypoxic load leads to a significantly stronger SNS activation. These findings are in line with previous studies [5, 6, 27–30] because longer respiratory events cause more severe desaturations [19]. Based on the present findings and previous literature [2, 3, 15, 31], OSA patients with more severe hypoxic load are at higher risk for cardiovascular diseases because sympathetic overdrive and severe desaturations are major risk factors for numerous other cardiovascular diseases. Long-term intermittent hypoxaemia increases the chemosensitivity of the carotid body and leads to increased SNS activity [32]. Sympathetic overdrive is also one predisposing factor for early morning cardiovascular events [29] and has been suggested to partly explain the increased propensity for arrhythmias in OSA patients [33]; the risk of arrhythmia is higher in OSA patients and is shown to significantly increase shortly after respiratory events [34, 35]. In the present study, patients with a history of arrhythmias had lower sympathetic tone than patients with no arrhythmia history (table 5) but this might have been due to possible medications stabilising heart function. Although these results may manifest the increased cardiovascular risk linked to more severe desaturations, more detailed research on respiratory events, desaturations and their effect on cardiovascular stress is warranted.
In agreement with our hypothesis, the time-domain HRV parameter values increased with the increasing Ev and DesSev. In Ev and DesSev groups, the mean RR interval decreased significantly and other time-domain parameters (SDNN, RMSSD and pRR50) increased towards the more severe Ev and DesSev groups (table 3). This finding is expected because the number, duration and depth of desaturations increased towards the more severe groups (table 4). Furthermore, longer respiratory events result in greater changes in RR intervals, higher ultra-short-term (<5 min) HRV and more severe desaturations [12, 19]. Some studies [6, 28] have also shown that more severe OSA can lead to higher long-term time-domain HRV, although OSA is commonly linked with decreased long-term HRV [5, 11] and autonomic nervous system regulation [36]. However, because we investigated the short-term HRV, these previous studies [6, 28] are not directly comparable to our present findings. Guilleminault et al. [23] have shown that respiratory events are accompanied by cyclical heart rate variation, and the vagally mediated respiratory sinus arrhythmia is known to contribute to e.g SDNN variation within short-term analyses [7]. The higher time-domain HRV parameter values could, therefore, be affected by the hyperpnoea following respiratory events. Another potential factor explaining these findings is that the body adapts to the recurrent intermittent hypoxaemia [37]. Thus, increased short-term HRV might be more harmful owing to increased beat-to-beat variation within a short time, although decreased long-term HRV generally indicates poor health [10, 13]. Based on both time- and frequency-domain HRV results, more severe desaturations and a higher Ev lead to stronger SNS activation. HRV analyses could thus be used alongside the conventional OSA severity parameters, e.g. the AHI, to assess the physiological consequences of OSA and their severity.
The main limitation of this study was not considering the sleep stages or arousals from sleep when analysing HRV within the 5-min segments. Transitions between sleep stages modulate HRV alongside arousals because heart rate decreases towards deeper sleep [26, 37–40] and, therefore, affected our results. We decided not to include sleep stages in our analysis because the 5-min segments consisted of ten 30-s sleep epochs: the determination of a single sleep stage for the individual segment is complicated and shorter segments would not have been suitable for frequency-domain HRV analyses [7]. However, to prevent the bias caused by the wake periods, we excluded segments consisting of ≥50% wake.
A second limitation was only using a short-term period for HRV analysis instead of other analysis periods. To study the immediate effects of the desaturations on HRV, we selected short-term analyses owing to their suitability for this purpose [7]. Third, including patients with multiple comorbidities and medications in this study was a limitation because several comorbidities and anthropometric factors can affect HRV results alongside OSA (table 1) [37]. The results of the covariate-adjusted regression analyses showed that sex, BMI and history of arrhythmia were the confounding factors with the greatest effect on HRV findings alongside desaturation severity (table 5). However, the list of certain comorbidities and medications was incomplete and, thus, they were not used as exclusion criteria. Fourth, the desaturations and respiratory events were counted into the segment in which they started. The last desaturation or respiratory event of the 5-min segment could have continued to the next segment and so the immediate physiological response was not fully reflected in the segment in which it was counted into. A more thorough investigation is warranted that studies the simultaneous effect of the desaturations, respiratory events, sleep stages and other above-mentioned aspects on short-term HRV.
In conclusion, our results show that within a short-term period, the HRV of the patients with OSA suspicion increases with an increasing rate of respiratory events and severity of desaturations. These findings indicate that patients with suspected OSA with more severe desaturations demonstrate a shift towards sympathetic predominance, which in turn increases the risk for cardiovascular diseases. Moreover, ECG is routinely recorded as a part of PSG but it is not used in the current OSA diagnostics. With short-term HRV measurements, the immediate physiological consequences of OSA can be assessed, which cannot be done with long-term analyses. Therefore, a more detailed analysis of the ECG, HRV and desaturation characteristics could provide valuable information on the cardiovascular stress alongside the AHI when diagnosing OSA.
Acknowledgements
We thank Brett Duce (Princess Alexandra Hospital and the Queensland University of Technology, Brisbane, Australia) for providing the polysomnographic data for our analyses.
Footnotes
Provenance: Submitted article, peer reviewed.
Data availability: The data include medical records and personal information, and therefore the data can only be shared within the confinements of the Australian legislation and ethical conventions. Reasonable requests considering data sharing will be individually assessed.
Author contributions: J. Töyräs, S. Myllymaa, T. Leppänen and S. Kainulainen devised the project and the main conceptual ideas for the analysis. S. Hietakoste, T. Karhu, S. Kainulainen, S. Sillanmäki, R. Bailón and T. Penzel carried out analyses and interpretation. S. Hietakoste drafted the manuscript and prepared all the tables and figures. All the authors have revised the manuscript critically, approved the version submitted for publication and agreed to be accountable for all aspects of the work.
Conflict of interest: S. Hietakoste reports grants from The Research Committee of the Kuopio University Hospital Catchment Area, the Research Foundation of the Pulmonary Diseases, the Foundation of the Finnish Anti-Tuberculosis Association, the Instrumentarium Science Foundation and the Päivikki and Sakari Sohlberg Foundation during the conduct of the study.
Conflict of interest: T. Karhu reports grants from the Research Committee of the Kuopio University Hospital Catchment Area for the State Research Funding, the Instrumentarium Science Foundation and the Tampere Tuberculosis Foundation during the conduct of the study.
Conflict of interest: S. Sillanmäki reports grants from Kuopio University Hospital Catchment Area State Research Funding (no 5041798) during the conduct of the study; and personal fees from the Finnish Society of Cardiology, the Finnish Society of Nuclear Medicine, the Finnish Society of Clinical Physiology, the Onni and Hilja Tuovinen Foundation, the Paavo Nurmi Foundation, and the Scandinavian Society of Clinical Physiology, outside the submitted work.
Conflict of interest: R. Bailón reports other support from the Aragon Government and European Social Fund (EU) during the conduct of the study.
Conflict of interest: T. Penzel reports grants from Cidelec and Novartis, grants and personal fees from Löwenstein Medical, other support from Bayer and MyCardio, and personal fees from Neuwirth, the National Sleep Foundation, Cerebra and Jazz Pharma, outside the submitted work; and is a shareholder in Advanced Sleep Research, The Siestagroup GmbH and Nukute.
Conflict of interest: J. Töyräs reports grants from Business Finland, the National Health and Medical Research Council, Australia, and European Commission Horizon 2020, during the conduct of the study.
Conflict of interest: T. Leppänen reports grants from European Union's Horizon 2020 Research and Innovation Programme under grant agreement number 965417, NordForsk (NordSleep project 90458) via Business Finland (5133/31/2018), the Academy of Finland (323536) and the Research Committee of the Kuopio University Hospital Catchment Area for the State Research Funding (5041767 and 5041794), and personal fees from Turku University Hospital, during the conduct of the study.
Conflict of interest: S. Myllymaa reports grants from the Research Committee of the Kuopio University Hospital Catchment Area for the State Research Funding (project numbers 5041770 and 5041797), the Research Foundation of the Pulmonary Diseases and the Foundation of the Finnish Anti-Tuberculosis Association during the conduct of the study.
Conflict of interest: S. Kainulainen reports grants from the Maud Kuistila Memorial Foundation, the Finnish Cultural Foundation central fund, the Tampere Tuberculosis Foundation and Kuopio University Hospital, and other support from the Regional Fund of Northern Savo, during the conduct of the study.
Support statement: The financial support for this study was provided by several sources including the European Union's Horizon 2020 Research and Innovation Programme (decision number 965417), the Research Committee of the Kuopio University Hospital Catchment Area for the State Research Funding (projects 5041767, 5041770, 5041787, 5041790, 5041794, 5041797, 5041798 and 5041804), the Academy of Finland (decision number 323536), Business Finland (decision number 5133/31/2018), Aragon Government and European Social Fund through BSICoS group (EU, T39-20R), the Finnish Cultural Foundation, Foundation of the Finnish Anti-Tuberculosis Association, Instrumentarium Science Foundation, Maud Kuistila Memorial Foundation, Päivikki and Sakari Sohlberg Foundation, Respiratory Foundation of Kuopio Region, Scientific Foundation of the Pulmonary Diseases, and Tampere Tuberculosis Foundation. Funding information for this article has been deposited with the Crossref Funder Registry.
- Received March 4, 2022.
- Accepted July 5, 2022.
- Copyright ©The authors 2022
This version is distributed under the terms of the Creative Commons Attribution Licence 4.0.