PLoS Genet plos plosgen PLOS Genetics 1553-7390 1553-7404 Public Library of Science San Francisco, CA USA 10.1371/journal.pgen.1009679 PGENETICS-D-20-01841 Research Article Medicine and health sciencesMedical conditionsCardiovascular diseasesCongenital heart defects Medicine and health sciencesCardiologyCardiovascular medicineCardiovascular diseasesCongenital heart defects Medicine and health sciencesMedical conditionsCongenital disordersBirth defectsCongenital heart defects Biology and life sciencesDevelopmental biologyMorphogenesisBirth defectsCongenital heart defects Biology and life sciencesAnatomyCardiovascular anatomyHeart Medicine and health sciencesAnatomyCardiovascular anatomyHeart Biology and life sciencesGeneticsGene expression Computer and information sciencesNetwork analysisProtein interaction networks Biology and life sciencesBiochemistryProteomicsProtein interaction networks Biology and life sciencesGeneticsHeredityHaploinsufficiency Biology and life sciencesGeneticsGenomics Biology and life sciencesGeneticsGenetic loci Physical sciencesMathematicsDiscrete mathematicsCombinatoricsPermutation Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease Integrative analysis of genomic variants in CHD https://orcid.org/0000-0002-9201-7840 Audain Enrique Conceptualization Data curation Formal analysis Writing – original draft Writing – review & editing 1 2 https://orcid.org/0000-0003-2686-1039 Wilsdon Anna Data curation Investigation Writing – original draft 3 Breckpot Jeroen Data curation Investigation Writing – review & editing 4 Izarzugaza Jose M. G. Data curation Formal analysis 5 https://orcid.org/0000-0002-2370-8496 Fitzgerald Tomas W. Methodology Writing – review & editing 6 Kahlert Anne-Karin Data curation Writing – review & editing 1 2 7 https://orcid.org/0000-0001-8247-4020 Sifrim Alejandro Formal analysis Methodology 8 9 https://orcid.org/0000-0002-0978-6076 Wünnemann Florian Methodology Writing – review & editing 10 https://orcid.org/0000-0001-6579-6941 Perez-Riverol Yasset Methodology Writing – review & editing 11 Abdul-Khaliq Hashim Data curation Writing – review & editing 12 https://orcid.org/0000-0003-2762-1002 Bak Mads Data curation Writing – review & editing 13 14 https://orcid.org/0000-0002-0681-7279 Bassett Anne S. Data curation Writing – review & editing 15 16 Benson Woodrow D. Data curation Writing – review & editing 17 https://orcid.org/0000-0001-7881-1557 Berger Felix Data curation Writing – review & editing 18 https://orcid.org/0000-0002-2773-689X Daehnert Ingo Data curation Writing – review & editing 19 https://orcid.org/0000-0002-3651-2548 Devriendt Koenraad Data curation Writing – review & editing 4 https://orcid.org/0000-0002-8455-3416 Dittrich Sven Data curation 20 https://orcid.org/0000-0001-7573-8235 Daubeney Piers EF Data curation 21 https://orcid.org/0000-0002-3778-5927 Garg Vidu Data curation Writing – review & editing 22 23 24 25 https://orcid.org/0000-0003-4181-2120 Hackmann Karl Data curation Writing – review & editing 7 https://orcid.org/0000-0002-3259-2536 Hoff Kirstin Data curation Writing – review & editing 1 2 Hofmann Philipp Data curation Writing – review & editing 1 2 https://orcid.org/0000-0003-2591-2682 Dombrowsky Gregor Data curation Writing – review & editing 1 2 https://orcid.org/0000-0003-0019-6325 Pickardt Thomas Data curation Resources 26 Bauer Ulrike Resources 26 https://orcid.org/0000-0001-9573-0812 Keavney Bernard D. Data curation Writing – review & editing 27 28 Klaassen Sabine Data curation Writing – review & editing 29 30 31 Kramer Hans-Heiner Funding acquisition Resources Writing – review & editing 1 2 Marshall Christian R. Data curation Writing – review & editing 32 33 https://orcid.org/0000-0002-7806-0068 Milewicz Dianna M. Data curation Writing – review & editing 34 Lemaire Scott Data curation Writing – review & editing 35 Coselli Joseph S. Data curation Writing – review & editing 36 Mitchell Michael E. Data curation Writing – review & editing 36 Tomita-Mitchell Aoy Data curation Writing – review & editing 36 https://orcid.org/0000-0001-6341-9624 Prakash Siddharth K. Data curation Writing – review & editing 34 https://orcid.org/0000-0002-7844-7326 Stamm Karl Data curation Writing – review & editing 36 https://orcid.org/0000-0003-2673-9164 Stewart Alexandre F. R. Data curation Writing – review & editing 37 Silversides Candice K. Data curation Writing – review & editing 15 Siebert Reiner Data curation Writing – review & editing 38 39 https://orcid.org/0000-0002-7657-4647 Stiller Brigitte Data curation Writing – review & editing 40 https://orcid.org/0000-0001-5664-7987 Rosenfeld Jill A. Data curation Writing – original draft 17 Vater Inga Data curation Writing – review & editing 39 Postma Alex V. Data curation Writing – review & editing 41 42 https://orcid.org/0000-0003-2157-425X Caliebe Almuth Data curation Writing – review & editing 39 Brook J. David Writing – review & editing 3 https://orcid.org/0000-0002-9047-0502 Andelfinger Gregor Methodology Writing – review & editing 43 https://orcid.org/0000-0002-2333-7015 Hurles Matthew E. Conceptualization Funding acquisition Writing – review & editing 44 https://orcid.org/0000-0002-8772-6845 Thienpont Bernard Conceptualization Writing – review & editing 4 45 https://orcid.org/0000-0002-7536-1172 Larsen Lars Allan Conceptualization Data curation Formal analysis Writing – review & editing 13 * Hitz Marc-Phillip Conceptualization Data curation Formal analysis Writing – review & editing 1 2 39 44 * Department of Congenital Heart Disease and Pediatric Cardiology, University Hospital of Schleswig-Holstein, Kiel, Germany German Center for Cardiovascular Research (DZHK), Kiel, Germany School of Life Sciences, University of Nottingham, University Park, Nottingham, United Kingdom Centre for Human Genetics, Katholieke Universiteit Leuven, Leuven, Belgium Department of Health Technology, Technical University of Denmark, Lyngby, Denmark European Bioinformatics Institute (EMBL-EBI), Wellcome Genome Campus, Cambridge, United Kingdom Institute for Clinical Genetics, Faculty of Medicine Carl Gustav Carus, TU Dresden, Dresden, Germany Department of Human Genetics, University of Leuven, KU Leuven, Leuven, Belgium Sanger Institute-EBI Single-Cell Genomics Centre, Wellcome Trust Sanger Institute, Hinxton, United Kingdom Montreal Heart Institute, Université de Montréal, Québec, Canada European Molecular Biology Laboratory, European Bioinformatics Institute (EMBL-EBI), Wellcome Trust Genome Campus, Hinxton, Cambridge, United Kingdom Clinic for Pediatric Cardiology—University Hospital of Saarland, Homburg (Saar), Germany Department of Cellular and Molecular Medicine, University of Copenhagen, Copenhagen, Denmark Department of Clinical Genetics, Rigshospitalet, Copenhagen University Hospital, Copenhagen, Denmark Toronto Congenital Cardiac Centre for Adults, and Division of Cardiology, Department of Medicine, University Health Network, Toronto, Canada Department of Psychiatry, University of Toronto, Toronto, Canada Department of Pediatrics, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America Department of Congenital Heart Disease—Pediatric Cardiology, German Heart Center Berlin, Berlin, Germany Department of Pediatric Cardiology and Congenital Heart Disease, Heart Center, University of Leipzig, Leipzig, Germany Department of Pediatric Cardiology, University Hospital Erlangen, Friedrich-Alexander-University Erlangen-Nürnberg (FAU), Erlangen, Germany Division of Paediatric Cardiology, Royal Brompton Hospital, London, United Kingdom The Heart Center, Nationwide Children’s Hospital, Columbus, Ohio, United States of America Department of Molecular Genetics, The Ohio State University, Columbus, Ohio, United States of America Center for Cardiovascular Research, Nationwide Children’s Hospital, Columbus, Ohio, United States of America Department of Pediatrics, The Ohio State University, Columbus, Ohio, United States of America Competence Network for Congenital Heart Defects, Berlin, Germany Division of Cardiovascular Sciences, School of Medical Sciences, Faculty of Biology, Medicine and Health, The University of Manchester, Manchester, United Kingdom Division of Evolution & Genomic Sciences, School of Biological Sciences, Faculty of Biology, Medicine and Health, The University of Manchester, Manchester, United Kingdom Experimental and Clinical Research Center (ECRC), a joint cooperation between the Charité Medical Faculty and the Max-Delbrück-Center for Molecular Medicine (MDC), Berlin, Germany Charité—Universitätsmedizin Berlin, corporate member of Freie Universität Berlin, Humboldt-Universität zu Berlin, and Berlin Institute of Health, Department of Pediatric Cardiology, Berlin, Germany DZHK (German Centre for Cardiovascular Research), partner site Berlin, Berlin, Germany The Centre for Applied Genomics, The Hospital for Sick Children, Toronto, Canada Genome Diagnostics, Department of Paediatric Laboratory Medicine, The Hospital for Sick Children, Toronto, Canada Department of Internal Medicine, McGovern Medical School, University of Texas Health Science Center at Houston, Houston, Texas, United States of America Michael E. DeBakey Department of Surgery, Baylor College of Medicine, Houston, Texas, United States of America Department of Surgery, Division of Cardiothoracic Surgery, Medical College of Wisconsin, Milwaukee, Wisconsin, United States of America Ruddy Canadian Cardiovascular Genetics Centre, University of Ottawa Heart Institute, Ottawa, Canada Institute of Human Genetics, University Hospital Ulm, Ulm, Germany Department of Human Genetics, University Medical Center Schleswig-Holstein (UKSH), Kiel, Germany Department of Congenital Heart Disease and Pediatric Cardiology, University Heart Center Freiburg—Bad Krozingen, Freiburg, Germany Department of Medical Biology, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands Department of Clinical Genetics, Amsterdam UMC, University of Amsterdam, Amsterdam, The Netherlands Cardiovascular Genetics, Department of Pediatrics, Centre Hospitalier Universitaire Saint-Justine Research Centre, Université de Montréal, Montreal, Canada Wellcome Sanger Institute, Wellcome Genome Campus, Hinxton, Cambridge, United Kingdom Laboratory of Translational Genetics, Department of Human Genetics, KU Leuven, Leuven, Belgium Firulli Anthony B. Editor Indiana University Purdue University at Indianapolis, UNITED STATES

I have read the journal’s policy and the authors of this manuscript have the following competing interests: The Department of Molecular and Human Genetics at Baylor College of Medicine receives revenue from clinical genetic testing conducted at Baylor Genetics Laboratories. M.E.H. is a co-founder of, consultant to and holds shares in Congenica, a genetics diagnostic company.

 * E-mail: larsal@sund.ku.dk (LAL); Marc-Phillip.Hitz@uksh.de (MPH) 29 7 2021 7 2021 17 7 e1009679 3 12 2020 23 6 2021 This is an open access article, free of all copyright, and may be freely reproduced, distributed, transmitted, modified, built upon, or otherwise used by anyone for any lawful purpose. The work is made available under the Creative Commons CC0 public domain dedication.

Numerous genetic studies have established a role for rare genomic variants in Congenital Heart Disease (CHD) at the copy number variation (CNV) and de novo variant (DNV) level. To identify novel haploinsufficient CHD disease genes, we performed an integrative analysis of CNVs and DNVs identified in probands with CHD including cases with sporadic thoracic aortic aneurysm. We assembled CNV data from 7,958 cases and 14,082 controls and performed a gene-wise analysis of the burden of rare genomic deletions in cases versus controls. In addition, we performed variation rate testing for DNVs identified in 2,489 parent-offspring trios. Our analysis revealed 21 genes which were significantly affected by rare CNVs and/or DNVs in probands. Fourteen of these genes have previously been associated with CHD while the remaining genes (FEZ1, MYO16, ARID1B, NALCN, WAC, KDM5B and WHSC1) have only been associated in small cases series or show new associations with CHD. In addition, a systems level analysis revealed affected protein-protein interaction networks involved in Notch signaling pathway, heart morphogenesis, DNA repair and cilia/centrosome function. Taken together, this approach highlights the importance of re-analyzing existing datasets to strengthen disease association and identify novel disease genes and pathways.

 Author summary

