- 1Department of ENT and Head and Neck Surgery, Antwerp University Hospital, Edegem, Belgium
- 2Faculty of Medicine and Health Sciences, University of Antwerp, Wilrijk, Belgium
- 3Division of Sleep and Circadian Disorders, Brigham and Women's Hospital and Harvard Medical School, Boston, MA, United States
- 4Multidisciplinary Sleep Disorders Center, Antwerp University Hospital, Edegem, Belgium
Study Objectives: Obstructive sleep apnea (OSA) is increasingly recognized as a complex and heterogenous disorder. As a result, a “one-size-fits-all” management approach should be avoided. Therefore, evaluation of pathophysiological endotyping in OSA patients is emphasized, with upper airway collapse during sleep as one of the main features. To assess the site(s) and pattern(s) of upper airway collapse, natural sleep endoscopy (NSE) is defined as the gold standard. As NSE is labor-intensive and time-consuming, it is not feasible in routine practice. Instead, drug-induced sleep endoscopy (DISE) is the most frequently used technique and can be considered as the clinical standard. Flow shape and snoring analysis are non-invasive measurement techniques, yet are still evolving. Although DISE is used as the clinical alternative to assess upper airway collapse, associations between DISE and NSE observations, and associated flow and snoring signals, have not been quantified satisfactorily. In the current project we aim to compare upper airway collapse identified in patients with OSA using endoscopic techniques as well as flow shape analysis and analysis of tracheal snoring sounds between natural and drug-induced sleep.
Methods: This study is a blinded prospective comparative multicenter cohort study. The study population will consist of adult patients with a recent diagnosis of OSA. Eligible patients will undergo a polysomnography (PSG) with NSE overnight and a DISE within 3 months. During DISE the upper airway is assessed under sedation by an experienced ear, nose, throat (ENT) surgeon using a flexible fiberoptic endoscope in the operating theater. In contrast to DISE, NSE is performed during natural sleep using a pediatric bronchoscope. During research DISE and NSE, the standard set-up is expanded with additional PSG measurements, including gold standard flow and analysis of tracheal snoring sounds.
Conclusions: This project will be one of the first studies to formally compare collapse patterns during natural and drug-induced sleep. Moreover, this will be, to the authors' best knowledge, the first comparative research in airflow shape and tracheal snoring sounds analysis between DISE and NSE. These novel and non-invasive diagnostic methods studying upper airway mechanics during sleep will be simultaneously validated against DISE and NSE.
Clinical Trial Registration: www.ClinicalTrials.gov, identifier: NCT04729478.
Introduction
Obstructive sleep apnea (OSA) is a highly prevalent disorder affecting almost one billion individuals worldwide (1). This entity can be best described as a heterogenous condition with recurrent upper airway collapse provoking nocturnal breathing cessation (2). The pathophysiology of OSA is multifactorial and varies considerably between individuals (3). One of the key determinants is an anatomically narrow and collapsible upper airway (4, 5). However, non-anatomical traits are also important contributors to OSA pathogenesis, including an ineffective pharyngeal muscle response during sleep, an oversensitive ventilatory control system (i.e., high loop gain), and a low arousal threshold (6). The PALM concept combines these pathophysiological traits on a three-point scale (passive critical closing pressure, arousal threshold, loop gain, and muscle responsiveness), enabling categorization of patients into three categories: PALM scale 1 consists of patients with a high collapsible upper airway, necessitating a major anatomic intervention. PALM scale 2 includes patients with a moderate collapsibility who are candidates for an anatomic or non-anatomic intervention or a combination of both according to the other OSA traits, and PALM scale 3 is constituted of patients with predominantly non-anatomic traits and overall a less severe OSA (6). This research project will focus on the anatomical OSA traits, with a particular focus on the site of upper airway collapse.
The standard treatment for OSA is to pneumatically splint the upper airway throughout the respiratory cycle using continuous positive airway pressure (CPAP) (7). However, CPAP-therapy is often poorly tolerated with non-adherence rates ranging between 46 and 83% (8). Other treatment modalities, such as mandibular advancement devices (9), pharyngeal surgery (10), and upper airway stimulation (11) produce beneficial results in many cases but are limited by a variable efficacy. Thus, in light of the complex and heterogenous pathophysiology of OSA, treatment should be individually tailored according to various phenotypic traits (12–14).
At present, endoscopy-based methods are the clinical mainstay of outcome prediction and patient selection, allowing to identify the sites and patterns of upper airway collapse (15). The gold standard method, natural sleep endoscopy (NSE), visualizes upper airway dynamics using a pediatric bronchoscope during normal sleep (16). However, due to the labor-intensive and time-consuming nature of NSE, including overnight measurements, drug-induced sleep endoscopy (DISE) was developed as an alternative method for routine clinical practice (17). Using sedation to induce a sleep-like state, DISE offers the advantage of being performed during daytime hours in a high-volume setting (18). In the last two decades, DISE has gained great popularity for assessing the sites and patterns of upper airway collapse in patients with OSA.
Collapse can occur at multiple levels along the pharynx, including the soft palate, oropharynx, tongue base, and epiglottis. The absence of complete concentric collapse at the level of the palate (CCCp) and the presence of tongue base collapse during DISE may be positive predictors for therapeutic success with hypoglossal nerve stimulation therapy and mandibular advancement devices (19–21). The same favorable association between tongue base collapse and mandibular advancement device outcome was found by a recent study using NSE (22). Conversely, epiglottis collapse during DISE may negatively affect the outcome of mandibular advancement devices, CPAP, and upper airway surgery (23). However, this collapse type may respond well to changes in body position from supine to lateral as seen during both DISE and NSE (23, 24).
A recent literature review by our research group concluded that velar collapse is seen during NSE in 58.8% of patients, tongue base collapse in 43.2%, lateral wall collapse in 29.9%, and epiglottis collapse in 22.4% (25). Drug-induced sleep endoscopy studies in the past reported comparable findings with the highest prevalence of collapse at the level of the soft palate (81.0%), followed by the tongue base (46.6%), hypopharynx (38.7%), and oropharynx (21.9%) (15).
Park et al. recently explored the differences in obstruction patterns between DISE using midazolam and NSE (26). They demonstrated significant concordance in epiglottis collapse (92.3%), lateral wall collapse (88.5%), and velar collapse (76.9%). However, the agreement in tongue base collapse was relatively low (69.2%). Compared to NSE, CCCp was observed slightly more often during DISE. Regarding the extent of obstruction, higher obstruction grades were present during DISE observations, with more sites in need of treatment, in comparison to NSE. Overall, these findings would suggest that DISE is a reliable test. However, this study only included a small patient sample in the supine position.
An evolving method to non-invasively determine the sites of upper airway collapse is airflow shape analysis. Specific airflow features have been previously validated by using DISE and NSE observations (27–31). Negative effort dependence (NED) as a flow limitation characteristic is quantified as the percentage reduction in inspiratory flow from peak to plateau associated with increasing respiratory effort ([flowpeak − flowplateau]/flowpeak). In a preliminary assessment of correlations between flow measurements and upper airway collapse sites during DISE, the authors observed that the presence of epiglottis collapse was associated with a higher NED compared to the absence of epiglottis collapse or collapse at other pharyngeal sites (30). This is in line with previous findings of Genta et al., who also found that tongue base collapse is associated with a low NED, and isolated palatal and lateral wall collapse with a moderate NED (27). Using a machine learning algorithm, Azarbarzin et al. expanded this work by developing a model that discriminated epiglottis from non-epiglottis collapse with 84% cross-validated accuracy (29). The main distinctive parameters of this model were the discontinuity index (slope of the steepest line fitted to the inspiratory airflow signal) and inspiratory jaggedness (extent of deviation from a flat line during inspiration). The same research group also identified expiratory flow limitation (recognizable as pinching of the airflow signal) as the hallmark of palatal prolapse (29). According to recent evidence, this pinching feature is associated with a negative response to mandibular advancement device treatment, indicating the potential of this non-invasive assessment (32).
Different parts of the upper airway may generate different snoring sounds. Therefore, acoustic analysis of tracheal snoring sounds could be a useful indicator of the site of obstruction in patients with OSA by determining the number of snore events, the intensity of snoring, the sound frequencies, and/or the spectrogram (33, 34). Liistro et al. were the first in 1991 who examined the relationship between snoring characteristics and the anatomical site of collapse by using a supraglottic pressure catheter (35). A study of Xu et al. showed a difference in sound spectrum between upper and lower-level obstructive apneas using esophageal and oropharyngeal pressure measurements. An obstruction level above the free margin of the soft palate produced a characteristic frequency in the low frequency domain and an obstruction level below the free margin of the soft palate in the high frequency domain (36). However, most studies have not used direct visualization of the upper airway with DISE. Recent studies using DISE have shown similar associations between the sites of collapse and certain frequency bands (37, 38). Palatal snoring was identified at low frequencies (137, 105–189 Hz), tongue base snoring at high frequencies (1,243, 1,215–1,277 Hz), and epiglottic snoring in a mid-frequency range (490, 331–510 Hz) (38). Multilevel obstruction during DISE might be associated with complex snoring sounds, expressing as an irregular pattern of snoring (37) or a high maximal intensity of low-frequency snoring sounds (≥60 dB) (33). Koo et al. stated that the spectrograms and frequencies of snoring sounds of drug-induced sleep did not differ significantly from those of natural sleep (39). Nevertheless, induced sleep might show a slightly higher intensity than natural sleep, predominantly on a retro-lingual obstruction level, and might be more irregular in comparison with natural sleep on spectrogram (34, 38–41). Sebastian et al. used an automated technique from audio signal recordings to identify the primary site of upper airway collapse. The automated system showed a high accuracy for identifying tongue and non-tongue related collapse (macro-average recall of 80%) but with only low unweighted average recalls for classifying all sites (34). A review of Penzel and Sabil supports the hypothesis of tracheal sound analysis to be a surrogate for respiratory flow (42). In a more recent study of the same research group with a small sample size (n = 32), the detection of apneas was compared using four different methods of airflow signals, concluding tracheal sounds to be a possible future substitute for oral thermistors (43). However, more studies utilizing endoscopic visualization of the upper airway during sleep, are necessary to validate this hypothesis. To conclude, tracheal snoring sounds analysis may be used as a useful screening tool to determine the site of obstruction during sleep.
In the current project, we aim to compare endoscopic findings of upper airway collapse as well as flow and tracheal snoring sounds analysis between natural and drug-induced sleep in patients with OSA.
Materials and Methods
Study Objectives
The study aims are as follows:
1. To compare upper airway collapse observed during NSE and DISE in patients with OSA.
2. To compare the airflow shapes caused by the different sites of collapse, between NSE and DISE.
3. To compare the acoustics of snores (i.e., tracheal snoring sounds) generated by the different sites of upper airway collapse between NSE and DISE.
4. To validate previously developed models for predicting the sites of collapse from airflow and snoring sounds and to compare prediction performance between NSE and DISE (27–31, 44, 45).
Patient Recruitment
Regarding all four aims, patients participating in this study must follow the inclusion and exclusion criteria as listed in Table 1.
To achieve the objectives, subjects will be prospectively recruited at the Department of Otorhinolaryngology in the Antwerp University Hospital and at the Department of Sleep and Circadian Disorders in the Brigham and Women's Hospital between March 1, 2021 and September 30, 2022. Every year, approximately 400 patients undergo a DISE at both departments combined. As a result, we expect to be able to include 40 participants during this time frame. Being a pilot study, no formal sample size calculations can be performed. Power analyses will be performed retrospectively to include additional patients if necessary.
Study Design
This study is a blinded prospective comparative multicenter cohort study. The study population will consist of adult patients with a recent diagnosis of OSA [baseline apnea-hypopnea index (AHI) ≥15/h based on polysomnography (PSG)]. After receiving informed consent, eligible patients will undergo a research PSG (with NSE) and a subsequent DISE procedure within 3 months (Figure 1). The research PSG will be identical to the clinical PSG (i.e., overnight sleep in a sleep laboratory with standard measurements), except that it includes NSE, as well as airflow measurements using a pneumotachometer, and acoustic measurements of snoring using a microphone placed over the trachea. For the DISE procedure, the patient is lying supine in a semi-dark and silent operating theater. An intravenous bolus injection of midazolam 1.5 mg will be used to induce sleep. Maintenance of sedated sleep will be obtained by a target-controlled infusion of propofol (2.0–3.0 μg/ml). Scoring of obstruction sites (for DISE and NSE) will be done using two standardized scoring systems as previously described by our research team (Figure 2) (25). Different ear, nose, throat (ENT) surgeons within the research team will perform and score the DISE and NSE of the same participating patient to ensure that all investigations are scored in a blinded fashion. During DISE, we will also use electroencephalography (EEG) and electromyography (EMG) measurements to enable scoring of sleep stages, and to control for them in our comparison of DISE and NSE.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/768973/fneur-12-768973-HTML-r1/image_m/fneur-12-768973-g001.jpg)
Figure 1. Generic study protocol. AHI, apnea-hypopnea index; BMI, body mass index; CPAP, continuous positive airway pressure; DISE, drug-induced sleep endoscopy; EEG, electroencephalography; EMG, electromyography; NSE, natural sleep endoscopy; OSA, obstructive sleep apnea; PSG, polysomnography.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/768973/fneur-12-768973-HTML-r1/image_m/fneur-12-768973-g002.jpg)
Figure 2. Endoscopic classification systems. HMS, Harvard Medical School. VOTE classification system (26, 46). HMS classification system (22, 27, 29).
Study Materials
A type I full-night PSG will be performed using standard sleep study equipment, following latest American Academy of Sleep Medicine (AASM) guidelines (47, 48). Classic PSG consists of EEG channels (prefrontal, central, and occipital leads with a reference lead to the mastoids), three surface EMG channels measuring the activity of the submentalis and the bilateral tibialis anterior muscles, two electro-oculography (EOG) channels, one electrocardiography (ECG) channel, oxygen saturation measurement, a position detector and measurement of thoraco-abdominal movements.
Airway visualization during natural sleep will be achieved by a 2.8 mm diameter pediatric bronchoscope (BF-XP-190, Olympus Europe, Hamburg, Germany) with attached video-processor (CV-190, Olympus Europe, Hamburg, Germany) inserted through an oronasal mask. Before insertion, topical application of a decongestant and topical anesthetic will be applied. The patient will be asked to sleep in the supine position during the measurements. Before sleep onset, the tip of the endoscope will be positioned just above the soft palate. Subsequently, collapse at different upper airway levels will be observed during the night by repositioning the endoscope repeatedly at the naso- and oropharynx. An overview of the set-up is given in Figure 3.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/768973/fneur-12-768973-HTML-r1/image_m/fneur-12-768973-g003.jpg)
Figure 3. Overview of the set-up to simultaneously measure airflow and capture endoscopy footage. ECG, electrocardiography; EEG, electroencephalography; EMG, electromyography; EOG, electrooculography.
Airflow will be assessed using a calibrated pneumotachograph with heater control (3,700 A pneumotachographs and 3,850 AF pneumotach heater control, Hans-Rudolph, Shawnee, KS, USA) attached to a calibrated pressure and volume transducer with matching carrier demodulator (Validyne Engineering, Northridge, CA, USA). Mask pressure will be measured using a pressure transducer referenced to atmosphere and arterial oxygen will be monitored at the fingertip. Moreover, an esophageal pressure catheter will be inserted into one of the nostrils through a sealed hole in the mask.
During DISE, the airway will be visualized with a flexible nasopharyngoscope (END-GP, 3.7 mm diameter, Olympus Europe, Hamburg, Germany). In addition to the endoscopic video recordings, we will also measure pneumotachograph flow (Hans-Rudolph), calibrated mask pressure, EEG, EOG, chin EMG, SaO2, and thoraco-abdominal movements. These additional measurements will be recorded using an integrated Alice 6 LDx system (Philips Respironics, Murrysville, Pennsylvania, USA), synchronized with the DISE video acquisition system (Figure 3).
Additional acoustic measurements will be performed during NSE and DISE to combine tracheal snoring sounds with flow patterns and examine associations with the site(s) of collapse. A calibrated omnidirectional zoom microphone (Sony ECM-77B) will be attached to the trachea. A custom-made holder will keep the microphone in place. Gain will be obtained through an M-track 2 × 2 (M-audio). Sound signals will be integrated into the DISE and NSE set-up to allow for time-coupling of all signals.
Signal Processing
Several features, including respiratory measures in time and frequency domains, will be extracted from the flow and acoustic signal. The endotypic traits—collapsibility, loop gain, arousal threshold, and muscle responsiveness—will be determined as described previously (49–52). First, ventilatory drive, defined as ventilation during unobstructed breaths, will be calculated as tidal volume × respiratory rate (L.min−1) using diaphragm EMG. Determination of the percentage of flow:drive will be done by dividing ventilation (L.min−1) by ventilatory drive (L.min−1). Collapsibility of the pharynx will be characterized by the ventilation at a normal ventilatory drive (Vpassive) during sleep. Ventilatory control stability is most commonly defined by loop gain, which will be calculated as the ratio of the ventilatory response and the associated ventilatory disturbance (50). The arousal threshold will be measured as the mean ventilatory drive preceding the start of an arousal (51). The increment in ventilation from Vpassive to ventilation at the arousal threshold will be computed, indicating the muscle responsiveness (52).
Individual flow limited breaths will be scored and flow shape analysis will be performed according to Mann et al. (44). A particular focus will be on NED, as this feature was commonly used in previous research to study upper airway mechanics during DISE (27, 30).
Acoustic features from tracheal snores will be estimated primarily in the frequency domains. Tracheal snoring sounds first undergo preprocessing with a high pass filter with a cutoff frequency of 5 Hz. Features are then extracted from the processed snore signal on a window-level with overlap (e.g., 100 ms window with 50% overlap). These features will include those developed for speech analysis that have been adapted for use in analyzing snoring sounds (33, 34, 37, 38). They include, but are not limited to, power in relevant frequency bands (45), formant frequencies, pitch, harmonicity, mel-frequency cepstral coefficients, and several other metrics describing the power spectral density of the snore signal. Time domain features will also be explored, such as zero-crossing rate, energy, snore duration, and within-breath timing of peak snore amplitude. Snores will be evaluated at the breath-level by averaging windows within the bounds of inspiratory start and stop times (for inspiratory snores) and expiratory start and stop times (for expiratory snores) for all breaths. All breath-level features, such as snore features, flow shape features, and ventilation, are stored in a table detailing all breaths for the night (along with event type, sleep state, etc.), which can be used for various analysis applications. These applications include models for predicting the site of collapse from snore sounds and/or flow shapes that can be applied to predict the site of collapse at the breath-level, which can then be utilized to assess likelihood of response to alternative non-CPAP therapies.
The primary site of collapse of each individual flow-limited breath will be labeled using a custom-made scoring tool. Patient level analysis requires summarizing flow shapes and acoustic features for the patient (e.g., mean values of each feature) and using patient level scoring of site of collapse as the outcome. The overall clinical scoring will be used for analysis on a patient level (Figure 2) (25). The comparison of flow shape and tracheal snore features of specific sites of collapse between DISE and NSE (Aim 3), will be performed on both a breath level and a patient level. Further, model validation for predicting site of collapse using flow shape and tracheal snore features (Aim 4) will also be done on both a breath level and a patient level.
Statistical Analysis
The chi squared test will be performed to compare the site, degree, and pattern of upper airway collapse on a patient level between NSE and DISE (Aim 1). Analysis, linking flow, and acoustic parameters to the OSA traits, will be done on a patient level and on a breath level using Matlab (The MathWorks, Inc., Natick, Massachusetts, USA) (Aim 2). To compare flow shape metrics and acoustic features from snores for each primary site of collapse between DISE and NSE (on a patient level and a breath level), unpaired comparison tests (e.g., t-test, Mann-Whitney U-test) will be performed (Aim 3). Lastly, for Aim 4 we will prospectively validate previously developed models (27–31, 44, 45) that predict the site of pharyngeal collapse on DISE and NSE data. Prediction will be compared to the true site of collapse labels generated by expert reviewers, as previously described. Prediction performance between DISE and NSE will be qualitatively compared.
Discussion
Obstructive sleep apnea is increasingly recognized as a complex and heterogenous disorder in terms of its causes, clinical expression and susceptibility to comorbidities. As a result, a “one-size-fits-all” management approach is not appropriate and should therefore be avoided. A more personalized patient selection through an endotyping approach would allow not only to increase true clinical effectiveness but also to avoid the current “trial-and-error” method which can significantly delay appropriate treatment and lead to a potential waste of resources. Therefore, there is a need to evaluate pathophysiological endotyping in OSA patients, with upper airway collapse at one or more upper airway levels during sleep as one of the main pathophysiological features. Knowing which trait(s) should be addressed in each patient is crucial to predict and achieve response to non-CPAP therapies (53, 54). To assess the site(s) and pattern(s) of upper airway collapse, NSE is defined as the gold standard. Flow shape and tracheal snoring analysis are evolving as non-invasive measurement techniques. Although DISE is used as the clinical alternative to assess the site(s) and pattern(s) of upper airway collapse, associations between DISE and NSE collapse, and associated flow and snoring patterns, have not been quantified satisfactorily.
Study Protocol
The proposed study protocol has different crucial aspects. A recent clinical PSG no more than 2 years old is needed to confirm OSA diagnosis. Hereby, we prevent patients being included in the study protocol while they do not suffer from moderate to severe OSA. The completion of the research PSG marks the official start for the study patient. During this PSG, baseline gold standard measurements will be performed and coupled with NSE video recordings. After the research PSG, a DISE with flow measurements and tracheal snoring sounds analysis will be scheduled. During both NSE and DISE, flow measurements and acoustic analysis of tracheal snoring sounds will be performed as the latter are assumed to be innovative methods to predict upper airway collapse, in accordance with previous research (33, 34, 42, 43). This comparison not only allows us to compare collapse patterns during NSE and DISE, but also to validate previously developed models for predicting the sites of collapse from airflow and snoring sounds and to compare prediction performance between NSE and DISE.
Study Aims
Our first aim focuses on the comparison in collapse patters between NSE and DISE. This project will be one of the first studies to formally compare both examinations in each subject. No restriction will be put on the desired treatment for these patients. This is crucial as we aim to collect data in OSA patients with different collapse sites. To allow proper prediction, sufficient data is needed for each site of obstruction. For example, patients with CCCp, as assessed during DISE, are currently excluded for upper airway stimulation therapy (19). Thus, we would exclude this collapse type by only including patients eligible for upper airway stimulation therapy.
Moreover, this will be, to the authors' best knowledge, the first comparative research in flow shape analysis and tracheal snoring sounds analysis between DISE and NSE. As shown by Vanderveken et al. (55), DISE requires an additional step in the clinical pathway. A future goal would be to omit this step in certain patients to achieve a faster screening process resulting in more time- and cost-effective treatment. Upfront response prediction using the combination of flow characteristics and tracheal snoring sounds has the potential to optimize and shorten the clinical pathway. By reducing the clinical pathway and as such reducing costs, OSA treatment will be accessible for a larger patient population, reducing the portion of patients remaining undertreated.
Limitations
The performance of NSE requires the presence of one or two investigators during a whole night to examine only one to two patients. Moreover, blinding will be attained by having different ENT surgeons performing the DISE and NSE of the same participating patient, which will necessitate another additional researcher. Furthermore, NSE is an invasive study for the patient, accompanied with numerous challenges for the clinicians, such as the risk of awakening the patient when moving the scope to a different upper airway level. However, in previous studies subjects seemed to tolerate scope insertion and were able to sleep during the examination with minimal sleep disruption during manipulation of the endoscope (16, 27–29, 31, 56, 57). Nonetheless, it is a highly time consuming and labor-intensive study, limiting the feasibility to examine many patients in a short time frame and making it a certain challenge to achieve a large sample size.
Natural sleep endoscopy will only be performed in the supine position, which makes it hard to evaluate upper airway collapse in different sleeping positions. However, upper airway collapse is often more pronounced in the supine position, depending on the responsible structure, which makes it a reliable body position for this study (56, 58). During DISE, the upper airway will additionally be evaluated after lateral head rotation. In this way, changes in upper airway collapse related to body position can still be assessed in each patient.
The effect of sedative drugs used during DISE (propofol and midazolam) on upper airway collapse is not fully known. Previous studies showed that respiratory parameters such as the AHI and the mean oxygen saturation remained unchanged during propofol sedation (59). Moreover, the use of propofol did not induce snoring in healthy subjects in earlier publications (59). However, according to previous studies, propofol might lead to a 40% greater decrease in genioglossal muscle activity compared to natural nonrapid eye movement sleep (60, 61) and increasing propofol concentrations were shown to be associated with an increase in collapsibility (62–65). Therefore, the comparison with an endoscopic evaluation during natural sleep might become even more interesting.
Normal sleep at home or during a clinical PSG differs from a research measurement during NSE because of the presence of a pediatric bronchoscope, a sealed mask, an esophageal pressure catheter and a pneumotachograph, which may influence the results. However, according to Maddison et al., a catheter into the pharynx might not have an influence on upper airway collapsibility (66). Therefore, we argue that NSE remains a relatively reliable examination that closely resembles natural sleep under normal conditions.
Study Set-Up Advantages
By including measurements during both natural and drug-induced sleep, we will be able to compare the results obtained during both conditions. This is of great importance as only limited formal comparison of collapse patterns during natural and drug-induced sleep has been done. Moreover, since CCCp, as assessed during DISE, is currently considered as a formal exclusion parameter for upper airway stimulation therapy, this comparison is of utmost importance to understand OSA pathophysiology.
Second, an esophageal pressure catheter will be added to the set-up, which may improve the accuracy of detecting flow limitation. Therefore, esophageal pressure measurements may play an important role in distinguishing obstructive from central apneas. However, the presence of endoscopy during drug-induced or natural sleep also provides certitude in this differentiation. If the quality of the endoscopic images during natural sleep would be disturbed at certain hours of the night, the esophageal pressure measurements might provide us a back-up in evaluating the presence or absence of upper airway obstruction.
Third, this is an investigator-blinded comparative cohort study with different ENT surgeons conducting DISE and NSE in each patient, thereby minimizing the risk of biased observations.
Fourth, both the NSE set-up and DISE set-up with additional measurements have been successful in preliminary testing and are operational. One of the primary advantages of this dataset is that it enables within-patient comparison of the differences in the site of collapse when under natural and drug-induced sleep. Further, it allows for comparison between DISE and NSE of non-invasive measures of the site of collapse derived from flow shapes and snoring sounds. This would then facilitate the validation of some earlier developed models, based on drug-induced sleep and imply that they are applicable to normal sleep. Furthermore, the finding that gold standard measures (e.g., site of pharyngeal collapse) from DISE and NSE are sufficiently similar, would imply that new non-invasive measures (e.g., flow shape/snoring sound derived site of collapse) could be effectively developed and validated from DISE studies. The key advantage being that DISE studies can be collected more rapidly than NSE studies, will make it feasible to collect large patient datasets required to develop more complex prediction models.
Conclusion
This project will be one of the first studies to formally compare collapse patterns during natural and drug-induced sleep. Moreover, this will be, to the authors' best knowledge, the first comparative research in airflow shape and tracheal snoring sounds analysis between DISE and NSE. These novel and non-invasive diagnostic methods studying upper airway mechanics during sleep will be simultaneously validated against DISE and NSE.
Ethics Statement
The studies involving human participants were reviewed and approved by the Institutional Ethics Committee of the University of Antwerp and Antwerp University Hospital (Protocol no. 21/02/009). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
OV, AW, DV, and SO worked on the conception. OV, AW, DV, SO, KV, EV, EK, MW, and JV designed the methodology. KV, EV, DV, and SO contributed to the manuscript first draft. KV drafted the final version of the article. All authors contributed to revision and final approval of the manuscript.
Funding
This study will be supported by grants from the National Institutes of Health (HL102321 and HL128658) to AW, a Senior Clinical Fellowship Grant (Fundamenteel Klinisch Mandaat) from the Research Foundation Flanders-Vlaanderen (1833517N) to OV, and a postdoctoral fellowship at the Research Foundation Flanders (FWO) (1299822N) to SO.
Conflict of Interest
AW reports grants and personal fees from SomniFix; grants from Sanofi and fees from Nox, Apnimed, Bayer, and Inspire Medical Systems outside the submitted work. JV reports grants and fees from SomnoMed, AirLiquide, Vivisol, Mediq Tefa, Medidis, OSG, Agfa-Gevaert, Accuramed, Bioprojet, Jazz Pharmaceutics, Desitin, Idorsia, Nightbalance, Inspire Medical Systems, Heinen and Löwenstein, Ectosense, Philips, and ResMed outside the submitted work. OV reports grants from Inspire Medical Systems, Nightbalance, GlaxoSmithKline, and LivaNova at the Antwerp University Hospital outside the submitted work. The funders were not involved in the study design, collection, the writing of this article or the decision to submit it for publication.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher's Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Abbreviations
AASM, American Academy of Sleep Medicine; AHI, apnea-hypopnea index; BMI, body mass index; CCCp, complete concentric collapse at the level of the palate; CPAP, continuous positive airway pressure; DISE, drug-induced sleep endoscopy; EEG, electroencephalography; ECG, electrocardiography; EMG, electromyography; EOG, electro-oculography; ENT, ear, nose, throat; HMS, Harvard Medical School; NED, negative effort dependence; NSE, natural sleep endoscopy; OSA, obstructive sleep apnea; PSG, polysomnography.
References
1. Benjafield AV, Ayas NT, Eastwood PR, Heinzer R, Ip MSM, Morrell MJ, et al. Estimation of the global prevalence and burden of obstructive sleep apnoea: a literature-based analysis. Lancet Respir Med. (2019) 7:687–98. doi: 10.1016/S2213-2600(19)30198-5
2. Malhotra A, White DP. Obstructive sleep apnoea. Lancet. (2002) 360:237–45. doi: 10.1016/S0140-6736(02)09464-3
3. Eckert DJ. Phenotypic approaches to obstructive sleep apnoea – new pathways for targeted therapy. Sleep Med Rev. (2018) 37:45–59. doi: 10.1016/j.smrv.2016.12.003
4. Schwab RJ, Gupta KB, Gefter WB, Metzger LJ, Hoffman EA, Pack AI. Upper airway and soft tissue anatomy in normal subjects and patients with sleep-disordered breathing. Significance of the lateral pharyngeal walls. Am J Respir Crit Care Med. (1995) 152:1673–89. doi: 10.1164/ajrccm.152.5.7582313
5. Schwab RJ, Pasirstein M, Pierson R, Mackley A, Hachadoorian R, Arens R, et al. Identification of upper airway anatomic risk factors for obstructive sleep apnea with volumetric magnetic resonance imaging. Am J Respir Crit Care Med. (2003) 168:522–30. doi: 10.1164/rccm.200208-866OC
6. Eckert DJ, White DP, Jordan AS, Malhotra A, Wellman A. Defining phenotypic causes of obstructive sleep apnea. Identification of novel therapeutic targets. Am J Respir Crit Care Med. (2013) 188:996–1004. doi: 10.1164/rccm.201303-0448OC
7. Kushida CA, Littner MR, Hirshkowitz M, Morgenthaler TI, Alessi CA, Bailey D, et al. Practice parameters for the use of continuous and bilevel positive airway pressure devices to treat adult patients with sleep-related breathing disorders. Sleep. (2006) 29:375–80. doi: 10.1093/sleep/29.3.375
8. Weaver TE, Grunstein RR. Adherence to continuous positive airway pressure therapy: the challenge to effective treatment. Proc Am Thorac Soc. (2008) 5:173–8. doi: 10.1513/pats.200708-119MG
9. Ramar K, Dort LC, Katz SG, Lettieri CJ, Harrod CG, Thomas SM, et al. Clinical practice guideline for the treatment of obstructive sleep apnea and snoring with oral appliance therapy: an update for 2015. J Clin Sleep Med. (2015) 11:773–827. doi: 10.5664/jcsm.4858
10. Smith DF, Cohen AP, Ishman SL. Surgical management of OSA in adults. Chest. (2015) 147:1681–90. doi: 10.1378/chest.14-2078
11. Strollo PJ. Jr., Soose RJ, Maurer JT, de Vries N, Cornelius J, Froymovich O, et al. Upper-airway stimulation for obstructive sleep apnea. N Engl J Med. (2014) 370:139–49. doi: 10.1056/NEJMoa1308659
12. Koutsourelakis I, Safiruddin F, Ravesloot M, Zakynthinos S, de Vries N. Surgery for obstructive sleep apnea: sleep endoscopy determinants of outcome. Laryngoscope. (2012) 122:2587–91. doi: 10.1002/lary.23462
13. Edwards BA, Andara C, Landry S, Sands SA, Joosten SA, Owens RL, et al. Upper-airway collapsibility and loop gain predict the response to oral appliance therapy in patients with obstructive sleep apnea. Am J Respir Crit Care Med. (2016) 194:1413–22. doi: 10.1164/rccm.201601-0099OC
14. Joosten SA, Leong P, Landry SA, Sands SA, Terrill PI, Mann D, et al. Loop gain predicts the response to upper airway surgery in patients with obstructive sleep apnea. Sleep. (2017). 40:1–10. doi: 10.1093/sleep/zsx094
15. Vroegop AV, Vanderveken OM, Boudewyns AN, Scholman J, Saldien V, Wouters K, et al. Drug-induced sleep endoscopy in sleep-disordered breathing: report on 1,249 cases. Laryngoscope. (2014) 124:797–802. doi: 10.1002/lary.24479
16. Borowiecki B, Pollak CP, Weitzman ED, Rakoff S, Imperato J. Fibro-optic study of pharyngeal airway during sleep in patients with hypersomnia obstructive sleep-apnea syndrome. Laryngoscope. (1978) 88:1310–3. doi: 10.1288/00005537-197808000-00010
17. Pringle MB, Croft CB. A grading system for patients with obstructive sleep apnoea–based on sleep nasendoscopy. Clin Otolaryngol Allied Sci. (1993) 18:480–4. doi: 10.1111/j.1365-2273.1993.tb00618.x
18. De Vito A, Carrasco Llatas M, Ravesloot MJ, Kotecha B, De Vries N, Hamans E, et al. European position paper on drug-induced sleep endoscopy: 2017 update. Clin Otolaryngol. (2018) 43:1541–52. doi: 10.1111/coa.13213
19. Vanderveken OM, Maurer JT, Hohenhorst W, Hamans E, Lin HS, Vroegop AV, et al. Evaluation of drug-induced sleep endoscopy as a patient selection tool for implanted upper airway stimulation for obstructive sleep apnea. J Clin Sleep Med. (2013) 9:433–8. doi: 10.5664/jcsm.2658
20. Huyett P, Kent DT, D'Agostino MA, Green KK, Soose RJ, Kaffenberger TM, et al. Drug-induced sleep endoscopy and hypoglossal nerve stimulation outcomes: a multicenter cohort study. Laryngoscope. (2021) 131:1676–82. doi: 10.1002/lary.29396
21. Op de Beeck S, Dieltjens M, Verbruggen AE, Vroegop AV, Wouters K, Hamans E, et al. Phenotypic labelling using drug-induced sleep endoscopy improves patient selection for mandibular advancement device outcome: a prospective study. J Clin Sleep Med. (2019) 15:1089–99. doi: 10.5664/jcsm.7796
22. Marques M, Genta PR, Azarbarzin A, Taranto-Montemurro L, Messineo L, Hess LB, et al. Structure and severity of pharyngeal obstruction determine oral appliance efficacy in sleep apnoea. J Physiol. (2019) 597:5399–410. doi: 10.1113/JP278164
23. Torre C, Camacho M, Liu SY, Huon LK, Capasso R. Epiglottis collapse in adult obstructive sleep apnea: A systematic review. Laryngoscope. (2016) 126:515–23. doi: 10.1002/lary.25589
24. Vonk PE, Ravesloot MJL, Kasius KM, van Maanen JP, de Vries N. Floppy epiglottis during drug-induced sleep endoscopy: an almost complete resolution by adopting the lateral posture. Sleep Breath. (2020) 24:103–9. doi: 10.1007/s11325-019-01847-x
25. Van den Bossche K, Van de Perck E, Kazemeini E, Willemen M, Van de Heyning PH, Verbraecken J, et al. Natural sleep endoscopy in obstructive sleep apnea: a systematic review. Sleep Med Rev. (2021) 2021:101534. doi: 10.1016/j.smrv.2021.101534
26. Park D, Kim JS, Heo SJ. Obstruction patterns during drug-induced sleep endoscopy vs natural sleep endoscopy in patients with obstructive sleep apnea. JAMA Otolaryngol Head Neck Surg. (2019) 145:730–4. doi: 10.1001/jamaoto.2019.1437
27. Genta PR, Sands SA, Butler JP, Loring SH, Katz ES, Demko BG, et al. Airflow shape is associated with the pharyngeal structure causing OSA. Chest. (2017) 152:537–46. doi: 10.1016/j.chest.2017.06.017
28. Azarbarzin A, Marques M, Sands SA, Op de Beeck S, Genta PR, Taranto-Montemurro L, et al. Predicting epiglottic collapse in patients with obstructive sleep apnoea. Eur Respir J. (2017) 50:1700345. doi: 10.1183/13993003.00345-2017
29. Azarbarzin A, Sands SA, Marques M, Genta PR, Taranto-Montemurro L, Messineo L, et al. Palatal prolapse as a signature of expiratory flow limitation and inspiratory palatal collapse in patients with obstructive sleep apnoea. Eur Respir J. (2018). 51:1701419. doi: 10.1183/13993003.01419-2017
30. Op de Beeck S, Van de Perck E, Vena D, Kazemeini E, Dieltjens M, Willemen M, et al. Flow identified site of collapse during drug-induced sleep endoscopy: feasibility and preliminary results. Chest. (2020) 159:828–32. doi: 10.1016/j.chest.2020.09.079
31. Marques M, Genta PR, Azarbarzin A, Sands SA, Taranto-Montemurro L, Messineo L, et al. Retropalatal and retroglossal airway compliance in patients with obstructive sleep apnea. Respir Physiol Neurobiol. (2018) 258:98–103. doi: 10.1016/j.resp.2018.06.008
32. Vena D, Azarbarzin A, Marques M, Op de Beeck S, Vanderveken OM, Edwards BA, et al. Predicting sleep apnea responses to oral appliance therapy using polysomnographic airflow. Sleep. (2020) 43:zsaa004. doi: 10.1093/sleep/zsaa004
33. Lee LA, Lo YL Yu JF, Lee GS Ni YL, Chen NH, et al. Snoring sounds predict obstruction sites and surgical response in patients with obstructive sleep apnea hypopnea syndrome. Sci Rep. (2016) 6:30629. doi: 10.1038/srep30629
34. Sebastian A, Cistulli PA, Cohen G, de Chazal P. Automated identification of the predominant site of upper airway collapse in obstructive sleep apnoea patients using snore signal. Physiol Meas. (2020) 41:095005. doi: 10.1088/1361-6579/abaa33
35. Liistro G, Stănescu D, Veriter C. Pattern of simulated snoring is different through mouth and nose. J Appl Physiol (1985). (1991). 70:2736–41. doi: 10.1152/jappl.1991.70.6.2736
36. Xu H, Huang W, Yu L, Chen L. Sound spectral analysis of snoring sound and site of obstruction in obstructive sleep apnea syndrome. Acta Otolaryngol. (2010) 130:1175–9. doi: 10.3109/00016481003694774
37. Gurpinar B, Salturk Z, Kumral TL, Civelek S, Izel O, Uyar Y. Analysis of snoring to determine the site of obstruction in obstructive sleep apnea syndrome. Sleep Breath. (2020) 25:1427–32. doi: 10.1007/s11325-020-02252-5
38. Agrawal S, Stone P, McGuinness K, Morris J, Camilleri AE. Sound frequency analysis and the site of snoring in natural and induced sleep. Clin Otolaryngol All Sci. (2002) 27:162–6. doi: 10.1046/j.1365-2273.2002.00554.x
39. Koo SK, Kwon SB, Moon JS, Lee SH, Lee HB, Lee SJ. Comparison of snoring sounds between natural and drug-induced sleep recorded using a smartphone. Auris Nasus Larynx. (2018) 45:777–82. doi: 10.1016/j.anl.2017.09.005
40. Pevernagie D, Aarts RM, De Meyer M. The acoustics of snoring. Sleep Med Rev. (2010) 14:131–44. doi: 10.1016/j.smrv.2009.06.002
41. Jones TM, Earis JE, Swift AC, Charters P. Acoustic parameters of snoring sound to compare natural snores with snores during 'steady-state' propofol sedation - response. Clin Otolaryngol. (2006). 31:341–2; author reply 342. doi: 10.1111/j.1749-4486.2006.01252.x
42. Penzel T, Sabil A. The use of tracheal sounds for the diagnosis of sleep apnoea. Breathe (Sheff). (2017) 13:e37–45. doi: 10.1183/20734735.008817
43. Sabil A, Glos M, Gunther A, Schobel C, Veauthier C, Fietze I, et al. Comparison of apnea detection using oronasal thermal airflow sensor, nasal pressure transducer, respiratory inductance plethysmography and tracheal sound sensor. J Clin Sleep Med. (2019) 15:285–92. doi: 10.5664/jcsm.7634
44. Mann DL, Terrill PI, Azarbarzin A, Mariani S, Franciosini A, Camassa A, et al. Quantifying the magnitude of pharyngeal obstruction during sleep using airflow shape. Eur Respir J. (2019). 54:1802262. doi: 10.1183/13993003.02262-2018
45. Huyett P, Calianese N, Azarbarzin A, Gell L, Taranto-Montemurro LT, White D, et al. Identifying site of pharyngeal collapse using snore sounds. Amer J Respir Crit Care Med. (2021) 203:A4700. doi: 10.1164/ajrccm-conference.2021.203.1_MeetingAbstracts.A4700
46. Kezirian EJ, Hohenhorst W, de Vries N. Drug-induced sleep endoscopy: the VOTE 565 classification. Eur Arch Otorhinolaryngol. (2011) 268:1233–6.
47. Berry RB, Budhiraja R, Gottlieb DJ, Gozal D, Iber C, Kapur VK, et al. Rules for scoring respiratory events in sleep: update of the 2007 AASM manual for the scoring of sleep and associated events. Deliberations of the sleep apnea definitions task force of the American Academy of Sleep Medicine. J Clin Sleep Med. (2012) 8:597–619. doi: 10.5664/jcsm.2172
48. Berry RB, Brooks R, Gamaldo C, Harding SM, Lloyd RM, Quan SF, et al. AASM scoring manual updates for 2017 (version 2.4). J Clin Sleep Med. (2017) 13:665–6. doi: 10.5664/jcsm.6576
49. Sands SA, Edwards BA, Terrill PI, Taranto-Montemurro L, Azarbarzin A, Marques M, et al. Phenotyping pharyngeal pathophysiology using polysomnography in patients with obstructive sleep apnea. Am J Respir Crit Care Med. (2018) 197:1187–97. doi: 10.1164/rccm.201707-1435OC
50. Terrill PI, Edwards BA, Nemati S, Butler JP, Owens RL, Eckert DJ, et al. Quantifying the ventilatory control contribution to sleep apnoea using polysomnography. Eur Respir J. (2015) 45:408–18. doi: 10.1183/09031936.00062914
51. Sands SA, Terrill PI, Edwards BA, Taranto Montemurro L, Azarbarzin A, Marques M, et al. Quantifying the arousal threshold using polysomnography in obstructive sleep apnea. Sleep. (2018). 41:zsx183. doi: 10.1093/sleep/zsx183
52. Sands SA, Edwards BA, Terrill PI, Butler JP, Owens RL, Taranto-Montemurro L, et al. Identifying obstructive sleep apnoea patients responsive to supplemental oxygen therapy. Eur Respir J. (2018) 52:1800674. doi: 10.1183/13993003.00674-2018
53. Owens RL, Edwards BA, Eckert DJ, Jordan AS, Sands SA, Malhotra A, et al. An integrative model of physiological traits can be used to predict obstructive sleep apnea and response to non positive airway pressure therapy. Sleep. (2015) 38:961–70. doi: 10.5665/sleep.4750
54. Wellman A, Edwards BA, Sands SA, Owens RL, Nemati S, Butler J, et al. A simplified method for determining phenotypic traits in patients with obstructive sleep apnea. J Appl Physiol (1985). (2013). 114:911–22. doi: 10.1152/japplphysiol.00747.2012
55. Vanderveken OM, Beyers J, Op de Beeck S, Dieltjens M, Willemen M, Verbraecken JA, et al. Development of a clinical pathway and technical aspects of upper airway stimulation therapy for obstructive sleep apnea. Front Neurosci. (2017) 11:523. doi: 10.3389/fnins.2017.00523
56. Marques M, Genta PR, Sands SA, Azarbazin A, de Melo C, Taranto-Montemurro L, et al. Effect of sleeping position on upper airway patency in obstructive sleep apnea is determined by the pharyngeal structure causing collapse. Sleep. (2017). 40:zsx005. doi: 10.1093/sleep/zsx005
57. Genta PR, Edwards BA, Sands SA, Owens RL, Butler JP, Loring SH, et al. Tube law of the pharyngeal airway in sleeping patients with obstructive sleep apnea. Sleep. (2016) 39:337–43. doi: 10.5665/sleep.5440
58. Menon A, Kumar M. Influence of body position on severity of obstructive sleep apnea: a systematic review. ISRN Otolaryngol. (2013) 2013:670381. doi: 10.1155/2013/670381
59. Rabelo FA, Braga A, Kupper DS, De Oliveira JA, Lopes FM, de Lima Mattos PL, et al. Propofol-induced sleep: polysomnographic evaluation of patients with obstructive sleep apnea and controls. Otolaryngol Head Neck Surg. (2010) 142:218–24. doi: 10.1016/j.otohns.2009.11.002
60. Fogel RB, Trinder J, White DP, Malhotra A, Raneri J, Schory K, et al. The effect of sleep onset on upper airway muscle activity in patients with sleep apnoea versus controls. J Physiol. (2005). 564(Pt 2):549–62. doi: 10.1113/jphysiol.2005.083659
61. Eckert DJ, Malhotra A, Lo YL, White DP, Jordan AS. The influence of obstructive sleep apnea and gender on genioglossus activity during rapid eye movement sleep. Chest. (2009) 135:957–64. doi: 10.1378/chest.08-2292
62. Hillman DR, Walsh JH, Maddison KJ, Platt PR, Kirkness JP, Noffsinger WJ, et al. Evolution of changes in upper airway collapsibility during slow induction of anesthesia with propofol. Anesthesiology. (2009) 111:63–71. doi: 10.1097/ALN.0b013e3181a7ec68
63. Eastwood PR, Platt PR, Shepherd K, Maddison K, Hillman DR. Collapsibility of the upper airway at different concentrations of propofol anesthesia. Anesthesiology. (2005) 103:470–7. doi: 10.1097/00000542-200509000-00007
64. Kellner P, Herzog B, Plossl S, Rohrmeier C, Kuhnel T, Wanzek R, et al. Depth-dependent changes of obstruction patterns under increasing sedation during drug-induced sedation endoscopy: results of a German monocentric clinical trial. Sleep Breath. (2016) 20:1035–43. doi: 10.1007/s11325-016-1348-6
65. Lo YL Ni YL, Wang TY, Lin TY Li HY, White DP, et al. Bispectral index in evaluating effects of sedation depth on drug-induced sleep endoscopy. J Clin Sleep Med. (2015) 11:1011–20. doi: 10.5664/jcsm.5016
Keywords: diagnosis, endotyping, OSA (obstructive sleep apnea), personalized medicine, upper airway obstruction
Citation: Van den Bossche K, Van de Perck E, Wellman A, Kazemeini E, Willemen M, Verbraecken J, Vanderveken OM, Vena D and Op de Beeck S (2021) Comparison of Drug-Induced Sleep Endoscopy and Natural Sleep Endoscopy in the Assessment of Upper Airway Pathophysiology During Sleep: Protocol and Study Design. Front. Neurol. 12:768973. doi: 10.3389/fneur.2021.768973
Received: 01 September 2021; Accepted: 15 November 2021;
Published: 07 December 2021.
Edited by:
Sairam Parthasarathy, University of Arizona, United StatesReviewed by:
Irma Rukhadze, Albany Medical College, United StatesThomas Penzel, Charité University Medicine Berlin, Germany
Copyright © 2021 Van den Bossche, Van de Perck, Wellman, Kazemeini, Willemen, Verbraecken, Vanderveken, Vena and Op de Beeck. 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: Karlien Van den Bossche, karlien.vandenbossche@student.uantwerpen.be
†These authors have contributed equally to this work