Association of congenital heart disease and neurodevelopmental disorders: an observational and Mendelian randomization study

Abstract

Background

This study aims to thoroughly study the connection between congenital heart disease (CHD) and neurodevelopmental disorders (NDDs) through observational and Mendelian randomization (MR) designs.

Methods

This observational study uses data from the National Survey of Children’s Health (2020–2021). Multivariable logistic regression and propensity score matching (PSM) were performed to analyze the association. PSM was used to minimize bias for covariates such as age, race, gender, maternal age, birth weight, concussion or brain injury, preterm birth, cerebral palsy, Down syndrome, and other inherited conditions. In MR analyses, inverse variance-weighted measures, weighted median, and MR-Egger were employed to calculate causal effects.

Results

A total of 85,314 children aged 0–17 were analyzed in this study. In regression analysis, CHD (p = 0.04), the current heart condition (p = 0.03), and the severity of current heart condition (p < 0.05) had a suggestive association with speech or language disorders. The severity of current heart condition (p = 0.08) has a potential statistically significant association with attention deficit hyperactivity disorder(ADHD). In PSM samples, ADHD(p = 0.003), intellectual disability(p = 0.012), and speech or language disorders(p < 0.001) were all significantly associated with CHD. The severity of current heart condition (p < 0.001) also had a significant association with autism. MR analysis did not find causality between genetically proxied congenital cardiac malformations and the risk of NDDs.

Conclusions

Our study shows that children with CHD have an increased risk of developing NDDs. Heart conditions currently and severity of current heart conditions were also significantly associated with these NDDs. In the future, we need to try more methods to clarify the causal relationship between CHD and NDDs.

Introduction

Congenital heart disease (CHD) is a common congenital structural malformation and a leading cause of infant mortality [1]. Advances in surgical and intensive care have significantly increased the survival rate of children with CHD, including those with critical CHD, to 82.5% [2]. With the prolongation of lifespan in children with CHD, the focus of current research has shifted, and neurodevelopmental disorders (NDDs), which have a major impact on their long-term quality of life have received increasing attention. In 2012, the American Heart Association issued a statement aimed at improving the assessment of and attention to neurodevelopmental outcomes in children with CHD [3]. Despite numerous studies on this topic, the exact conclusion for the neurodevelopmental outcomes of CHD is still limited to date. Even the causal association between CHD and NDDs remains unclear.

Due to the low incidence rate of CHD, the sample size of many studies is small [4, 5]. Besides, some studies only focused on children with complex CHD, such as single and left ventricular hypoplasia syndromes [6, 7]. The generalizability of conclusions from these studies may be limited or not applicable to children with more common mild or moderate cardiac malformations. Recently, Tsao [8] and Sigman et al. [9] studies with a large number of samples have found a link between CHD and NDDs, but the populations of children in these studies were limited to the medical system, which may cause selection bias. Similarly, ascertainment bias cannot be ignored when children are exposed to hospitals or treatments frequently. Additionally, many confounding factors, such as the type of cardiac malformation [10], surgical methods [11], postoperative complications [12], and gene mutations [13], can influence developmental disorders in children with CHD, which also reduces the reliability and stability of the study conclusion to varying degrees.

Therefore, we conducted a large nationwide observational study to investigate the association between CHD and NDDs. To control for confounding, we utilized propensity score matching (PSM) and regression analyses. Additionally, we performed a two-sample Mendelian randomization (MR) study, which is a statistical genetics method that draws causal inferences between exposure and outcome by utilizing genetic variants from the whole genome [14]. Since genetic variants are assigned randomly at meiosis and before disease onset, MR studies can minimize confounding biases and reverse causation [15]. To our knowledge, this is the first MR study investigating the causal relationship between CHD and NDDs. We hypothesized that there is an association and causality between CHD and NDDs in children, and congenital heart malformations increase the risk of NDDs, including autism spectrum disorder (ASD), attention deficit hyperactivity disorder, intellectual disability, and speech or language disorders. Through this comprehensive method, we hope to provide a strong scientific basis for the relationship between CHD and NDD.

Methods

Study design and participants

