Introduction
Children born with critical congenital heart disease (CHD) are at increased risk for neurodevelopmental impairments that influence educational, occupational, and quality of life outcomes. Reference Marino, Lipkin and Newburger1,Reference Sood, Newburger and Anixt2 With decreasing CHD mortality related to advances in care, attention has shifted towards identifying modifiable risk factors that may optimise neurodevelopmental outcomes for this growing population. Reference Gilboa, Devine and Kucik3 While numerous clinical risk factors for neurodevelopmental impairment have been identified, most are medically necessary components of care and not easily modified. Reference Wernovsky and Licht4,Reference Marelli, Miller, Marino, Jefferson and Newburger5 Family-centred care, recognised as a core value and approach to care in paediatrics, 6–8 acknowledges the importance of providing care and support not only to the paediatric patient, but also to the child’s parents and caregivers. Furthermore, family-centred care includes both health care and social care of paediatric patients and their families. 9–Reference Ho, Shaul, Chapman and Ford-Jones12
Within this context of family-centred care, maternal mental health Reference Lisanti13–Reference Woolf-King, Anger, Arnold, Weiss and Teitel16 and social determinants of health Reference Lopez, Baker-Smith and Flores17 are potential modifiable factors influencing child neurodevelopment in CHD that have yet to be studied in combination using longitudinal datasets. Undergoing a foetal echocardiogram and receiving a prenatal diagnosis of CHD is a significant stressor that impacts maternal mental health, Reference Bravo-Valenzuela, Peixoto and Araujo Junior18 resulting in the development of enduring maternal mental health symptoms. Reference Lisanti13 While increased stress is expected in this context, prenatal stress can contribute to the development of maternal mental health pathologies such as major depressive disorder and post-traumatic stress disorder. Reference Woolf-King, Anger, Arnold, Weiss and Teitel16,Reference Wray, Cassedy and Ernst19,Reference Brosig, Mussatto, Kuhn and Tweddell20 Maternal mental health risks of depressive and stress symptoms have also been shown to be associated with increased risk for worse pregnancy outcomes such as preterm birth and low birth weight infants. Reference Brown, Yelland, Sutherland, Baghurst and Robinson21,Reference Dunkel Schetter22 Mothers experiencing additional social determinants of health adversities may be at greater risk and these risks may ultimately negatively influence child hospitalisation and neurodevelopmental outcomes in children with CHD. Reference Lisanti13
Social determinants of health contribute to population health outcomes, including mortality and quality of life. Reference Hood, Gennuso, Swain and Catlin23 Social determinants of health disparities have been shown to adversely influence neurodevelopmental outcomes and to be associated with higher mortality and worse quality of life. Reference Magnan24–Reference Wang, Liu, Druschel and Kirby35 These disparities have been demonstrated by social determinants of health factors, including maternal age and education, Reference Tran, Forman, Mossialos, Nasir and Kulkarni27 insurance type, Reference Kucik, Cassell and Alverson34 community-level socioeconomic risk, Reference Kucik, Nembhard and Donohue33 and race and ethnicity. Reference Tran, Forman, Mossialos, Nasir and Kulkarni27,Reference Wang, Liu, Druschel and Kirby35 However, the extent to which these disparities, combined with prenatal maternal mental health symptoms, affect children with CHD is unknown. Using clustering methods of social determinants of health and mental health variables may provide new insights into the differing levels vulnerability of these paediatric patients and their mothers. Reference Zhang, Clay, Lee and Howell36
We utilised the Parental Stress and Resilience in CHD Model to study the association between maternal mental health symptoms and social determinants of health risks on outcomes in children with CHD. The Parental Stress and Resilience in CHD Model (Figure 1) was specifically developed to provide a framework of parental stress in CHD through a health equity lens and identifies parent demographic, mental health, and social determinants of health variables influencing parental and child outcomes, including parent sex, race and ethnicity, socioeconomic status, and neighbourhood environment. Reference Lisanti13 These relationships have been posited for the general population, as a recent review of literature highlighted the critical need for interventions targeting social determinants of health and psychosocial conditions in early pregnancy. Reference Simoncic, Deguen, Enaux, Vandentorren and Kihal-Talantikite37 This study explored the combined influence of prenatal maternal mental health symptoms and maternal social determinants of health risks on pregnancy outcomes and hospitalisation and neurodevelopmental outcomes in children with CHD.