Congenital heart disease (CHD) is the most common congenital anomaly and represents a major global health burden. Multiple studies have identified a key genetic component contributing to the aetiology of CHD. However, despite the advances in the field of CHD within the last three decades, the genetic causes underlying CHD are still not fully understood. Herein we have assembled a large patient CHD cohort and performed a data-driven meta-analysis of genomic variants in CHD. This analysis has allowed us to strengthen the disease association of known CHD genes, as well as identifying novel haploinsufficient CHD candidate genes.

 german center for cardiovascular research (dzhk) https://orcid.org/0000-0002-9201-7840 Audain Enrique german center for cardiovascular research (dzhk) Hitz Marc-Phillip national register for congenital heart defects via Competence Network for Congenital Heart Defects https://orcid.org/0000-0002-9201-7840 Audain Enrique national register for congenital heart defects via Competence Network for Congenital Heart Defects Hitz Marc-Phillip kinderherz https://orcid.org/0000-0002-9201-7840 Audain Enrique kinderherz Hitz Marc-Phillip British Heart Foundation CH/13/2/30154 https://orcid.org/0000-0001-9573-0812 Keavney Bernard D. This study was supported by the German Center for Cardiovascular Research (DZHK) partner sites Berlin, Kiel; the Competence Network for Congenital Heart Defects, National Register for Congenital Heart Defects and KinderHerz (E.A. and M.P.H). BK is supported by a British Heart Foundation Personal Chair (CH/13/2/30154). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. PLOS Publication Stage vor-update-to-uncorrected-proof Publication Update 2021-08-10 Data Availability All relevant data are within the manuscript and its Supporting Information files. The data used in this study have been already published. The reference for each individual study is shown in S2 and S3 Tables. The assembled DNV dataset used in this study is provided in the S14 Table. In addition, we have provided a BED file with the CNV dataset (S15 Table). Introduction

Congenital Heart Disease (CHD) accounts for a large fraction of foetal and infant deaths, with incidence rates ranging from 7–9 per 1000 live births [1]. Within the last 30 years, survival rates have substantially increased due to improvements in surgical, interventional and clinical intensive care resulting in a rapidly growing number of CHD survivors reaching adulthood [2]. Nevertheless, there is still increased morbidity and mortality in individuals with CHD, resource utilization is high especially among severely affected patients, and importantly, the underlying etiology remains unclear for the majority of cases.

CHD is multifactorial, with both environmental and genetic risk factors [3,4]. Familial aggregation of CHD including Thoracic aortic aneurysm (TAA), as well as a large proportion of genomic copy number variants (CNVs) and de novo intragenic variations (DNVs) in probands with CHD suggest a strong genetic component. An estimated 4–20% of CHD cases are due to rare CNVs, suggesting that a significant part of CHD is caused by gene-dosage defects [5]. Recently, exome sequencing in large cohorts has been used to identify novel disease genes and strengthen known disease associations through the demonstration of an excess of de novo protein truncating variants (PTV) and rare inherited loss-of-function (LOF) variants in probands with CHD [6,7].

Overlaying both CNVs and PTVs has been used to define novel CHD relevant disease genes in contiguous gene disorders [8,9]. Following this principle, we have performed a genome-wide integrative meta-analysis of published and publicly available datasets of CNVs and DNVs identified in probands with CHD. This analysis, which is one of the larger meta-analyses of genomic variants in CHD so far, strengthens the disease association of known CHD genes and identifies novel haploinsufficient CHD candidate genes.

 Results Cohort description and workflow

We assembled a cohort with 7,958 cases (comprising both non-syndromic CHD, syndromic CHD and TAA cases) and 14,082 controls (S1 Table). Of the total of cases, 777 (~10%) were diagnosed with Thoracic Aortic Aneurysm (TAA). An overview of the sources used to assemble the present cohort is listed in S2 Table (for CHD cases) and S3 Table (for controls). We applied a set of quality control filters to our assembled CNV data before performing case-control association tests (Materials and Methods). In addition, common CNVs (minor allele frequency (MAF) in controls > 0.01) were excluded from the analysis. After filtering, 6,746 cases and 14,024 controls remained for further downstream analysis. Furthermore, we built a dataset of de novo variations (DNVs) identified in 2,489 probands with CHD from parent-offspring trios [6,7].

 CNV burden test of known CHD genes

Haploinsufficiency has been shown to cause a reasonable proportion of CHD [5]. Thus, genes known to be associated with CHD and genes which are intolerant for LOF variations should be deleted more often in probands with CHD than in controls. To test this hypothesis, we performed a CNV burden test using sets of genes known to be involved in CHD. In addition, we included genes known to be associated with developmental disorders, a curated list of known haploinsufficient disease genes, autosomal recessive disease genes and genes predicted to be intolerant to LOF variations (based on the observed/expected LOF ratio from gnomAD [10]). The burden test was performed using a logistic regression framework [11] (implemented in PLINK v1.7). Fig 1 and S4 Table summarize the results from the burden test on the different gene sets: known CHD genes (grouped in syndromic, non-syndromic, monoallelic and biallelic), developmental disorder genes, haploinsufficiency disease genes, autosomal recessive genes and all protein-coding genes. We tested all protein-coding genes to address the possibility that the analyses could be biased by differences in the CNV rate within the case and control groups, since we have assembled our cohort from different datasets. We did not observe genome-wide (all tested protein-coding genes) enrichment (P = 0.39, OR = 0.99) nor enrichment in the autosomal recessive gene set (P = 0.52, OR = 1.03) when comparing rare CNV deletions in cases vs controls. In contrast, the analysis revealed significant differences in the burden of CNV deletions between cases and controls for the set of haploinsufficiency genes (P = 8.29 x 10−13, OR = 2.27). As expected, our analysis revealed significant enrichment for the set of known CHD genes, which is mainly explained by the contribution of monoallelic CHD genes (P = 2.04 x 10−31, OR = 4.13) and syndromic CHD gene set (P = 1.66 x 10−33, OR = 4.06). Unlike the monoallelic and syndromic CHD gene sets, no significant enrichment was found for the nonsyndromic (P = 0.75, OR = 1.16) and biallelic (P = 0.08, OR = 1.87) CHD gene sets. Our analysis revealed a moderate enrichment of rare CNVs in the developmental disorder gene set (P = 6.90 x 10−11, OR = 1.75).

10.1371/journal.pgen.1009679.g001 CNV burden test on known gene sets.

The forest plot shows the odds ratio (dots), the 95% confidence intervals indicating the certainty about the OR (interrupted line) and the P-value in the indicated gene sets.

When the regression-based analysis was performed at different levels of the observed/expected LOF ratio (oeLOF) constraint metric (Fig 2 and S5 Table), we observed the higher enrichment toward the most LOF constrained genes (oeLOF < 0.01, P = 9.55 x 10−18, OR = 1.40) and still a moderate enrichment for genes with oeLOF < 0.1 (P = 0.002, OR = 1.09). No enrichment was observed in the set with oeLOF ratio > = 0.1 (P = 0.03, OR = 0.99). Based on these results we conclude that haploinsufficiency causes a significant component of CHD.

10.1371/journal.pgen.1009679.g002 CNV burden test on constraint LOF genes at different observed/expected LOF ratio thresholds.

The forest plot shows the odds ratio (dots), the 95% confidence intervals indicating the certainty about the OR (interrupted line) and the P-value in the indicated gene sets.

 Genome-wide identification of haploinsufficiency candidate disease genes for CHD

To perform a systematic, genome-wide identification of potential haploinsufficient CHD disease genes and loci, we analysed the CNV burden of 19,969 protein-coding genes (GENCODE v19). To this end, we compared the number of rare CNV deletions (MAF < 0.01) among cases and controls for each gene, and identified genes with significant CNV burden using a permutation test (significance level of adjusted P < 0.05, see Materials and Methods). If a CNV spanned two or more genes, all affected protein-coding genes were considered in the analysis. The distributions of rare CNV deletions in CHD cases across all 22 human autosomes is shown in Fig 3. Significant candidate genes had a median number of 12 overlapping CNVs in cases, compared to a median of 0 overlapping CNVs in controls (S1A Fig). Because CNVs can be large chromosomal aberrations, multiple genes were affected by some of the CNVs. In total, 528 genes (Sheet A in S6 Table) reached significance (Permutation test, P adjusted < 0.05). These 528 genes encompass a total of 63 loci (Sheet B in S6 Table, highlighted in magenta in Fig 3). The sizes of these loci range from 558 bp to 10.5 Mbp, with a median value of 243 Kbp (S1B Fig). The number of genes per locus ranged from 1 to 48, with a median value of 3 (S1C Fig). Only 16 loci contained a single gene (Sheet B in S6 Table).

10.1371/journal.pgen.1009679.g003 CNV deletion distribution across the 22 autosomes.

The plot shows the distribution of rare CNV deletions (green track) in CHD cases, the differences between the overlapping CNV deletions in cases and controls (black track) and highlight the location of the 63 significant loci discovered (in magenta).