The present study used cross-sectional data from the 2020–2021 National Survey of Children’s Health (NSCH), which is funded and directed by the Health Resources and Services Administration’s Maternal and Child Health Bureau [16]. The NSCH is designed to generate national and state-level data on the physical and emotional health of children in the United States aged 0–17 years. The NSCH (2020–2021) had a total of 93,669 completed surveys with an overall weighted response rate of 42.4% for 2020 (N = 42,777) and 40.3% for 2021 (N = 50,892). Participants with missing values of study variables were excluded, leaving a final analysis sample of 85,314 children. The flow of participant screening and the frequency of missing values are illustrated in Fig. 1.

Fig. 1

Flow chart showing the patient selection process and analysis in the observational study

Sampling methods and procedures

In this survey, randomly selected households were contacted by mail, web, or paper (offered in both English and Spanish). In each household with children aged < 18 years old, one child was randomly chosen for participation in the survey. The survey oversampled children with special health care needs and children aged 0–5 years. Additional information regarding sampling, administration, and methodology of the NSCH can be found on the Data Resource Center website [17].

Variables and data sources

The interest variable, congenital heart condition, was assessed based on the question, “Was this child born with the condition?” If the respondents answered “yes” to this question, the child was categorized as “yes” and “no” otherwise. In addition, respondents who answered “yes” were further asked whether the child currently has the condition, and, if so, the severity level was assessed with the question, “If yes, is it mild, moderate, or severe?” The variables of the heart condition currently and severity of current heart condition were based on these two questions, respectively. Children with current heart conditions were classified as mild, moderate, or severe based on their responses.

Autism, the outcome variable of the study, was assessed based on the question, “Has a doctor or other health care provider ever told you that this child has autism or ASD? Include diagnoses of Asperger’s disorder or pervasive developmental disorder.” ADHD, another outcome variable, was assessed based on the question, “Has a doctor or other health care provider ever told you that this child has attention deficit disorder or attention deficit/hyperactivity disorder, that is, ADD or ADHD?” Intellectual disability, another outcome variable, was assessed based on the question, “Has a doctor, other health care provider, or educator ever told you that this child has an intellectual disability (formerly known as mental retardation)?” Speech or language disorder, another outcome variable, was assessed based on the question, “Has a doctor, other health care provider, or educator ever told you that this child has a speech or other language disorder?” The child was categorized as “yes” or “no” based on the responses to all the above-mentioned questions.

Based on the risk factors for NDDs, we examined the following covariates as potential confounding variables: maternal age, birth weight (in ounces), concussion or brain injury, preterm birth, cerebral palsy, Down syndrome, and other genetic or inherited conditions. Demographic variables including age (in years), race (classified as Hispanic; White, non-Hispanic; Black, non-Hispanic; Asian, non-Hispanic; American Indian or Alaska Native Non-Hispanic; Multi-Race Non-Hispanic), and gender were also included as covariates.

MR analysis

We conducted a two-sample MR analysis using summary-level genome-wide association study statistics data. The datasets for genetic associations with congenital heart conditions were obtained from FinnGen, which included 3,42,499 individuals of European ancestry (3,459 cases and 3,39,040 controls). Instrumental variables were selected based on the following criteria: (i) associated at the level of genome-wide significance (p < 5 × 10^− 6); (ii) linkage disequilibrium r² < 0.001; (iii) < 5000 KB from the index variant. A total of 14 single-nucleotide polymorphisms (SNPs) were identified and listed in Supplementary Table 1.

The outcome data for NDDs were obtained from different sources, including the Integrative Psychiatric Research-Psychiatric Genomics Consortium (18,382 cases and 27,969 controls) for ASD [18], the Psychiatric Genomics Consortium (20,183 cases and 35,191 controls) for ADHD [19], the Social Science Genetic Association Consortium (12,441 individuals) for childhood intelligence [20], and FinnGen (2,288 cases and 25,4976 controls) for speech and linguistic disorders (Supplementary Table 2). European pedigree population data were restricted to reduce the bias caused by population level.

Statistical analysis

