Patient identification
We retrospectively reviewed the initial neonatal admission of all patients with PAIVS admitted in the cardiac ICU to one of the 19 participating CoRe-PCICS pediatric centers in the United States between 2009 and 2019. A list of participating institutions can be found in Supplemental Table 1. The study was approved by the institutional review boards of all participating study centers and the data coordinating center at Cincinnati Children’s Hospital Medical Center (ID: 2019 − 0708; 24 May 2019) and given the retrospective nature of the data collected, the need for informed consent was waived (IRB Protocol H-45938).
We used demographic and anatomic data (obtained from echocardiogram, cardiac catheterization, and computed tomography angiography) during the first week of life to determine association with ECMO and mortality. Inclusion criteria were patients with ECMO use at any point during the initial neonatal hospitalization, including before first intervention, intraoperatively, postoperatively and/or prior to hospital discharge. Hospital mortality was defined as death prior to hospital discharge. Exclusion criteria were presence of ventricular septal defect, evidence of blood flow across the pulmonary valve, and Ebstein’s anomaly or severely dysplastic tricuspid valve. Exclusion criteria were determined based on the initial echocardiogram report.
Variables of interest
Study relevant descriptive information was collected. Also, as one of the aims was to model the use of ECMO and associated inpatient mortality using demographic, echocardiographic, and catheterization data routinely obtained in the first week of life, variables of interest were selected accordingly. The following demographic information was obtained: presence of prenatal diagnosis, gestational age (weeks), birth weight (kg), birth height (cm), gender, ethnicity, presence of chromosomal abnormality, presence of a noncardiac anomaly.
Based on previously reported data, the following echocardiographic data were obtained from the initial echocardiogram: tricuspid valve z-score, right ventricular morphology, right ventricular size, estimated right ventricular pressure, type of pulmonary atresia, pulmonary valve annulus z-score, right pulmonary artery z-score, left pulmonary artery z-score, left ventricular ejection fraction, presence of normal right coronary artery, presence of normal left coronary artery, evidence of to-from flow in the coronary circulation, and presence of ventriculocoronary connections. All participating institutions used echocardiography-derived calculations for TV Z-score determination with the majority (17 out of 19 centers) utilizing the Boston Children’s Hospital dataset. No data imputation was done.
The following cardiac catheterization data was obtained from the initial catheterization: presence of ventriculocoronary connections, presence of coronary artery stenosis, and presence of coronary ostial atresia.
Statistical analyses
Descriptive variables were described as absolute frequency and percentage. Continuous variables were described as mean and standard deviation due to data for most variables being normally distributed.
Demographic characteristics, echocardiographic data, and cardiac catheterization were compared between those who received extracorporeal membrane oxygenation and those who did not, using chi-square or Fisher exact tests for descriptive variables and independent t-tests or Mann-Whitney-U tests for continuous variables. Variables with a p-value of less than 0.2 were entered into a backward logistic regression with need for extracorporeal membrane oxygenation as the dependent variable. A model was developed using the beta-coefficients for independent variables that were found to have significant association with the dependent variable. The developed model was used to calculate the estimated risk of extracorporeal membrane oxygenation. A receiver operator curve analysis was then conducted to determine the accuracy of this risk score to predict extracorporeal membrane oxygenation use.
Among ECMO patients, analyses compared those that did and did not experience inpatient mortality during the admission. Demographic characteristics, echocardiographic data, and cardiac catheterization data were compared using chi-square of Fisher exact analyses for descriptive variables and independent t-tests of Mann-Whitney-U tests for continuous variables. Variables with a p-value of less than 0.2 were entered into a backward logistic regression with inpatient mortality as the dependent variable. A model was developed using the beta-coefficients for independent variables that were found to have significant association with the dependent variable. The developed model was used to calculate the estimated risk of inpatient mortality on extracorporeal membrane oxygenation. A receiver operator curve analysis was then conducted to determine the accuracy of this risk score to predict inpatient mortality in patients requiring extracorporeal membrane oxygenation.
All statistical analyses were conducted using SPSS Version 23.0. A p-value of less than 0.05 was considered to be statistically significant. Any use of the word “significant” or “significantly” in the manuscript refers to statistical significance unless explicitly specified.