Abstract
The use of lung-protective ventilation (LPV) strategies may minimize iatrogenic lung injury in surgical patients. However, the identification of an ideal LPV strategy, particularly during one-lung ventilation (OLV), remains elusive. This study examines the role of ventilator management during OLV and its impact on clinical outcomes.
Data were retrospectively collected from the hospital electronic medical record and the Society of Thoracic Surgery database for subjects undergoing thoracic surgery with OLV between 2012 and 2014. Mean tidal volume (VT) during two-lung ventilation and OLV and ventilator driving pressure (ΔP) (plateau pressure − positive end-expiratory pressure [PEEP]) were analyzed for the 1,019 cases that met the inclusion criteria. Associations between ventilator parameters and clinical outcomes were examined by multivariate linear regression.
After the initiation of OLV, 73.3, 43.3, 18.8, and 7.2% of patients received VT greater than 5, 6, 7, and 8 ml/kg predicted body weight, respectively. One hundred and eighty-four primary and 288 secondary outcome events were recorded. In multivariate logistic regression modeling, VT was inversely related to the incidence of respiratory complications (odds ratio, 0.837; 95% CI, 0.729 to 0.958), while ΔP predicted the development of major morbidity when modeled with VT (odds ratio, 1.034; 95% CI, 1.001 to 1.068).
Low VTper se (i.e., in the absence of sufficient PEEP) has not been unambiguously demonstrated to be beneficial. The authors found that a large proportion of patients continue to receive high VT during OLV and that VT was inversely related to the incidence of respiratory complications and major postoperative morbidity. While low (physiologically appropriate) VT is an important component of an LPV strategy for surgical patients during OLV, current evidence suggests that, without adequate PEEP, low VT does not prevent postoperative respiratory complications. Thus, use of physiologic VT may represent a necessary, but not independently sufficient, component of LPV.
Abstract
Analysis from 1,019 patients undergoing one-lung ventilation indicated that low tidal volume in the presence of low positive end-expiratory pressure is associated with increased pulmonary complications. This suggests that low tidal volume during one-lung ventilation is protective only when accompanied by adequate positive end-expiratory pressure.
Low tidal volume is an important component of protective ventilation and may minimize lung injury during surgery, but the optimal combination of tidal volume and positive end-expiratory pressure (especially during one-lung ventilation) is unknown.
Analysis from 1,019 patients undergoing one-lung ventilation indicated that low tidal volume in the presence of low positive end-expiratory pressure is associated with increased pulmonary complications. This suggests that low tidal volume during one-lung ventilation is protective only when accompanied by adequate positive end-expiratory pressure.
MECHANICAL ventilation is a necessary supportive therapy for critically ill patients and those undergoing major surgeries. However, phasic lung expansion under positive pressure subjects the lungs to a variety of potentially injurious stimuli, which can ultimately result in clinically significant ventilator-induced lung injury (VILI). Historically, approaches to intraoperative mechanical ventilation focused primarily on preventing intraoperative atelectasis and thus endorsed the use of high tidal volumes (VTs).1 More recent experimental and clinical studies have demonstrated that a high VT approach to ventilator management can be injurious. The demonstration that high VT ventilation resulted in significantly higher mortality in patients with acute respiratory distress syndrome (ARDS)2 led to the concept of “protective” ventilation strategies—that is, limiting alveolar overdistension through the application of smaller physiologic VT, and this approach represents the standard of care for critically ill patients with lung injury. Subsequent studies of high VT ventilation in surgical patients at risk for lung injury confirm that a similar approach may decrease systemic and pulmonary inflammation3,4 and improve postoperative pulmonary function3,5 and clinical outcomes including pulmonary complications6 and hospital stay.6
Patients undergoing thoracic surgery may be at increased risk for complications as a result of preexisting disease processes, the nature of the planned surgery, loss of functional lung parenchyma (pulmonary resection), and the detrimental effects of mechanical ventilation, particularly one-lung ventilation (OLV). Thus, these patients could potentially derive even greater benefit from the application of protective ventilation principles. Despite significant advances in our understanding of protective ventilation in patients subjected to two-lung ventilation (TLV), considerably less evidence is available to guide management of OLV, a technique commonly used to optimize operating conditions for thoracic surgery. In this study, we hypothesized that large VTs and higher driving pressures during OLV are associated with an increased risk of postoperative pulmonary complications and overall morbidity after thoracic surgery.
Materials and Methods
This study was approved by the University of Virginia Institutional Review Board for Health Sciences Research. Data were retrospectively collected from the hospital electronic medical record and the Society of Thoracic Surgeons (STS) database for all patients undergoing thoracic surgery with OLV between January 1, 2012, and June 30, 2014. Cases were excluded for the following reasons: age less than 18 yr and incomplete height and/or weight data. A total of 1,232 cases met the initial inclusion criteria. Additionally, we excluded cases due to reoperation, incomplete data, or erroneous duplication of case records. The final study cohort consisted of 1,019 cases. Patient demographics and characteristics included age, sex, height, weight, and medical diagnoses. Ventilator parameters recorded during TLV and OLV included VTs, PEEP, plateau ventilator pressures (Pplat), inspired oxygen fraction, expiratory carbon dioxide, and respiratory rate. Ventilation parameters for TLV were recorded for a 5-min epoch beginning 10 min before initiation of OLV. For OLV, the 5-min epoch began 10 min after the initiation of OLV. Standard practice during the study period included initiation of OLV well before chest wall incision. Thus, the transition from TLV to OLV in the current study generally reflects a transition to OLV well before opening of the pleura. Also, since the 5-min epoch started 10 min after the initiation of OLV, the nonventilated lung can be considered to be at least partially collapsed. VT was calculated on the basis of actual body weight (ABW) and predicted body weight (PBW). PBW was calculated as follows: PBW for males = 50 kg + 2.3 kg × (Height [in] − 60); PBW for females = 45.5 kg + 2.3 kg × (Height [in] − 60).
The STS database was used to obtain information for risk prediction based on previously published thoracic surgery risk models7–9 and for specific information on postoperative outcomes. Definitions of risk predictors and specific outcome events are as specified by the STS (STS GTSD Version 2.3, updated January 2015) and are available via the following link: http://www.sts.org/sites/default/files/documents/STSThoracicDataSpecsV2_3.pdf (accessed March 21, 2016). The candidate risk predictor—major preoperative morbidity—was defined to include any patient with a preoperative diagnosis of coronary artery disease, congestive heart failure, peripheral vascular disease, and/or diabetes mellitus.
The primary outcome was respiratory complications, including tracheostomy, empyema requiring treatment, pneumonia, reintubation, initial ventilator support greater than 48 h, ARDS, bronchopleural fistula, pulmonary embolism, air leak greater than 5 days, atelectasis requiring bronchoscopy, and respiratory failure. The secondary complication was overall postoperative morbidity, including all respiratory complications listed above and major nonrespiratory complications such as unexpected return to the operating room, atrial or ventricular dysrhythmias requiring treatment, myocardial infarction, sepsis, renal failure, central neurologic event, unexpected intensive care unit admission, and anastomotic leak.
Driving pressure (ΔP) and static compliance (Cs) were defined and calculated as follows: ΔP = Pplat − PEEP; Cs = VT/(Pplat − PEEP).
Statistical Analysis
Summary and descriptive statistics were obtained for all clinical data. VT during TLV and OLV were calculated relative to ABW and PBW.
Categorical variables were compared using chi-square tests, and continuous variables were compared across patient groups using means and the two-sample t test. Changes within patients from TLV to OLV were compared using paired t tests for continuous variables and McNemar test for discrete variables. Associations between continuous variables were assessed using Pearson correlation.
The association of primary and secondary outcomes with patient and surgical characteristics was assessed using multivariate logistic regression. The regression models were constructed with variables identified a priori, based on previous studies,7–9 as significant predictors of adverse outcomes after major thoracic surgeries. Since there were fewer primary than secondary outcome events, adherence to the statistical rule allowing 10 outcome events per predictor mandated a reduced set of predictors in the primary outcome models.10 Forced expiratory volume 1-s (FEV1) data were missing in 37.4% of cohort patients. Since spirometric testing data, including FEV1, are more likely to be available in patients with advanced lung disease and/or in those undergoing major pulmonary resection surgery, the absence of such data may itself represent a meaningful predictor. To control for potential confounding, we added a “FEV1 missing” indicator variable and replaced the missing FEV1 data with zeros. This approach leaves the FEV1 coefficient identical to the condition in which data were restricted to only patients with nonmissing FEV1 and separately models risk for those patients with missing FEV1 data.
Area under the receiver operating characteristic curve was used to evaluate the regression models’ sensitivity and specificity for predicting outcome. Differences were considered significant at P values less than 0.05. All statistical analyses were performed in R (version 3.1.3; R Foundation for Statistical Computing, Austria. Available at: https://www.R-project.org/. Accessed March 21, 2016).
Results
Baseline patient characteristics and procedural frequencies are shown in table 1. Figure 1 depicts the derivation of the final study cohort. Table 2 depicts ventilation parameters during TLV and OLV periods. Mean VT during TLV was 6.0 ml/kg ABW and 7.3 ml/kg PBW. After initiation of OLV, the mean VT decreased from 6.0 to 4.9 ml/kg ABW (P < 0.001) and from 7.3 to 5.9 ml/kg PBW (P < 0.001). Pplat increased from 17.7 to 21.7 cm H2O (P < 0.001), end-tidal carbon dioxide increased from 37.4 to 39.5 mmHg (P < 0.001), and the frequency of cases utilizing PEEP greater than or equal to 5 cm H2O increased from 32.5 to 46.6% (P < 0.001). Static compliance decreased during this transition from 35.6 to 23.2 ml/cm H2O (P < 0.001).
Mean PEEP values during OLV for the cohort was 4.2 cm H2O (SD, 1.6 cm H2O), and this did not significantly differ between patients who developed complications (4.2 cm H2O) and those who did not (4.2 cm H2O).
As a function of ABW, 42.4% of this cohort received VT greater than 5 ml/kg during OLV, 19.4% received greater than 6 ml/kg, 6.6% greater than 7 ml/kg, and 2.3% greater than 8 ml/kg. As a function of PBW, the percentage of patients receiving VT greater than 5, 6, 7, and 8 ml/kg during OLV were 73.3, 43.3, 18.8, and 7.2%, respectively (data not shown).
The frequency and types of major morbidities are shown in table 3. Multivariate logistic regression models for the primary and secondary outcomes are shown in table 4. Generally, the patient and procedural risk predictors were consistent with previously published models in the surgical literature. Risk predictors for the primary outcome included Zubrod score, preoperative FEV1, thoracotomy incision, and segmentectomy/lobectomy procedures. Significant predictors for the secondary outcome include patient age, Zubrod score, induction chemotherapy and/or radiation, FEV1, thoracotomy incision, esophageal surgery, segmentectomy and lobectomy procedures, and blood product transfusion.
VT was inversely related to the incidence of the primary outcome—respiratory complications (odds ratio [OR], 0.837; 97.5% CI, 0.729 to 0.958). Thus, an increase in VT of 1 ml/kg PBW was associated with approximately a 16% reduction in the risk of respiratory complications. The relationship between VT and major morbidity was qualitatively similar and approached, but did not reach, statistical significance (P = 0.08). The inverse relationship between VT and the primary outcome was also seen in model 2 (OR, 0.854; 97.5% CI, 0.748 to 0.973), which excluded ΔP and included PEEP, although PEEP itself was not predictive for development of either outcome. Ventilator ΔP was identified as a risk factor for the development of the secondary outcome (OR, 1.034; 97.5% CI, 1.001 to 1.068) in model 1. Each unit of driving pressure (1 cm H2O) was associated with a 3.4% increase in the risk of major morbidity. When modeled without VT (model 3), the effect of ΔP approached, but did not reach, significance (P = 0.084). Receiver operating characteristic curve analysis of the models for the primary and secondary outcomes (model 1) yielded area under the curve of 0.78 and 0.79, respectively, indicating moderate ability to predict complications.
Derivation of the study cohort and exclusions are shown in figure 1. Static compliance during OLV and delivered VT were significantly (P < 0.001) but not strongly correlated (fig. 2), yielding a correlation coefficient of 0.467. To assess whether the relationship between low VT and increased risk of primary and secondary outcome events was independent of compliance, multivariate models were recalculated including compliance. Compliance was not itself significantly predictive of outcomes and did not affect the association of either VT with the primary outcome or ΔP with the secondary outcome (data not shown). The correlation between driving pressure and VT (r = 0.126) is depicted in figure 3. The relationships between VT or ΔP and the log odds of developing respiratory complications or major postoperative complications are depicted in figures 4 and 5.
Discussion
In this study cohort, we found that, during OLV, patients received higher VT when calculated on the basis of PBW as compared to ABW. In multivariate logistic regression models, VT was inversely related to the incidence of the primary outcome—respiratory complications. After controlling for other candidate risk predictors, the odds of developing respiratory complications after thoracic surgery were lower in patients ventilated with higher VT. Ventilator ΔP was identified as a risk factor for the development of major morbid complications in a regression model inclusive of VT. However, when modeled without VT, the effect of ΔP approached, but did not reach, statistical significance.
Protective Ventilation in TLV and OLV
Studies of general surgery patients at high risk for postoperative respiratory complications have reported that lung-protective ventilation (LPV) strategies (lower VTs and higher PEEP) improve clinical outcomes in patients with preoperatively healthy lungs.5,6 However, these studies have been criticized for the use of a control group with high VTs without the use of PEEP (zero end-expiratory pressure [ZEEP]). Although these studies may contribute to the identification of a protective ventilation strategy, they do not permit the elucidation of specific factors (VT, PEEP, airway pressure, driving pressure, and transpulmonary pressure) responsible for the observed clinical effect. More recent trials have evaluated the role of VT and PEEP in lung protection via controlled modification of a single ventilator variable. Treschan et al.11 demonstrated no differences in postoperative lung function in upper abdominal surgery patients ventilated with a high (VT 12 ml/kg PBW) or low (VT 6 ml/kg PBW) VT strategy that included comparable levels of PEEP (5 cm H2O) in both groups.
In thoracic surgery patients, OLV contributes to the development of lung injury and the development of serious complications. Very little data exist to specifically support a particular approach to management of OLV with regard to clinical outcomes. Some but not all12 prospective studies examining putative protective OLV (reduced VT and moderate PEEP) have demonstrated a reduction in pulmonary13 or systemic inflammation,3 extravascular lung water,14,15 or pulmonary complications.13,16 A retrospective study following institution of a protective ventilation protocol incorporating reduced VT, increased PEEP, limited ventilator pressures, and recruitment maneuvers during OLV for lung cancer surgery is also consistent with a reduced risk of acute lung injury.17 In retrospective reviews of pneumonectomy, VTs,17,18 ventilation pressures,19,20 and duration of OLV19 have been identified as risk factors for the development of lung injury.
It is important to note that the results of the current study complement rather than contradict those of recent high-quality prospective trials. No studies have yet unambiguously demonstrated a specific advantage of low VT ventilation in the absence of other ventilatory strategies (PEEP, airway pressure limitation, and recruitment maneuvers), and it is not yet clear which ventilator parameters, if any, are most likely to predict adverse outcomes. That a low VT regimen per se is not inherently protective is supported by a high-quality study of surgical patients, which demonstrated no difference in postoperative pulmonary function between high and low VT regimens, both with equivalent moderate levels of PEEP11 ; a large retrospective study of surgical patients, which demonstrated an association between low intraoperative VT with minimal PEEP and the subsequent risk of mortality21 ; and the results of the present study, in which we found a similar inverse relationship between VT during OLV and the risk of complications after thoracic surgery. In our study cohort, as in the study by Levin et al.,21 patients received low PEEP (4.2 cm H2O), which may have been insufficient to stabilize alveoli, reduce alveolar strain, and prevent atelectasis. While atelectasis is a significant consideration in all anesthetized surgical patients, it may be of greater importance during OLV, due to the use of higher inspiratory oxygen fractions (absorption atelectasis) and the greater potential for dependent lung compression (compression atelectasis). Without adequate PEEP, low VT during OLV may predispose to atelectasis and thus contribute to an increased risk of morbidity. Perioperative atelectasis results in a profoundly injurious inflammatory state and has been extensively reviewed.22 The presence of nonlobar atelectasis appears to act as a “stress concentrator” and can even cause inflammation and alveolar injury in adjacent healthy lung parenchyma.23 Furthermore, the presence of atelectatic regions may promote bacterial translocation24 and increase the risk of pneumonia.24–26
That atelectasis contributes to the increased risk of morbidity in surgical patients ventilated with both low VT and low PEEP is supported by studies that demonstrate that the addition of significant PEEP is required to promote lung protection, improve lung function, and reduce the risk of postoperative complications during both TLV5,6 and OLV.3,13,16 The administration of either high PEEP27 or low VT11 as a sole variable between groups has not been shown to be protective. The synergistic interaction between VT and PEEP has been described in experimental models that demonstrate an interdependence between these variables.28,29 Moreover, although the use of all studied ventilation regimens resulted in substantial release of cytokines in an animal model of lung injury, the combination of low VT and ZEEP resulted in a much higher mortality relative to groups that received either higher VT or PEEP,30 suggesting that the nature, and perhaps severity, of lung injury resulting from atelectasis versus overdistension may be qualitatively different. The ideal amount of PEEP and the ideal approach to titrating PEEP in surgical patients have not yet been elucidated. However, the use of pulmonary mechanical parameters to identify intratidal recruitment during TLV suggests that 5 cm H2O PEEP may be insufficient to prevent this injurious process.31 Likewise, titrating PEEP on the basis of static compliance during OLV resulted in the selection of a much higher PEEP (10 cm H2O) than is conventionally used by most practitioners.32 Collectively, these studies indicate that an effective LPV strategy results from the use of a physiologic (low) VT and sufficient PEEP to minimize both overdistension and atelectasis and suggest that the level of PEEP utilized in the current OLV study cohort and that of a previously reported TLV cohort21 were insufficient to prevent the pathophysiologic processes and sequelae of atelectasis and tidal recruitment/derecruitment.
An alternative explanation for the results of both studies is that unidentified patient factors, which may lead anesthesia care providers to use lower VT, might be independently linked to adverse outcomes. That the use of lower VT in their cohort might be caused by reduced compliance and higher airway pressures was excluded by Levin et al.21 based on their inclusion of compliance and inspiratory pressure in the multivariate model. Our results paralleled those of Levin et al. VT and static compliance (VT/Pplat − PEEP) were significantly (P < 0.001) but not strongly correlated (r = 0.467). However, inclusion of compliance in regression models for both the primary and secondary outcomes did not alter the primary findings. This indicates that, despite a correlation between VT and compliance, a low VT is an independent risk factor for the development of these adverse outcomes.
Practice Patterns of OLV
In the present study examining single tertiary center practice patterns of OLV for thoracic surgery, we found that relative to current recommendations,33 practitioners continue to use high VT and low levels of PEEP. Moreover, high patient body mass index was associated with the use of larger VT upon normalization to PBW. After normalization to PBW, 73.3, 43.3, 18.8, and 7.2% of patients in the cohort were ventilated with VT greater than 5, 6, 7, and 8 ml/kg, respectively. These findings mirror those of large database studies of surgical patients, which demonstrated that patients of short stature, high body mass index, and female sex received higher VT than the general population when normalized to predicted or ideal body weight.21,34 This occurs presumably because anesthesia ventilators in general do not receive input regarding patient height or weight and because the selection of VT by the anesthesia care provider tends to inadequately compensate for the discrepancy between ABW and PBW. Mean PEEP levels during OLV in our cohort were 4.2 cm H2O and did not significantly differ between patients with and without major postoperative morbidity. During OLV, fewer than half (47%) of the patient cohort received PEEP greater than or equal to 5 cm H2O.
Pathophysiologic Determinants of VILI
It is important to point out that much of the contradictory results in studies of VILI in critically ill and surgical patients may arise from the imprecision inherent in the studied ventilator variables. Neither VT nor PEEP data contain information inherently important from a pathophysiologic standpoint. That is, the pathophysiologic effect of a delivered positive pressure tidal breath, if any, derives from the generated transpulmonary pressure (PL) and its subsequent impact on tissue deformation (stress and strain). Dynamic alveolar strain (VT/functional residual capacity) defines the degree of alveolar tissue deformation between inspiration and expiration. Dynamic but not static alveolar strain is injurious in healthy pigs ventilated at total lung capacity35 and appears to be the more important determinant of VILI.36
Although ΔP has been identified as a risk factor for the development of ARDS in a general surgical population, it has not been previously studied in thoracic surgery.37 In a large retrospective study of ARDS patients, Amato et al.38 identified ΔP but not VT as the ventilation variable that best stratified risk of mortality. Using double-stratification analysis, these authors were able to identify ΔP as the only ventilator variable that predicted mortality in ARDS. Thus, ΔP may serve as a surrogate for dynamic alveolar strain. PEEP is important to the extent that it reduces ΔP and dynamic strain. That PEEP stabilizes the alveolus and can be used to reduce dynamic strain is well established from animal28 and clinical models.35,36 The effect of ΔP on major morbidity after OLV is not yet entirely clear. While ΔP predicts major morbidity after thoracic surgery with OLV in a model including VT, the effect of ΔP in a model without VT approaches, but does not reach, statistical significance.
Limitations
There are several limitations to this study. First, as a consequence of its retrospective design, the relationship between VT, driving pressure, and subsequent complications cannot yet be construed as causal. Second, in an effort to reduce the risk of inaccurate sampling of anesthetic/ventilator data (due to undocumented termination of OLV), we evaluated brief periods (5-min epochs) early during both TLV and OLV. This approach was mandated by the absence of a clearly documented OLV end-time in approximately 40% of these cases. We believe that this approach is likely to result in a sampling of ventilator management data most representative of overall OLV since standard OLV management at our institution is established early during this period and early sampling of OLV data is most likely to avoid the undocumented reinstitution of TLV. Also, this is a single-center study and as such may not be reflective of practice patterns at other institutions.
Finally, in this study as in that of Amato et al.,38 statistical analysis of ΔP and VT is limited by the mechanistic linkage between these variables. Although the correlation between these variables in the current study is low, there is nonetheless a concern about physiologic interdependence. That is, for a given respiratory system compliance, ΔP should be proportional to VT. Thus, one would expect a qualitatively similar relationship between both variables and the studied outcome measures. However, we found that, while VT was inversely related to the primary outcome, ΔP predicted major morbidity after thoracic surgery with OLV in a regression model including VT and approached, but did not reach, statistical significance in a model without VT. One explanation for this discrepancy relates to the variation in ΔP resulting from a given VT based on the “size” of functional lung parenchyma in a manner analogous to that of the “baby lung” concept in ARDS.39 Significant variation in nature and severity of lung diseases and resultant effects on functional lung volume and elastance is expected in a large cohort of thoracic surgical patients. Severely diseased lungs, particularly those with restrictive disease, would be expected to exhibit a larger ΔP for a given VT, while normal lungs would be expected to exhibit a lesser increase in ΔP. A larger VT delivered to normal lung parenchyma (with a resultant lower ΔP relative to that of a diseased lung) could conceivably contribute to a reduced likelihood of atelectasis development, while a low VT in such a patient would be more likely to promote atelectasis and related morbidity. The discrepant ability of VT and ΔP to predict postoperative morbidity in this cohort might then be explained by a differential potential of a given VT to cause overdistension or contribute to atelectasis based on functional lung volume and respiratory system elastance. Additionally, since composite outcome measures used in this study are not likely to reflect a single pathophysiologic mechanism, it is also possible that lower VT and higher ΔP might differentially promote distinct pathophysiologic pathways—atelectasis versus overdistension. Such a possibility is supported by experimental evidence that while high VT with ZEEP produces a similar cytokine response, low VT with ZEEP results in higher mortality.30
Conclusions
Mechanical ventilation has the potential to adversely affect outcomes in surgical patients. These risks clearly extend to patients subjected to OLV. Protective ventilation with low VT and PEEP has been recommended by a number of experts,33,40,41 but no standardized guidelines exist. While it is clear that a high VT regimen without PEEP is injurious, especially in high-risk patients and surgeries, low VT without adequate PEEP may also be injurious.21,30 It is important for the practitioner to note that we are not advocating the use of high or supraphysiologic VT. High-quality prospective data strongly support the use of LPV in TLV and, to a lesser extent, in OLV. While the ideal LPV regimen is yet to be completely elucidated, the current evidence strongly indicates that protection from iatrogenic lung injury follows from the concomitant use of both physiologic (low) VT and sufficient PEEP to prevent overdistension, atelectasis, and tidal recruitment/derecruitment phenomena. Ultimately, advances in our understanding of protective ventilation during OLV are likely to derive from well-designed randomized trials controlling for variables of inherent pathophysiologic significance. The control of PL, ΔP, and alveolar strain in such trials of protective ventilation will be important for the identification of optimal ventilation strategies.
Acknowledgments
Support was provided solely from institutional and/or departmental sources.
Competing Interests
The authors declare no competing interests.