Figure 1. The parental stress and resilience in CHD model.
Materials and methods
This retrospective cohort study examined 189 mother–child dyads of children born with CHD who underwent neonatal cardiac surgery, linking child electronic health record data with their mothers’ prenatal demographics and maternal mental health screening. Inclusion criteria for our sample were all children with CHD who had neonatal cardiac surgery, who completed neurodevelopmental assessment at 13 to 29 months of age between 1 December, 2012, and 30 November, 2019, and whose mothers completed prenatal mental health screening for depression and post-traumatic stress at one of their early prenatal visits to our centre for foetal diagnosis and treatment. This window of time reflects the time period after we began a prenatal mental health screening programme at our institution as a part of standard clinical practice Reference Cole, Olkkola, Zarrin, Berger and Moldenhauer38,Reference Cole, Moldenhauer and Berger39 and before the start of the COVID-19 pandemic, which may confound variables in our study, including maternal mental health and child neurodevelopment. Children who underwent cardiac surgery for CHD >30 days of age were excluded, as well as those children with mothers who did not have prenatal maternal mental health screening.
Study data were collected and managed using Arcus resources hosted at Children’s Hospital in Philadelphia. Arcus is a suite of tools and services developed to enhance research efforts at Children’s Hospital in Philadelphia by helping researchers to explore available data, see overlaps among datasets, build new cohorts, and determine if there are data or samples available for additional research projects. Incubated within the Department of Biomedical and Health Informatics, Arcus connects the hospital’s clinical and research data to enable biomedical researchers to conduct highly innovative, data-driven, reproducible research within a managed, scalable framework. This framework includes 1) user access controls; 2) patient privacy and confidentiality protections through regulatory review; 3) electronic honest-brokered data de-identification and re-identification; and 4) data retention, management, sharing, and destruction services in an auditable computational environment.
Measures
Child demographic and clinical data
Data were abstracted from our electronic health record and local registries on the following demographic and clinical variables: sex, extracardiac conditions (i.e. chromosomal anomaly, extracardiac anomaly, syndrome), pre-operational complications (i.e. preoperative cardiac arrest, extracorporeal membrane oxygenation, mechanical ventilation, seizures, or inotropes), STAT category (i.e. Society of Thoracic Surgeons and European Association for Cardiothoracic Surgery Congenital Heart Surgery Mortality Score), Reference O’Brien, Clarke and Jacobs40 and cardiopulmonary bypass times (i.e. total minutes for the first STAT category operation in the hospitalisation).
Prenatal maternal mental health symptoms
Prenatal maternal mental health screening provided data on depressive and post-traumatic stress symptoms of mothers of children diagnosed with CHD using the Postpartum Depression Screening Scale and the Impact of Events Scale—Revised. The Postpartum Depression Screening Scale is a 35-item self-report instrument that assesses for depressive risk in postpartum women. Reference Beck and Gable41 The Postpartum Depression Screening Scale is a reliable and valid measure that has been used with women across the perinatal period. Reference Beck and Gable41,Reference Pereira, Bos and Marques42 Scores of at least 60 signify significant depressive symptoms, and scores of at least 80 screen positive for likely postpartum depression. The Impact of Events Scale—Revised is a 22-item self-report instrument that assesses subjective distress caused by traumatic events. Reference Weiss, Marmar, Wilson and Keane43 Pregnant mothers were asked to identify a specific stressful life event (in this case, at the time of the foetal CHD diagnosis) and then indicate their level of distress for each “difficulty” listed during the past seven days. The Impact of Events Scale—Revised yields a total score (ranging from 0 to 88), and subscale scores can also be calculated for the Intrusion, Avoidance, and Hyperarousal subscales. Reference Pillai, Aigalikar, Vishwasrao and Husainy44 Scores of at least 23 signify significant post-traumatic stress symptoms and of at least 33 screen positive for likely post-traumatic stress disorder. Although both instruments have cut-off scores for screening, we used continuous Postpartum Depression Screening Scale and Impact of Events Scale—Revised scores to capture the extent of depressive and post-traumatic stress symptoms, since categories with low frequencies are difficult to fit into a latent model due to limited distributional information.
Maternal prenatal social determinants of health
Due to the lack of standardised data elements measuring social determinants of health in the electronic health record, Reference Li, Mowery and Ma45 we incorporated individual-level social determinants of health and neighbourhood-level social determinants of health variables that were available in the electronic health record to provide a more comprehensive social determinants of health risk profile. Reference Berg, Doktorchik, Quan and Saini46
Individual-level social determinants of health variables
Demographic data from the child’s mothers were pulled from our local data registry of prenatally diagnosed mothers cared for in our foetal care centre. We pulled data available from the mothers’ first prenatal visit, including maternal age, race and ethnicity, insurance type, education level, and home address.
Neighbourhood-level social determinants of health variable
We geocoded patients’ home addresses at the time of ICU admission and defined their residential census tracts to represent their neighbourhood-level socioeconomic risk. We identified the neighbourhood socioeconomic characteristics using the US Centre for Disease Control and Prevention’s 2014, 2016, and 2018 social vulnerability index. Reference Flanagan, Gregory, Hallisey, Heitgerd and Lewis47 The social vulnerability index was matched with the year of the mother’s prenatal visit: 2014–2015 visits with social vulnerability index 2014, 2016–2017 visits with social vulnerability index 2016 and 2018 visits with social vulnerability index 2018. The social vulnerability index is a ranking system that evaluates a community’s resilience following public health disasters across all US census tracts as percentiles. It is based on 15 social factors in four domains, including socioeconomic status, household composition-disability, minority status-language, and housing type-transportation. A higher score indicates that a community is more vulnerable. We classified the overall social vulnerability index score as low (< 25th percentile), moderate (25–49th percentile), substantial (50–74th percentile), and high (≥ 75th percentile).
Pregnancy outcomes
Child gestational age and weight at birth were abstracted from our electronic health record to examine as pregnancy outcomes.
Child hospitalisation outcome
Hospital length of stay (HLOS) was calculated in days using the hospital admission and hospital discharge dates abstracted from our electronic health record.
Child neurodevelopment outcome
Neurodevelopment was assessed using the Bayley Scales of Infant and Toddler Development, 3rd Edition. Reference Bayley48 This is a standardised measure that evaluates cognitive, language, and motor skills through age-appropriate materials and is administered by a trained developmental practitioner. Items are administered until a basal and ceiling are obtained, and standard scores are calculated for each domain based on age-grouped normative samples. Neurodevelopmental assessment of children in our sample was completed in the recommended window between 13 and 29 months of age. Reference Ware, Butcher and Latal49
Analysis
We described mother and infant characteristics using descriptive statistics. We examined possible selection bias by comparing infant characteristics between our study sample and the entire infant cohort using Chi-squared test, t-test, and Mann–Whitney U test. To identify maternal classes, first, we selected maternal demographic and social determinants of health characteristics as age, race and ethnicity, insurance type, education, and the social vulnerability index overall score based on their significant associations with any infant outcomes (i.e. gestational age at birth, birth weight, HLOS, neurodevelopment) in bivariate regression analysis (p < 0.05). Using latent profile analysis, including the above five characteristics and Postpartum Depression Screening Scale and Impact of Events Scale—Revised scores, we identified the two distinct maternal classes. Latent profile analysis is a probabilistic modelling algorithm that allows clustering of data and statistical inference. Once the clusters are mathematically identified, the classes are homogeneous within but distinct from each other. Regarding the optimal number of classes, we consider both model fit indices and clinical interpretability. The separation of the classes should be meaningful from a clinical standpoint. Reference Sinha, Calfee and Delucchi50 As model complexity increases, the number of observations in each class inevitably gets smaller, potentially decreasing the generalizability of clustering. We selected two classes (Bayesian Information Criterion = 6964.15) since the three classes (Bayesian Information Criterion = 6875.64) ended up with low generalizability (less than 10% in the smallest group; n = 17, 72, and 136). Reference Sinha, Calfee and Delucchi50 We conducted a sensitivity analysis with three maternal classes. We also duplicated the latent profile analysis with depressive and post-traumatic stress cut-off scores instead of continuous Postpartum Depression Screening Scale and Impact of Events Scale—Revised scores, and these models yielded trivial differences in effect estimates.
To examine the association between maternal classes and pregnancy, hospitalisation, and neurodevelopmental outcomes, we used continuous scores for all outcome measures. We used regression models after adjusting for infant sex as a biological variable for all outcomes. Due to the small sample size, we refrained from adjusting for several covariates. For the hospitalisation and neurodevelopmental outcomes, we chose two covariates in addition to infant sex that provided context on the infant’s pre- and post-operative trajectory: pre-operational complications and STAT category. We did not adjust for cardiopulmonary bypass time due to its collinearity with STAT categories (p < 0.001). STAT category and HLOS (<14 days, 14–29 days, and ≥30 days) showed a dose-response relationship (e.g. proportions of HLOS <14d across STAT categories 1 through 5: 43, 59, 21, 18%, and 7%, p < 0.001). Finally, we did not include extracardiac conditions (chromosomal anomaly, extracardiac anomaly, syndrome) in the model because the high-risk mothers were more likely to have extracardiac conditions than the low-risk mothers (59 vs. 39%, p = 0.02). Given the temporal relationship of prenatally diagnosed fetal complications and maternal mental health, maternal classes inherently capture the impact of previously observed foetal extracardiac conditions. Therefore, adjusting for foetal extracardiac conditions in the model could result in over-adjustment when analysing neurodevelopmental outcomes. However, we did perform a sensitivity analysis in the subset of mother-infant dyads without foetal extracardiac conditions to determine whether observed relationships remained significant in that subgroup (n = 133).
Results
Among 189 mothers, the majority were non-Hispanic White (n = 135, 71%), and 79% of them were privately insured (Table 1). The mean gestational age at the maternal mental health screening was 29.3 weeks. Mean scores on the Postpartum Depression Screening Scale were 63.3 (Standard Deviation, SD = 21.7) and for the Impact of Events Scale—Revised were 18.1 (SD = 16.4). Mothers represented a full range of social vulnerability index categories, with approximately one-third of mothers in the substantial or high social vulnerability categories, one-third in the moderate social vulnerability category, and one-third in the low social vulnerability category. Infants in our cohort underwent surgeries representing a full range of STAT categories, Reference O’Brien, Clarke and Jacobs40 44% had extracardiac conditions related to chromosome anomalies, extracardiac anomalies, or syndromes, and 31% had pre-operational complications, such as preoperative ventilation, stroke, inotropes, and sepsis (Table 1). There was no selection bias in infant demographics and clinical characteristics when comparing our sample (n = 189) to the entire infant cohort (n = 647).
Table 1. Maternal and infant characteristics

