Oxygenation Threshold Derived from Near-Infrared Spectroscopy: Reliability and Its Relationship with the First Ventilatory Threshold

Background Near-infrared spectroscopy (NIRS) measurements of oxygenation reflect O2 delivery and utilization in exercising muscle and may improve detection of a critical exercise threshold. Purpose First, to detect an oxygenation breakpoint (Δ[O2HbMb-HHbMb]-BP) and compare this breakpoint to ventilatory thresholds during a maximal incremental test across sexes and training status. Second, to assess reproducibility of NIRS signals and exercise thresholds and investigate confounding effects of adipose tissue thickness on NIRS measurements. Methods Forty subjects (10 trained male cyclists, 10 trained female cyclists, 11 endurance trained males and 9 recreationally trained males) performed maximal incremental cycling exercise to determine Δ[O2HbMb-HHbMb]-BP and ventilatory thresholds (VT1 and VT2). Muscle haemoglobin and myoglobin O2 oxygenation ([HHbMb], [O2HbMb], SmO2) was determined in m. vastus lateralis. Δ[O2HbMb-HHbMb]-BP was determined by double linear regression. Trained cyclists performed the maximal incremental test twice to assess reproducibility. Adipose tissue thickness (ATT) was determined by skinfold measurements. Results Δ[O2HbMb-HHbMb]-BP was not different from VT1, but only moderately related (r = 0.58–0.63, p<0.001). VT1 was different across sexes and training status, whereas Δ[O2HbMb-HHbMb]-BP differed only across sexes. Reproducibility was high for SmO2 (ICC = 0.69–0.97), Δ[O2HbMb-HHbMb]-BP (ICC = 0.80–0.88) and ventilatory thresholds (ICC = 0.96–0.99). SmO2 at peak exercise and at occlusion were strongly related to adipose tissue thickness (r2 = 0.81, p<0.001; r2 = 0.79, p<0.001). Moreover, ATT was related to asymmetric changes in Δ[HHbMb] and Δ[O2HbMb] during incremental exercise (r = -0.64, p<0.001) and during occlusion (r = -0.50, p<0.05). Conclusion Although the oxygenation threshold is reproducible and potentially a suitable exercise threshold, VT1 discriminates better across sexes and training status during maximal stepwise incremental exercise. Continuous-wave NIRS measurements are reproducible, but strongly affected by adipose tissue thickness.