Observational analyses were conducted based on NSCH 2020–2021 data. Continuous data were expressed as mean ± standard deviation or median (interquartile range) (for skew distributional data). The t-test or nonparametric rank-sum test, was performed to compare the differences between the groups with and without NDDs. Categorical data are described by frequency, and the differences between groups were compared using χ2 or Fisher exact test. To estimate the association between each factor and outcome, univariate and multivariable logistic regression analysis was used, with NDDs as the outcome variable and study and confounding variables (inclusion of NDDs other than outcomes) as explanatory variables. Multivariate adjusted odds ratios (ORs), 95% confidence intervals (Cis), and p-values were calculated. To reduce the impact of confounding variables and selection bias on results, PSM was also used to match case groups to control groups. Demographic variables (age, race, and gender), confounding variables (maternal age, birth weight, concussion or brain injury, preterm birth, cerebral palsy, Down syndrome, and other genetic or inherited conditions), and NDDs (other than outcomes) were used as covariates in the PSM model, using 1:1 nearest-neighbor matching techniques with a 0.01 caliper level. In the matched samples, NDDs were the outcome variable, and heart condition-related variables (congenital heart condition, current heart condition, and its severity) were the explanatory variables. A p-value of < 0.05 indicates a statistically significant difference.

In the MR analyses, we used inverse variance-weighted measures (IVW-MR), weighted median, and MR-Egger to estimate causal effects, and results are presented as ORs with 95% CIs per standard deviation increment in exposures. The F-statistic was employed to assess the strength of the instrument, and a value of at least 10 indicates weak instrument bias [21]. Additionally, we conducted several sensitivity analyses, including MR-Egger regression, MR Pleiotropy Residual Sum and Outlier (MR-PRESSO) method, to investigate the possibility of horizontal pleiotropy (p < 0.05) [22, 23]. The Cochran Q test in the IVW-MR method was performed to determine heterogeneity (p < 0.05) [24]. We further conducted a leave-one-out analysis and made a scatter plot to visually examine possible outliers. All analyses were performed using R (version 4.2.2) and the following packages: dplyr, MatchIt, cobalt, gmodels, ggplot2, knitr, kableExtra, tableone, survey, reshape2, and TwoSampleMR.

Research restrictions

We recognize that there are some restrictions in this study. First of all, due to the relatively low incidence of severe CHD, the sample volume of some study variables may be small. Secondly, the data sources adopted by this research may have problems with insufficient information in some aspects, such as surgical methods for CHD. In the end, although MR analysis can reduce confusion, there are still some assumptions and limitations. These restrictions will be fully discussed in the conclusion and provide inspiration for future research.

Ethical statement

This study was approved by the Ethics Committee of Anhui children’s Hospital (approval no. EYLL-2022-020). This study involved secondary analysis using data from genome-wide association studies and the National Survey of Children’s Health (2020–2021). The participants had granted written informed consent before either study started. We confirm that all methods were performed in accordance with the ethical standards as laid down in the Declaration of Helsinki and its later amendments or comparable ethical standards.

Results

Observational study

A total of 85,314 children and adolescents aged 0–17 were included in the analysis, of whom 1,977 participants had CHD (2.3%), 2,581 had ASD (3.0%), 8,198 had ADHD (9.6%), 877 had intellectual disability (1%), and 7,141 had speech or language disorders (8.3%). The baseline characteristics of children with and without NDDs are presented in Table 1. The Kruskal-Wallis rank sum test and χ2 test indicated significant differences in covariates between groups with and without NDDs (p < 0.001). Then we conducted a multivariable logistic regression analysis to explore the association between congenital heart conditions and the risk of NDDs (Table 2 and Supplementary Table 3). For ADHD, the severity of current heart condition (OR, 2.4; 95% CI, 1.4–3.4; p = 0.08) has a potential statistically significant association with it. For speech or language disorders, regression analysis showed that congenital heart disease (OR, 1.61; 95% CI, 1.14–2.08; p = 0.04), the current heart condition (OR, 2.45; 95% CI, 1.65–3.25; p = 0.03), and the severity of current heart condition (p = 0.01 for mild, p = 0.02 for moderate) had a suggestive association with it. However, these variables were not associated with ASD or intellectual disability (p > 0.05) in multivariable regression analyses. Notably, the association between current severe heart conditions and NDDs was not available due to the small number of cases in groups.

Table 1 Baseline characteristics of children according to the presence of neurodevelopmental disorders

Table 2 The association between congenital heart condition and the risk of neurodevelopmental disorders in multivariable regression analyses and propensity score matching