In addition, we tested previously described CNV deletion syndrome regions (https://decipher.sanger.ac.uk/disorders#syndromes/overview) associated with developmental disorders and/or CHD for enrichment in our analysis (Materials and Methods). We found eight of these regions enriched in the dataset (S7 Table), with the 16p11.2-p12.2 locus being the region with the largest number of deletions in cases (n = 230).

 Shared genetic architecture of CHD and TAA

We independently performed a genome-wide test without the TAA cases to evaluate its impact on CHD. As expected, most of the genes (447 out of 528) remained significant after removing the contribution of TAA cases, since ~90% of the cases in the analyzed CNV cohort were CHD. Ten genes were significantly enriched independently when analyzing CHD and TAA cases, while 61 were significantly enriched only in TAA cases (S8 Table).

 De novo variation analysis

To identify an independent set of haploinsufficient CHD candidate genes, we combined de novo variations identified in two large-scale CHD case-control studies [6,7] and performed a gene-based de novo variation (DNV) burden test [12]. We analysed a total of 4,195 rare DNVs within 2,534 genes in the patient cohort. After classifying every variant into functional groups (Materials and Methods) 526 of these variants were predicted to be protein-truncating and 2,647 were missense. We evaluated for potential differences of the DNV rates between cohorts (see Materials and Methods). Comparison of the rate of each variant type across the groups was non-significant (P > 0.05, Poisson test, S9 Table).

We used two available statistical methods, Mupit [12] and DeNovoWEST [13], which test the significance of observed DNV at gene level, by comparing the number of observed variations with the number of expected variations (based on a sequence-dependent variation recurrence rate, see Materials and Methods). While Mupit focuses on enrichment of protein-truncating DNVs specifically, the DeNovoWEST test incorporates missense constraint information at variant level and applies a unified severity scale at variant level based on the empirically-estimated positive predictive value of being pathogenic. Based on the complementary results of both tests [13], we reported the minimal observed DNV p-value (Pdnv) per gene.

We identified 14 genes significantly enriched in the DNV analysis (P < 0.05 after Bonferroni correction for multiple testing, S10 Table). All of these genes were affected by at least two constrained non-synonymous DNV (nsDNV) and show significant overlap with 11/14 (78.6%) of the genes being known CHD disease genes. CHD7 (OMIM 214800) was the most significant haploinsufficient gene (P = 2.84 x 10−26) with 18 nsDNVs identified in the patient cohort. Other highly enriched genes for nsDNV—KMT2D (OMIM 147920), KMT2A (OMIM 605130), NSD1 (OMIM 117550), TAB2 (OMIM 614980), and ADNP (OMIM 615873)—have been previously associated with different types of neurodevelopmental disorders with co-occurrence of CHD. In the case of KDM5B (OMIM 618109), it has only been described in the context of a recessive neurodevelopmental phenotype with cases presenting ASD (Atrial septal defects) [14,15].

We next evaluated the distribution of o/e LOF ratio at different levels of DNV enrichment (genes were split based on Pdnv). Since the o/e ratio of LOF variation in each gene is strongly affected by its length, we instead used the 90% upper bound of its confidence interval (termed LOEUF), which keeps the direct estimate of the o/e ratio and allows to distinguish small genes from large genes, as suggested by Karczewski et al [10]. We observed that the genes with higher enrichment for nsDNV (lower Pdnv) show a significant decreased LOEUF compared to the mean of all protein-coding genes (Fig 4).

10.1371/journal.pgen.1009679.g004 Comparison of the distribution of LOEUF metric at different level of significance of nsDNV-enriched genes.

X-axis denotes the P-values from the DNV analysis (binned). Y-axis denotes the o/e LOF ratio upper bound fraction (LOEUF). All groups were compared against the LOEUF distribution of all protein-coding genes (purple). Differences between the distributions were tested using a two-sided Wilcoxon rank sum test. ****: P<0.0001, ns: non-significant.

 Integration of DNV and CNV results

To identify high confidence haploinsufficient CHD disease genes, we performed a joint analysis integrating the results from the CNV and the DNV analysis. We combined the results from both analyses (Pdnv and Pcnv) using the Fisher combine method. We demonstrated that both enriched genes for DNV and CNV deletions are significantly represented among LOF constraint genes (measured by the o/e LOF ratio). Therefore, we applied a Bonferroni multiple testing correction using independent hypothesis weighting [16] (IHW) by incorporating the gene o/e LOF ratio, as a measure of haploinsufficiency (S2 Fig). Our analysis revealed 21 genes that were significantly enriched for CNV deletions and/or non-synonymous DNV (Table 1). A gene was included in the final set of haploinsufficient CHD disease genes if it reached a significant corrected metaP < 0.05 (after Bonferroni adjustment with IHW).

10.1371/journal.pgen.1009679.t001 Top 21 significant genes arising from both the permutation-based test and the DNV rate-based test.

Cases/Controls: Number of cases and controls carrying CNV deletions overlapping the gene in the CNV analysis. Pcnv: p-value from the CNV permutation test. nsDNV: Number of constrained non-synonymous variations in the de novo analysis. Pdnv: p-value from the DNV analysis. Significant: The analysis where the gene was significant (dnv: DNV analysis, cnv: CNV analysis, both: Both analysis, none: Non-significant neither DNV nor CNV analysis). metaP: combined p-value (Pdnv and Pcnv) using the Fisher method. Pihw: Bonferroni corrected p-value using independent hypothesis weighting (IHW) and LOEUF metric as covariate. LOEUF: o/e LOF ratio upper bound fraction from gnomAD. *All the 21 genes were significant after combining their p-values and applying Bonferroni correction. 1Evidence is from mouse models [24,62].

CNV DNV Combined
Gene case control Pcnv nsDNV Pdnv *Significant metaP Pihw LOEUF Known CHD
CHD7 6 1 6.80E-03 18 2.84E-26 dnv 1.25E-26 8.05E-23 0.076 Yes
KMT2D 0 0 1.00E+00 18 1.32E-25 dnv 7.67E-24 4.93E-20 0.103 Yes
NSD1 1 1 5.63E-01 12 1.00E-14 dnv 1.90E-13 2.14E-09 0.095 Yes
KMT2A 0 0 1.00E+00 7 1.00E-14 dnv 3.32E-13 1.86E-09 0.065 Yes
NOTCH1 10 24 1.00E+00 7 1.00E-14 dnv 3.32E-13 2.14E-09 0.097 Yes
TAB2 12 0 1.00E-04 5 3.46E-09 both 1.03E-11 5.75E-08 0.098 Yes
ANKRD11 13 0 1.00E-04 3 2.32E-05 cnv 4.85E-08 2.72E-04 0.107 Yes
WHSC1 11 0 1.00E-04 3 8.73E-05 cnv 1.71E-07 9.96E-04 0.119 No1
ADNP 0 0 1.00E+00 4 9.94E-09 dnv 1.93E-07 1.13E-03 0.123 Yes
DYRK1A 4 0 1.43E-02 4 9.46E-07 dnv 2.59E-07 1.64E-03 0.214 Yes
NALCN 10 1 1.00E-04 3 1.76E-04 cnv 3.32E-07 6.83E-03 0.522 No
ELN 30 0 1.00E-04 2 1.77E-04 cnv 3.34E-07 7.50E-03 0.871 Yes
WAC 7 0 4.00E-04 3 1.31E-04 none 9.33E-07 5.44E-03 0.084 No
RBFOX2 1 0 3.45E-01 4 1.59E-07 dnv 9.72E-07 6.25E-03 0.194 Yes
KANSL1 94 110 2.00E-04 2 3.38E-04 none 1.19E-06 6.92E-03 0.238 Yes
MYO16 13 2 1.00E-04 2 8.38E-04 cnv 1.45E-06 1.74E-02 0.272 No
MED13L 2 1 2.70E-01 4 4.58E-07 dnv 2.09E-06 1.22E-02 0.064 Yes
KDM5B 0 0 1.00E+00 4 1.45E-07 dnv 2.43E-06 4.97E-02 0.572 No
GATA6 0 0 1.00E+00 5 1.80E-07 dnv 2.98E-06 1.92E-02 0.174 Yes
ARID1B 4 0 1.31E-02 4 1.48E-05 none 3.20E-06 1.87E-02 0.102 No
FEZ1 10 0 1.00E-04 2 2.22E-03 cnv 3.62E-06 4.30E-02 0.414 No
Subclassification of CHD phenotypes

We performed a further analysis based on specific CHD subtypes in addition to the collective analysis of all CHD phenotypes. The analysis focused on simplex CHD cases, within two main categories: Conotruncal (NCNV = 873, NDNV = 234) and LVOTO (NCNV = 594, NDNV = 351). We only included cases with a clear phenotypic description and without any overlapping phenotypic features between the two categories (LVOTO/Conotruncal). The conotruncal group consisted mostly of Tetralogy of Fallot (TOF), Truncus arteriosus and Transposition of great arteries (TGA), whereas LVOTO mainly constituted Aortic stenosis (AS) including Bicuspid aortic valve disease (BAV), Coarctation and Hypoplastic left heart syndrome (HLHS).

As described above, we performed the same integrative approach for the LVOTO and conotruncal groups to identify CHD subtype-specific genes. Our analysis showed four significant genes (S11 Table). Three are observed in LVOTO (KMT2D, KMT2A and TAB2) and a single gene showed significant enrichment in the conotruncal subtype (NSD1).

 Significant CHD genes are highly and/or differentially expressed in the heart

We next evaluated the expression pattern of the 21 significant genes (Bonferroni corrected metaP < 0.05) in the heart using RNA-Seq data from human tissues at different developmental time points [17]. We stratified the analysis based on stages of heart development (see Materials and Methods). Our analysis revealed that the most significant genes (metaP < 1 x 10e-5) show significantly increased mean expression in the heart (P < 0.0001, Wilcox test) at different developmental stages (development, maturation and infant/adult), compared to all protein-coding genes (Fig 5). Moreover, 18 out of 21 genes fall in the in the top quartile of heart expression in both developmental and maturation stages (S3 Fig). To complement our expression analysis, we compared gene expression during human heart development with expression in two other mesodermal organs: kidney and liver. This allowed identification of genes with significant changes in its expression levels during crucial heart developmental stages, which would have not been possible when focusing on expression levels alone (Materials and Methods). We found that 17 out of 21 CHD candidate genes are differentially expressed in the heart (R2 > 0.50, Bonferroni corrected P < 0.01) when compared to its expression levels in kidney and/or liver. Interestingly, the three genes (FEZ1, NALCN and MYO16) which are not among the highly expressed genes, were found to be significantly differentially expressed during heart development compared to kidney and/or liver (S12 Table).

10.1371/journal.pgen.1009679.g005 Comparison of the mean expression (heart) distribution at different <italic>metaP</italic> cut-offs.

Panels show three different heart development stages: early development, maturation and infant/adult. X-axis denotes the combined p-value from DNV and CNV analysis (metaP, at different cut-offs). Y-axis denotes the genes’ mean expression in the heart (log scale). The 21 significant candidate CHD genes (Table 1) are contained in the fraction with the higher expression (red box). Differences between the distributions were tested using a two-sided Wilcoxon rank sum test (reference group: all genes). ****: P<0.0001, ***: P<0.001, **: P<0.01, *: P<0.05, ns: non-significant.

 CNV/DNVs burden of specific protein complexes

CNVs and DNVs can affect heart development either through haploinsufficiency of a single gene, or through its combined impact on the function of several genes. Indeed, oligogenic models have been implicated in CHD, and proteins acting in the same complex or pathway are known to be encoded in genomic clusters [18,19]. We therefore conducted a systems-level analysis to identify global mechanisms by which haploinsufficiency might promote CHD. In particular, we assessed the combined effect of CNVs and DNVs with respect to human protein-protein interactions (PPIs). The InWeb and ConsesusPathDB databases provides ranked information about experimentally determined physical interactions and, therefore, serves as a proxy to understand the functional effects of CNV/DNVs on human protein complexes (Materials and Methods). The genes with Benjamini–Hochberg adjusted metaP < 0.05 (n = 492 genes) were used as seeds to build a PPI network from the data available in InWEb and ConsessusPathDB. No additional interections were considered. The final network consisted of 164 proteins and 290 interactions (S4 Fig). A total of 10 overlapping sub-clusters within this network were identified using the in-built clustering algorithm implemented in GeNets [20] (Materials and Methods). Gene-ontology (GO) enrichment analysis suggested that four out of these ten sub-clusters are enriched for genes involved in Notch signaling pathway, cardiocyte differentiation, DNA repair and centrosome function (Fig 6). All the four clusters accommodate more CNV deletions in CHD cases compared to controls. Six out of the ten sub-clusters did not show significant enrichment for any particular biological process.

10.1371/journal.pgen.1009679.g006 Identification of functional networks enriched for proteins encoded by genes affected by CNVs and/or DNVs associated with CHD.

The protein-protein interaction networks (a-d, for clusters 1, 3, 8 and 9 respectively) were identified using GeNets (S4 Fig). Proteins are shown as nodes, interactions as edges. Enrichment for CNVs (blue) and DNVs (green) are highlighted. Proteins with no specific enrichment for CNV and/or DNVs but with B-H adjusted metaP < 0.05 are highlighted in red. The size of the circles denotes if the genes was found significantly highly and/or differentially expressed in the heart (large circles: significant expression; small circles: non-significant). The distribution of CHD case-CNVs and control-CNVs are shown for each cluster. Significant difference in the CNV distribution was calculated using a Wilcox rank sum test. The horizontal bar plots show the top ten GO enriched terms for each cluster (output from Enrichr tool). X-axis in the horizontal bar plot denotes the combined score from Enrichr, which is computed by multiplying the log-transformed p-value and the z-score. Bar color encoded the GO biological process significant level (dark blue: FDR < 5%, light blue: FDR 5–10%, grey: FDR > 10%).

 Discussion

We performed a meta-analysis of rare genomic variants in a cohort of 10,447 CHD probands, which provides a useful resource for interpreting CNVs and DNVs identified in patients with CHD. We implemented a statistical approach which allows the integration of different types of genomic variants to discover novel genes associated with CHD. Our data-driven integrative analysis took into account three major criteria at the genomic level: a) gene enrichment for DNVs, b) gene enrichment for CNV deletions and c) gene intolerance for LOF variations. Our analysis identified 21 significant haploinsufficient CHD genes. Fourteen of these are known CHD genes, and the remaining seven genes have not previously been associated with CHD (Table 1).

To further strengthen associations, we made use of a newly published human transcriptome atlas covering different developmental, maturation and adult stages in numerous organs [17]. Similar to previous results [7], our analysis highlights that the majority of the 21 significant genes are highly expressed during critical stages of heart development. Unlike earlier studies [7,21] which did not address the importance of expression changes over time, we evaluated the differential expression patterns of genes by comparing levels of expression in the heart, kidney and liver at different time points in development. This analysis allowed us to strengthen disease association for genes not falling under the high expression group and highlight the critical importance of all 21 genes independently of the genomic approach. This aspect is complemented by the fact that the majority of genes (14/21) were already known to cause CHD. To further strengthen disease associations, spatiotemporal expression at single-cell resolution during critical cardiac developmental timepoints and analyses of animal models with targeted mutation in the candidate disease genes is warranted. This could strengthen disease association further and provide pathophysiological information.

Among the 21 likely haploinsufficient disease genes for which the combined analyses showed enrichment (Bonferroni corrected metaP < 0.05), 14 genes (CHD7, KMT2D, KMT2A, NOTCH1, NSD1, TAB2, ANKRD11, ADNP, DYRK1A, RBFOX2, KANSL1, ELN, MED13L and GATA6) are well-established CHD genes, and our data confirms this association. To the best of our knowledge, association between CHD and seven genes (KDM5B, WHSC1, WAC, NALCN, ARID1B, FEZ1 and MYO16) had either not been established, or had been reported in small cases studies or a single individual only.