* PDSS = Postpartum Depression Screening Scale; IESR = Impact of Events Scale—Revised; SVI = Social Vulnerability Index; CPB = cardiopulmonary bypass.** No one had pre-operational ECMO.***Bayley measures had 3% missing in motor domain and 1% missing in language domain.
A latent profile analysis identified two distinct groups that differed in maternal mental health measures and demographics: a high-risk group (n = 46) and a low-risk group (n = 143; Figure 2). Mothers in the high-risk group had significantly higher Impact of Events Scale—Revised scores (mean (SE): 23.2 (2.7) vs. 16.3 (1.4)) and Postpartum Depression Screening Scale scores (mean (SE): 69.7 (3.6) vs. 61.1 (1.9)) than those in the low-risk group. They were also younger (24.5 (0.7) vs. 32.6 (0.4)), had a higher social vulnerability index (47.9 (4.4) vs. 37.7 (2.5)), lower education levels (less than college: 91% vs. 14%), and were more likely to receive Medicaid (63% vs. 6%) and belong to a racial/ethnic minority group (60% vs. 17%).

Figure 2. Two defined maternal classes by latent profile analysis. PDSS = postpartum depression screening scale; IESR = impact of events scalerevised; SVI = social vulnerability index.
Infants born to high-risk mothers had a shorter gestational age (Beta (SE) = –0.45 (0.21) days) and lower weight (–0.21 (0.09) kg) at birth. They also had a longer HLOS (24.5 (6.9) days), and lower cognitive (–4.3 (2.1)), language (–8.6 (2.8)), and motor scale scores (–5.4 (2.5)) than infants in the low-risk group (all p < 0.05; Table 2). A sensitivity analysis, which compared three maternal classes, also showed consistent results, although some significant differences disappeared (Supplementary Table 1). We also performed a sensitivity analysis in the subgroup of mother infant dyads without foetal extracardiac conditions and found that among this group of dyads (n = 133), infant born to high-risk mothers had a shorter gestational age (Beta (SE) = –0.62 (.27); p = 0.02) and lower language scale scores (Beta (SE) = –9.9 (3.4); p = 0.004).
Table 2. Differences in child health outcomes between high-risk mothers (n = 46) and low-risk mothers (n = 143)

