key: cord-0812173-lorcfe5r
authors: Weber, Jonas; Straka, Leonie; Borgmann, Silke; Schmidt, Johannes; Wirth, Steffen; Schumann, Stefan
title: Flow-controlled ventilation (FCV) improves regional ventilation in obese patients – a randomized controlled crossover trial
date: 2020-01-28
journal: BMC Anesthesiol
DOI: 10.1186/s12871-020-0944-y
sha: 9d9a9b7b8ca606bdcc199a39cd19512690fd26d5
doc_id: 812173
cord_uid: lorcfe5r

BACKGROUND: In obese patients, high closing capacity and low functional residual capacity increase the risk for expiratory alveolar collapse. Constant expiratory flow, as provided by the new flow-controlled ventilation (FCV) mode, was shown to improve lung recruitment. We hypothesized that lung aeration and respiratory mechanics improve in obese patients during FCV. METHODS: We compared FCV and volume-controlled (VCV) ventilation in 23 obese patients in a randomized crossover setting. Starting with baseline measurements, ventilation settings were kept identical except for the ventilation mode related differences (VCV: inspiration to expiration ratio 1:2 with passive expiration, FCV: inspiration to expiration ratio 1:1 with active, linearized expiration). Primary endpoint of the study was the change of end-expiratory lung volume compared to baseline ventilation. Secondary endpoints were the change of mean lung volume, respiratory mechanics and hemodynamic variables. RESULTS: The loss of end-expiratory lung volume and mean lung volume compared to baseline was lower during FCV compared to VCV (end-expiratory lung volume: FCV, − 126 ± 207 ml; VCV, − 316 ± 254 ml; p < 0.001, mean lung volume: FCV, − 108.2 ± 198.6 ml; VCV, − 315.8 ± 252.1 ml; p < 0.001) and at comparable plateau pressure (baseline, 19.6 ± 3.7; VCV, 20.2 ± 3.4; FCV, 20.2 ± 3.8 cmH(2)O; p = 0.441), mean tracheal pressure was higher (baseline, 13.1 ± 1.1; VCV, 12.9 ± 1.2; FCV, 14.8 ± 2.2 cmH(2)O; p < 0.001). All other respiratory and hemodynamic variables were comparable between the ventilation modes. CONCLUSIONS: This study demonstrates that, compared to VCV, FCV improves regional ventilation distribution of the lung at comparable PEEP, tidal volume, P(Plat) and ventilation frequency. The increase in end-expiratory lung volume during FCV was probably caused by the increased mean tracheal pressure which can be attributed to the linearized expiratory pressure decline. TRIAL REGISTRATION: German Clinical Trials Register: DRKS00014925. Registered 12 July 2018.

In obese patients, the excessive adipose tissue around the thorax and the visceral organs reduce the functional residual capacity and expiratory reserve volume [1] . Obesity also leads to a low respiratory system compliance, early expiratory alveolar collapse with consecutive atelectasis, increased airway resistance [2] and increased risk for airway closure [3] . All these changes make mechanical ventilation in obese patients prone to respiratory complications [4, 5] .

An emerging ventilation technique to linearize expiratory flow is flow-controlled ventilation (FCV), provided by the new ventilator Evone (Ventinova Medical B.V., Eindhoven, the Netherlands). This device provides a constant positive flow during inspiration and a constant negative flow during expiration. Thereby pressure increases linearly during inspiration [comparable to volume-controlled ventilation (VCV)] and decreases linearly during expiration. Recently, we demonstrated that linearizing the expiratory flow improved lung recruitment, the homogeneity of lung aeration [6, 7] , gas exchange [8] and further attenuated experimental lung injury [9] . Since FCV is a new emerging technique comparative clinical studies in humans, particularly in patients with impaired respiratory system mechanics, are lacking.

We hypothesized that FCV improves regional ventilation distribution of the lung and respiratory system mechanics in obese patients. Therefore, we compared regional ventilation using electrical impedance tomography (EIT) and respiratory system mechanics during FCV and VCV in obese patients in a randomized controlled crossover trial.

The study was approved by the Ethics Committee of the University Medical Center of Freiburg (Engelbergstr. 21, 79106 Freiburg, Germany, Ethical Committee N°179/ 18) on 29th March 2018 (Chairperson Prof. Dr. R. Korinthenberg) and registered at the German Register for Clinical Trials (DRKS00014925). Please note that this study adheres to the CONSORT guidelines.