Introduction
Performance benefits from adequate training. Even though training intensity can be quantified in many different ways, endurance training intensity is often quantified by the lactate thresholds (LT1 and LT2) obtained from blood sampling or the ventilatory thresholds (VT1 and VT2) obtained from gas exchange data [1]. Above LT1 and VT1, blood lactate concentration is elevated [1] and CO 2 production exceeds oxygen uptake due to anaerobic energy production in the muscle cells as bicarbonate buffers the produced lactic acid [2]. Above LT2 and VT2, blood lactate accumulates as lactate production exceeds removal rates [1] and respiratory compensation starts by an excessive increase in minute ventilation [3]. There is consensus that one can obtain ventilatory or lactate thresholds reproducibly from an incremental exercise test [1,4,5] and can use these thresholds for training purposes [1]. However, ventilatory and lactate thresholds are derived from measurements at the mouth or (closer to the muscle) in the blood. Therefore, the first ventilatory and lactate threshold rather indirectly and with a delay reflect an increase of anaerobic ATP resynthesis. Both exercise intensity and oxygen availability contribute to this transition. Increases in the former lead to recruitment of additional (larger) motor units with lower oxidative capacity, while the latter depends on O 2 consumption by the contracting muscle fibres and O 2 supply by the O 2 cascade from the air to the mitochondria [6]. Near-infrared spectroscopy (NIRS) provides non-invasive measures of the balance between oxygen delivery and consumption in the muscle [7], and thereby NIRS is potentially a very suitable technique to detect a critical exercise threshold directly in the exercising muscle.
In active muscle, the matching of oxygen supply ( _ QO 2 ) and oxygen consumption (m _ V O 2 ) is crucial for O 2 diffusion and metabolic control [8]. This _ QO 2 =m _ V O 2 matching can theoretically be assessed non-invasively by NIRS [8]. Hereto, concentration changes of deoxygenated haemoglobin and myoglobin ([HHbMb]), oxygenated haemoglobin and myoglobin ([O 2 HbMb]) and tissue saturation (SmO 2 ) are measured. These NIRS signals present multiple breakpoints during incremental exercise. For instance, the [HHbMb] breakpoint during ramp exercise is suggested to resemble an upper limit of fractional oxygen extraction in the exercising leg muscles during incremental exercise [9][10][11] and approximates VT2 [10,12,11,[13][14][15], maximal lactate steady state [14,16] and critical power [14]. At somewhat lower exercise intensities, a breakpoint in Hb difference ([O 2 HbMb-HHbMb]) occurs during maximal stepwise incremental exercise and shows high correlations with VT1 and the lactate threshold (r>0.88) [17][18][19]. At this [O 2 HbMb-HHbMb] breakpoint, the fractional oxygen extraction supposedly has not yet reached its upper limit in the exercising leg muscles during incremental exercise. Therefore, this breakpoint could reflect a transition in oxygen extraction that is likely attributed to capillary-venular PO 2 reaching its critical value (PO 2crit ) and lactic acidosis inducing a rightward shift of the O 2 Hb dissociation curve (i.e. Bohr effect) [17]. Hence, the [O 2 HbMb-HHbMb] breakpoint may relate to the occurrence of oxygen limitation in the skeletal muscle fibres.
Although reproducibility of ventilatory thresholds [4] and lactate threshold [5] are known to be high, data on reproducibility of NIRS signals [5,[20][21][22] and NIRS derived breakpoints [23] is limited. Only in a single study with sedentary subjects, reproducibility of breakpoints in [O 2 HbMb] has been determined (r = 0.67-0.85) [23]. However, it is yet unknown whether reproducibility of these [O 2 HbMb] breakpoints also applies to trained subjects and to other NIRS derived breakpoints.
At present, the [O 2 HbMb-HHbMb] breakpoint has been determined during cycling exercise in mountain climbers [17], fin swimmers [18] and college students [19]. However, the [O 2 HbMb-HHbMb] breakpoint has not been assessed with subjects differing in training status and sex in one study. It is important to note that adipose tissue thickness (ATT) may vary widely between male and female subjects and across subjects with heterogeneous training status. Also, it is well-known that scattering and absorbance due to adipose tissue affect NIRS signals, reducing absorbance by underlying muscle tissue [24][25][26]. One may reduce these effects of ATT by normalizing oxygenation changes to the oxygenation at peak exercise [10,27], maximal voluntary contraction [26,28] or cuff occlusion [21,26], where the effect of ATT on the amplitude of separate [O 2 HbMb] and [HHbMb] signals can be assessed in the absence of blood volume changes.
The main purpose of this study was to investigate how the Δ[O 2 HbMb-HHbMb] breakpoint relates to ventilatory thresholds during incremental step exercise across sexes and training status. Secondary purposes of the study were to determine the reproducibility of NIRS signals and exercise thresholds and to assess the confounding effect of adipose tissue thickness on NIRS signals. We expected that exercise intensity is similar at the [O 2 HbMb-HHbMb] breakpoint and VT1. In addition, we hypothesized that both the [O 2 HbMb-HHbMb] breakpoint and VT1 differ across sexes and training status due to differences in oxygen supply [29][30][31][32] and oxygen utilization [29,30] within the exercising muscles.

Methods Subjects
Forty subjects participated in the present study: 10 trained male cyclists (CM), 10 trained female cyclists (CF), 11 endurance trained males (EM) and 9 recreationally trained males (RM) as classified by their _ VO 2max [33] (characteristics are summarized in Table 1). Female and male cyclists (CF, CM) were recruited from Dutch cycling teams, and RM and EM were recruited based on training volume (EM >5 hours endurance training per week and RM <3 hours training per week). Note that EM contained 7 cyclists and 4 non-cyclists, since many endurance athletes perform cycling exercise in the Netherlands. Reliability measurements were performed in CM and CF. Prior to participation, experimental procedures, risks and aims of the study were explained and all subjects provided written informed consent. The study was conducted according to the principles of the Declaration of Helsinki and was approved by the ethics committee of the Department of Human Movement Sciences, Vrije Universiteit, Amsterdam, the Netherlands.
After three minutes of seated rest, subjects started at a workload of 1.5WÁkg -1 (85-145W), which was increased with 0.5WÁkg -1 every three minutes (30-50WÁ3min -1 ) until voluntary exhaustion. Trained cyclists performed two maximal incremental tests on separate days to assess reproducibility of NIRS signals and exercise thresholds. The test was terminated if cadence dropped below 80 rpm, despite verbal encouragement. Active cool-down consisted of five minutes cycling at 1 WÁkg -1 (60-100W). Thereafter, an arterial occlusion was performed to assess maximal desaturation of the m. vastus lateralis. Hereto, a pressure cuff was inflated to approximately~300mmHg at the most proximal part of the thigh until maximal desaturation was reached (duration 314 ± 44 s). The cuff was secured by an extra strap that prevented it from unwrapping [34].
Prior to the maximal incremental test, subjects were instructed to avoid strenuous exercise and alcohol consumption within 24 hours before the test and to consume their last meal and caffeinated beverages at least three hours prior to testing. Subsequent measurements were separated by at least one day of rest. Saddle and handle bar position were adjusted to individual preferences and were recorded for subsequent measurements. Environmental circumstances were controlled in a climate-controlled environment (temperature 16.3±0.9°C, relative humidity 46.4±6.4%).

Data collection
Pulmonary measures. Ventilation and gas exchange were analysed breath-by-breath using open circuit spirometry (Cosmed Quark CPET, Cosmed S. R. L., Rome, Italy). Prior to every test, the gas analyser and volume transducer were calibrated according to the manufacturer's instructions.
Near-Infrared Spectroscopy. Relative changes with respect to the start of maximal incremental test were measured for haemoglobin + myoglobin ([HHbMb], [O 2 HbMb]) and muscle saturation (SmO 2 ) by a continuous-wavelength portable NIRS device using spatially resolved spectroscopy (Portamon, Artinis Medical Systems, Arnhem, the Netherlands). Concentration changes of the right leg were determined at 10Hz at the vastus lateralis (VL) muscle belly and calculated from light absorbance at 758nm and 847nm using the modified Lambert-Beer law. The NIRS device was positioned above the upper patella border at 1/3 distance of the patella and greater trochanter, parallel to the longitudinal femur axis. Light emitting optodes were positioned at 30, 35 and 40mm distal to the receiver, sending near-infrared light into the thigh with a corresponding penetration depth of 15-20mm [24]. How much light the muscle tissue can absorb depends on light absorbance and scattering by skin and subcutaneous adipose tissue (e.g. influenced by adipose tissue thickness, temperature, and pigment) and on secured contact of the NIRS device with the skin. Motion artefacts were minimized by fixating the optode position using adhesive tape (Fixomull1) and elastic bandages. Additionally, placement of a light absorbing cloth over the NIRS device minimized the influence from ambient light. Finally, probe position was marked for replication in subsequent measurements. Skinfold. Skinfold thickness of the m. vastus lateralis at the site of the NIRS device was measured (median of three measurements) in seated position using Harpenden skinfold callipers (British Indicators Ltd, Burgess Hill, UK). Adipose tissue thickness was calculated by dividing skinfold thickness by two [25], resembling subcutaneous fat and skin. Subjects showed a wide range in adipose tissue thickness (2.0-11.7 mm), yet participants were excluded if adipose tissue thickness was above 12.5 mm.

Data analysis
Gas exchange, power output (PO) and NIRS data were converted to 1Hz samples using linear interpolation, filtered for extreme values (data points deviating > 2 SD from a local window of 6 data points were replaced by the local average) and time-assembled. Critical exercise thresholds and maximal physiological parameters were derived from 1Hz data using smoothed 30-second moving averages and peak PO was defined as average PO over the last 180 seconds prior to termination of the test. VT1 was derived from respiratory exchange measurements using the Vslope method and ventilatory equivalents [35] and VT2 was derived using the minute ventilation versus _ VCO 2 graph and ventilatory equivalents [3]. Both thresholds were detected by automatic regression and subsequently (re)evaluated by two independent observers. In addition, the difference between concentration changes in [HHbMb] and [O 2 HbMb] were used to derive an oxygenation breakpoint (Δ[O 2 HbMb-HHbMb]-BP) using a double linear regression method [17][18][19]. Both linear regressions were fitted through Δ[O 2 HbMb-HHbMb] versus relative PO (%PO) plots. Δ[O 2 HbMb-HHbMb]-BP was defined as the intercept of two congregating regression lines with the least combined residuals sum of squares (Eq 1). Warm-up stages were excluded and regression lengths were set at !10% of relative workload.
In addition, the confounding effect of ATT on NIRS signals was determined for SmO 2 at the end of maximal exercise and arterial occlusion. The amplitudes of [HHbMb] and [O 2 HbMb] were assessed during rest and at the end of maximal exercise and arterial occlusion to determine the symmetry of these NIRS signals in relation to ATT. Moreover, next to the SmO 2 values provided by the NIRS apparatus, SmO 2 was also normalized to maximal changes during exercise (Eq 2). All data processing procedures were performed using custom written software in Matlab (Mathworks Co, Natick, Massachusetts, USA).
Statistics All data are presented as individual values or as mean±SD, unless indicated otherwise. Differences in _ VO 2 and PO between Δ[O 2 HbMb-HHbMb]-BP and ventilatory thresholds were assessed by paired samples t-test and Pearson's product-moment correlation analyses. A priori sample size calculations have been performed in GPower 3.1.2 (University Kiel, Germany) using the t-tests means (difference between two dependent means, matched pairs, a priori) with significance level 0.05 and power 80%. These calculations showed that differences in the Δ[O 2 HbMb-HHbMb]-BP and VT1 could be detected with a sample size of 9 subjects for PO (W) and 4-8 subjects for _ VO 2 (L/min) [18,19]. Differences between groups were assessed by one-way ANOVAs (between factor group: CM-CF-EM-RM). Moreover, reproducibility was determined by intra class-correlations (i.e. single measures ICC 3,1 ) and within-subject coefficient of variation (CV) after logarithmic transformation of the data. Qualification of correlation coefficients was performed according to Evans [36]. Relationships between ATT and SmO 2 were assessed by hyperbolic regressions and relationships between ATT and the symmetry of changes in [HHbMb] and [O 2 HbMb] amplitude were assessed by linear regressions. Presented R-squared values were adjusted for the number of predictors (i.e. coefficients) in the regression model, being a more conservative measure of explained variance. Differences were considered to be significant if p < 0.05.

Adipose tissue
SmO 2 values during exercise and occlusion were higher in groups with larger adipose tissue thickness (Table 2). Accordingly, ATT explained~80% of SmO 2 variance at peak exercise intensity and at the end of arterial occlusion (Fig 3). However, when SmO 2 was normalized to correct for ATT, at VT1 (75±15.3%) and VT2 (89±8.5%) group differences disappeared. During maximal incremental exercise and occlusion, we observed asymmetrical changes in the amplitude of [ as well as for maximal incremental exercise this asymmetry was significantly related to ATT (r = -0.50, p<0.05 and r = -0.64, p<0.01 respectively; Fig 3).  [18,19] and LT1 [18] or that Δ[O 2 HbMb-HHbMb]-BP was not different from LT1 [17,19]. In addition, we showed that the correlation between Δ[O 2 HbMb-HHbMb]-BP and ventilatory thresholds was rather moderate (r = 0.58-0.68). This deviates from previous findings in studies using similar methodology, demonstrating an excellent relationship between Δ[O 2 HbMb-HHbMb]-BP and VT1 or lactate threshold during stepwise incremental cycling exercise (r>0.88) [17][18][19]. A possible explanation for these contrasting findings is that step increments (30-50W) in the present study were larger than increments in previous studies (30W) [17][18][19]. Even though the regression model quality is high (r 2 = 0.97), this could have led to more stepwise kinetics affecting Δ[O 2 HbMb-HHbMb]-BP's goodness of fit (illustrated in Fig 1). Moreover, measurement  [17][18][19], which may be due to a more heterogeneous subject group in the present study. This heterogeneous group could have resulted in larger inter-individual variation in NIRS patterns during exercise. Note that we observed better relationships between Δ[O 2 HbMb-HHbMb]-BP and ventilatory thresholds in a more homogeneous group of trained cyclists (r = 0.68-0.84) as opposed to a more heterogeneous group of endurance and recreationally trained males (r = 0.48-0.50). Physiological differences between sexes are reflected by differences in Δ[O 2 HbMb-HHbMb]-BP and ventilatory thresholds. The absolute exercise intensity at VT1 and Δ [O 2 HbMb-HHbMb]-BP (reflecting when aerobic ATP resynthesis can no longer match ATP use in the working muscles) was lower in females compared to males. The lower absolute exercise intensity in females may be explained by a lower oxygen carrying capacity and lower cardiac output [29,30,32] and possibly by smaller arteriovenous oxygen differences [29] in females. Expressed in percentage peak power output, Δ[O 2 HbMb-HHbMb]-BP and VT1 were found to be similar for trained male and female cyclists, indicating that relative to PO peak the _ QO 2 =m _ V O 2 matching and consequent occurrence of increased anaerobic energy production in the muscle are similar between sexes. Previously, VT1 (in % _ V O 2max ) has also shown to be similar between recreationally trained males and females [10]. However, at a higher intensity exercise, the Δ[HHbMb]-BP (in % _ V O 2max ) has demonstrated to be lower in females compared to males, but is likely explained by less sufficient blood distribution towards active muscle in females [10].
Physiological differences between males with different training status did become apparent by differences in ventilatory thresholds, but not by differences in Δ[O 2 HbMb-HHbMb]-BP. In contrast, it was expected that exercise intensity (absolute and in %PO peak ) at both ventilatory thresholds and Δ[O 2 HbMb-HHbMb]-BP would be higher in trained males, because of a higher gross efficiency in trained subjects [37] and since endurance training enhances oxidative metabolism (e.g. increases in mitochondrial density, oxidative enzymes and percentage slowtwitch myosin heavy chain types; [29,30]) and oxygen supply to the mitochondria (e.g. stroke volume, cardiac output, capillary density, haematocrit, Ca-vO 2 and O 2 diffusion; [29,30]). Also, subjects with a higher aerobic fitness level have been reported to demonstrate a higher percentage of slow-twitch fibres, which has been associated with more effective _ QO 2 =m _ V O 2 matching [31]. However, our results on Δ[O 2 HbMb-HHbMb]-BP did not confirm these expectations. Possibly, in the present study, physiological differences in training status between male groups were too small and/or measurements of Δ[O 2 HbMb-HHbMb]-BP were not sensitive enough to quantify differences in exercise intensity at Δ[O 2 HbMb-HHbMb]-BP. Although the Δ[O 2 HbMb-HHbMb] breakpoint is potentially a suitable exercise threshold revealing when anaerobic energy production starts to increase in the muscle, VT1 being a rather indirect measure of these changes in energy status of the muscle, discriminates better across sexes and training status in the present study.

Reproducibility
Saturation measured by our NIRS device was obtained in trained cyclists and showed high reproducibility in rest, during exercise and at the end of arterial occlusion. These findings are similar to previously reported reproducibility values of SmO 2 measurements during maximal incremental exercise (ICC = 0.81-0.95) [20], at _ V O 2max (r = 0.99) [5] and following arterial arm or leg occlusion (ICC = 0.95-0.96) [21,22]. Thus, these results imply that in trained cyclists one is able to measure oxygenation reproducibly with a continuous-wave NIRS device.
Previous results [23] showed that in sedentary subjects, breakpoints in [O 2 HbMb] are detected reproducibly during ramp exercise (r = 0.67-0.85, p<0.05). The present study shows that also in a group of trained cyclists the reproducibility was high for the Δ[O 2 HbMb-HHbMb] breakpoint (ICC = 0.80-0.88). Note that in contrast to other studies that assessed an oxygenation threshold from NIRS signals [27,[38][39][40][41], the present study did not exclude any of the participating subjects from analysis. In addition, our assessment of ventilatory thresholds demonstrates excellent reproducibility (ICC = 0.96-0.99). It is well-known that ventilatory thresholds can be obtained with good reproducibility (r = 0.91-0.98) and can accurately predict endurance time trial performance [4,5]. Our reproducibility results showed that systematic measurement errors related to visual inspection of ventilatory thresholds are smaller than the measurement errors related to detection of Δ[O 2 HbMb-HHbMb]-BP, which may be due to effects of stepwise kinetics on the double linear regression model or effects of motion artefacts on the combined Δ[O 2 HbMb-HHbMb] signal. Hence, our results indicated that measurement errors related to the reproducibility of NIRS and respiratory measurements only marginally account for the unexplained variance in the relation between Δ[O 2 HbMb-HHbMb]-BP and VT1.

Adipose tissue thickness
ATT is known to affect the amplitude of NIRS measurements [7,[24][25][26]42]. Essentially, the NIRS device assumes one homogeneous medium for extending the Lambert-Beer law of light attenuation [43]. This assumption is violated by light scattering and absorption by adipose tissue [7]. Even though the VL adipose tissue thickness was below 12.5mm in all subjects of the present study, ATT comprised 11.0-64.7% of the measurement depth [24,44]. Our results confirmed that SmO 2 values were substantially affected by ATT, explaining~80% of the variance in SmO 2 at peak exercise and arterial occlusion. These results indicated that NIRS saturation measurements were largely determined by ATT. Although frequently assessed by linear regressions, the relationship between NIRS signals and ATT is better described hyperbolically [22,25,26]. Clearly, the hyperbola (Fig 3) showed that SmO 2 cannot be calculated for extremely low ATT values (i.e. below~2 mm). This finding corresponds with our observation that the NIRS receiver detects no light in extremely lean subjects as the light is fully absorbed or scattered by underlying muscle tissue. Note that in extremely lean subjects the skinfold thickness predominantly consists of skin tissue. In addition to adipose tissue, cutaneous dilation in response to increased temperature will likely affect NIRS measurements as well [45]. Furthermore, lean subjects (trained males) showed SmO 2 values that were roughly in line with mixed venous saturation values reported in literature at rest (~70%), VT1 (~35%) and peak exercise (15-30%) [46][47][48]. Thus, NIRS saturation measurements (SmO 2 ) that are obtained from a combination of [HHbMb] and [O 2 HbMb] signals are clearly affected by ATT.
One may reduce effects of ATT on the amplitude of [HHbMb] and [O 2 HbMb] signals by normalizing oxygenation changes to the oxygenation at peak exercise [10,27], maximal voluntary contraction [26,28] or cuff occlusion [21,26]. Still, this correction does not account for the evident asymmetry in Δ[O 2 HbMb] and Δ[HHbMb] amplitude. Asymmetrical changes during exercise may partially be explained by blood volume changes [26,49], increased haemoglobin concentration (i.e. haemoconcentration) [50] or scattering alterations [43]. For example, continuous-wave NIRS device incorrectly assumes fixed scattering coefficients and therefore can overestimate changes in [HHbMb] (i.e. short wavelengths <750nm) during exercise [43,51]. However, the asymmetry during arterial occlusion may not be explained by changes in blood volume and haemoglobin concentration (both assumed rather constant during occlusion) or scattering coefficients (contrastingly Δ[O 2 HbMb] was generally larger than Δ[HHbMb]). Recent findings in canine muscle clearly showed that [O 2 HbMb] and [HHbMb] signals change symmetrically when no skin and adipose tissue layer is present [52]. In the present study, the asymmetry during incremental exercise as well as occlusion was significantly related to ATT.

Conclusions
During maximal incremental step exercise, the Δ[O 2 HbMb-HHbMb] breakpoint is reproducible and coincides with VT1, but is only moderately related to the ventilatory thresholds. Although the Δ[O 2 HbMb-HHbMb] breakpoint is potentially a suitable exercise threshold revealing when anaerobic energy production starts to increase in the muscle, the VT1 (being a more indirect measure of these changes in muscle energy status) shows higher reproducibility and discriminates better across sexes and training status. Continuous-wave NIRS measurements are reproducible, but adipose tissue thickness strongly affects the amplitude of [O 2 HbMb] and [HHbMb] signals to a different extent and thereby may also affect the kinetics of combined [O 2 HbMb] and [HHbMb] signals. The large inter-individual variability observed at the Δ[O 2 HbMb-HHbMb] breakpoint and its underlying nature is yet poorly understood, although some of this variability may be explained by differences in adipose tissue thickness. Correcting the Δ[O 2 HbMb-HHbMb] breakpoint for ATT however is difficult, since it remains to be established how ATT affects the absorbance by O 2 HbMb and HHbMb at the measured wavelengths. Therefore, the first ventilatory threshold is favoured over the Δ[O 2 HbMb-HHbMb] breakpoint in context of determining an exercise threshold that reflects the metabolic demands of active muscle in maximal stepwise incremental exercise.
Supporting Information S1 File. Supporting excel-file with data of the present study. (XLSX)