KDM5B is not an established CHD gene thus far, although one patient with compound heterozygous frameshift variants had an ASD [14]. While some have argued against the haploinsufficiency of the gene [22], our analysis suggests KDM5B as a plausible haploinsufficient CHD gene. Additional functional studies are warranted to confirm its role in CHD.

A recent CNV meta-analysis [23] based on non-syndromic CHD patients found that duplication of WHSC1 (also known as NSD2) is a possible cause of CHD. However, haploinsufficiency of WHSC1 has not previously been associated with CHD. In support of its role in CHD, Whsc1 has been reported to cause heart malformations in mouse models [24]. In addition, WHSC1 is known to interact with NKX2.5 [24]. In spite of this, the low incidence of CHD in individuals with Wolf-Hirschhorn syndrome suggests that haploinsufficiency of WHSC1 alone does not cause CHD.

Heterozygous truncating variations in WAC, as well as CNV deletions involving this gene, have been recently associated with the DeSanto-Shinawi syndrome, a rare neurodevelopmental disorder characterized by global developmental delay [25,26]. Furthermore, in two non-consanguineous unrelated individuals with heart malformations, among other disorders [27], microdeletions at 10p11.23-p12.1 (overlapping ARMC4, MPP7, BAMBI and WAC) were identified. Despite these isolated reports, no definite association between WAC and CHD has been established.

DNVs in NALCN have been reported to cause a dominant condition characterized by multiple features including developmental delay, congenital contractures of the limbs and face and hypotonia [28,29]. However, among the phenotypes observed, CHD have been not described thus far.

Heterozygous variation of ARID1B is a frequent cause of intellectual disability [30,31]. A recent analysis of 143 patients with ARID1B variations showed that individuals display a spectrum of clinical characteristics. Congenital heart defects were observed in 19.5% of the patients [32].

FEZ1 is a neurodevelopmental gene, which has been associated with schizophrenia [33]. Fez1 has been reported to be regulated by Nkx2-5 in heart progenitors in mice, suggesting a possible role in heart development [34].

MYO16 (NYAP3) encodes an unconventional myosin protein, involved in regulation of neuronal morphogenesis [35]. We have not found an association between MYO16 and heart development in the literature.

Although, several genes have been shown to be altered in syndromic and non-syndromic cases with CHD and TAA (e.g. HEY2 [36], MYH11 and NOTCH1 [37]), among the 10 genes significant in our analysis for TAA, CHD and the combined scenario, none has been reported previously to be associated with either CHD or TAA. Given the limited data size and only accessing CNV calls from TAA cases, future studies looking at CNV and DNV in both phenotypes are required to establish stronger genotype-phenotype correlation to better understand a possibly shared genetic architecture for the two disease entities.

Also, our study did not identify strong signals associated with specific CHD subtypes. We identified four genes with significant enrichment (adjusted metaP < 0.05) within the two evaluated CHD subtypes (LVOTO and conotruncal defects). Three were associated with LVOTO, and the contribution was mainly from DNV. KMT2D was enriched in LVOTO, which is consistent with the reported spectrum of CHD in patients with Kabuki syndrome, where a large proportion of individuals have LVOTO type CHD [38,39]. PDA and septal defects predominate in Wiedemann-Steiner Syndrome (KMT2A). However, aortic insufficiency and BAV have all been reported [40] suggesting that LVOTO might form part of the phenotypic spectrum. TAB2 was also enriched in LVOTO. Mutations in TAB2 are associated with a wide range of cardiac phenotypes [41–43]. Deletions at 6q24 causing haploinsufficiency of TAB2 have been associated with outflow tract abnormalities, including LVOTO [44,45], which might point to a role of TAB2 LVOTO pathogenesis. NSD1 was the only gene identified as significant among conotruncal defects. The significance of this is unclear as septal defects predominate in Sotos Syndrome. However, patients in previous studies were ascertained focusing on Sotos syndrome (OMIM 117550) rather than CHD. Given the current size of each CHD subgroup and the low number of CNV and DNV events in these genes, these results cannot be considered conclusive. More precise phenotypic descriptions for each individual are necessary to increase the genotype-phenotype correlation for individual genes and improve the identification of genetic subnetworks, along with larger sample sizes.

In addition to the gene-centered analysis, we also applied a systems-level analysis in order to identify potential novel pathophysiological mechanisms affected by haploinsufficiency. In this approach, we took advantage of GeNets [20], a computational framework for the analysis of protein-protein interactions, developed for the interpretation of genomic data. Our analysis allowed us to identify PPI clusters enriched for genes affected by CNVs and/or DNV in patients. Furthermore, GO enrichment analysis suggested distinct biological functions for four of these clusters.

Cluster 1 (Fig 6A) contains proteins involved in the Notch signaling pathway. Our data corroborate previous studies that confirm the central role of Notch pathway in the pathophysiology of CHD [46] and highlights the shared contribution of CNVs and DNVs within the cluster. Cluster 3 (Fig 6B) contains proteins driving essentials processes in the development of the heart such as atrial septum and cardiac right ventricle morphogenesis as well as proteins playing significant role in the positive regulation of gene expression. These mechanisms has been well studied elsewhere [47]. Interestingly, three out of the seven candidate novel CHD genes (WAC, ARID1B and KDMB5) were found to be contributing to these two clusters. Cluster 8 (Fig 6C) showed enrichment for processes related with chromosome organization and DNA repair. Association between DNA repair and CHD is not well established thus far. Cluster 9 (Fig 6D) was found to be associated with microtubule organizing function. This biological process has been not described in the context of CHD, although an earlier report [48] describes complex CHD among the phenotypes in individuals with 15q11.2 deletion syndrome, which involves the tubulin gamma complex protein 5.

Given the heterogenous data sources and the complex inheritance patterns often observed in patients with CHD, our study has limitations. Firstly, the patient data was collected from almost 200 different published sources, and in many cases it was only possible to obtain data from CNV calls which had already been suggested to be pathogenic. Thus, we are aware that our patient data are incomplete because genome-wide CNV data are missing from a large part of the patient cohort and re-emphasizes already established associations. This is not the case for controls, for which genome-wide data was used. As a direct consequence, even though the difference between the rates of CNV deletions in controls and cases decreased dramatically after applying a quality control filtering step, a slight difference remained between both cohorts. The lack of collected CNV data spanning sex chromosomes limited the analysis only to autosomal chromosomes. In addition, the distribution of CNVs that overlap known microdeletion syndromes such as DiGeorge syndrome and Williams syndrome is overrepresented in the dataset. Similarly, the degree of phenotyping varied across the different studies, and often only basic phenotypic terms relating to CHD were available. This made it impossible to refine the diagnosis to a precise phenotypic class of CHD in many individuals.

Previous research has shown that the chances of finding a genetic cause of CHD is higher in syndromic, rather than non-syndromic CHD [6,7]. Therefore, it is not surprising that known CHD genes were enriched in this cohort. To identify non-syndromic causes of CHD, it is important to take into account previous findings, which have shown a significant excess of apparently deleterious inherited PTVs in unaffected parents [6]. To help address the challenges of identifying non-syndromic CHD genes, we have used an integrative approach, which has allowed us to look for novel CHD associations in a larger sample size, in a binary fashion. This will help to facilitate future studies.

Variable expression and reduced penetrance are common features in CHD, including in even well-established conditions such as Noonan Syndrome [49]. Most cases of CHD occur as a one off in the family and when recurrences do happen, CHD subtypes are more likely to be discordant than concordant [50]. This suggests that other modifying factors may be at play, including genetic, or environmental factors, or both. This presents a further challenge in identifying genetic causes of CHD. Moving forward, detailed phenotypic descriptions using a standardized system, and better data sharing strategies (including primary data), will facilitate further gene discovery and improved genotype-phenotype correlation in CHD subgroups.

In summary, we have performed an integrative analysis of CNVs, DNVs, o/e LOF ratio and expression during heart development amongst more than 10,000 CHD patients. Our analyses identify seven potential disease genes and mechanisms with novel association with CHD and strengthen previously reported associations.

 Materials and methods Ethics statement

This project was based on unidentifiable data and did not require approval by science ethics committees in Denmark or Christian-Albrechts-Universität in Kiel.

 Cohort description

Our cohort contains 7,958 CHD cases and 14,082 controls (see summary at S1 Table). Data from both affected and unaffected individuals were collected from 190 different CNV studies (S2 Table). Most of the CNV data included in the present study were assembled from public repositories, data available from literature as well as clinical data (see S2 and S3 Tables for a more detailed description). We sampled all available and accessible studies as of February 2018 cited in PubMed. We have focused on studies incorporating Caucasian samples (the larger sampled population) to decrease population effects. Given no primary data were available for most of the studies, we cannot exclude the possibility that non-Caucasian samples were included. Both non-syndromic and syndromic CHD cases were included and CNVs were mapped to the human genome build NCBI37/hg19. Phenotype information was reviewed and if possible standardised across studies, to ensure consistency and accuracy. We used as controls re-analysed samples from the Wellcome Trust Case Control Consortium 2, the Genetic Association Network (GAIN) and the Ottawa Heart Institute (S3 Table). In addition, we built a dataset from the two largest DNV studies in CHD published thus far, which include a total of 2,489 parent-offspring trios [6,7].

 Determining the gene sets

The compiled highly confident CHD gene list consisted of genes with plausible disease-causing mutations in an interpretable functional region of a gene reported in association with CHD in three or more unrelated individuals. Genes, where fewer than three reports of CHD exist, were also considered if there was available solid functional evidence, such as a mouse model that displays CHD or in silico evidence. CHD genes were then coded as either non-syndromic (isolated CHD), or syndromic based on published phenotypes. The list of genes associated with developmental disorders was derived from the Developmental Disorder Genotype to Phenotype (DDG2P) list maintained by Decipher and the European bioinformatics Institute [51]. All genes were annotated with either monoallelic or biallelic as appropriate, based on the published literature (S13 Table).

 CNV analysis