Reference group in the regression model: Low-risk mothers. Adjusted models included for infant sex for all outcomes. For hospital length of stay (HLOS) and Bayley measures, we adjusted for pre-operational complications and STAT category additionally.
Discussion
Despite decades of research describing neurodevelopmental outcomes of children born with CHD, Reference Marino, Lipkin and Newburger1 only recently have psychosocial factors such as maternal mental health and social determinants of health been recognised as important contributors. Reference Sood, Newburger and Anixt2,Reference Sood, Jacobs and Marino51 Grounded in the Parental Stress and Resilience in CHD model, this exploratory latent profile analysis is the first to describe an at-risk group of mother–child dyads based on prenatal maternal mental health symptoms in combination with social determinants of health factors. We found worse outcomes in children with CHD born to pregnant mothers experiencing maternal mental health and social determinants of health adversities, including higher depressive and traumatic stress symptoms, lower age, higher social vulnerability index, lower education, and more likely to have Medicaid insurance and represent a minority group. These latent profiles of high- versus low-risk mother groups were differentially associated with worse pregnancy outcomes, including lower gestational age and birth weight, worse child hospitalisation outcome of longer HLOS, and worse child neurodevelopmental outcomes by lower cognitive, language, and motor scales. Significant findings persisted for gestational age and language scale scores in the sensitivity analysis of the subgroup without foetal extracardiac conditions, indicating that maternal mental health and social determinants of health risks have an influence on outcomes that should be explored in future research.
We acknowledge that multiple factors contribute to the presence of neurodevelopmental impairment in children with CHD, including prenatal abnormalities of the central nervous system (e.g. periventricular leucomalacia, smaller head circumference), Reference Marino, Lipkin and Newburger1,Reference Licht, Shera and Clancy52 the complexity of CHD, Reference Wernovsky53 underlying genetic abnormalities, Reference Newburger, Sleeper and Bellinger54–Reference Brosig, Bear and Allen58 low birth weight, Reference Newburger, Sleeper and Bellinger54,Reference Fuller, Nord and Gerdes56,Reference Gaynor, Stopp and Wypij57 and longer hospital lengths of stay. Reference Newburger, Sleeper and Bellinger54,Reference Fuller, Nord and Gerdes56,Reference Gaynor, Stopp and Wypij57,Reference Bellinger, Watson and Rivkin59–Reference Sananes, Manlhiot and Kelly61 However, these factors have been shown to explain only one-third of the variance contributing to worse neurodevelopmental outcomes in CHD. Reference Gaynor, Wernovsky and Jarvik62–64 Thus, it is important to identify new, modifiable risk factors to optimise neurodevelopmental outcomes for the growing population of CHD survivors. The use of latent profile analysis provides a greater understanding of the combination of maternal mental health and social determinants of health risks that may influence pregnancy and child outcomes in CHD. More research is needed to understand whether any causal pathways or direct associations exist between the variables examined. Although there are likely complex interactions that can be studied further, of clinical importance is noting that at-risk mothers and their children with CHD need supportive intervention to mitigate maternal mental health and social determinants of health risks. Mothers of children born with CHD are known to be at risk for maternal mental health adversity and profound psychological distress. Reference Woolf-King, Anger, Arnold, Weiss and Teitel16 Yet, routine mental health screening and social care screening are not the universal standard of care in paediatrics. Reference Kasparian, Kan and Sood15,Reference Sokol, Ammer and Stein65,Reference LeLaurin, De La Cruz and Theis66 Mothers experiencing mental health or social determinants of health adversities need to be identified and supported during routine, family-centred, paediatric practice.
Future research with larger sample sizes may uncover the complexities of the relationships and underlying mechanisms influencing child clinical and neurodevelopmental outcomes. For example, prenatal maternal mental health and social determinants of health adversities may negatively influence the growth and development of the foetus in utero, thereby placing the child at greater risk for negative outcomes after birth. Maternal mental health and social determinants of health adversities experienced by mothers prenatally may have persisted postnatally, which can further impact both the mother and the child. The negative impact of untreated parental mental health problems on long-term child cognitive and behavioural development is well established. Reference Lyons-Ruth, Wolfe and Lyubchik67–Reference Tronick and Reck69 Negative child outcomes may be the result of impaired parent–child interactions and changes in typical parenting behaviour. Reference Abidin70 Preliminary evidence suggests maternal mental health has greater influence on infant and child behavioural outcomes than the child’s severity of illness in CHD. Reference Re, Dean, Mullaert, Guedeney and Menahem71–Reference Landolt, Ystrom, Stene-Larsen, Holmstrom and Vollrath73 Interestingly, a meta-analysis of early interventions involving parents of premature infants has demonstrated positive impacts on infant neurodevelopment. Reference Vanderveen, Bassler, Robertson and Kirpalani74 Limited interventional studies have targeted maternal mental health in CHD; Reference Kasparian, Kan and Sood15 however, none have targeted both maternal mental health and social determinants of health. Additional research is needed to demonstrate effective strategies to support outcomes in both mothers and their children born with CHD.
Limitations
We acknowledge the limitations of retrospective studies and that inherent biases may exist in the final analytic sample. Studies have found that patients with social determinants of health adversity have lower rates of CHD prenatal detection and higher incidence of loss to follow up. Reference Davey, Sinha, Lee, Gauthier and Flores26,Reference Krishnan, Jacobs and Morris75 Future studies with larger, prospective samples should explore whether the relationships found in this study are confirmed. We also acknowledge that there may be an intervention effect from the psychosocial services provided to pregnant mothers in our foetal care centre that we could not capture or measure. We did not have access to the mothers’ electronic health record and were limited to maternal data available in the foetal care centre’s data registry. Therefore, we could not include other clinical information on the mothers, such as pregnancy complications like hypertension or gestational diabetes, that may influence the outcomes we examined in this study. There may also be other maternal mental health symptoms that influence outcomes, including anxiety, stress, or other mental health diagnoses that were not covered in this study. We did not explore trajectories of maternal mental health symptoms and whether maternal mental health symptoms preceding the pregnancy or those developing during or in response to pregnancy or during child hospitalisation influenced outcomes. Future longitudinal studies are recommended that can provide a more complete examination of maternal mental health symptoms over time.
Conclusion
Our findings provide preliminary evidence for identifying potential disparities that can inform the design of interventions, improving the identification of at-risk mothers, and targeting individualised treatment approaches for maternal mental health symptoms and social determinants of health risks. Ongoing efforts to strengthen community and legislative engagement focused on supporting maternal mental health and reducing social determinants of health risk to address health care disparities may improve outcomes for children born with CHD and their families.
Supplementary material
The supplementary material for this article can be found at https://doi.org/10.1017/S1047951125100772.
Acknowledgements
None.
Financial support
This work was supported by the Cardiac Centre at Children’s Hospital in Philadelphia. Study data were collected and managed using Arcus resources hosted at the Children’s Hospital in Philadelphia.
Competing interests
The authors declare none.
Ethics standards
The authors assert that all procedures contributing to this work comply with the ethical standards of the relevant United States of America national guidelines on human experimentation and with the Helsinki Declaration of 1975, as revised in 2008, and has been approved by the Children’s Hospital in Philadelphia’s institutional review board.