Using pulse oximetry waveforms to detect coarctation of the aorta
BioMedical Engineering OnLine volume 19, Article number: 31 (2020)
Coarctation of the aorta is a common form of critical congenital heart disease that remains challenging to diagnose prior to clinical deterioration. Despite current screening methods, infants with coarctation may present with life-threatening cardiogenic shock requiring urgent hospitalization and intervention. We sought to improve critical congenital heart disease screening by using a novel pulse oximetry waveform analysis, specifically focused on detection of coarctation of the aorta.
Methods and results
Over a 2-year period, we obtained pulse oximetry waveform data on 18 neonates with coarctation of the aorta and 18 age-matched controls hospitalized in the cardiac intensive care unit at Children’s Healthcare of Atlanta. Patients with coarctation were receiving prostaglandin E1 and had a patent ductus arteriosus. By analyzing discrete features in the waveforms, we identified statistically significant differences in the maximum rate of fall between patients with and without coarctation. This was accentuated when comparing the difference between the upper and lower extremities, with the lower extremities having a shallow slope angle when a coarctation was present (p-value 0.001). Postoperatively, there were still differences in the maximum rate of fall between the repaired coarctation patients and controls; however, these differences normalized when compared with the same individual’s upper vs. lower extremities. Coarctation patients compared to themselves (preoperatively and postoperatively), demonstrated waveform differences between upper and lower extremities that were significantly reduced after successful surgery (p-value 0.028). This screening algorithm had an accuracy of detection of 72% with 0.61 sensitivity and 0.94 specificity.
We were able to identify specific features in pulse oximetry waveforms that were able to accurately identify patients with coarctation and further demonstrated that these changes normalized after surgical repair. Pulse oximetry screening for congenital heart disease in neonates may thus be improved by including waveform analysis, aiming to identify coarctation of the aorta prior to critical illness. Further large-scale testing is required to validate this screening model among patients in a newborn nursery setting who are low risk for having coarctation.
Coarctation of the aorta (COA) is one of the most common lesions of congenital heart disease, accounting for approximately 7% of all cases , 4 per 10,000 births  or approximately 1600 newborns per year. Despite this prevalence, it frequently eludes detection in both the prenatal and neonatal periods . Current neonatal screening detects hypoxemia in newborns using pulse oximetry analysis from different extremities in an algorithmic approach (see Fig. 1) in order to identify critical congenital heart disease (CCHD)  before hospital discharge. It is important to evaluate the PPG waveform from extremities supplied by branch vessels both proximal and distal to the patent ductus arteriosus (PDA) and potential coarctation (see Fig. 2). This method of screening does well to detect cyanotic mixing lesions using only numeric values of oxygen saturation. Of the 12 CCHDs, COA is the most common, yet has the worst false-positive rate and only 46% sensitivity with current screening methods [5,6,7]. Even with prenatal ultrasound and pulse oximetry screening the majority (53–62%) of COA cases are late diagnoses [6, 7] and many such infants present to medical care only once they are in life-threatening extremis.
Current newborn screens are inadequate to detect COA, but repurposing existing technology may allow for more timely diagnosis by analyzing waveform data in place of solely using integer values of arterial oxygen saturation. Photoplethysmography (PPG) is routinely seen as a live waveform on pulse oximetry machines that is a composite of multiple waves which have shown promise in detecting changes in peripheral perfusion, blood pressure and local vasomotor tone . Comparing waveform analysis with current non-invasive blood pressure monitoring, pulse oximetry technology may be equal to if not better than current non-invasive techniques . The clinical utility of analyzing PPG waveforms has been demonstrated in newborns with COA in order to quantify physical exam characteristics such as diminished lower extremity pulses and pulse arrival time between extremities . Additional studies have shown promise in PPG signals identifying the presence of significant patent ductus arteriosus in premature infants  as well as PPG may be a useful marker to trend cardiac output and stroke volume . The purpose of this study was to incorporate PPG waveform data to develop an algorithm that can increase the sensitivity of detection of coarctation of the aorta as compared to existing methods. We hypothesized that an algorithmic approach to analyze PPG waveforms would increase the sensitivity of detection of coarctation of the aorta compared to existing methods.
We identified a cohort of 18 patients with COA requiring surgical intervention as neonates and matched them with 18 control patients who met inclusion and exclusion criteria listed in the Methods section. See Table 1 for patient group characteristics. There were no significant differences in gender or ethnicity. Gender did not affect the upper or lower extremity PPG waveforms as the underlying vascular anatomy is the same for monitoring PPG on the digits or distal extremities. Of the control group, 9 patients (50%) had some degree of anomalous pulmonary venous return, 4 (22%) had neonatal arrhythmias, 2 (11%) had valvar abnormalities and 3 others had cardiomyopathy, pulmonary hypertension and a pericardial tumor (representing 6% each).
Preoperative COA vs. unaffected controls
There was a near-linear association between the average rate of rise and fall compared to heart rate (Additional file 1: Figure S1). Presurgical waveform characteristics between COA and control patients were calculated and shown in Table 2. The two characteristics that were significant (p-value < 0.01) were the maximum rate of fall (MRF) in the lower extremity and the difference between the MRF of the upper extremity minus the lower extremity (∆MRF on Table 2 and Fig. 3).
Postoperative COA vs. unaffected controls
Of the 18 patients with COA, one patient did not have upper extremity PPG data recorded postoperatively, therefore they were withheld from the postop statistical analysis. Table 3 shows the same waveform characteristics evaluated postoperatively. Here, we note statistical significance in the MRF in both the upper and lower extremities but not in the ∆MRF, indicating that the difference normalized after surgical correction of the COA (see Fig. 4).
Preoperative COA vs. postoperative COA
Table 4 shows the COA group compared to themselves preoperatively vs postoperatively relative to the COA repair.
In order to test the validity of our algorithm to identify and appropriately characterize patients with COA, we implemented leave-one-out cross-validation (LOOCV) stratified by patient as the number of subjects is too few to split the dataset into a train and test set. The area under receiver operating characteristic curve (AUROC) is generated for the LOOCV performance of the classifier. The MRF in the upper and lower extremity are used to classify subjects in the COA or control groups. A linear discriminant analysis (LDA) classifier achieved an average accuracy of 72%. The receiver operating characteristic curve (ROC) is given in Fig. 5. The area under the curve (AUROC) is 0.78. The suggested operating point is indicated by a red circle with a sensitivity of 0.61 and a specificity of 0.94, providing approximately double the sensitivity of detection of COA compared to current screening.
Our data show that it is indeed possible to correctly identify COA by identifying slope discrepancies of PPG waveforms between the upper and lower extremities of neonates. This was most notable when comparing the maximum rate of fall of the PPG signal in both the upper and lower extremities and obtaining the difference between them. In our data, the MRF remains significant in the lower extremity for the COA group when compared to controls and the ∆MRF does not, emphasizing the importance of needing an upper extremity as a reference for measuring the rate of flow in the lower extremity. This observed difference between COA and controls resolves following surgical correction suggesting that the fixed obstruction of COA and that the COA itself limits the maximum rate of blood flow in the lower extremities. This corresponds to a well-documented finding on echocardiographic assessment of coarctation of the aorta, known as diastolic tailing [12, 13]. When choosing a threshold for this new screening tool, a lower false-positive rate would be preferred, even if the sensitivity is sacrificed. Our data suggest that this detection algorithm could improve the sensitivity of screening from 40 to 60% which would be a significant improvement from the standard of care.
Our study is limited by its retrospective nature and small cohort of patients. Additionally, the signal analysis between upper and lower extremities did not occur simultaneously as was done in the studies of neonates with COA  and premature neonates with PDAs  with a custom-built photodetector for PPG analysis. Simultaneous acquisition reduces error creeping into the data as the clinical situation of a critically ill neonate fluctuates from moment to moment, especially in terms of heart rate and blood pressure. However, we attempted to eliminate these confounders in the data accordingly by normalizing for heart rate and timing changes. Our study may also be limited by the general practice that any neonate with suspected COA is started on a prostaglandin infusion prior to transfer to our quaternary children’s hospital to prevent the patent ductus arteriosus from closing as an attempt to prolong fetal circulation and buy time for further clinical evaluation. For this reason, we excluded infants from our control group that had a significant ductus arteriosus. All of our coarctation group had PDAs at the time of our study, and it is not yet determined what waveform differences there may be from the PDA in our dataset, although studies have shown promise in using PPG to detect the presence and closure of PDA . Our study’s strength lies in the number of patients with usable waveforms obtained that had definitive COA requiring surgical repair in the neonatal period.
Future directions for this research include validating this algorithm with simultaneously acquired pulse oximetry waveforms from upper and lower extremities to eliminate signal noise acquired due to changing patient condition, validating the screening utility in a larger study in a newborn nursery among patients without suspected coarctation as well as exploring additional clinical conditions for this non-invasive screening technology. One potential area of additional clinical utility would be to screen for vascular anomalies (interrupted aortic arch, aberrant subclavian insertion, aortic dissection, aneurysm and perhaps even obstruction due to atherosclerosis or thrombosis) in a new way with pulse oximetry that is readily available at the patient bedside without waiting for advanced imaging.
We identified specific differences in waveform characteristics of PPG signals in patients with COA using beat-by-beat analysis of multiple extremities in the same patients. This provides a valuable clinical application for data analytics in healthcare by comparing waveform morphology to improve the sensitivity of COA diagnosis. We hope to eventually implement this potentially life-saving screening method to identify neonates with COA early enough to prevent cardiogenic shock commonly associated with acute COA. Future directions include validating this screening method on a cohort of prospective, unscreened newborns.
We retrospectively identified a cohort of neonates less than 30 days of age at time of admission to the cardiac intensive care unit who required surgical repair of COA over a 2-year period beginning in 10/2016 through 9/2018 at a quaternary children’s hospital. IRB approval was obtained for a retrospective chart review. Each of these patients was started on a prostaglandin infusion to maintain patency of the ductus arteriosus prior to surgery. These patients were matched 1:1 with controls who were the next closest admission meeting inclusion and exclusion parameters. Controls were excluded for age greater than 30 days of age at time of admission, with concurrent left heart obstructive lesions or a PDA. Single-ventricle patients and those with significant anatomic abnormality of the great vessels were also excluded.
PPG waveforms were digitally recorded from multiple patient sites at different times (upper and lower extremities) for the duration of their intensive care stay. Patients with inadequate recordings (due to poor signal quality or lack of recording) were omitted from analysis. The PPG waveform was recorded at least at 125 Hz to maintain waveform fidelity (waveforms that were recorded at 250 Hz were resampled at 125 Hz to maintain consistency and for ease of use due to file size). Nursing notes documented pulse oximetry probe acquisition site, and we allowed a 30-min grace period to minimize risk of any imprecise documentation of the timing of site change. If there was at least 1 h of data from a single acquisition site, that continuous waveform was saved for further analysis.
These signals were then processed through a bandpass filter (0.5–8 Hz) in order to remove baseline drift and high-frequency noise . The continuous waveform for each site was traversed incrementally and analyzed in Matlab to find a 1-min segment that met criteria for a clean waveform defined as a signal quality index greater than or equal to 90% (as described in previously published work ). Clean segments were normalized by amplitude and by heart rate to eliminate confounders of different pulse oximetry acquisition with different devices and at different times. Each of the waveform characteristics was normalized by heart rate to reduce physiologic differences in patient condition from moment to moment.
We evaluated four main waveform features which are illustrated in Fig. 6. The slope of a line connecting points at 20% and 80% of the max PPG wave amplitude on the upslope produced the average rate of rise (ARR) and on the down slope produced the average rate of fall (ARF). This slope is representative of the average rate of rise and fall of the blood flow rate based on the method demonstrated by Itu et al. in order to estimate the blood flow rate in the descending aorta . Maximum rate of rise (MRR) and maximum rate of fall (MRF) were calculated over short 40-ms windows incrementally along the waveform and the greatest value was reported. This method allowed estimation of the max rate of rise while filtering against noise, compared to evaluation of the max rate from the first derivative of the PPG wave .
Waveform features for each site were computed by taking the median of all the indexed median values for that feature for each 1-min segment. This was repeated for each PPG site change of sufficient duration and signal quality (as defined above). We also calculated the difference between the upper and lower extremities to express the difference between those two site locations for each patient. In this way, we could evaluate the PPG waveform from extremities supplied by branch vessels both proximal and distal to the site of COA (see Fig. 2).
Each of these waveform features was evaluated for each patient’s PPG signal. Those resultant median wave features were analyzed in two main group analyses. First, we compared patients with pre-surgical COA to the unaffected controls, and second we compared pre-surgical COA to themselves as post-surgical COA following surgical repair. Statistical methods to evaluate the differences in waveform characteristics between groups utilized the two-sided Wilcoxon rank sum. We identified the features with the greatest differences that had statistically significant p-values (p < 0.05) as features that classify a PPG waveform as either belonging to a COA or control patient. Then, LOOCV was used to evaluate an LDA classifier trained on the resultant features that are useful in discriminating between the control and COA populations.
Availability of data and materials
The data that support the findings of this study are available from Children’s Healthcare of Atlanta, but restrictions apply to the availability of these data, which were used under license for the current study, and so are not publicly available. Data are however available from the authors upon reasonable request and with permission of Children’s Healthcare of Atlanta.
Coarctation of the aorta
Critical congenital heart disease
Patent ductus arteriosus
Average rate of rise
Average rate of fall
Maximum rate of rise
Maximum rate of fall
Linear discriminant analysis
Rosenthal E. Coarctation of the aorta from fetus to adult: curable condition or life long disease process? Heart. 2005;91(11):1495–502.
Bjornard K, et al. Patterns in the prevalence of congenital heart defects, metropolitan Atlanta, 1978 to 2005. Birth Defects Res A Clin Mol Teratol. 2013;97(2):87–94.
Chang RK, Gurvitz M, Rodriguez S. Missed diagnosis of critical congenital heart disease. Arch Pediatr Adolesc Med. 2008;162(10):969–74.
Diller CL, et al. A modified algorithm for critical congenital heart disease screening using pulse oximetry. Pediatrics. 2018. https://doi.org/10.1542/peds.2017-4065.
Ailes EC, et al. Estimated number of infants detected and missed by critical congenital heart defect screening. Pediatrics. 2015;135(6):1000–8.
Peterson C, et al. Late detection of critical congenital heart disease among US infants: estimation of the potential impact of proposed universal screening using pulse oximetry. JAMA Pediatr. 2014;168(4):361–70.
Lannering K, Bartos M, Mellander M. Late diagnosis of coarctation despite prenatal ultrasound and postnatal pulse oximetry. Pediatrics. 2015;136(2):e406–12.
Bartels K, Thiele RH. Advances in photoplethysmography: beyond arterial oxygen saturation. Can J Anaesth. 2015;62(12):1313–28.
Palmeri L, et al. Photoplethysmographic waveform characteristics of newborns with coarctation of the aorta. J Perinatol. 2017;37(1):77–80.
Goudjil S, et al. Noninvasive technique for the diagnosis of patent ductus arteriosus in premature infants by analyzing pulse wave phases on photoplethysmography signals measured in the right hand and the left foot. PLoS ONE. 2014;9(6):e98763.
Amirtharaj CR, et al. Photoplethysmographic assessment of pulse transit time correlates with echocardiographic measurement of stroke volume in preterm infants with patent ductus arteriosus. J Perinatol. 2018;38(9):1220–6.
Lai WW. Echocardiography in pediatric and congenital heart disease: from fetus to adult. 2nd ed. New York: Wiley; 2016. p. 928.
Thiele RH, Durieux ME. Arterial waveform analysis for the anesthesiologist: past, present, and future concepts. Anesth Analg. 2011;113(4):766–76.
Zhang Q, Zhou D, Zeng X. Highly wearable cuff-less blood pressure and heart rate monitoring with single-arm electrocardiogram and photoplethysmogram signals. Biomed Eng Online. 2017;16(1):23.
Li Q, Clifford GD. Dynamic time warping and machine learning for signal quality assessment of pulsatile signals. Physiol Meas. 2012;33(9):1491–501.
Itu L, et al. Non-invasive hemodynamic assessment of aortic coarctation: validation with in vivo measurements. Ann Biomed Eng. 2013;41(4):669–81.
Shin H, Min SD. Feasibility study for the non-invasive blood pressure estimation based on ppg morphology: normotensive subject study. Biomed Eng Online. 2017;16(1):10.
CDC. Congenital heart defects information for healthcare providers. https://www.cdc.gov/ncbddd/heartdefects/hcp.html. Accessed 2 Nov 2018, 17 Oct 2019.
We would like to thank the Information Systems Department at Children’s Healthcare of Atlanta and the Bioinformatics Department at Emory University for their assistance in converting large amounts of data into useable formats pertinent to this project for our analysis.
Supported by the Orefice Family Foundation.
Ethics approval and consent to participate
The Institutional Review Board of Children’s Healthcare of Atlanta authorized the data review for this study (IRB number 17-196). The need for informed consent was waived due to the retrospective nature of the study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Normalizing Slope by Heart Rate. Using a control patient, we observed a near linear association between the average rate of rise and fall compared to heart rate. The faster the heart rate the greater the slope of the rise and fall of the waveforms, in order to accommodate the rapid cardiac cycle. (A) Shows the strong negative correlation between heart rate and rate of fall. (B) Shows the strong positive correlation between heart rate and rate of rise.
About this article
Cite this article
Sorensen, M.W., Sadiq, I., Clifford, G.D. et al. Using pulse oximetry waveforms to detect coarctation of the aorta. BioMed Eng OnLine 19, 31 (2020). https://doi.org/10.1186/s12938-020-00775-2