Only autosomal CNVs were included in the analysis. All the CNV boundaries were determined using genome build NCBI37/hg19. For the CNVs provided in hg18, we used the Assembly Converter (https://www.ensembl.org/Homo_sapiens/Tools/AssemblyConverter) build on CrossMap (http://crossmap.sourceforge.net/) to convert samples to NCBI37/hg19. Also, an extra validation step of all CNV boundaries was performed using the R-package BSgenome.Hsapiens.UCSC.hg19. Smaller and longer CNVs were filtered out by applying a size cut-off of 5 Kb and 20 Mb as lower and upper limit, respectively. It has been demonstrated before that smaller and larger CNV calls tend to have a high rate of false positives [52,53]. We removed CNVs overlapping more than 50% of telomeres, centromeres and segmental duplication regions. In addition, we computed the internal CNV frequencies by counting the number of relative overlaps (>50% reciprocal overlap) on the CNVs control subset divided by the total number of controls. The internal MAF was computed for deletions and duplications subsets separately. Only CNVs with a minor allele frequency (MAF) < 0.01 in controls and overlapping ten or more CNV platform calling probes (Affymetrix Array 6.0 and Illumina Human660W-Quad) were considered for downstream analysis. Our analysis was focused only on CNV deletions. The distributions of the number of CNV deletions per individual within the case and control groups were compared (two-sided Wilcoxon rank sum test) to evaluate the impact of the quality control filtering step (S5 Fig). After filtering, 6,746 cases (3, 929 harbouring CNV deletions) and 14,024 controls (12,585 harbouring CNV deletions) remained for further analysis. A region-based permutation test (using PLINK version 1.07, test ‘—cnv-test-region—mperm 10000’) was used on the filtered set to perform a case-control association analysis. For the gene-based permutation analysis, we reported both the ‘point-wise’ empirical p-value (EMP1) and the empirical adjusted p-value (EMP2), which controls the family-wise error rate (FWER) (http://zzz.bwh.harvard.edu/plink/perm.shtml). In addition to the gene-centered permutation testing, a similar region-based permutation analysis was performed to access enrichment in known CNV deletion syndromes. All CNVs deletions passing QC filtering overlapping these regions were considered in the analysis. The region genomic coordinates and syndrome descriptions were downloaded from the Database of genomic variation and phenotype in humans using Ensembl resources (Decipher, https://decipher.sanger.ac.uk/disorders#syndromes/overview).

 CNV burden test on gene sets

A logistic regression-based burden test (‘cnv-enrichment-test’ in PLINK v1.7) [11] was performed on different gene sets (known CHD genes (non-syndromic/syndromic/biallelic/monoallelic), developmental disorder genes, known haploinsufficient disease genes [54], autosomal recessive disease genes [55,56] and low observed/expected LOF ratio genes, S13 Table) using the rare CNV deletions passing the quality control and filtering stage. For every gene set examined, the binary phenotype (CHD case or control) was regressed on the number of genes disrupted by one or more CNVs. The averaged CNV size and the number of segments per individual were used as covariates into the model to control for potential differences between cases and controls as suggested by Raychaudhuri et al [11]. In addition, the PLINK implementation of this test was slightly modified by including a third (categorical) covariate, the sample study ID, since we have assembled the CNV data from different sources.

 DNV analysis

The assembled DNV dataset (Sheet A in S14 Table) was re-annotated using the Variant Effect Predictor (VEP version 90) tool. All the DNVs included in this study were validated with the VariantValidator tool [57] (Sheet B in S14 Table). Based on the VEP annotation, we classified every variation into three major functional groups as follows: a) Protein truncation variant (stop_gained, splice_acceptor, splice_donor, frameshift, initiator_codon, start_lost, conserved_exon_terminus), b) missense variant (stop_lost, missense, inframe_deletion, inframe_insertion, coding_sequence, protein_altering) and c) silent variant (synonymous). Variants with minor allelic frequency (MAF) > 0.01 in gnomAD database were excluded from the analysis. The rates of rare DNVs (MAF < 0.01) in both DNV studies [6,7] were compared (Poisson test) for different variant consequence groups (PTV, missense and synonymous). No significant differences were found between the DNV rates for any of the evaluated groups (S9 Table). De novo variation recurrence significance testing was performed to evaluate the impact of DNVs at gene level using the Mupit tool [12]. By default, Mupit uses the sequence-specific variation rate published by Samocha et al [58]. A second test, DeNovoWEST [13], was used to assess gene-wise de novo variation enrichment. DeNovoWEST assigns a variation severity score (based on the variant consequences and the CADD score) to all classes of variants as a proxy of its deleteriousness. For each tested gene, the minimal p-value obtained from Mupit and DeNovoWEST was reported (Pdnv). The corrected P value was computed using the Bonferroni method with n = 18,272.

 Inferring differentially and highly expressed genes

Differentially expressed genes (DEGs) were identified by comparing the gene expression profile in heart to kidney and liver at matched time points. We used maSigPro R-package [59] for inferring genes with dynamic temporal profiles from time-course transcriptomic data as previously described by Cardoso-Moreira et al [17]. As the input for maSigPro, we used the count per million matrix (CPM, output from EdgeR package) hosted in ArrayExpress (E-MTAB-6814). Genes which did not reach a CPM > 0.5 in at least five samples were excluded from the analysis. We ran maSigPro on the time scale measured in days post-conception using defaults parameters and only included time points with at least two biological replicates. A gene was selected as DEG if the R2 (goodness-of-fit) parameter was higher than 0.50 and Bonferroni corrected P < 0.01. The R2 parameter distinguish genes with clear expression trends from genes with ‘flat’ expression profile. S12 Table lists the final DEGs identified in the heart (R2 > 0.50). To assess the gene expression levels in the heart, the RPKM matrix was used. Gene expression levels were averaged among samples in the different development stages of the heart as follow: early development (4wpc-8wpc), maturation (9wpc-20wpc), infant/adult (newborn-adulthood). Genes were ranked based on the computed mean expression.

 Identification of CNV/DNV enriched PPI sub-clusters

A protein-protein interaction network was constructed using the GeNets framework [20] and the information from InWeb [60] and ConsensusPathDB [61] protein-protein interaction databases. Nodes in the network correspond to proteins whereas edges represent their physical interactions. The network was strictly seeded with 492 candidate genes, those with a significant adjusted metaP < 0.05 (Benjamini-Hochberg’s false discovery rate, FDR). The PPI network was partitioned into overlapping sub-clusters using the in-built clustering method described in GeNets [20]. Only statistically significant sub-clusters (p-value < 0.05, permutation test) with at least 5 proteins were considered for further analysis. Finally, Gene Ontology enrichment analysis (Biological Process database 2018) of each identified sub-cluster was performed using the enrichr tool (https://maayanlab.cloud/Enrichr/).

 Supporting information Distribution of overlapping CNVs, size and genes in 63 CHD loci.

A) CNVs per gene. Overlapping CNVs for each of the 528 significant candidate genes are shown as box-and-whiskers plots. Statistically significant difference was observed between the two distributions (Mann-Whitney test, ***: P<0.001). B) Size of loci in kilobase-pairs (kbp). C) Number of genes per locus. Median values are shown above each box.

(TIF)

 Statistical framework to discover novel candidate CHD genes by integrating DNV and CNV deletions.

The workflow follows four major steps: 1) Data aggregation and quality control of both DNV data and CNV data, 2) DNV rate-based enrichment testing and CNV deletions case/control association analysis at gene level are performed independently, 3) the results are combined using the Fisher method and 4) P-values are Bonferroni corrected using the Independent Hypothesis Weighting method (IHW). As independent covariate for the IHW method, the o/e LOF ratio upper bound fraction (LOEUF) was used.

(TIF)

 Heart expression pattern of the 21 significant genes at different heart development stages.

Panels show three different heart development stages: early development (red), maturation (green) and infant/adult (blue). The x-axis denotes the percentile rank of heart expression in the heart. The y-axis denotes the o/e LOF ratio upper bound fraction (LOEUF) from gnomAD. Dashed lines denote the threshold for highly expressed genes (expression rank > = 0.75) and highly LOF constrained genes (LOEUF < = 0.30).

(TIF)

 The functional network enriched for proteins encoded by genes affected by CNVs and/or DNVs associated with CHD.

Ten sub-clusters were identified using GeNets. Proteins are shown as nodes, interactions as edges. Enrichment for CNVs (blue), DNVs (green) or both independently (purple) are highlighted. Proteins with no specific enrichment for CNV and/or DNVs but with B-H adjusted metaP < 0.05 are highlighted in red. The size of the circles denotes if the gene was found significantly highly and/or differentially expressed in the heart (large circles: significant expression; small circles: non-significant).

(TIF)

 <p>Distribution of the number of CNV deletions per individual in both control and CHD case cohorts before (A) and after (B) applying the quality control filtering approach. Differences between the distributions were tested using a two-sided Wilcoxon rank sum test. ****: P<0.0001.</p> <p>(TIF)</p> </caption> </supplementary-material> <supplementary-material id="pgen.1009679.s006" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" position="float" xlink:href="info:doi/10.1371/journal.pgen.1009679.s006" xlink:type="simple"> <label>S1 Table</label> <caption> <title>Number of probands in the CNV cohort.

Stratified by CHD cases/controls and CNV type (deletion/duplication).

(XLSX)

 Sources of the CNV-case cohorts used in this study.

(XLSX)

 Sources of the CNV-control cohorts used in this study.

(XLSX)

 Gene set-based logistical regression enrichment CNV analysis.

(XLSX)

 Gene set-based logistical regression enrichment CNV analysis stratified by observed/expected LOF ratio (from gnomAD).

(XLSX)

 <p>Sheet A) Gene-based case/control CNV-deletions permutation testing (PLINK results). Sheet B) Significant locus (Locus ID and contributing genes).</p> <p>(XLSX)</p> </caption> </supplementary-material> <supplementary-material id="pgen.1009679.s012" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" position="float" xlink:href="info:doi/10.1371/journal.pgen.1009679.s012" xlink:type="simple"> <label>S7 Table</label> <caption> <title>Case/control CNV permutation testing on known deletion syndrome (PLINK results).

(XLSX)

 CNV case/control permuatation testing output from PLINK.

The table shows the 528 significant genes combining both CHD and TAA cases (CHD+TAA), the contribution of only CHD cases (only CHD) and the contribution of only TAA cases (only TAA).

(XLSX)

 Comparision of the DNV rates, stratified by variant consequences, between two independent cohorts.

(XLSX)

 Gene-based DNV analysis.

(XLSX)

 Meta-analysis of CNV/DNV stratified by CHD sub-types (Conotruncal and LVOTO).

Table shows the four genes with Bonferroni corrected metaP < 0.05. Cases/Controls: Number of cases and controls carrying CNV deletions overlapping the gene in the CNV analysis. p_cnv: p-value from the CNV permutation test. nsDNV: Number of constrained non-synonymous variations in the de novo analysis. p_dnv: p-value from the DNV analysis. metaP: combined p-value (Pdnv and Pcnv) using the Fisher method. adj metaP: Bonferroni corrected p-value using independent hypothesis weighting (IHW) and LOEUF metric from gnomad as covariate.

(XLSX)

 Differentially expressed genes in the heart compared to kidney and liver.

(XLSX)

 List of gene sets used in the CNV enrichment analysis.

(XLSX)

 <p>Sheet A) List of the <italic>de novo</italic> variants analized in this study. Sheet B) Validation results of the DNVs from the VarinatValidator tool.</p> <p>(XLSX)</p> </caption> </supplementary-material> <supplementary-material id="pgen.1009679.s020" mimetype="application/vnd.openxmlformats-officedocument.spreadsheetml.sheet" position="float" xlink:href="info:doi/10.1371/journal.pgen.1009679.s020" xlink:type="simple"> <label>S15 Table</label> <caption> <title>Rare CNV deletions (MAF 0.01) used in this study.

CNVs bounderies are determinated in genome build hg19.

(XLSX)