In order to cope with potential interindividual variability, the study was designed as a randomized controlled interventional crossover trial. After obtaining written informed consent, we studied twenty-three obese patients with body mass index (BMI) ≥ 30 kg•m − 2 . Patients eligible for enrolment were patients with physical status ASA ≤ III undergoing elective bariatric surgery. Exclusion criteria were ASA physical status >III, age < 18 years, pregnancy, emergency procedure, cardiac pacemaker and other active implants, chronic obstructive pulmonary disease classified as GOLD stage > II or refusal to participate. The trial was conducted at the University Medical Center Freiburg, Germany. Participants were enrolled and assigned by a study related anesthesiologist. Data were collected at the University Medical Center of Freiburg, Germany.

After obtaining written informed consent, 23 patients were included in the study. After primary recruitment and preoperative evaluation, the patients received routine monitoring (electrocardiography, SpO 2 , noninvasive blood pressure measurement; Infinity Delta XL, Dräger Medical, Lübeck, Germany) and a 18-20-G intravenous catheter was established. After preoxygenation to an fraction of expired oxygen of 0.8, anesthesia was induced with 0.3-0.5 μg•kg − 1 predicted body weight [10] iv sufentanil (Janssen-Cilag, Neuss, Germany) and 2-3 mg•kg − 1 actual body weight iv propofol (Fresenius Kabi, Bad Homburg vor der Höhe, Germany). Tracheal intubation was facilitated with 0.6 mg•kg − 1 predicted body weight iv rocuronium (Fresenius Kabi). If the patient required a rapid sequence induction, neuromuscular blockage was performed by the administration of 1.0 mg•kg − 1 predicted body weight iv rocuronium. Neuromuscular blockage was monitored with a mechanomyograph (TOFscan; Dräger Medical). For tracheal intubation, we used tracheal tubes with low pressure cuffs (internal diameter of 7.0-7.5 mm for women and 8.0 mm for men; Mallinckrodt Hallo-Contour; Covidien, Neustadt an der Donau, Germany). After adequate placement of the tracheal tube, iv propofol was administered continuously (110-150 μg•kg − 1 •min − 1 ). Potential hypotension (defined as mean arterial pressure < 65 mmHg) was treated with a continuous infusion of iv noradrenaline (0.03-0.2 μg•kg − 1 •min − 1 ). Perioperative volume requirements were addressed with a crystalloid solution (8 ml•kg − 1 •h − 1 , Jonosteril; Fresenius Kabi). According to our local standard, mechanical ventilation was started as volume-controlled baseline ventilation (Fabius Tiro, Dräger Medical) with a tidal volume of 7 ml•kg − 1 predicted body weight, inspiration-to-expiration ratio of 1:2, a positive end-expiratory-pressure (PEEP) of 9 cmH 2 O and ventilation frequency set to maintain an end-tidal carbon dioxide partial pressure between 4.7 and 5.1 kPa. These ventilation settings were based on our study protocol and in accordance with our clinical routine in obese patients. After 7 min of baseline ventilation, all patients were randomly allocated to one of two crossover groups to receive ventilation sequences either VCV-FCV or FCV-VCV for 7 min per ventilation mode. To avoid irritations due to the surgical procedure (e.g. impaired respiratory mechanics by the capnoperitoneum and electrical irritations of the measurement of Electrical Impedance Tomography), our study was performed prior to the surgical intervention. For adequate allocation, a computer generated randomization in blocks was used. Disclosure of the randomization was requested right after induction of anesthesia. A study related anesthesiologist conducted the randomization in blocks, enrolled participants and assigned participants to the interventions. During the study protocol, ventilation variables were kept constant as set during the baseline measurements. To prevent from the risks of extubation and reintubation, FCV was performed by introducing the narrow-bore tracheal tube (Tribute, Ventinova Medical B.V.) into the standard tracheal tube. Blocking the cuff of the Tritube in the lumen of the tracheal tube provided a sufficient seal. By controlling both tube's markings, placement of the Tritube's tip exceeding that of the standard tracheal tube by 2-5 mm was ensured, and the potential risk of bronchial intubation was avoided. Respiratory data were collected from both ventilators via the respective serial communication interface and analyzed offline. Electrical impedance tomography (EIT) was performed with PulmoVista 500 (Dräger Medical) in all patients to measure regional ventilation, changes in relative thoracic electrical impedance during the different ventilation phases, relative end-expiratory lung volume (ΔEELV) and to compare the expiratory decrease in intrapulmonary air [11] [12] [13] .

