Conotruncal heart defects (CTDs) are among the most common and severe groups of congenital heart defects. Despite evidence of an inherited genetic contribution to CTDs, little is known about the specific genes that contribute to the development of CTDs. We performed gene-based genome-wide analyses using microarray-genotyped and imputed common and rare variants data from two large studies of CTDs in the United States. We performed two case-parent trio analyses (N = 640 and 317 trios), using an extension of the family-based multi-marker association test, and two case-control analyses (N = 482 and 406 patients and comparable numbers of controls), using a sequence kernel association test. We also undertook two meta-analyses to combine the results from the analyses that used the same approach (i.e. family-based or case-control). To our knowledge, these analyses are the first reported gene-based, genome-wide association studies of CTDs. Based on our findings, we propose eight CTD candidate genes (ARF5, EIF4E, KPNA1, MAP4K3, MBNL1, NCAPG, NDFUS1 and PSMG3). Four of these genes (ARF5, KPNA1, NDUFS1 and PSMG3) have not been previously associated with normal or abnormal heart development. In addition, our analyses provide additional evidence that genes involved in chromatin-modification and in ribonucleic acid splicing are associated with congenital heart defects.

The 22q11.2 deletion syndrome (22q11.2DS; DiGeorge syndrome/velocardiofacial syndrome) occurs in 1 of 4000 live births, and 60% to 70% of affected individuals have congenital heart disease, ranging from mild to severe. In our cohort of 1472 subjects with 22q11.2DS, a total of 62% (n=906) have congenital heart disease and 36% (n=326) of these have tetralogy of Fallot (TOF), comprising the largest subset of severe congenital heart disease in the cohort. To identify common genetic variants associated with TOF in individuals with 22q11.2DS, we performed a genome-wide association study using Affymetrix 6.0 array and imputed genotype data. In our cohort, TOF was significantly associated with a genotyped single-nucleotide polymorphism (rs12519770, P=2.98×10-8) in an intron of the adhesion GPR98 (G-protein-coupled receptor V1) gene on chromosome 5q14.3. There was also suggestive evidence of association between TOF and several additional single-nucleotide polymorphisms in this region. Some genome-wide significant loci in introns or noncoding regions could affect regulation of genes nearby or at a distance. On the basis of this possibility, we examined existing Hi-C chromatin conformation data to identify genes that might be under shared transcriptional regulation within the region on 5q14.3. There are 6 genes in a topologically associated domain of chromatin with GPR98, including MEF2C (Myocyte-specific enhancer factor 2C). MEF2C is the only gene that is known to affect heart development in mammals and might be of interest with respect to 22q11.2DS. In conclusion, common variants may contribute to TOF in 22q11.2DS and may function in cardiac outflow tract development.

The 22q11.2 deletion syndrome (22q11.2DS) results from non-allelic homologous recombination between low-copy repeats termed LCR22. About 60%-70% of individuals with the typical 3 megabase (Mb) deletion from LCR22A-D have congenital heart disease, mostly of the conotruncal type (CTD), whereas others have normal cardiac anatomy. In this study, we tested whether variants in the hemizygous LCR22A-D region are associated with risk for CTDs on the basis of the sequence of the 22q11.2 region from 1,053 22q11.2DS individuals. We found a significant association (FDR p < 0.05) of the CTD subset with 62 common variants in a single linkage disequilibrium (LD) block in a 350 kb interval harboring CRKL. A total of 45 of the 62 variants were associated with increased risk for CTDs (odds ratio [OR) ranges: 1.64-4.75). Associations of four variants were replicated in a meta-analysis of three genome-wide association studies of CTDs in affected individuals without 22q11.2DS. One of the replicated variants, rs178252, is located in an open chromatin region and resides in the double-elite enhancer, GH22J020947, that is predicted to regulate CRKL (CRK-like proto-oncogene, cytoplasmic adaptor) expression. Approximately 23% of patients with nested LCR22C-D deletions have CTDs, and inactivation of Crkl in mice causes CTDs, thus implicating this gene as a modifier. Rs178252 and rs6004160 are expression quantitative trait loci (eQTLs) of CRKL. Furthermore, set-based tests identified an enhancer that is predicted to target CRKL and is significantly associated with CTD risk (GH22J020946, sequence kernal association test (SKAT) p = 7.21 × 10^-5) in the 22q11.2DS cohort. These findings suggest that variance in CTD penetrance in the 22q11.2DS population can be explained in part by variants affecting CRKL expression.