We express our gratitude to the patients and their families for their participation in the analysed studies. We would like to thank the Genetic Association Information Network (GAIN) and dbGAP for making the data available. We would like to thank the Wellcome Trust Case Control Consortium (WTCCC) for making the data accessible https://www.wtccc.org.uk/info/access_to_data_samples.html. We used data from the the Deciphering Developmental Disorders (DDD) study. The DDD study presents independent research commissioned by the Health Innovation Challenge Fund, a parallel funding partnership between the Wellcome Trust and the UK Department of Health, and the Wellcome Trust Sanger Institute. The views expressed in this publication are those of the author(s) and not necessarily those of the Wellcome Trust or the UK Department of Health. We would like to thank the Pediatric Cardiac Genomics Consortium (PCGC) and dbGAP for making the data publicly available. This study makes use of data generated by the DECIPHER community. A full list of centres who contributed to the generation of the data is available from http://decipher.sanger.ac.uk and via email from decipher@sanger.ac.uk. We thanks to Margarida C. Moreira and the Kaessmann Lab by making accessible the human RNA-Seq data and the support for the data analysis. We thanks Rasmus Wernersson and Federico de Masi for facilitating the use of the protein-protein interaction database, InWeb, in this work. We thanks the Lage Lab and Taibo Li for their support with GeNets. And finally, we would like to thank all data submitters and collaborators for their contributions.

 References van der Linde D, Konings EEM, Slager MA, Witsenburg M, Helbing WA, Takkenberg JJM, et al. Birth Prevalence of Congenital Heart Disease Worldwide. J Am Coll Cardiol. 2011;58: 2241–2247. doi: 10.1016/j.jacc.2011.08.025 22078432 van der Bom T, Zomer AC, Zwinderman AH, Meijboom FJ, Bouma BJ, Mulder BJM. The changing epidemiology of congenital heart disease. Nat Rev Cardiol. 2011;8: 50–60. doi: 10.1038/nrcardio.2010.166 21045784 Warnes CA, Botto L, Correa A, Britt AE, Elixson M, Jenkins KJ, et al. Noninherited Risk Factors and Congenital Cardiovascular Defects: Current Knowledge. Circulation. 2007;115: 2995–3014. doi: 10.1161/CIRCULATIONAHA.106.183216 17519397 Zaidi S, Brueckner M. Genetics and Genomics of Congenital Heart Disease. Circ Res. 2017;120: 923–940. doi: 10.1161/CIRCRESAHA.116.309140 28302740 Andersen TA, Troelsen KDLL, Larsen LA. Of mice and men: Molecular genetics of congenital heart disease. Cell Mol Life Sci. 2014;71: 1327–1352. doi: 10.1007/s00018-013-1430-1 23934094 Sifrim A, Hitz M-P, Wilsdon A, Breckpot J, Turki SH Al, Thienpont B, et al. Distinct genetic architectures for syndromic and nonsyndromic congenital heart defects identified by exome sequencing. Nat Genet. 2016;48: 1060–5. doi: 10.1038/ng.3627 27479907 Jin SC, Homsy J, Zaidi S, Lu Q, Morton S, DePalma SR, et al. Contribution of rare inherited and de novo variants in 2,871 congenital heart disease probands. Nat Genet. 2017. doi: 10.1038/ng.3970 28991257 Tan TY, Gonzaga-Jauregui C, Bhoj EJ, Strauss KA, Brigatti K, Puffenberger E, et al. Monoallelic BMP2 Variants Predicted to Result in Haploinsufficiency Cause Craniofacial, Skeletal, and Cardiac Features Overlapping Those of 20p12 Deletions. Am J Hum Genet. 2017;101: 985–994. doi: 10.1016/j.ajhg.2017.10.006 29198724 Glessner JT, Bick AG, Ito K, Homsy J, Rodriguez-Murillo L, Fromer M, et al. Increased frequency of de novo copy number variants in congenital heart disease by integrative analysis of single nucleotide polymorphism array and exome sequence data. Circ Res. 2014;115: 884–896. doi: 10.1161/CIRCRESAHA.115.304458 25205790 Karczewski KJ, Francioli LC, Tiao G, Cummings BB, Alföldi J, Wang Q, et al. Variation across 141,456 human exomes and genomes reveals the spectrum of loss-of-function intolerance across human protein-coding genes. bioRxiv. 2019; 531210. doi: 10.1101/531210 Raychaudhuri S, Korn JM, McCarroll SA, International Schizophrenia Consortium, Altshuler D, Sklar P, et al. Accurately assessing the risk of schizophrenia conferred by rare copy-number variation affecting genes with brain function. Allison DB, editor. PLoS Genet. 2010;6: e1001097. doi: 10.1371/journal.pgen.1001097 20838587 McRae JF, Clayton S, Fitzgerald TW, Kaplanis J, Prigmore E, Rajan D, et al. Prevalence and architecture of de novo mutations in developmental disorders. Nature. 2017;542: 433–438. doi: 10.1038/nature21062 28135719 Kaplanis J, Samocha KE, Wiel L, Zhang Z, Kevin J, Eberhardt RY, et al. Integrating healthcare and research genetic data empowers the discovery of 49 novel developmental disorders. bioRxiv. 2019; 797787. doi: 10.1101/797787 Faundes V, Newman WG, Bernardini L, Canham N, Clayton-Smith J, Dallapiccola B, et al. Histone Lysine Methylases and Demethylases in the Landscape of Human Developmental Disorders. Am J Hum Genet. 2018;102: 175–187. doi: 10.1016/j.ajhg.2017.11.013 29276005 Martin HC, Jones WD, McIntyre R, Sanchez-Andrade G, Sanderson M, Stephenson JD, et al. Quantifying the contribution of recessive coding variation to developmental disorders. Science. 2018;362: 1161–1164. doi: 10.1126/science.aar6731 30409806 Ignatiadis N, Klaus B, Zaugg JB, Huber W. Data-driven hypothesis weighting increases detection power in genome-scale multiple testing. Nat Methods. 2016;13: 577–80. doi: 10.1038/nmeth.3885 27240256 Cardoso-Moreira M, Halbert J, Valloton D, Velten B, Chen C, Shao Y, et al. Gene expression across mammalian organ development. Nature. 2019;571: 505–509. doi: 10.1038/s41586-019-1338-5 31243369 Priest JR, Osoegawa K, Mohammed N, Nanda V, Kundu R, Schultz K, et al. De Novo and Rare Variants at Multiple Loci Support the Oligogenic Origins of Atrioventricular Septal Heart Defects. Jiao K, editor. PLoS Genet. 2016;12: e1005963. doi: 10.1371/journal.pgen.1005963 27058611 Reuter MS, Jobling R, Chaturvedi RR, Manshaei R, Costain G, Heung T, et al. Haploinsufficiency of vascular endothelial growth factor related signaling genes is associated with tetralogy of Fallot. Genet Med. 2019;21: 1001–1007. doi: 10.1038/s41436-018-0260-9 30232381 Li T, Kim A, Rosenbluh J, Horn H, Greenfeld L, An D, et al. GeNets: a unified web platform for network-based genomic analyses. Nat Methods. 2018;15: 543–546. doi: 10.1038/s41592-018-0039-6 29915188 Sevim Bayrak C, Zhang P, Tristani-Firouzi M, Gelb BD, Itan Y. De novo variants in exomes of congenital heart disease patients identify risk genes and pathways. Genome Med. 2020;12: 9. doi: 10.1186/s13073-019-0709-8 31941532 Lebrun N, Mehler-Jacob C, Poirier K, Zordan C, Lacombe D, Carion N, et al. Novel KDM5B splice variants identified in patients with developmental disorders: Functional consequences. Gene. 2018;679: 305–313. doi: 10.1016/j.gene.2018.09.016 30217758 Fotiou E, Williams S, Martin-Geary A, Robertson DL, Tenin G, Hentges KE, et al. Integration of Large-Scale Genomic Data Sources With Evolutionary History Reveals Novel Genetic Loci for Congenital Heart Disease. Circ Genomic Precis Med. 2019;12: 442–451. doi: 10.1161/CIRCGEN.119.002694 31613678 Nimura K, Ura K, Shiratori H, Ikawa M, Okabe M, Schwartz RJ, et al. A histone H3 lysine 36 trimethyltransferase links Nkx2-5 to Wolf-Hirschhorn syndrome. Nature. 2009;460: 287–291. doi: 10.1038/nature08086 19483677 DeSanto C, K D’Aco, Araujo GC, Shannon N, DDD Study, Vernon H, et al. WAC loss-of-function mutations cause a recognisable syndrome characterised by dysmorphic features, developmental delay and hypotonia and recapitulate 10p11.23 microdeletion syndrome. J Med Genet. 2015;52: 754–61. doi: 10.1136/jmedgenet-2015-103069 26264232 Wentzel C, Rajcan-Separovic E, Ruivenkamp CAL, Chantot-Bastaraud S, Metay C, Andrieux J, et al. Genomic and clinical characteristics of six patients with partially overlapping interstitial deletions at 10p12p11. Eur J Hum Genet. 2011;19: 959–964. doi: 10.1038/ejhg.2011.71 21522184 Okamoto N, Hayashi S, Masui A, Kosaki R, Oguri I, Hasegawa T, et al. Deletion at chromosome 10p11.23-p12.1 defines characteristic phenotypes with marked midface retrusion. J Hum Genet. 2012;57: 191–6. doi: 10.1038/jhg.2011.154 22258158 Chong JX, McMillin MJ, Shively KM, Beck AE, Marvin CT, Armenteros JR, et al. De novo mutations in NALCN cause a syndrome characterized by congenital contractures of the limbs and face, hypotonia, and developmental delay. Am J Hum Genet. 2015;96: 462–73. doi: 10.1016/j.ajhg.2015.01.003 25683120 Fukai R, Saitsu H, Okamoto N, Sakai Y, Fattal-Valevski A, Masaaki S, et al. De novo missense mutations in NALCN cause developmental and intellectual impairment with hypotonia. J Hum Genet. 2016;61: 451–5. doi: 10.1038/jhg.2015.163 26763878 Hoyer J, Ekici AB, Endele S, Popp B, Zweier C, Wiesener A, et al. Haploinsufficiency of ARID1B, a member of the SWI/SNF-A chromatin-remodeling complex, is a frequent cause of intellectual disability. Am J Hum Genet. 2012;90: 565–572. doi: 10.1016/j.ajhg.2012.02.007 22405089 Fitzgerald TW, Gerety SS, Jones WD, Van Kogelenberg M, King DA, McRae J, et al. Large-scale discovery of novel genetic causes of developmental disorders. Nature. Nature Publishing Group; 2015. pp. 223–228. doi: 10.1038/nature14135 25533962 van der Sluijs PJ, Jansen S, Vergano SA, Adachi-Fukuda M, Alanay Y, AlKindy A, et al. The ARID1B spectrum in 143 patients: from nonsyndromic intellectual disability to Coffin-Siris syndrome. Genet Med. 2019;21: 1295–1307. doi: 10.1038/s41436-018-0330-z 30349098 Kang E, Burdick KE, Kim JY, Duan X, Guo JU, Sailor KA, et al. Interaction between FEZ1 and DISC1 in regulation of neuronal development and risk for schizophrenia. Neuron. 2011;72: 559–71. doi: 10.1016/j.neuron.2011.09.032 22099459 Barth JL, Clark CD, Fresco VM, Knoll EP, Lee B, Argraves WS, et al. Jarid2 is among a set of genes differentially regulated by Nkx2.5 during outflow tract morphogenesis. Dev Dyn. 2010;239: 2024–33. doi: 10.1002/dvdy.22341 20549724 Yokoyama K, Tezuka T, Kotani M, Nakazawa T, Hoshina N, Shimoda Y, et al. NYAP: a phosphoprotein family that links PI3K to WAVE1 signalling in neurons. EMBO J. 2011;30: 4739–54. doi: 10.1038/emboj.2011.348 21946561 van Walree ES, Dombrowsky G, Jansen IE, Mirkov MU, Zwart R, Ilgun A, et al. Germline variants in HEY2 functional domains lead to congenital heart defects and thoracic aortic aneurysms. Genet Med. 2021;23: 103–110. doi: 10.1038/s41436-020-00939-4 32820247 Lindsay ME, Dietz HC. The genetic basis of aortic aneurysm. Cold Spring Harb Perspect Med. 2014;4. doi: 10.1101/cshperspect.a015909 25183854 Digilio MC, Gnazzo M, Lepri F, Dentici ML, Pisaneschi E, Baban A, et al. Congenital heart defects in molecularly proven Kabuki syndrome patients. Am J Med Genet Part A. 2017;173: 2912–2922. doi: 10.1002/ajmg.a.38417 28884922 Yuan SM. Congenital heart defects in Kabuki syndrome. Cardiology Journal. Cardiol J; 2013. pp. 121–124. doi: 10.5603/CJ.2013.0023 23558868 Baer S, Afenjar A, Smol T, Piton A, Gérard B, Alembik Y, et al. Wiedemann-Steiner syndrome as a major cause of syndromic intellectual disability: A study of 33 French cases. Clin Genet. 2018;94: 141–152. doi: 10.1111/cge.13254 29574747 Permanyer E, Laurie S, Blasco-Lucas A, Maldonado G, Amador-Catalan A, Ferrer-Curriu G, et al. A single nucleotide deletion resulting in a frameshift in exon 4 of TAB2 is associated with a polyvalular syndrome. Eur J Med Genet. 2020;63. doi: 10.1016/j.ejmg.2020.103854 31981616 Ackerman JP, Smestad JA, Tester DJ, Qureshi MY, Crabb BA, Mendelsohn NJ, et al. Whole Exome Sequencing, Familial Genomic Triangulation, and Systems Biology Converge to Identify a Novel Nonsense Mutation in TAB2-encoded TGF-beta Activated Kinase 1 in a Child with Polyvalvular Syndrome. Congenit Heart Dis. 2016;11: 452–461. doi: 10.1111/chd.12400 27452334 Ritelli M, Morlino S, Giacopuzzi E, Bernardini L, Torres B, Santoro G, et al. A recognizable systemic connective tissue disorder with polyvalvular heart dystrophy and dysmorphism associated with TAB2 mutations. Clin Genet. 2018;93: 126–133. doi: 10.1111/cge.13032 28386937 Thienpont B, Zhang L, Postma A V, Breckpot J, Tranchevent LC, Van Loo P, et al. Haploinsufficiency of TAB2 Causes Congenital Heart Defects in Humans. Am J Hum Genet. 2010;86: 839–849. doi: 10.1016/j.ajhg.2010.04.011 20493459 Cheng A, Neufeld-Kaiser W, Byers PH, Liu YJ. 6q25.1 (TAB2) microdeletion is a risk factor for hypoplastic left heart: A case report that expands the phenotype. BMC Cardiovasc Disord. 2020;20. doi: 10.1186/s12872-020-01328-0 31952508 Chapman G, Moreau JLM, I P E, Szot JO, Iyer KR, Shi H, et al. Functional genomics and gene-environment interaction highlight the complexity of congenital heart disease caused by Notch pathway variants. Hum Mol Genet. 2020;29: 566–579. doi: 10.1093/hmg/ddz270 31813956 Tomita-Mitchell A, Maslen CL, Morris CD, Garg V, Goldmuntz E. GATA4 sequence variants in patients with congenital heart disease. J Med Genet. 2007;44: 779–83. doi: 10.1136/jmg.2007.052183 18055909 Doornbos M, Sikkema-Raddatz B, Ruijvenkamp CAL, Dijkhuizen T, Bijlsma EK, Gijsbers ACJ, et al. Nine patients with a microdeletion 15q11.2 between breakpoints 1 and 2 of the Prader-Willi critical region, possibly associated with behavioural disturbances. Eur J Med Genet. 2009;52: 108–115. doi: 10.1016/j.ejmg.2009.03.010 19328872 Roberts AE, Allanson JE, Tartaglia M, Gelb BD. Noonan syndrome. The Lancet. Elsevier B.V.; 2013. pp. 333–342. doi: 10.1016/S2213-8587(13)70153-0 24703051 Øyen N, Poulsen G, Boyd HA, Wohlfahrt J, Jensen PKA, Melbye M. Recurrence of congenital heart defects in families. Circulation. 2009;120: 295–301. doi: 10.1161/CIRCULATIONAHA.109.857987 19597048 Wright CF, Fitzgerald TW, Jones WD, Clayton S, McRae JF, Van Kogelenberg M, et al. Genetic diagnosis of developmental disorders in the DDD study: A scalable analysis of genome-wide research data. Lancet. 2015;385: 1305–1314. doi: 10.1016/S0140-6736(14)61705-0 25529582 Zhang D, Qian Y, Akula N, Alliey-Rodriguez N, Tang J, Bipolar Genome Study, et al. Accuracy of CNV Detection from GWAS Data. PLoS One. 2011;6: e14511. doi: 10.1371/journal.pone.0014511 21249187 Singh T, Walters JTR, Johnstone M, Curtis D, Suvisaari J, Torniainen M, et al. The contribution of rare variants to risk of schizophrenia in individuals with and without intellectual disability. Nat Genet. 2017;49: 1167–1173. doi: 10.1038/ng.3903 28650482 Rehm HL, Berg JS, Brooks LD, Bustamante CD, Evans JP, Landrum MJ, et al. ClinGen—The clinical genome resource. N Engl J Med. 2015;372: 2235–2242. doi: 10.1056/NEJMsr1406261 26014595 Blekhman R, Man O, Herrmann L, Boyko AR, Indap A, Kosiol C, et al. Natural selection on genes that underlie human disease susceptibility. Curr Biol. 2008;18: 883–9. doi: 10.1016/j.cub.2008.04.074 18571414 Berg JS, Adams M, Nassar N, Bizon C, Lee K, Schmitt CP, et al. An informatics approach to analyzing the incidentalome. Genet Med. 2013;15: 36–44. doi: 10.1038/gim.2012.112 22995991 Freeman PJ, Hart RK, Gretton LJ, Brookes AJ, Dalgleish R. VariantValidator: Accurate validation, mapping, and formatting of sequence variation descriptions. Hum Mutat. 2018;39: 61–68. doi: 10.1002/humu.23348 28967166 Samocha KE, Robinson EB, Sanders SJ, Stevens C, Sabo A, McGrath LM, et al. A framework for the interpretation of de novo mutation in human disease. Nat Genet. 2014;46: 944–950. doi: 10.1038/ng.3050 25086666 Nueda MJ, Tarazona S, Conesa A. Next maSigPro: updating maSigPro bioconductor package for RNA-seq time series. Bioinformatics. 2014;30: 2598–602. doi: 10.1093/bioinformatics/btu333 24894503 Li T, Wernersson R, Hansen RB, Horn H, Mercer J, Slodkowicz G, et al. A scored human protein-protein interaction network to catalyze genomic interpretation. Nat Methods. 2017;14: 61–64. doi: 10.1038/nmeth.4083 27892958 Herwig R, Hardt C, Lienhard M, Kamburov A. Analyzing and interpreting genome data at the network level with ConsensusPathDB. Nat Protoc. 2016;11: 1889–907. doi: 10.1038/nprot.2016.117 27606777 Dilg D, Saleh RNM, Phelps SEL, Rose Y, Dupays L, Murphy C, et al. HIRA is required for heart development and directly regulates Tnni2 and Tnnt3. Fraidenraich D, editor. PLoS One. 2016;11: e0161096. doi: 10.1371/journal.pone.0161096 27518902 10.1371/journal.pgen.1009679.r001 Decision Letter 0 Barsh Gregory S. Editor-in-Chief Firulli Anthony B Associate Editor 2021 Barsh, Firulli This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Submission Version 0