Ventilation settings during baseline measurements and VCV were identical. In each patient, baseline measurements were performed prior to the intervention. During FCV, patients were ventilated with a constant positive flow during inspiration and a constant negative flow during expiration (Fig. 1 ). To avoid intrinsic PEEP, the intratracheal pressure is monitored continuously via a dedicated pressure measurement lumen of the Tritube. During FCV, the operator is able to adjust the inspiratory flow rate, inspiration to expiration ratio, peak inspiratory pressure, end-expiratory pressure and the inspiratory concentration of oxygen. In this special ventilation mode, there is no direct way to control minute volume via tidal volumes and/or respiratory rate. However, the respiratory rate depends on the peak inspiratory pressure, the set (positive) end-expiratory pressure, the set inspiratory flow rate, the inspiration to expiration ratio and the patient's lung compliance [14] . The (end) expiratory pressure was kept constant in all conditions during the study procedure.

ΔEELV was the primary endpoint of this study. EIT recordings were analyzed using software developed in Matlab (R2014, The MathWorks Inc.). We derived ΔEELV from adjusting end-expiratory impedance changes by tidal volume and tidal impedance changes as described before [7, 11] . As a first step, the lung area estimation method was applied to all EIT recordings to estimate the relevant lung area [15] . Afterwards, global tidal impedance curves were calculated. These curves represent the sum of impedance of all pixels per frame over time. To scale the absolute impedance values to milliliters, the relation between tidal impedance change and tidal volume was used. Changes of the baseline of these tidal impedance curves were determined as estimates for changes of the end-expiratory lung volume. ΔEELV was then calculated as the difference of endexpiratory lung volume during the different ventilation phases [11] . Secondary endpoints were the respiratory system variables: plateau pressure (P Plat ), mean tracheal pressure (P mean ), mean tracheal pressure during expiration (P mean expiration ), peripheral oxygen saturation (SpO 2 ), fraction of inspired oxygen (FiO 2 ) and quasistatic respiratory system compliance (C RS ). To calculate C RS during FCV, the plateau pressure was determined from a short (approximately 0.1 s) end-inspiratory pause. This pause is performed automatically by the Evone ventilator (Ventinova Medical B.V.) with every ten breaths and used to calculate C RS . Non-invasively collected hemodynamic variables included mean systolic blood pressure, mean diastolic blood pressure, mean arterial pressure and heart rate. To compare relative intrapulmonary air distribution, baseline tidal impedance curves for ventral and dorsal lung areas were determined and compared as described before [7, 12] . The differences in mean lung volume (ΔMLV) between baseline ventilation and VCV and FCV were calculated, respectively. Further, the decrease in global thoracic electrical impedance during each ventilation mode was separated into four equal sections (ΔEI 25 , ΔEI 50 , ΔEI 75 and ΔEI 100 ), then matched with the correlating decrease in tidal volume and compared successively.

Pressure data from the Evone are based on direct tracheal pressure measurement via a dedicated lumen of the Tritube. To allow for comparability of pressure data from both ventilators and to calculate quasistatic compliance of the respiratory system, airway pressure data from the Dräger Fabius Tiro were generally converted into tracheal pressure data by calculating the flow dependent pressure drop across the respective tracheal tube and pointwise subtracting this Baseline_VCV_FCV, randomization that was characterized by baseline measurements, followed by volume-controlled ventilation (VCV), followed by flow-controlled ventilation (FCV); baseline_FCV_VCV, randomization that was characterized by baseline measurements, followed by flow-controlled ventilation (FCV), followed by volume-controlled ventilation (VCV); ASA American Society of Anesthesiologists, PBW Predicted body weight, ABW Actual body weight, BMI Body mass index value from airway pressure [16] . Thus all pressure data in the following refer to the respective tracheal pressure.

The datasets used and analyzed during the current study are available from the corresponding author on request. Please note that EIT data files require large memory.