To address possible bias caused by confounding variables, we constructed a 1:1 PSM model. After including these covariates in the model, we successfully matched 2,510 pairs for the ASD group, 8,063 pairs for the ADHD group, 780 pairs for the intellectual disability group, and 6,492 pairs for the speech or language disorder group. As shown in Supplementary Table 4, most baseline characteristics did not have statistically significant differences between the groups (p > 0.05), indicating that the PS-matched samples achieved balance. For ASD, the severity of current heart condition (p < 0.001) had a significant association with it in PS-matched samples (Table 2 and Supplementary Table 4). For ADHD, the variables of congenital heart disease (p = 0.003), heart condition currently (p = 0.002), and severity of current heart condition (p < 0.001) were all associated with it. Similarly, intellectual disability was significantly related to congenital heart disease (p = 0.012), heart condition currently (p = 0.033), and severity of current heart condition(p = 0.033). For speech or language disorders, congenital heart disease (p < 0.001), heart condition currently (p < 0.001), and severity of current heart condition (p = 0.001) were also significantly associated with it, consistent with the result of multivariable regression analysis.

MR analysis

IVW-MR analysis did not find a causal relationship between genetically proxied congenital cardiac malformations and the risk of NDDs. This result was consistent with the findings from the weighted median and MR-Egger methods (Fig. 2 and Supplementary Table 5). The F-statistics of all SNPs used as instrument variants were > 20, indicating no significant weak instrument bias (Supplementary Table 1). Additionally, sensitivity and heterogeneity analyses demonstrated that there was no significant evidence of pleiotropy or heterogeneity (Supplementary Table 6), indicating the stability and reliability of the MR results. The leave-one‐out analysis and scatter plots also did not identify any outlying SNPs (Supplementary Figs. 1–4).

Fig. 2

Association of genetically proxied congenital cardiac malformations with four neurodevelopmental disorders. Odds ratios are per standard deviation increment in the exposure. IVW, inverse-variance-weighted; MR, Mendelian randomization

Discussion

In this large, nationally representative study, our results provided positive evidence of an association between CHD and NDDs through regression analysis and PSM. In addition, our results also found that heart condition currently and the severity of the current heart condition are significantly associated with these neurodevelopmental disorders.

Our results are consistent with some of the previous studies [4,5,6,7,8,9,10,11,12,13]. Previous small sample studies have already found that children with left ventricular hypoplasia syndromes have a significantly higher risk of being diagnosed with autism than children without CHD [4,5,6]. Specifically, 4 out of 58 children with left ventricular hypoplasia syndrome were diagnosed with ASD, while none of the control group were diagnosed. However, Wier et al. [25] a retrospective case–control study (ASD n = 417; control n = 2,067) yielded a contradictory report. They did not find a significant association between ASD and cardiac defects (OR = 1.5, 95% CI = 0.7–2.8). It is possible that this study were not well remove the influence of confounding factors on the results, especially they include all congenital malformations.

To address this potential reporting bias, Tsao et al. [8] conducted a retrospective case-control study with 3,552 CHD cases and 14,208 Non-CHD controls. Multivariate Cox regression analysis showed that children with CHD had an increased risk of developing ADHD (hazard ratio [HR] 2.52, 95% CI 1.96–3.25) and ASD (HR 1.97, 95% CI 1.11–3.52). Recently, using a large Military Health System (MHS) administrative database with 8,760 autism spectrum disorders cases and 26,280 controls, Sigman et al. [9] also demonstrated an increased odds of ASD in patients with CHD (OR 1.32; 95% CI 1.10–1.59). However, these previous studies only focused on the association between developmental delays and CHD but did not clarify the causal relationship between them.

To further obtain causal associations, we conducted a two-sample MR study. However, MR analysis did not find causality between genetically proxied congenital cardiac malformations and the risk of NDDs. This result of MR may seem counterintuitive at first sight when compared to our prior studies. Therefore, some interpretations of this contradictory finding should be listed before drawing our final conclusions.

First, in MR studies, we did not include the impact of external intervention measures (such as surgery). This is necessary for the survival of children, but it may also complicate neurodevelopment. These interventions generate various factors, including preoperative and postoperative brain injury and intraoperative hemodynamic factors, among others [11, 12, 26]. Indeed, previous studies have largely identified the significant impact of surgery on neurodevelopmental outcomes in children. However, because the genetic effect is lifelong, the MR study cannot estimate the impact of intervention measures at a specific time in life. This may result in the effect sizes of exposure inferred by genetic analysis being different from those in reality. In this MR study, we only explored the impact of congenital heart malformation itself on children, which may imply that the impact of external interventions such as surgery on neurological development may not be significantly reflected.