4 Feb 2021

Dear Dr Hitz,

Thank you very much for submitting your Research Article entitled 'Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease' to PLOS Genetics.

The manuscript was fully evaluated at the editorial level and by independent peer reviewers. The reviewers appreciated the attention to an important problem, but raised some substantial concerns about the current manuscript. Based on the reviews, we will not be able to accept this version of the manuscript, but we would be willing to review a much-revised version. We cannot, of course, promise publication at that time.

Should you decide to revise the manuscript for further consideration here, your revisions should address the specific points made by each reviewer. We will also require a detailed list of your responses to the review comments and a description of the changes you have made in the manuscript.

If you decide to revise the manuscript for further consideration at PLOS Genetics, please aim to resubmit within the next 60 days, unless it will take extra time to address the concerns of the reviewers, in which case we would appreciate an expected resubmission date by email to plosgenetics@plos.org.

If present, accompanying reviewer attachments are included with this email; please notify the journal office if any appear to be missing. They will also be available for download from the link below. You can use this link to log into the system when you are ready to submit a revised version, having first consulted our Submission Checklist.

To enhance the reproducibility of your results, we recommend that you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. For instructions see our guidelines.

Please be aware that our data availability policy requires that all numerical data underlying graphs or summary statistics are included with the submission, and you will need to provide this upon resubmission if not already present. In addition, we do not permit the inclusion of phrases such as "data not shown" or "unpublished results" in manuscripts. All points should be backed up by data provided with the submission.

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool.  PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org.

PLOS has incorporated Similarity Check, powered by iThenticate, into its journal-wide submission system in order to screen submitted content for originality before publication. Each PLOS journal undertakes screening on a proportion of submitted articles. You will be contacted if needed following the screening process.

To resubmit, use the link below and 'Revise Submission' in the 'Submissions Needing Revision' folder.

[LINK]

We are sorry that we cannot be more positive about your manuscript at this stage. Please do not hesitate to contact us if you have any concerns or questions.

Yours sincerely,

Anthony B Firulli

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

Reviewer's Responses to Questions

Comments to the Authors:

Please note here if the review is uploaded as an attachment.

Reviewer #1: This paper by Audain et al is a massive study collecting CNVs/DNVs data across multiple sources in an attempt to increase the power to determine involvement in CHD, and then explore these CNVs and DNVs to find individual genes and networks. They focused only on solidly pathogenic changes due to the nature of the case data, a trade off to using all CNVs that could identify new CNVs worth consideration as novel pathogenic variants. A large control group was selected across several datasets. The authors have used orthogonal analyses for several parts of the study (such as using Mupit and DeNovoWEST) when no one gold standard exists, and then considered the genes in the CNVs in different ways (straight association, as part of heart tissue expression, and then as differential expression, and as a network analysis). Overall, the study is quite well executed, and the paper provides nice confirmation of known information (specific CNVs, known genes and involvement of specific pathways). In addition, they identify several novel genes that warrant investigation.

The methods are scant in the description of the study design. The authors should state their search strategy for the literature and the public datasets chosen. Of the studies identified, how were they selected for inclusion or exclusion? Were the CNV boundaries determined using the same build of the genome? How did the authors determine the controls and cases were of similar ancestral groups? While this is likely less of a problem for CNVs, it would be helpful to note any differences between groups. Did the authors include all genes within a CNV, or was there some selection?

A definition of the different gene sets used in the analysis is needed. For example, it is not clear what difference there is between monoallelic and non-syndromic sets.

Was the phenotype information clean enough and were there enough individuals to explore the CHDs by subgroups? The genetic architecture is quite different for instance between conotruncal defects and left sided defects.

The authors should also acknowledge an additional bias of the case data collected. The CNVs are highly selected from only the cases published due the presence of a CNV so are not representative of all individuals with CHD.

Reviewer #2: In this manuscript, the authors present their study on an integrative analysis to identify copy number variants (CNVs) and de novo intragenic variants (DNVs) in patients with congenital heart disease (CHD) with sporadic thoracic aortic aneurysm. The authors used published and publicly available datasets of CNVs and DNVs identified in probands with CHD (10,447) to perform a genome-wide integrative meta-analysis. To achieve their goal the authors implemented a statistical approach which allows the integration of different types of genomic variants. Based on their study they identified 21 significant haploinsufficient CHD genes. Among these genes 14 are already known as CHD genes and 7 have not previously been associated with CHD. The authors evaluated expression pattern of these genes in the heart using available data from the recent human atlas study. In summary, this study used astutely available datasets to identify novel CHD genes.

Although the basic insight regarding the identification of novel genes involved in CHD is very interesting, there are some weaknesses in the data reported in this study.

Major concerns:

- My main concern is related to diversity of CHD used in this study. The authors used CHD as a generic denomination. There is no detail regarding the type and class of CHD used in this study. For instance, conotruncal and non-conotruncal groups would be a good indicator to better classify these CHDs in order to improve the analysis.

- I wonder if there is a bias in using these publicly available datasets since their integrative analysis revealed significant enrichment for the known CHD genes, which is mainly explained by the contribution of monoallelic CHD genes and syndromic CHD gene set. This should be better discussed and the authors should propose alternative analysis to identify non-syndromic CHD genes.

- It is very surprising to see that one of the most frequent syndromes (22q11/DiGeorge syndrome) is not among the CNV deletion syndrome regions identified in this study. Particularly, the authors mentioned that the distribution of CNVs that overlap known microdeletion syndrome such as DiGeorge syndrome is overrepresented in the dataset. I wonder about the quality of the analysis and the filters applied to identify CHD genes.

- The authors used bulk-RNAseq dataset to determine if the 21 significant genes are expressed in the heart at different developmental time points. However, there is no information regarding the region of the heart where these genes are expressed. This information could be useful in relation to the type of CHD they are associated with.

Minor points:

- S6 Table: It is not evident to find the 528 genes. The authors should better highlight them).

- The authors should better explain why it is important to known that genes remain highly expressed in adulthood.

- It is not clear why the authors have conducted a systems-level analysis to identify global mechanisms if there is no mechanism proposed in this study.

Reviewer #3: The manuscript by Audain et al. performs a meta-analysis of copy number variants and de novo variants identified by exome sequencing from previously published studies. The case and control cohort sizes are large and the authors have established an extensive cohort. Using burden analyses the authors identified 21 genes associated with CHD including 14 previously known genes and 7 “novel” genes. In many ways this manuscript provides proof of principle that meta-analyses from a large number of studies, including single case studies, can be performed and platform results can be harmonized for both DNV and CNVs and this is an important contribution. However, the novel findings are limited and the goal of using integrative analyses to identify CHD disease genes in multifactorial inheritance is not accomplished here but it does provide a path forward. It would improve the impact of the manuscript to apply the expression and/or network analyses to a somewhat broader group of candidate genes.

Major comments:

1. The study focuses on CNV deletions but it would be useful to perform analyses for duplications as well particularly for overlaps with the deletions as this would further support dosage sensitivity and mechanisms.

2. The control population is not described well. How carefully phenotyped are these controls?

3. Please clarify how genes were assigned to syndromic, non-syndromic, monoallelic, biallelic gene sets and other gene sets. What proportion of the genes are shared between lists and how is this accounted for in analyses?

4. Previous studies from some of the authors of this study have demonstrated enrichment of CNVs in specific subphenotypes of CHD. Additional information and analyses performed on subphenotypes would be useful to include as well as whether any previous findings can be expanded.

5. For TAA, previous studies from some of the authors indicate a higher burden of CNVs in familial TAA. What is the representation of familial vs sporadic in this meta-analysis and can any additional comments be made about CNV burden?

6. In Table 7 the method for identifying CNVs in controls needs explanation or better descriptors. Why are the CNV counts so high in well described genomic disorders such as 22q11.2/VCFS/DiGeorge (986 controls), 1p36 deletion syndrome (886 controls) or Cri du Chat (1488 controls). These values should also have established a MAF > .01 prompting filtering. More broadly, it is surprising that TBX1 was not identified in the burden analyses. Could the authors comment particularly as it was noted in the discussion that 22q11.2 cases were overrepresented and so TBX1 seemingly should have shown significant association by burden testing.

