Abstract
Rationale Pulmonary hypertension (PH) is associated with significant perioperative morbidity and mortality. We hypothesised that pulmonary arterial hypertension (PAH) composite risk assessment scores could estimate perioperative risk for PH patients when adjusted for inherent procedural risk.
Methods We identified patients in the Johns Hopkins PH Center Registry that had noncardiac surgery (including endoscopies) between September 2015 and January 2020. We collected information on preoperative patient-level and procedural variables and used logistic regression to evaluate associations with a composite outcome of death within 30 days or serious postoperative complication. We generated composite patient-level risk assessment scores for each subject and used logistic regression to estimate the association with adverse surgical outcomes. We adjusted multivariable models for inherent procedural risk of major cardiovascular events and used these models to generate a numerical PH perioperative risk (PHPR) score.
Results Among 150 subjects, 19 (12.7%) reached the primary outcome, including 7 deaths (4.7%). Individual patient-level and procedural variables were associated with the primary outcome (all p<0.05). A composite patient-level risk assessment score built on three noninvasive parameters was strongly associated with reduced risk for poor outcomes (OR=0.4, p=0.03). This association was strengthened after adjusting the model for procedural risk. A PHPR score derived from the multivariable model stratified patients into low (0%), intermediate (≤10%), or high (>10%) risk of reaching the primary outcome.
Conclusion Composite PAH risk assessment scores can predict perioperative risk for PH patients after accounting for inherent procedural risk. Validation of the PHPR score in a multicentre, prospective cohort is warranted.
Abstract
For patients with pulmonary hypertension undergoing noncardiac surgery, perioperative risk can be estimated using a model that combines inherent procedural risk with composite PAH risk assessment scores https://bit.ly/3uTCibP
Introduction
Perioperative risk assessment presents a formidable challenge for providers managing patients with pulmonary hypertension (PH). Multiple studies have demonstrated increased risk of mortality and serious morbidity among PH patients undergoing noncardiac surgery [1–6], and some have identified PH as an independent risk for perioperative complications and death [6–8]. This clinical challenge becomes increasingly relevant for patients with pulmonary arterial hypertension (PAH) as advances in medical therapy expand the population of patients requiring invasive procedures [9]. Current PH guidelines [10, 11] provide limited guidance for assessing perioperative risk. Better risk assessment tools may assist providers and patients in complex perioperative decision-making.
Previous studies have identified preoperative variables associated with poor surgical outcomes for patients with PH [1, 3, 5, 8]. The clinical utility of these associations has been limited, as they have been inconsistent among studies with differing patient populations, surgical procedures, and small samples sizes. In parallel, there is a growing effort to address longitudinal clinical risk assessment for PAH patients, unrelated to surgery. Several approaches have been used to combine multiple clinical variables into composite risk assessment scores [11–16]. These scoring systems have been validated in predicting clinical [12–14, 17, 18] and patient-centred [19] outcomes in diverse PAH groups. We hypothesised that composite risk assessment scores could estimate perioperative risk for PH patients when adjusted for procedural risk. Moreover, we present “proof-of-concept” data showing that this approach can be used to develop a clinically pragmatic perioperative risk assessment tool.
Methods
Study population and data extraction
We included participants in the Johns Hopkins Pulmonary Hypertension Center Registry (Institutional Review Board no. NA_00027124), which enrolls consecutive patients from all World Symposium on Pulmonary Hypertension (WSPH) groups [20]. From September 2015 to January 2020, we identified upcoming surgical procedures among Registry participants during weekly meetings, including all noncardiac surgical procedures (excluding lung transplantation) requiring at least moderate sedation. Some patients had multiple procedures during the study period; only the last procedure was analysed.
Preoperative patient-level risk variables of interest
Individual patient-level risk variables were selected a priori from previous studies evaluating operative risk in PH [1–4, 7, 8] and studies evaluating predictors of overall survival in PAH [11, 13–16, 21, 22]. These included American Society of Anesthesiologists’ (ASA) classification, World Health Organization functional class (WHO FC), serum brain natriuretic peptide (BNP) or N-terminal pro-hormone BNP (NT-proBNP) level, 6-min walking distance (6MWD), echocardiographic parameters (right atrial and ventricular dilation, tricuspid annular plane systolic excursion (TAPSE), maximum Doppler tricuspid regurgitant velocity), and right heart catheterisation parameters (right atrial pressure (RAP), cardiac index, mixed venous oxygen saturation (SvO2)). We used multiple approaches to generate composite risk assessment scores [12–14] from patient-level risk variables as defined in current PAH guidelines [11]. In one approach (“low-risk focused”), we quantified the number of “low-risk” features present among three noninvasive parameters (WHO FC I or II, 6MWD >440 m, BNP <50 ng·L−1 or NT-proBNP <300 ng·L−1) or a combination of four invasive and noninvasive parameters (WHO FC I or II, 6MWD >440 m, RAP <8 mmHg, cardiac index ≥2.5 L·min−1·m−2) to generate composite scores [14]. In another approach (“score and average”), we assigned a value of 1 (low risk), 2 (intermediate risk), or 3 (high risk) for each of six clinical parameters (supplemental table S1), then averaged the values to generate a composite score [13]. In all cases, we used the closest measurement before the surgical procedure for analyses or composite risk assessment score calculation.
Procedural risk variables of interest
We identified procedural risk variables of interest using American College of Cardiology/American Heart Association (ACC/AHA) perioperative guidelines [23] and previous studies [3, 4, 7, 8, 24]. These included inherent risk of major cardiovascular events (“elevated” or ≥1% risk versus “low” or <1% risk), emergency procedure (patient requires access to the operating room within 24 h of the decision to operate [25]), length of procedure (>3 h), and location of patient immediately prior to procedure (inpatient versus outpatient).
Outcome of interest
The primary outcome was the composite of death within 30 days or serious postoperative complication. Serious postoperative complications required intensive care unit (ICU) admission and included any of the following: 1) haemodynamic instability requiring inotropes and/or vasopressors; 2) respiratory failure requiring invasive mechanical ventilation ≥48 h; 3) initiation of inhaled pulmonary vasodilators; 4) acute coronary syndrome; 5) cerebrovascular accident; 6) arrhythmia; 7) renal failure (serum creatinine >2 mg·dL−1 for patients with normal preoperative levels, or a 50% increase in serum creatinine, or initiation of dialysis); 8) hepatic injury (elevation of aspartate aminotransferase and/or alanine aminotransferase by 50% or more); and 9) sepsis [1–3, 8, 26]. Secondary outcomes included postoperative hospital length of stay (LOS), ICU LOS, and hospital readmission for a medical deterioration within 30 days of the procedure.
Statistical analysis
For descriptive analyses, we used the independent t-test for normally distributed continuous variables, the Wilcoxon rank-sum test for skewed continuous variables, and the Chi-squared test for categorical variables. We used univariable logistic regression to evaluate associations between preoperative variables and the primary outcome. For secondary outcomes, we used linear regression to estimate associations between preoperative variables and hospital or ICU LOS and logistic regression to estimate the association with readmission within 30 days.
We used logistic regression to estimate associations between composite risk assessment scores (“low-risk focused” and “score and average”) and the primary outcome. We used multivariable regression to account for differences in procedural risk based on inherent risk of major cardiovascular events (elevated versus low) [23]. We selected this variable to capture procedural risk because it had the strongest association with the primary outcome, potentially incorporates the other procedural risk variables, and is known to the perioperative physician prior to surgery. We generated receiver-operating characteristic (ROC) curves to quantify the discriminatory power of multivariable models, and used the DeLong test to compare the area under ROC curves. In all analyses, we considered a p-value (two-tailed) <0.05 to represent statistical significance. We used Stata version 15 (StataCorp LLC, College Station, TX, USA) for all analyses.
Risk assessment tool development
We used the multivariable model including the three-parameter noninvasive “low-risk focused” score to develop a risk assessment tool (PH perioperative risk or PHPR score) that can be calculated before surgery. We used this model to develop the PHPR score because of its superior performance characteristics and because additional invasive procedures are not required. The PHPR score (range 0–7) is calculated by adding a patient-level component and a procedural component. The patient-level component is the number of features in the composite score that do not reach “low-risk” thresholds (0–3). The procedural component is binary. “Low” risk procedures are assigned a value of 0 (no added risk). “Elevated” risk procedures are assigned a value of 4, derived from the ratio of odds ratios for “elevated” procedural risk and number of composite score features that do not reach “low-risk” thresholds in the multivariable regression model (e2.7/e1.3=e1.4=4). For example, a patient with WHO FC II symptoms, NT-proBNP=600 ng·L−1, and 6MWD=250 m undergoing an endoscopy (“low” risk) would have a PHPR score of (2+0)=2, while the same patient undergoing a total knee replacement (“elevated” risk) would have a PHPR score of (2+4)=6. We used the distribution of PHPR scores to assign three risk categories (0–2=low, 3=intermediate, 4–7=high). Thresholds were intentionally set to include all patients having an “elevated” procedural risk procedure in the highest category.
Results
Study population, surgical procedures, and preoperative evaluation
During the study period, we prospectively identified 150 subjects undergoing noncardiac surgery (table 1). The cohort was predominantly female (76%) with an average age of 59.7 years. Approximately half (52%) had PAH, and 47.4% of these had connective tissue disease-associated PAH (supplemental table S2). Most subjects were receiving PAH-specific therapy at the time of surgery (72%), and nearly all were deemed ASA classification of 3 or higher.
Preoperative variables
Most procedures (82%) were performed at the Johns Hopkins Hospital, and 73.3% were classified as “low” risk using ACC/AHA guidelines [23]. Of 110 procedures classified as “low” risk, 54 (49%) were endoscopies, primarily of the gastrointestinal tract (supplemental table S3). Endoscopic procedures are routinely staffed by cardiac anaesthesiology providers for PH patients at our institution [27]. Forty procedures (26.7%) classified as “elevated” risk included primarily intraperitoneal, orthopaedic, or intrathoracic surgeries. 6% of all procedures were emergency procedures, 32% were performed during inpatient hospitalisation, and 28% lasted longer than 3 h.
Outpatient clinical assessment by a PH specialist (67.3%) or anaesthesiologist (21.3%) was completed within 90 days of surgery for most patients. Among procedures performed at the Johns Hopkins Hospital, 75.6% were referred to the cardiac anaesthesiology liaison for case review. Most patients had WHO FC assessment (63%), natriuretic peptide (NT-proBNP or BNP) measurement (55%), and echocardiography (53%) within 3 months before surgery, while fewer had 6MWD testing (23%) or right heart catheterisation (21%) in this interval (fig 1). Average values for the closest measurement before surgery are shown in supplemental table S4.
Percentage of patients having preoperative assessments during specified time intervals before surgery. FC: World Health Organization functional class; NP: natriuretic peptide; Echo: echocardiography; 6MWT: 6-min walk test; RHC: right heart catheterisation.
Surgical outcomes
The primary outcome was observed in 19 patients (12.7%; table 2), including 7 deaths (4.7%). Median time from surgery until death was 13 days (interquartile range (IQR) 5–23 days), with causes listed in supplemental table S5. Seventeen patients (11.3%) developed serious postoperative complications. Among these, haemodynamic instability and respiratory failure were most common. Median hospital LOS was 1 day (IQR 0–6 days). Of the 37 patients (24.7%) admitted to the ICU after surgery, 20 were planned (enhanced monitoring) and 17 were unplanned. ICU LOS was longer for unplanned admissions (median 4 versus 1 day; p<0.001, supplemental figure S1). Nineteen patients (12.7%) were readmitted to the hospital within 30 days of surgery for clinical deterioration (supplemental table S6).
Outcomes
Associations between individual risk variables and surgical outcomes
Associations between individual preoperative variables and the primary outcome are shown in table 3. Among patient-level variables, we observed no associations between patient age, sex, presence of systemic sclerosis, WSPH group, or use of PAH-specific therapy and the primary outcome. Patients with the highest ASA classification (4 or 5) had increased risk compared to those with ASA classification 2 or 3. We did not find associations between preoperative evaluation by a PH specialist, anaesthesiologist, or case review by the cardiac anaesthesiology liaison and the primary outcome. WHO FC was associated with the primary outcome; subjects with WHO FC 3 or 4 limitations had greater risk than those with WHO FC 1 or 2 limitations (OR 5.87; p=0.02). Natriuretic peptide levels were also associated with the primary outcome; subjects with NT-proBNP >1400 ng·L−1 (or BNP >300 ng·L−1) had higher odds of adverse events compared with those with NT-proBNP <300 ng·L−1 (or BNP <50 ng·L−1) (OR 7.08; p=0.01). We did not find associations between right atrial dilation, TAPSE, or pericardial effusion and the primary outcome. However, severe right ventricular (RV) dilation (OR 6.27; p=0.01) and maximum Doppler tricuspid regurgitant velocity (OR 1.5 for every increase of 0.5 m·s−1; p=0.02) were associated with increased risk. We found no associations between 6MWD or invasive haemodynamic measurements and the primary outcome.
Associations between preoperative variables and the primary outcome
All procedural variables examined were associated with the primary outcome. The strongest association was observed with procedures designated as “elevated” risk for major cardiovascular events, with an OR of 8.35 (p<0.001) compared to “low” risk procedures. Emergency procedures, those lasting >3 h, and procedures performed on inpatients were also associated with increased odds of the primary outcome (table 3).
Associations between composite PAH risk assessment scores and surgical outcomes
We used three approaches to generate composite PAH risk assessment scores, then used univariable logistic regression to estimate associations with the primary outcome. Only the three-parameter noninvasive “low-risk focused” score was significantly associated with the primary outcome (OR 0.4, p=0.03) (table 4). When the models were adjusted for procedural risk, all three composite PAH risk assessment scores were significantly associated with the primary outcome. In multivariable models, the strongest association was with the three-parameter noninvasive “low-risk focused” score where each additional “low-risk” feature was associated with a 70% reduction in the estimated odds of an adverse surgical outcome (p=0.01). We generated ROC curves to assess the discriminatory power of the multivariable models. Each had robust discriminatory power for the primary outcome, with a c-statistic at or above 0.8 (supplemental figure S2). However, when these multivariable ROC curves were compared to ROC curves for their component patient-level and procedural risk scores, only the three-parameter noninvasive “low-risk focused”-based multivariable model had better discrimination than both univariable models (figure 2 and supplemental figure S3). Finally, we used regression analyses to evaluate associations between the composite PAH risk assessment scores and secondary outcomes (table 5). In univariable analyses, all approaches were associated with 30-day hospital readmission and none were associated with ICU LOS. The three-parameter noninvasive “low-risk focused” and “score and average” methods were also associated with hospital LOS in univariable models. After adjustment for procedural risk, all three composite PAH risk assessment scores were significantly associated with hospital LOS and odds of hospital readmission, while no score was associated with ICU LOS.
Association between composite patient-level risk scores and the primary outcome
Receiver-operating characteristic (ROC) curves for logistic regression models including only patient-level (Pt-level) risk, only procedural (Proc) risk, or both (Pt-level+Proc). The area under the curve (AUC) for each model is shown in the legend. p-values are for DeLong test comparing ROC curves for each univariable model with the multivariable model.
Association between composite patient-level risk scores and secondary outcomes
Development of a perioperative PH risk assessment tool
We used the three-parameter noninvasive “low-risk focused” multivariable model to develop a risk assessment tool (PHPR), as outlined in the methods. We used the distribution of patients to stratify PHPR scores into discrete categories that assign low (0%), intermediate (≤10%), or high (>10%) perioperative risk (figure 3a). As anticipated, the percentage of subjects reaching the primary outcome increased steadily across PHPR categories (p<0.001). ROC curve analysis showed that the discriminatory power of the PHPR categories approximated that of the original multivariable model (figure 3b). As an exploratory analysis, we used ROC curves to compare the performance of the PHPR categories among patients with PAH (WSPH group I) and those without PAH (WSPH group II–V). We did not observe a difference in the discriminatory power of the PHPR score categories between the groups (p=0.96; figure 3c).
a) Percentage of patients reaching the primary end-point in each pulmonary hypertension perioperative risk (PHPR) category is shown. The number of patients in each category is shown in parenthesis. p-value is for Chi-squared test. b) Receiver-operating characteristic (ROC) curve for the PHPR categories. c) PHPR category ROC curves for patients with pulmonary arterial hypertension (PAH) and other forms of pulmonary hypertension are compared. The area under the curve (AUC) for each group is shown in the legend. p-value is for DeLong test comparing ROC curves.
Discussion
We describe a novel approach to perioperative risk assessment for PH patients undergoing noncardiac surgery. We showed that a multivariable model combining patient-level risk (composite PAH risk assessment scores) and procedural risk strongly predicts adverse surgical outcomes with improved discriminatory power compared to either approach alone. Importantly, all data required to evaluate risk is routinely collected in clinical practice, does not require additional invasive assessment, and is known to the provider before the procedure. We also provide “proof of principle” that the model can be converted to a clinically pragmatic risk score that may be useful for patients with all forms of PH.
Several studies have described associations between individual patient-level variables and perioperative risk, but findings have been heterogeneous across studies and difficult to operationalise clinically. In a retrospective study of 145 PH patients (excluding those with left heart disease) undergoing noncardiac procedures under general anaesthesia, Ramakrishna et al. [1] showed that advanced WHO FC, RV hypertrophy, elevated RV systolic pressure relative to systemic systolic blood pressure and a history of pulmonary embolism were associated with perioperative morbidity and mortality. In another retrospective study of 196 PH patients (also excluding those with left heart disease) undergoing noncardiac procedures, Deljou et al. [5] showed that WHO FC and natriuretic peptide level are associated with postoperative complications. A prospective multicentre study conducted on 114 patients with PAH undergoing noncardiac nonobstetric surgeries identified elevated RAP and low 6MWD (<399 m) as patient-level factors associated with adverse perioperative outcome [3]. In a retrospective analysis of 173 patients undergoing right heart catheterisation and noncardiac surgery (96 with PH), mean pulmonary artery pressure, ASA classification, and chronic renal insufficiency were identified as independent risk factors for postoperative morbidity [8]. We confirmed some of these relationships in our cohort (ASA classification, WHO FC, natriuretic peptide), but did not observe statistically significant associations for others. We also identified novel associations with two echocardiographic parameters (severe RV dilation, maximum Doppler tricuspid regurgitant velocity). This variability likely reflects the small sample sizes and limited statistical power for all of these studies, an inherent limitation to analyses of subjects with a rare disease who are often advised to avoid surgery [10, 11].
To address these limitations, we hypothesised that a composite “score” of variables routinely obtained in clinical practice would capture patient-level risk. Initially developed to assess prognosis in PAH patients [12–16], this approach combines clinical, functional, exercise, RV function and haemodynamic parameters to stratify patients into low, intermediate and high risk for 1-year mortality. Several scoring methods have been validated in estimating 1-year survival at the time of diagnosis [13], after treatment for 1 year [12, 14], and in several PAH subgroups [17, 18]. More recently, these composite scores have been associated with health-related quality of life and rate of hospitalisations in PAH patients [19]. Our findings now extend this approach to perioperative risk assessment. In this domain, we found the strongest associations with the three-component noninvasive “low-risk focused” composite risk assessment score. However, we may have lacked sufficient power for statistical significance in univariable models with the other approaches, as there appeared to be meaningful trends in associations with the primary outcome and all approaches showed concordant associations with secondary outcomes. It is also possible that the three-component noninvasive model had superior performance in our analysis because the other approaches include invasive haemodynamic parameters, which were often obtained more than 1 year prior to surgery in our cohort (not reflective of perioperative clinical status). We believe this reflects clinical practice, as there are currently no strong recommendations for preoperative haemodynamic testing in patients with PAH.
Procedural factors intuitively inform perioperative risk, and procedural variables previously associated with adverse surgical outcomes among PH patients include inherent risk of major adverse cardiovascular event [1, 2, 5], emergency surgery [2–4, 7], and procedures lasting longer than 3 h [2, 24]. We confirmed these relationships in our cohort, with the inherent procedural risk (elevated versus low) having the strongest association with the primary outcome. As this variable potentially incorporates other procedural risk factors (e.g. procedure length) and is known to the perioperative physician before the procedure, we elected to use it to adjust for procedural risk in multivariable models. This adjustment improved the discriminatory power of all three composite PAH risk assessment scores, but was significantly better than procedural risk alone at predicting adverse outcomes only when combined with the three-component noninvasive “low-risk focused” score. Although we lacked sufficient statistical power to determine whether procedural risk is an effect modifier or a confounder of the patient-level risk, we note that four of the seven deaths in our cohort occurred among patients with low procedural risk and high patient-level risk.
Clinically pragmatic risk assessment tools should aid providers in perioperative decision-making. We used the three-parameter noninvasive “low-risk focused” multivariable model to develop PHPR scores and stratify patients into risk categories (low, intermediate, high). This approach requires only clinical data routinely collected at regular intervals [11] and knowledge of the type of upcoming surgery. It avoids the need for repeat invasive haemodynamic assessment, which could delay surgery or generate additional risk. Furthermore, previous analyses have suggested that the noninvasive “low-risk focused” approach may be more accurate than the “score and average” approach for identifying patients with the lowest risk (highest long-term survival) [28]. As devised, we believe the PHPR score is most useful in identifying patients at the lowest perioperative risk that may be able to avoid additional workup or substantial adjustments to the perioperative plan. Because category thresholds intentionally restrict “elevated” risk procedures to the “high” PHPR category, the PHPR score incorporates patient-level risk to stratify patients undergoing “low” risk procedures as low or intermediate risk. Notably, a low PHPR score differentiated the 0% of patients who reached the primary outcome from the 10% of patients with an intermediate PHPR score who reached the primary outcome. Future studies in additional cohorts will be useful for further calibration of the PHPR score, potentially allowing additional discrimination between intermediate and high-risk patients.
Limitations
The cohort is limited to a single PH referral centre and included mostly low-risk procedures, which may limit generalisability. Other factors that may limit generalisability include the enrichment of our cohort with connective tissue disease-associated PH patients and routine incorporation of cardiac anaesthesiology in perioperative management. However, our observed rates of perioperative mortality and morbidity (4.7% and 11.3%) are consistent with previous reports (mortality 1–9.7% [1–8]; morbidity 6.1–42% [1, 3, 5]). The small sample size limited some analyses (as outlined above), but is also consistent with prior studies [1–3, 5, 7, 8]. There was heterogeneity among subjects in the time between some preoperative patient-level variables and surgery, which could lead to ascertainment bias. This heterogeneity reflects the observational nature of our study and highlights a key knowledge gap in the field, namely the lack of evidence-based recommendations for specific preoperative testing protocols. Future prospective trials could evaluate the predictive power of a defined preoperative testing regimen (e.g. WHO FC, 6MWD, and NT-proBNP within 3 months of surgery). We were unable to assess the relationship between the US Registry to Evaluate Early and Long-term PAH Disease Management (REVEAL) 2.0 risk score, another composite PAH risk assessment tool [15, 16], and perioperative outcomes given the limited number of patients having sufficient information to calculate REVEAL 2.0 scores. Finally, we recognise that the PHPR score will require additional prospective validation, but present it here as “proof-of-principle” that a clinically pragmatic scoring system can be used to inform perioperative risk.
Conclusion
For PH patients undergoing noncardiac surgery, a composite PAH risk assessment score can be combined with procedural risk to improve perioperative risk assessment. When converted to a PHPR score, this approach provides a simple and clinically pragmatic tool for perioperative risk assessment for patients with PH.
Supplementary material
Supplementary Material
Please note: supplementary material is not edited by the Editorial Office, and is uploaded as it has been supplied by the author.
Supplementary material 00257-2021.SUPPLEMENT
Figure S1 00257-2021.FIGURES1
Figure S2 00257-2021.FIGURES2
Figure S3 00257-2021.FIGURES3
Acknowledgements
The authors would like to acknowledge David N. Hager (Division of Pulmonary and Critical Care Medicine, Johns Hopkins University, Baltimore, MD, USA) for his insightful comments during the design of this study.
Footnotes
This article has supplementary material available from openres.ersjournals.com
Conflict of interest: H.J. Hassan has nothing to disclose.
Conflict of interest: T. Housten has nothing to disclose.
Conflict of interest: A. Balasubramanian has nothing to disclose.
Conflict of interest: C.E. Simpson has nothing to disclose.
Conflict of interest: R.L. Damico has nothing to disclose.
Conflict of interest: S.C. Mathai reports personal fees from Actelion, Liquidia and United Therapeutics outside the submitted work, and is a member of the Rare Disease Advisory Panel at the Patient Centered Outcomes Research Institute (PCORI).
Conflict of interest: P.M. Hassoun serves on a scientific advisory board for MSD outside the submitted work.
Conflict of interest: J. Steppan has nothing to disclose.
Conflict of interest: P.J. Leary reports personal fees from Bayer, research support as a site PI for United Therapeutics and Actelion, and research grants from the American Heart Association, Pulmonary Hypertension Association, Chest Foundation, and NIH/NHLBI. All are outside the submitted work.
Conflict of interest: T.M. Kolb has nothing to disclose.
- Received April 12, 2021.
- Accepted May 2, 2021.
- Copyright ©The authors 2021
This version is distributed under the terms of the Creative Commons Attribution Non-Commercial Licence 4.0. For commercial reproduction rights and permissions contact permissions{at}ersnet.org