Second, pleiotropy in genetic loci may contribute to the risk of both disorders, but it needs to be avoided in MR analysis. During early morphogenesis, the brain shares the same genetic pathway with the heart, and the genetic factors that cause heart malformations may also disrupt the normal expression of gene fragments for brain development, thereby triggering congenital brain injury and NDDs [27]. In 2015, a meta-analysis showed that the incidence of CHARGE syndrome with ASD was 30%, and the incidence of 22q11.2 deletion syndrome with ASD was 11% [28]. Similarly, a review by Vorstman et al. [29] summarizes that at least 10 chromosomal abnormalities can lead to heart defects in patients with ASD. Genetic factors are an important cause of NDDs in children with CHD. However, in the MR study, all genetic variations with pleiotropic effects on the outcome cannot be used as instrumental variables for proxy exposure [30]. Instrument variables can only affect the outcome through exposure instead of other pathways [31]. This exclusivity assumption is essential for the MR design to ensure the validity of the causal assessment. Therefore, our results of MR did not preclude an association between both disorders, as cardiac malformations may increase the risk of NDDs via genetic factors.

Thirdly, it should be emphasized that genetic changes are also important in CHD and NDDs. Recent studies have shown that genes related to autism with damaging de novo mutations have been identified in patients with CHD, which highlighting the overlap between the development of the heart and neurodevelopment [13, 32]. However, the penetrance of these genetic mutations (include point mutations, copy number variations, and aneuploidy) varies. And their contribution to NDDs in patients with CHD remains unclear. Therefore, it is difficult to accurately evaluate this factor or conduct layered analysis in MR study, which may partly explain the contradiction results between our research. In summary, although the MR study did not fully elucidate causality between CHD and NDDs, it does not mean denying the impact of CHD on children’s neurodevelopment. In the future, controlled trials and longitudinal follow-up cohorts are needed to infer causality between them.

As discussed previously, our study shows that children with CHD have an increased risk of developing NDDs. The heart conditions currently and severity of current heart conditions were also significantly associated with these NDDs. This means that compared with normal children, the quality of life and social participation ability of CHD population may be lower than expected, and the burden of mental health is higher, especially in children with severe cardiac malformations. Therefore, pediatricians should not only pay attention to the circulation function of children, but also pay attention to the impact of the complexity of heart disease on their neurodevelopment. For children with severe cardiac malformations, comprehensive neurodevelopmental evaluation or treatment should be provided as soon as possible. In addition, CHD children should also undergo routine neurodevelopmental screening and regular follow-up, so that each CHD child with NDDs can benefit from timely and effective individualized rehabilitation treatment or special education according to their condition, and ultimately improve the neurodevelopmental prognosis of CHD children and the quality of life.

The advantages of our study include the use of a large, nationally representative dataset and the sample was selected from nationwide households instead of the medical system. As a result, our results are less biased and can be generalized to a larger population of children and adolescents. Additionally, we conducted an MR study using summary-level genome-wide association study data. MR analysis can overcome some limitations of observational studies, such as measurement or confounding errors, and facilitate evaluating the long-term effects of cardiac malformation.

There are some limitations of this study. First of all, we cannot rule out that we may have overlooked some associations due to inadequate statistical power. Although our research samples are large, the number of children with complex heart malformations is still not enough to perform some statistical analysis. In addition, the lack of causality in MR research does not exclude the biological effect of CHD, but rather indicates a lack of evidence that changes in the expression of genes are related to NDDs. Additionally, the measurement of some covariates was based on parental reports in our observational study. However, parental reports may be prejudiced or may not accurately reflect the child’s health status. Finally, the genetic summary data we used for the MR study were derived from populations with European ancestry, which may limit the universality of our results in other races.

Conclusions

Our study shows that children with CHD have an increased risk of developing NDDs. Heart conditions currently and severity of current heart conditions were also significantly associated with these neurodevelopmental disorders. In the future, we need to try more methods to clarify the causal relationship between CHD and NDDs.