7. Some of the “novel” genes cause syndromes that do have well described CHD such as Wolf-Hirschhorn or ARID1B. The manuscript could be improved by a better discussion of modifiers or oligogenic mechanisms and reduced penetrance of CHD as the rule even in genetic syndromic conditions.

8. The data on TAA do not contribute significantly to the manuscript.

9. The manuscript is difficult to follow from an analysis standpoint at times. It would benefit from a flow diagram of analyses and/or including more of the methods within the manuscript.

Minor:

1. The size cutoff of 5 kb and 20 Mb is relevant and should be moved to the methods of the main manuscript.

2. Fig. S4 – what does purple denote?

3. Fig. 5 - the x axis (Score) needs an explanation in the legend

4. Heading p. 12 CNV/DNVs hinder the function of specific protein complexes should be reworded since this has not been proven by this study.

**********

Have all data underlying the figures and results presented in the manuscript been provided?

Large-scale datasets should be made available via a public repository as described in the PLOS Genetics data availability policy, and numerical data that underlies graphs or summary statistics should be provided in spreadsheet form as supporting information.

Reviewer #1: Yes

Reviewer #2: Yes

Reviewer #3: Yes

**********

PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #1: Yes: Kim L McBride

Reviewer #2: No

Reviewer #3: No

 10.1371/journal.pgen.1009679.r002 Author response to Decision Letter 0 Submission Version 1

23 Apr 2021

Submitted filename: Response_to_Reviewers_PGEN.docx

 10.1371/journal.pgen.1009679.r003 Decision Letter 1 Barsh Gregory S. Editor-in-Chief Firulli Anthony B Associate Editor 2021 Barsh, Firulli This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Submission Version 1

24 May 2021

Dear Dr Hitz,

Thank you very much for submitting your Research Article entitled 'Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease' to PLOS Genetics.

The manuscript was fully evaluated at the editorial level and by independent peer reviewers. The reviewers appreciated the attention to an important topic but identified some concerns that we ask you address in a revised manuscript

We therefore ask you to modify the manuscript according to the review recommendations. Your revisions should address the specific points made by each reviewer.

In addition we ask that you:

1) Provide a detailed list of your responses to the review comments and a description of the changes you have made in the manuscript.

2) Upload a Striking Image with a corresponding caption to accompany your manuscript if one is available (either a new image or an existing one from within your manuscript). If this image is judged to be suitable, it may be featured on our website. Images should ideally be high resolution, eye-catching, single panel square images. For examples, please browse our archive. If your image is from someone other than yourself, please ensure that the artist has read and agreed to the terms and conditions of the Creative Commons Attribution License. Note: we cannot publish copyrighted images.

We hope to receive your revised manuscript within the next 30 days. If you anticipate any delay in its return, we would ask you to let us know the expected resubmission date by email to plosgenetics@plos.org.

If present, accompanying reviewer attachments should be included with this email; please notify the journal office if any appear to be missing. They will also be available for download from the link below. You can use this link to log into the system when you are ready to submit a revised version, having first consulted our Submission Checklist.

While revising your submission, please upload your figure files to the Preflight Analysis and Conversion Engine (PACE) digital diagnostic tool. PACE helps ensure that figures meet PLOS requirements. To use PACE, you must first register as a user. Then, login and navigate to the UPLOAD tab, where you will find detailed instructions on how to use the tool. If you encounter any issues or have any questions when using PACE, please email us at figures@plos.org.

Please be aware that our data availability policy requires that all numerical data underlying graphs or summary statistics are included with the submission, and you will need to provide this upon resubmission if not already present. In addition, we do not permit the inclusion of phrases such as "data not shown" or "unpublished results" in manuscripts. All points should be backed up by data provided with the submission.

To enhance the reproducibility of your results, we recommend that you deposit your laboratory protocols in protocols.io, where a protocol can be assigned its own identifier (DOI) such that it can be cited independently in the future. Additionally, PLOS ONE offers an option to publish peer-reviewed clinical study protocols. Read more information on sharing protocols at https://plos.org/protocols?utm_medium=editorial-email&utm_source=authorletters&utm_campaign=protocols

Please review your reference list to ensure that it is complete and correct. If you have cited papers that have been retracted, please include the rationale for doing so in the manuscript text, or remove these references and replace them with relevant current references. Any changes to the reference list should be mentioned in the rebuttal letter that accompanies your revised manuscript. If you need to cite a retracted article, indicate the article’s retracted status in the References list and also include a citation and full reference for the retraction notice.

PLOS has incorporated Similarity Check, powered by iThenticate, into its journal-wide submission system in order to screen submitted content for originality before publication. Each PLOS journal undertakes screening on a proportion of submitted articles. You will be contacted if needed following the screening process.

To resubmit, you will need to go to the link below and 'Revise Submission' in the 'Submissions Needing Revision' folder.

[LINK]

Please let us know if you have any questions while making these revisions.

Yours sincerely,

Anthony B Firulli

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

Reviewer's Responses to Questions

Comments to the Authors:

Please note here if the review is uploaded as an attachment.

Reviewer #2: The manuscript by Audain et al. has been well revised. The authors have completed significant revisions based on my comments. I find the current version of the manuscript considerably improved. The authors have replied to my main concern regarding the diagnosis of specific CHD type. They have performed a sub-analysis of the data based on specific groups of heart defects. They have also replied to my question concerning the possibly bias of their analysis since they used previously published datasets. The current version included a discussion on this specific point. The authors have satisfactorily answered my question regarding the 22q11 locus. Absence of detailed information about spatiotemporal expression of identified genes with CHD does not confirm the link these genes in the different part of the heart and the type of CHD in which they are involved. The authors should have performed a better analysis to reinforce this specific issue.

Reviewer #3: The authors have been responsive to the reviewers’ comments within the limitations of available data. The addition of the CHD phenotype subclassifications is useful and the discussion/interpretation of the findings is appropriate. The requested points for clarification have mostly been addressed. I have only minor comments.

I did not note in the original submission that CNVs were compiled for autosomes only. I would suggest the discussion/limitations make more explicit the fact that X-linked candidate genes will not be identified by the integrative approach.

The incorporation of the detailed methods is helpful. For the older platforms in the study, was the minimum CNV interval or maximum CNV interval (< 20 Mb) used? Are platform differences the explanation for the persistent differences between cases and controls even after filtering (FigS5)?

**********

Have all data underlying the figures and results presented in the manuscript been provided?

Large-scale datasets should be made available via a public repository as described in the PLOS Genetics data availability policy, and numerical data that underlies graphs or summary statistics should be provided in spreadsheet form as supporting information.

Reviewer #2: Yes

Reviewer #3: Yes

**********

PLOS authors have the option to publish the peer review history of their article (what does this mean?). If published, this will include your full peer review and any attached files.

If you choose “no”, your identity will remain anonymous but your review may still be made public.

Do you want your identity to be public for this peer review? For information about this choice, including consent withdrawal, please see our Privacy Policy.

Reviewer #2: Yes: Stéphane ZAFFRAN

Reviewer #3: No

 10.1371/journal.pgen.1009679.r004 Author response to Decision Letter 1 Submission Version 2

18 Jun 2021

Submitted filename: Response_to_Reviewers_PGEN.docx

 10.1371/journal.pgen.1009679.r005 Decision Letter 2 Barsh Gregory S. Editor-in-Chief Firulli Anthony B Associate Editor 2021 Barsh, Firulli This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Submission Version 2

23 Jun 2021

Dear Dr Hitz,

We are pleased to inform you that your manuscript entitled "Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease" has been editorially accepted for publication in PLOS Genetics. Congratulations!

Before your submission can be formally accepted and sent to production you will need to complete our formatting changes, which you will receive in a follow up email. Please be aware that it may take several days for you to receive this email; during this time no action is required by you. Please note: the accept date on your published article will reflect the date of this provisional acceptance, but your manuscript will not be scheduled for publication until the required changes have been made.

Once your paper is formally accepted, an uncorrected proof of your manuscript will be published online ahead of the final version, unless you’ve already opted out via the online submission form. If, for any reason, you do not want an earlier version of your manuscript published online or are unsure if you have already indicated as such, please let the journal staff know immediately at plosgenetics@plos.org.

In the meantime, please log into Editorial Manager at https://www.editorialmanager.com/pgenetics/, click the "Update My Information" link at the top of the page, and update your user information to ensure an efficient production and billing process. Note that PLOS requires an ORCID iD for all corresponding authors. Therefore, please ensure that you have an ORCID iD and that it is validated in Editorial Manager. To do this, go to ‘Update my Information’ (in the upper left-hand corner of the main menu), and click on the Fetch/Validate link next to the ORCID field.  This will take you to the ORCID site and allow you to create a new iD or authenticate a pre-existing iD in Editorial Manager.

If you have a press-related query, or would like to know about making your underlying data available (as you will be aware, this is required for publication), please see the end of this email. If your institution or institutions have a press office, please notify them about your upcoming article at this point, to enable them to help maximise its impact. Inform journal staff as soon as possible if you are preparing a press release for your article and need a publication date.

Thank you again for supporting open-access publishing; we are looking forward to publishing your work in PLOS Genetics!

Yours sincerely,

Anthony B Firulli

Associate Editor

PLOS Genetics

Gregory Barsh

Editor-in-Chief

PLOS Genetics

www.plosgenetics.org

Twitter: @PLOSGenetics

----------------------------------------------------

Comments from the reviewers (if applicable):

----------------------------------------------------

Data Deposition

If you have submitted a Research Article or Front Matter that has associated data that are not suitable for deposition in a subject-specific public repository (such as GenBank or ArrayExpress), one way to make that data available is to deposit it in the Dryad Digital Repository. As you may recall, we ask all authors to agree to make data available; this is one way to achieve that. A full list of recommended repositories can be found on our website.

The following link will take you to the Dryad record for your article, so you won't have to re‐enter its bibliographic information, and can upload your files directly: 

http://datadryad.org/submit?journalID=pgenetics&manu=PGENETICS-D-20-01841R2

More information about depositing data in Dryad is available at http://www.datadryad.org/depositing. If you experience any difficulties in submitting your data, please contact help@datadryad.org for support.

Additionally, please be aware that our data availability policy requires that all numerical data underlying display items are included with the submission, and you will need to provide this before we can formally accept your manuscript, if not already present.

----------------------------------------------------

Press Queries

If you or your institution will be preparing press materials for this manuscript, or if you need to know your paper's publication date for media purposes, please inform the journal staff as soon as possible so that your submission can be scheduled accordingly. Your manuscript will remain under a strict press embargo until the publication date and time. This means an early version of your manuscript will not be published ahead of your final version. PLOS Genetics may also choose to issue a press release for your article. If there's anything the journal should know or you'd like more information, please get in touch via plosgenetics@plos.org.

 10.1371/journal.pgen.1009679.r006 Acceptance letter Barsh Gregory S. Editor-in-Chief Firulli Anthony B Associate Editor 2021 Barsh, Firulli This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

26 Jul 2021

PGENETICS-D-20-01841R2

Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease

Dear Dr Hitz,

We are pleased to inform you that your manuscript entitled "Integrative analysis of genomic variants reveals new associations of candidate haploinsufficient genes with congenital heart disease" has been formally accepted for publication in PLOS Genetics! Your manuscript is now with our production department and you will be notified of the publication date in due course.

The corresponding author will soon be receiving a typeset proof for review, to ensure errors have not been introduced during production. Please review the PDF proof of your manuscript carefully, as this is the last chance to correct any errors. Please note that major changes, or those which affect the scientific understanding of the work, will likely cause delays to the publication date of your manuscript.

Soon after your final files are uploaded, unless you have opted out or your manuscript is a front-matter piece, the early version of your manuscript will be published online. The date of the early version will be your article's publication date. The final article will be published to the same URL, and all versions of the paper will be accessible to readers.

Thank you again for supporting PLOS Genetics and open-access publishing. We are looking forward to publishing your work!

With kind regards,

Agota Szep

PLOS Genetics

On behalf of:

The PLOS Genetics Team

Carlyle House, Carlyle Road, Cambridge CB4 3DN | United Kingdom

plosgenetics@plos.org | +44 (0) 1223-442823

plosgenetics.org | Twitter: @PLOSGenetics