Research

Congenital heart disease (CHD) is a leading cause of infant mortality. We analyzed de novo mutations (DNMs) and very rare transmitted/unphased damaging variants in 248 prespecified genes in 11,555 CHD probands. The results identified 60 genes with a significant burden of heterozygous damaging variants. Variants in these genes accounted for CHD in 10.1% of probands with similar contributions from de novo and transmitted variants in parent-offspring trios that showed incomplete penetrance. DNMs in these genes accounted for 58% of the signal from DNMs. Thirty-three genes were linked to a single CHD subtype while 12 genes were associated with 2 to 4 subtypes. Seven genes were only associated with isolated CHD, while 37 were associated with 1 or more extracardiac abnormalities. Genes selectively expressed in the cardiomyocyte lineage were associated with isolated CHD, while those widely expressed in the brain were also associated with neurodevelopmental delay (NDD). Missense variants introducing or removing cysteines in epidermal growth factor (EGF)-like domains of NOTCH1 were enriched in tetralogy of Fallot and conotruncal defects, unlike the broader CHD spectrum seen with loss of function variants. Transmitted damaging missense variants in MYH6 were enriched in multiple CHD phenotypes and account for  1% of all probands. Probands with characteristic mutations causing syndromic CHD were frequently not diagnosed clinically, often due to missing cardinal phenotypes. CHD genes that were positively or negatively associated with development of NDD suggest clinical value of genetic testing. These findings expand the understanding of CHD genetics and support the use of molecular diagnostics in CHD.

Examining the altered arrangement and patterning of sulcal folds offers insights into the mechanisms of neurodevelopmental differences in psychiatric and neurological disorders. Previous sulcal pattern analysis used spectral graph matching of sulcal pit-based graph structures to assess deviations from normative sulcal patterns. However, challenges exist, including the absence of a standard criterion for defining a typical reference set, time-consuming cost of graph matching, user-defined feature weight sets, and assumptions about uniform node distribution. We developed a deep learning-based sulcal pattern analysis to address these challenges by adapting prototype-based graph neural networks to sulcal pattern graphs. Additionally, we proposed a prototype inverse-projection for better interpretability. Unlike other prototype-based models, our approach inversely projects prototypes onto individual node representations to calculate the inverse-projection weights, enabling efficient visualization of prototypes and focusing the model on selective regions. We evaluated our method through a classification task between healthy controls (n = 174, age = 15.4 ±1.9 [mean ± standard deviation, years]) and patients with congenital heart disease (n = 345, age = 15.8 ±4.7) from four cohort studies and a public dataset. Our approach demonstrated superior classification performance compared to other state-of-the-art models, supported by extensive ablative studies. Furthermore, we visualized and examined the learned prototypes to enhance understanding. We believe our method has the potential to be a sensitive and understandable tool for sulcal pattern analysis.

Variants with large effect contribute to congenital heart disease (CHD). To date, recessive genotypes (RGs) have commonly been implicated through anecdotal ascertainment of consanguineous families and candidate gene-based analysis; the recessive contribution to the broad range of CHD phenotypes has been limited. We analyzed whole exome sequences of 5,424 CHD probands. Rare damaging RGs were estimated to contribute to at least 2.2% of CHD, with greater enrichment among laterality phenotypes (5.4%) versus other subsets (1.4%). Among 108 curated human recessive CHD genes, there were 66 RGs, with 54 in 11 genes with >1 RG, 12 genes with 1 RG, and 85 genes with zero. RGs were more prevalent among offspring of consanguineous union (4.7%, 32/675) than among nonconsanguineous probands (0.7%, 34/4749). Founder variants in GDF1 and PLD1 accounted for 74% of the contribution of RGs among 410 Ashkenazi Jewish probands. We identified genome-wide significant enrichment of RGs in C1orf127, encoding a likely secreted protein expressed in embryonic mouse notochord and associated with laterality defects. Single-cell transcriptomes from gastrulation-stage mouse embryos revealed enrichment of RGs in genes highly expressed in the cardiomyocyte lineage, including contractility-related genes MYH6, UNC45B, MYO18B, and MYBPC3 in probands with left-sided CHD, consistent with abnormal contractile function contributing to these malformations. Genes with significant RG burden account for 1.3% of probands, more than half the inferred total. These results reveal the recessive contribution to CHD, and indicate that many genes remain to be discovered, with each likely accounting for a very small fraction of the total.

BACKGROUND: SMAD2 is a coregulator that binds a variety of transcription factors in human development. Heterozygous SMAD2 loss-of-function and missense variants are identified in patients with congenital heart disease (CHD) or arterial aneurysms. Mechanisms that cause distinct cardiovascular phenotypes remain unknown. We aimed to define transcriptional and epigenetic effects of SMAD2 variants and their role in CHD. We also assessed the function of SMAD2 missense variants of uncertain significance.

METHODS AND RESULTS: Rare SMAD2 variants (minor allele frequency ≤10-5) were identified in exome sequencing of 11 336 participants with CHD. We constructed isogenic induced pluripotent stem cells with heterozygous or homozygous loss-of-function and missense SMAD2 variants identified in CHD probands. Wild-type and mutant induced pluripotent stem cells were analyzed using bulk RNA sequencing, chromatin accessibility (Assay for Transposase-Accessible Chromatin With Sequencing), and integrated with published SMAD2/3 chromatin immunoprecipitation data. Cardiomyocyte differentiation and contractility were evaluated. Thirty participants with CHD had heterozygous loss-of-function or missense SMAD2 variants. SMAD2 haploinsufficiency altered chromatin accessibility at promoters and dysregulated expression of 385 SMAD regulated genes, including 10 CHD-associated genes. Motifs enriched in differential Assay for Transposase-Accessible Chromatin peaks predicted that SMAD2 haploinsufficiency disrupts interactions with transcription factors NANOG (homeobox protein NANOG), ETS, TEAD3/4 (transcriptional enhanced associate domain 3/4), CREB1 (cAMP response element binding protein 1), and AP1 (activator protein 1). Compared with SMAD2-haploinsufficient cells, induced pluripotent stem cells with R114C or W274C variants exhibited distinct and shared chromatin accessibility and transcription factor binding changes.

CONCLUSIONS: SMAD2 haploinsufficiency disrupts transcription factor binding and chromatin interactions critical for cardiovascular development. Differences between the molecular consequences of loss-of-function and missense variants likely contribute to phenotypic heterogeneity. These findings indicate opportunities for molecular analyses to improve reclassification of SMAD2 variants of uncertain clinical significance.