In regard to previous investigations on gas exchange during FCV in a porcine model of ARDS [17] and the crossover design (paired test conditions) we assumed a standardized effect size of the primary endpoint of 0.7 (being the quotient of differences in means and SD). To reach a test power of 0.8 and a desired level of significance of 0.05, 19 patients were required. To compensate for potential incomplete data sets, 23 patients were included in the study. Lilliefors tests were used to confirm that the assumed normal distribution cannot be rejected.

Values are presented as mean ± standard deviation, unless indicated otherwise. Statistical analysis was done using Matlab (R2014, The MathWorks Inc., Natick, MA, USA). Linear mixed effects model analyses were performed to check for differences between respiratory variables and variables resulting from EIT measurements during the ventilation phases using R based software (jamovi project (2018), jamovi (Version 0.9.2.3), retrieved from https://www.jamovi.org). For each measured primary and secondary endpoint (dependent variable), the influence of the ventilation mode (baseline ventilation, VCV and FCV) and the ventilation sequence (baseline-VCV-FCV, baseline-FCV-VCV) (factors) was investigated. P < 0.05 was considered statistically significant.

In total, 23 consecutive patients presenting for elective bariatric surgery were included and 19 complete data sets could be recorded. Patients were recruited from Fig. 3 Exemplary relative global thoracic electrical impedance (EI Th ) of one patient during the study protocol. BL, baseline (volume-controlled) ventilation; VCV, volume-controlled ventilation; FCV, flow-controlled ventilation. The first slope represents the insertion of the Tritube® into the standard tracheal tube. The second slope represents the remove of the Tritube and re-connecting to the Dräger Fabius Tiro ventilator. Arrows indicate the switch between the respective ventilation modes 30th July 2018 to 23rd October 2018. One patient had to be excluded due to limited size of the EIT belt, three other patients due to incomplete data collection (Fig. 2 ). There were no adverse events during the study procedure. The study was ended regularly after the last subject was included. Age, gender, ASA physical status, predicted and actual body weight and BMI were comparable between the two interventional groups (Table 1) .

During mechanical ventilation, end-expiratory lung volume decreased generally (Fig. 3) . ΔEELV between baseline ventilation and FCV (− 126 ± 207 ml) was lower than between baseline and VCV (− 316 ± 254 ml, p < 0.001). ΔMLV between baseline and FCV (− 108 ± 198 ml) was lower than between baseline and VCV (− 315 ± 252 ml, p < 0.001) (Fig. 4) . P mean and P mean expiration was higher during FCV. No significant differences in tidal volume, ventilation frequency, P Plat , SpO 2 and C RS were found between FCV and VCV. All hemodynamic variables were comparable during FCV and VCV (Table 2) .

FCV was characterized by a more even decay of impedance throughout the expiration phase (Fig. 5) . ΔEI 25 , ΔEI 50 , ΔEI 75 and ΔEI 100 showed a more even decrease during FCV compared to VCV (Fig. 6) . ΔEI 25 decreases about 45% during baseline ventilation and VCV and 25% during FCV. ΔEI 50 showed no differences between the ventilation modes. ΔEI 75 and ΔEI 100 showed a lower decrease in global thoracic electrical impedance during baseline ventilation and VCV compared to FCV (Fig. 6) . 

In this study, we compared respiratory system mechanics and regional ventilation in obese patients during short application of FCV and VCV. The main findings of our study are that in obese patients, ΔEELV and mean lung volume decreased less during FCV than during VCVeven with identical respiratory and hemodynamic variables.

These effects were comparable to the effects one would expect from a PEEP increase and/or a tidal volume increase. However, minimal and maximal airway pressure and tidal volume remained unchanged. Our results are consistent with and enlarge upon earlier findings on the implications of a linearized expiratory pressure decrease in lung-healthy patients, lung healthy pigs and a porcine lung-injury model [6] [7] [8] [9] .

We observed the changes of respiratory mechanics during the VCV and FCV phases. Since baseline measurements were performed prior to the following ventilation sequence, we attribute the observed differences between baseline and VCV to the general tendency of the respiratory system to continued derecruitment during mechanical ventilation [18] which may be more pronounced in obese patients. The implications of obesity on respiratory system mechanics are well known: chest wall mechanics are impaired, and respiratory system compliance is reduced. Obese patients have an increased risk for early expiratory alveolar collapse and potential consecutive atelectrauma [3, 5, 17, 19, 20] and thus for decreased functional residual capacity and expiratory reserve volume [5, 13, 21, 22] . Values are stated as mean ± SD. Baseline, baseline measurements (consisting of volume-controlled ventilation); VCV Volume-controlled ventilation, FCV Flowcontrolled ventilation, V T tidal volume, V T PBW tidal volume per predicted body weight, VF Ventilation frequency, P Plat Plateau pressure, P mean mean tracheal pressure, P mean expiration mean tracheal pressure during expiration, C RS Respiratory system compliance, MAP Mean arterial pressure, SpO 2 Peripheral oxygen saturation, FiO 2 Fraction of inspired oxygen, p vent p for ventilation mode* = p < 0.001 for FCV vs. baseline and for # = p < 0.001 for FCV vs. VCV. Linear mixed effect model analyses were used to check for differences between the ventilation phases using R based software (jamovi project 2018, version 0.9.2.3). The randomization had no significant effect on the measured respiratory and hemodynamic variables Fig. 6 ).

Therefore, besides low tidal volume and optional recruitment maneuvers, lung-protective ventilation strategies include the application of adequate PEEP in these patients. However, the ideal adjustment of applied tidal volume, and PEEPwith respect to the potential injurious effects of alveolar overdistensionin obese patients still remain obscure [20] . In this regard, FCV improved lung recruitment without modifying PEEP or tidal volume. Further, it should be noted that because of the controlled endexpiratory pressure during FCV, intrinsic PEEP is nearly excluded. Caused by the active control of the expiration phase, FCV is able to overcome the usually passively driven expiration during conventional ventilation.

The mechanisms behind this recruiting effects may be time dependent: when the lung volume falls below the closing capacity airway closure can occur within the expiration [5, 23] . In obese patients, this airway closure can be observed frequently [3] . The overall delayed expiration during FCV delays the time point at which the lung volume falls below the closing capacity. Consequently, the time until the lung volume exceeds closing capacity within the next inspiration is reduced and thus the risk of airway closure may be lowered [22] . The characterization and correlation between the expiratory decrease in global electrical impedance and expiratory decrease in intrapulmonary air and the increased mean airway pressure during expiration support this conjecture.

Theoretical and clinical observations predict that the linearized decrease in expiratory airway pressure has a beneficial impact on the intrapulmonary inhomogeneity [6, 7, 9, 17, 24] . However, the comparison of tidal impedance variation revealed no differences in intrapulmonary gas distribution during the different ventilation phases. The reduced accessibility of EIT images in obese and morbidly obese patients was described earlier and may be caused by the excessive volume of fat tissue around the chest wall. In horizontal supine position, this fat tissue moves laterally and may create potential shortcuts for the electrical currents of the EIT [25] . Therefore, the resolution of the EIT is limited, which may have masked differences in intrapulmonary inhomogeneity in our patients.

C RS did not differ significantly between the investigated ventilation conditions. Reduced C RS in obese patients may be caused mainly by excess adipose tissue around the chest wall and poor posture caused by thoracic kyphosis and lumbar hyperlordosis, aggravated through excessive abdominal fat tissue [19] . Under these conditions, the recruiting effect of FCV may have influenced C RS only to a minor extent. This hypothesis is supported by earlier investigations that demonstrated that the intratidal recruitment state might differ without affecting C RS [26, 27] . To investigate potential effects of FCV on the C RS in obese patients, longer application of FCV may be necessary.

We did not perform arterial blood gas analyses to examine the effects of FCV on gas exchange in our patients. In preclinical [17] and clinical studies [8] , the controlled expiration improved oxygenation and CO 2 elimination. However, in contrast to other centers, placing an arterial line is not part of our standard treatment in this patient group. Therefore, we felt that such invasive approach was not justified for our study. Further, it should be stated that duration of ventilation of 7 min in each ventilation mode is too short to fully evaluate the effects of FCV on regional ventilation. It follows that further studies are required to investigate the long-term effects of FCV on measurements of regional ventilation, respiratory and hemodynamic variables in obese patients.

This is the first study to investigate the influence of FCV on respiratory mechanics and regional ventilation and in obese and morbidly obese patients. Utilizing measurement of regional ventilation, we could demonstrate that the linearized expiratory flow during FCV provided better maintenance of lung aeration with comparable tidal volume, P Plat and PEEP, compared to VCV. The recruiting effect caused by the linearized expiratory air flow and the elevated P mean during FCV may help prevent atelectasis and hypoxemia during mechanical ventilation in obese patients.

Abbreviations ASA: American Society of Anesthesiologists; BMI: Body mass index; C RS : Quasi-static compliance of the respiratory system; EELV: End-expiratory lung volume; EIT: Electrical impedance tomography; FCV: Flow-controlled ventilation; MLV: Mean lung volume; PEEP: Positive end-expiratory pressure; P mean : Mean airway pressure; P Plat : Plateau pressure; SpO 2 : Peripheral oxygen saturation (pulse oximetry); VCV: Volume-controlled ventilation

The effects of body mass index on lung volumes

Effects of obesity on respiratory resistance

Airway closure during surgical Pneumoperitoneum in obese patients

Perioperative lung protective ventilation in obese patients

The effects of body mass on lung volumes, respiratory mechanics, and gas exchange during general anesthesia

Application of the novel ventilation mode FLow-controlled EXpiration (FLEX): a crossover proof-of-principle study in lung-healthy patients

Dorsal recruitment with flow-controlled expiration (FLEX): an experimental study in mechanically ventilated lung-healthy and lung-injured pigs

Flow-controlled ventilation improves gas exchange in lung-healthy patients-a randomized interventional cross-over study

Flow-controlled expiration: a novel ventilation mode to attenuate experimental porcine lung injury

Gentamicin therapy

Positive end-expiratory pressure-induced changes in end-expiratory lung volume measured by spirometry and electric impedance tomography

Spontaneous breathing during general anesthesia prevents the ventral redistribution of ventilation as detected by electrical impedance tomography: a randomized trial

Physiology of obesity and effects on lung function

Flowcontrolled ventilation during ear, nose and throat surgery: a prospective observational study

Chest electrical impedance tomography examination, data analysis, terminology, clinical use and recommendations: consensus statement of the TRanslational EIT developmeNt stuDy group

Continuous calculation of intratracheal pressure in tracheally intubated patients

Improved lung recruitment and oxygenation during mandatory ventilation with a new expiratory ventilation assistance device: a controlled interventional trial in healthy pigs

Intraoperative positive end-expiratory pressure evaluation using the intratidal compliance-volume profile

Total respiratory system, lung, and chest wall mechanics in sedated-paralyzed postoperative morbidly obese patients

Recruitment maneuvers and positive end-expiratory pressure titration in morbidly obese ICU patients

Observational study of the effect of obesity on lung volumes

Altered respiratory physiology in obesity

Breathing mechanics, dead space and gas exchange in the extremely obese, breathing spontaneously and during anaesthesia with intermittent positive pressure ventilation

Determination of respiratory system mechanics during inspiration and expiration by FLow-controlled EXpiration (FLEX): a pilot study in anesthetized pigs

Positive end-expiratory pressure optimization using electric impedance tomography in morbidly obese patients during laparoscopic gastric bypass surgery

Intraoperative compliance profiles and regional lung ventilation improve with increasing positive endexpiratory pressure

Increasing positive endexpiratory pressure (re-)improves intraoperative respiratory mechanics and lung ventilation after prone positioning

Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations

We would like to thank Dietmar Enk, MD, PhD, inventor of Ventrain®, Evone® and Tritube®, consultant anesthetist (University Hospital Münster (UKM), Dept.of Anesthesiology and Intensive Care Medicine, Münster, Germany), for sharing his ideas on small lumen ventilation with us.

This project has received funding from the European Union's Horizon 2020 research and innovation program under grant agreement No. 691519. This funding did not influence the design, data collection, data analysis, interpretation of data and writing the manuscript of this study. The article processing charge was funded by the German Research Foundation (DFG) and the University of Freiburg in the funding program Open Access Publishing.

The datasets used and analyzed during the current study are available from the corresponding author on request. Please note that EIT data files require large memory. A separate data transfer service will be used to transfer EIT data files. 

Not applicable.Competing interests J.W., L.S., S.B., J.S. and S.W declare no conflicts of interest. S.S. has a consulting contract with Gründler GmbH, Freudenstadt (no relationship to this study).Received: 2 July 2019 Accepted: 16 January 2020