Impact Factor 3.258 | CiteScore 2.7
More on impact ›

Original Research ARTICLE

Front. Genet., 10 September 2019 | https://doi.org/10.3389/fgene.2019.00819

Association of Rare Recurrent Copy Number Variants With Congenital Heart Defects Based on Next-Generation Sequencing Data From Family Trios

Yichuan Liu1†, Xiao Chang1†, Joseph Glessner1†, Huiqi Qu1, Lifeng Tian1, Dong Li1, Kenny Nguyen1, Patrick M. A. Sleiman1,2 and Hakon Hakonarson1,2,3*
  • 1Center for Applied Genomics, Children’s Hospital of Philadelphia, Philadelphia, PA, United States
  • 2Division of Human Genetics, Department of Pediatrics, The Perelman School of Medicine, University of Pennsylvania, Philadelphia, PA, United States
  • 3Department of Human Genetics, Children’s Hospital of Philadelphia, Philadelphia, PA, United States

Congenital heart defects (CHDs) are a common birth defect, affecting approximately 1% of newborn children in the United States. As previously reported, a significant number of CHDs are potentially attributed to altered copy number variants (CNVs). However, as many genomic variants are rare, a large-scale CNV triad study is necessary to characterize the genetic architecture of CHD. We used whole-exome sequencing (WES) data generated by the Pediatric Cardiac Genomics Consortium (PCGC), including a discovery dataset of 2,103 individuals from 760 nuclear family trios and an independent replication set of 4,808 individuals from 1,712 trios. The candidate targets uncovered were further validated through different platforms, including the Omni single-nucleotide polymorphism (SNP) array chip in 1,860 individuals and the whole-genome sequencing (WGS) data in 33 trios. The genes harboring CNVs of interest were then investigated for expression alternations based on cardiac tissue RNA-Seq data. We identified multiple CNVs in the WES data that associated with specific sub-phenotypes of CHD in approximately 2,400 families, including 98 de novo CNV regions. We identified five CNV loci harboring LIMS1, GCC2, RANBP2, TTC3, and MAP3K7CL, respectively, where those genes are highly expressed in human heart and/or mouse embryo heart at 15 days. Five novel CNV loci were uncovered, demonstrating altered expression of the respective candidate genes involved. To our knowledge, this is the largest trio-based WES study of CHD and, in addition to uncovering novel CHD targets, presents an extensive resource with the potential to provide important insights to the architecture and impact of CNVs in CHD.

Introduction

Congenital heart defect (CHD) is the most common class of major congenital anomalies in humans and a major source of morbidity and pediatric mortality around the world (van der Linde et al., 2011; van der Bom et al., 2011). The incidence estimates range from 4 to 10 in 1,000 live births (Hoffman and Kaplan, 2002). Genome-wide rare copy number variants (CNVs) with a minor allele frequency (MAF) < 1% are now recognized as an important contributor to CHD (Silversides et al., 2012; Soemedi et al., 2012; Greenway et al., 2009; Fakhro et al., 2011). Burden of large CNVs has been observed across CHD subtypes. The canonical 3-Mb 22q11.2 deletion has been observed as the most common recurrent de novo CNV associated with syndromic CHD. Recurrent de novo CNVs in patients with CHD reported in multiple studies also occur at chromosomes 1q21.1, 3p25.1, 7q11.13, 8p23.1, 11q24–25, and 16p13.11. Using exome paired with dense arrays, causal genes in these intervals have been identified, including ELN (Williams syndrome), RAI1 (Smith–Magenis syndrome), TBX1 (22q11 deletion), GATA4 (8p23.1 deletion), GJA5 (1q21.1 duplication), and NKX2.5 (5q35.1 deletion). Recently, several large-scale studies have been conducted to explore the role of CNVs in CHD (Sifrim et al., 2016; Homsy et al., 2015; Glessner et al., 2014). Strikingly, recurrent de novo CNVs at 15q11.2 encompassing CYFIP1, NIPA1, and NIPA2 were identified. Genes that interact with established CHD proteins NKX2-5 and GATA4 had singular de novo CNVs encompassing DUSP1, JUN, JUP, MED15, MED9, PTPRE SREBF1, TOP2A, and ZEB2. However, due to limited sample sizes, platform limitations, and inadequate control sets, significant potential to identify new CHD-related genes exists, especially for de novo genomic CHD CNVs. In this study, we used the largest trio dataset available for CHD to search for rare CNVs associated with CHD. The dataset included a discovery set of 2,103 individuals from 760 nuclear family trios, all whole exome sequenced (WES), and a mutually exclusive and independent replication set of 4,808 individuals from 1,712 trios with WES data. Other independent replication datasets available included an Omni single-nucleotide polymorphism (SNP) array dataset of 1,860 individuals (922 cases and 938 controls), and a whole-genome sequencing (WGS) dataset of 99 individuals from 33 trios; all individuals in the replication sets are mutually exclusive from the discovery set. The CNV target genes were further validated based on a differential gene expression test in 55 cardiac biopsy samples from 55 CHD patients from the replication dataset. As CHD is associated with other diseases, including developmental delays, we also consider other disease-associated CNVs beside de novo CNVs. Collectively, we uncovered five CNV loci that associated with CHD. These loci corresponded to five genes, LIMS1, GCC2, RANBP2, TTC3, and MAP3K7CL, all of which were altered in expression due to the presence of CNVs we uncovered.

Materials and Methods

Patient Cohorts

The sequence data and phenotypes were downloaded from the Database of Genotypes and Phenotypes (dbGaP) (accession phs001194.v2.p2 and phs001194.v2.p2.c1) released by the Pediatric Cardiac Genomics Consortium (PCGC). CHD probands and parents were recruited into the CHD Genes Study of the PCGC (CHD genes: ClinicalTrials.gov identifier NCT01196182) as previously described (Pediatric Cardiac Genomics C et al., 2013), using protocols approved by the institutional review boards (IRBs) of each institution. The CHD trios selected for this study had no history of CHD in first-degree relatives, for example, parents and direct siblings. CHD diagnoses were obtained from echocardiograms, catheterization, and operative reports; extra-cardiac findings were extracted from medical records and included dysmorphic features, major anomalies, noncardiac medical problems, and deficiencies in growth or developmental delay. The etiologies for CHD were unknown; patients with previously identified cytogenetic anomalies or pathogenic CNVs identified through routine clinical evaluation were excluded. Whole-blood samples were collected, and genomic DNA was extracted. Cardiac tissue samples were collected during surgery for RNA extraction.

Patients recruited were divided into several datasets (Table 1), including a discovery dataset for WES (2,103 individuals from 760 nuclear family trios, including 652 families with at least one parent/sibling); a replication dataset for WES (4,808 individuals from 1,712 family trios, including 1,570 families with at least one parent/sibling); SNP array replication dataset (1,860 unrelated individuals); and WGS replication set (99 individuals from 33 family trios), and individuals in all replication datasets are mutually exclusive to the discovery set. An RNA-sequencing dataset was applied as the functional validation dataset for gene expression alternations (55 cardiac tissue samples). All the controls in these datasets are the first-degree relatives (parents or siblings) of the CHD patients except for the SNP array replication set. The discovery set is exclusive to the replication and validation datasets.

TABLE 1
www.frontiersin.org

Table 1 Summary of the patient cohorts and genomic technologies employed in this study.

Generating CNV Targets

In order to identify causal CNVs in the CHD samples, we set up a discovery pipeline based on several criteria (Figure 1). Methods applied to identify CNVs in different platforms were described in the following paragraphs. In general, for the CNVs found in the discovery set, we required that target CNVs were not inherited nor present in any independent family controls. The resulting CNVs were subsequently checked for individual replication in the independent datasets, and the CNVs were removed if they were found in any of the replication control datasets. In other words, the filtered discovery set (Table S1) required the CNVs be de novo while not identified in any controls and at least replicated in one set. After this step, additional filters were applied based on the prior biological knowledge of CHD, including the following: 1) the CNV corresponding genes are highly expressed in human or mouse embryo during development at 14.5 days (Zaidi et al., 2013); 2) the CNV is not physically overlapped with any previously reported CHD-associated/causal CNVs based on genomic locus (Glessner et al., 2014); 3) the CNV has not been previously reported as common CNVs in healthy individuals through the Database of Genomic Variations (DGV) (MacDonald et al., 2014). CNVs in Table S1 that fit at least one of these requirements were selected as a higher confidence target list (Table S2).

FIGURE 1
www.frontiersin.org

Figure 1 Copy number variant (CNV) identification and functional validation pipeline. This figure describes the conceptual pipeline for CNV target selection based on large sequencing/chip data of different platforms. The CNVs for the whole-exome sequencing (WES) data in the discovery set from 760 trios were detected based on the XHMM pipeline. Any inherited CNVs or CNVs identified in independent healthy controls were filtered at the first level; the remaining CNVs were further replicated in three exclusive datasets at different platforms and checked through previously reported results for further confirmation; the expressions for CNV-impacted genes were compared between the patients with/without the CNV based on RNA-Seq of cardiac tissues. CNVs associated with altered gene expression are listed in Table 2 as the targets for CHD.

Statistical Analysis

CNV calling software with different algorithms has varied reproducibility for NGS data. A previous study showed that CNV callers with read depth algorithm (XHMM, ExomeDepth, and CODEX) are beyond performance compared with others, especially in rare CNV detection (D O’Fallon et al., 2018). XHMM is less sensitive than ExomeDepth and CODEX, but much more precise than them (93% vs 54% and 72%) (Sadedin et al., 2018). In this study, CNV calling was performed using the standard XHMM pipeline consisting of six steps (Fromer et al., 2012). 1) The depth of coverage for all targets and all CHD samples used for subtype comparison was performed using GATK. 2) Target regions with extreme GC content (< 10% or >90%), and low-complexity regions were filtered out from further analysis. 3) Principal component analysis (PCA) normalization of read depth was performed for all samples to remove inherent biases in sample preparation and sequencing. 4) Samples with extreme variability in normalized read depth were removed. 5) Per-sample CNV detection with a hidden Markov model (HMM) was performed. 6) Quality metrics were assigned to all samples for discovered CNVs. One thousand eight hundred sixty individuals were genotyped on the Omni1M and Omni2.5M arrays. CNV calls were generated by the HMM-based software PennCNV as we described before (Wang et al., 2008). WGS CNVs were called by standard pipelines of BreakDancer (Fan et al., 2014), and we required at least 50% of reads to support the CNV identification. For RNA-Seq data of the corresponding CNV protein coding genes, we compared the estimated expression levels that were measured as fragments per kilobase of transcript per million (FPKM) in CHD individuals with/without CNVs by Cuffdiff and Cuffnorm (Trapnell et al., 2010).

Results

CNV Targets

Due to the high prevalence of CNVs across human genomes, the most significant difficulty of CNV analysis is to identify and remove the CNVs that are unrelated to a disease, such as CHD (Stankiewicz and Lupski, 2010). This study used three critical steps to remove CNV hits irrelevant to CHD. A recent study based on a large WES dataset showed that CNVs inherited from unaffected parents were usually benign in comparison with CNVs that were de novo and unique to the proband (Sifrim et al., 2016). Therefore, at the first step, we removed all the CNV hits which either were inherited or were found to overlap with other independent controls of our dataset before we check the replication dataset. Considering the limitation of this study that different genomic platforms with different sensitivities and preferential detection of different types of CNVs were used for CNV calling, at the second step, we check the CNV hits through three mutually exclusive independent replication datasets from different genomic platforms, including WES, SNP array chips, and WGS, in order to determine support for CNV causality of CHD. We identified 122 CNV regions which are supported by at least one replication dataset (Table S1). For biological relevance, we screened through previously published databases and reduced the candidate CNV targets group to 39 (Table S2). In addition, a list of 98 de novo CNVs never reported in any publication before based on the DGV database, including 35 recurrent CNVs, was generated for the CHD probands (Table S3).

Functional Validation Based on Cardiac Tissues

As DNA structural variations, such as CNVs, can lead to expression alternations for corresponding genes by copy number change or duplication/deletion of certain exons or regulatory elements, we examined the corresponding gene expression alternations between probands harboring the target CNVs and probands that do not. Due to the difficulty in collecting healthy people’s cardiac tissues, the validation was done in probands who had cardiac biopsy available through surgery. Four CNV regions were validated through this process (Tables 2 and 3, Figure 2).

TABLE 2
www.frontiersin.org

Table 2 Identification of CNV targets in CHD.

TABLE 3
www.frontiersin.org

Table 3 RNA-Seq validation for target CNVs.

FIGURE 2
www.frontiersin.org

Figure 2 Corresponding gene expression for target copy number variants (CNVs). These figures represent the corresponding CNV gene expression [Y axis of normalized fragments per kilobase of transcript per million (FPKM) value ranged from 0 to 100] for CHD probands who contain the CNV (case) versus CHD probands who do not contain the CNV (controls). Each bar is a CHD proband, and each cluster of bars indicates a certain gene. Housekeeping genes G6PD and ACTB were selected as the baseline control to normalize gene expression levels. (A) Deletion for LIMS1 and GCC2 genes at chr2:109113426-109287329 in aorta. (B) Deletion for RANBP2 gene at chr2:109363161-109371723 in aorta. (C) Deletion for TTC3 gene at chr21:38461093-38523202 in right atrium. (D) Deletion for MAP3K7CL AT chr21:30400216-30547213 in aorta.

We uncovered five CNV targets as shown in Table 2. A CNV residing at chr21:30400216-30547213 is a de novo deletion identified in a single CHD proband (1–01968) and was replicated in multiple dataset including the WES, SNP array, and WGS datasets as shown in Table 2. The corresponding gene MAP3K7CL is overexpressed in arteries, including coronary arteries and the aorta, and associated with coronary artery disease (Miller et al., 2016), and a recent genome-wide association study (GWAS) for heart failure suggests it is the top risk gene for coronary artery disease (Aragam et al., 2018). Aorta biopsy was selected, and the expression of MAP3K7CL was found to be reduced significantly by the CNV (p value = 0.014). In contrast, in atrial septum biopsy samples, the expression of MAP3K7CL increased about threefold in CHD patients compared to controls (Figure S1). These findings may imply tissue specificity for the CNV’s impact on gene expression and association with tissue-specific CHDs, which warrants further study.

The deletion CNV at chr2:109113426-109287320, corresponding to the LIMS1 and GCC2 genes, was detected in three CHD probands in the discovery dataset and replicated in two independent CHD probands in the WES replication dataset (permutation test p = 5.18 × 10−3). This region has been reported as a CHD-associated region (Glessner et al., 2014), and based on the Database of Genomic Variants (DGV), the region has been reported as a locus for developmental delay (Coe et al., 2014). LIMS1 is highly expressed in mice at day 14.5 and in human heart (Zaidi et al., 2013) and is associated with conotruncal heart malformations in multiple studies (Saraclar et al., 1996). Both LIMS1 and GCC2 show significant reduction of expression in the aorta tissue in CHD probands (subjects 1-02486 and 1-04724) who have this CNV, in comparison with probands (n = 7) who do not harbor this CNV (p value = 0.015 and 0.12).

The deletion CNV at chr2:109363161-109371723, corresponding to the gene RANBP2, existed in four CHD probands from the discovery dataset and was replicated in two independent individuals from the replication WES dataset without any occurrences in controls (permutation test p = 2.13 × 10−3). The deletion CNV is a de novo variation, and the RANBP2 gene is highly expressed in mouse at day 14.5 and in the human heart (Zaidi et al., 2013; Kim et al., 2012). Previous studies showed that the gene is the target for small ubiquitin-related modifier (SUMO) conjugation and could lead to CHDs and cardiac dysfunction in murine hearts (Kim et al., 2012; Sakin et al., 2015). The aorta tissue biopsy was collected for the CHD proband, 1-02846, who has the target CNV deletion and seven other CHD probands who do not have the CNV. The expression of RANBP2 was significantly reduced (p value = 0.05) while two housekeeping genes remained stable, indicating the expression alternations in RANBP2 as a consequence of the targeted CNV.

The deletion CNV at chr21:38461093-38523202, corresponding to the gene TTC3, existed in a single CHD proband, 1-04198, in the discovery dataset and was replicated in seven independent individuals from the WES replication dataset without any occurrences in controls. The deletion CNV overlaps with CNVs previously reported in CHD studies (Glessner et al., 2014) and associated with developmental delay (Coe et al., 2014). TTC3 is highly expressed in mouse at 14.5 days and in human heart tissue (Zaidi et al., 2013); a previous study showed that the gene correlates with cardiac defects in Down syndrome patients (Ripoll et al., 2012), and TTC3 is a novel candidate gene based on array comparative genomic hybridization (aCGH) for 316 nonsyndromic CHD children (Sanchez-Castro et al., 2016). The right ventricle and right atrium tissue biopsies were collected for CHD probands, 1-02727 and 1-03784, who have the target CNV deletions as well as for six other CHD probands who do not harbor the CNV. The expression of TTC3 was significantly reduced in both tissues while two housekeeping genes remained stable, indicating that the expression alternation is attributed to the targeted CNV. Mechanistic reasoning regarding the role of TTC3 in CHD warrants reconciliation with the observed increased gene dosage in Down syndrome and the deletion observed in our study.

Discussion

New CNV Loci Identified in This Study

CHD is the most common birth defect, affecting approximately 1% of newborns, and research has shown that genomic variations contribute significantly to CHD. Recent studies from large CHD populations demonstrate that CHD is associated with other birth conditions such as developmental delay (Homsy et al., 2015). CNVs have been reported previously in association with CHD where rare recurrent de novo CNVs are enriched in CHD cases (Glessner et al., 2014). Compared with the previous study where 65 de novo CNVs identified in CHD cases, our study identified 122 CNVs which are supported by at least one replication dataset, but not seen in controls We also performed a family-based data analysis from over 2,500 families with CHD probands to explore the impact of non-inherited CNVs in the pathogenesis of CHD. We used data from different platforms including WES, SNP arrays, WGS, and RNA sequencing. Four layers of filters based on independent replication sets were applied to remove CNVs unrelated to CHD, in order to reduce false-positive hits.

Beside the independent replications, we also performed evaluation based on various biological studies from existing knowledge bases, and functional studies at mRNA levels were used to lend further pathogenic support to variants of interest. We required target genes to be highly expressed in human heart, blood vessel tissues, or mouse embryo heart, an approach previously established in relation with CHD-associated CNVs (Glessner et al., 2014; Zaidi et al., 2013). In the filtered list of 39 CNVs (Table S2), 29 are novel CNVs while 10 CNVs overlapped with developmental delay regions. These results are consistent with previous reports which demonstrated association between CHD and developmental delay, as many of our targeted CNVs are associated with developmental delays and other developmental disease loci, which brings more confidence to the results.

The targeted CNVs were further validated through cardiac biopsies obtained from individuals harboring the targeted CNVs. Indeed, the corresponding genes impacted by targeted CNVs showed significantly altered gene expression in comparison with individuals without targeted CNVs, while the housekeeping genes are stable (Table 3, Figure 2). Ideally, affected probands with CNVs would be compared to unaffected controls; however, due to the natural difficulties in obtaining healthy cardiac tissues, we had to compare the gene expressions among CHD probands. These results demonstrate that the targeted genomic CNV does impact gene expression and are potential pathogenic culprits in CHD.

RANBP2: A Link Between Neurodevelopmental Disorders (NDDs) and CHD

Among recurrent de novo CNV examples, a deletion at chr2:109363161-109371723, which deletes exons 8 to 17 of the RANBP2 gene, was observed recurrently in four independent probands, but not in any of the control samples including their parents or siblings. RANBP2 is a small guanosine-5′-triphosphate (GTP)-binding protein of the RAS superfamily that is associated with the nuclear membrane and is thought to control a variety of cellular functions through its interactions with other proteins. The encoded protein directly interacts with the E2 enzyme UBC9 and strongly enhances SUMO1 transfer function. A previous study showed that both heterozygous and homozygous SUMO1 knockout mice exhibit atrial septal defect (ASD)/ventricular septal defect (atrial septal defect) and suffer from high mortality rates; this was rescued by cardiac re-expression of the SUMO1 transgene (Wang et al., 2011). Our results indicate that a large deletion of RANBP2 results in the reduction of the SUMO1 transfer function, which may explain the association with CHD in humans. The CNV region was further validated in two more CHD probands from an independent trio data. Aorta tissues were also obtained from one of the proband and from several CHD cases without this CNV deletion. As show in Figure 2B, the expression of RANBP2 (blue bar) is significantly reduced in the CNV probands, while the housekeeping genes G6PD and ACTB remain stable in all subjects tested. Thus, it is exceedingly likely that the deletion is responsible for the gene expression reduction and could further impact the SUMO1 transfer pathway in CHD cases with a high mortality rate.

NDD is commonly seen in CHD patients (Marino et al., 2012). Shared genetic etiology of NDD and CHD has been identified in genes involved in morphogenesis, chromatin modification, and transcriptional regulation (Homsy et al., 2015). As shown previously, RanBP2 works as a chaperone with the mitochondrial metallochaperone Cox11 (encoded by the cytochrome c oxidase copper chaperone COX11 gene) and plays critical roles in the modulation of neuronal hexokinase type I (HKI), which is the pacemaker of glycolysis (Aslanukov et al., 2006). Cox11 inhibits HKI activity, and RanBP2 suppresses this inhibition. Haploinsufficiency of RANBP2 in mice induces the downregulation of HKI and ATP levels selectively in the central nervous system (Aslanukov et al., 2006). Targeting at this mechanism, e.g., to activate HKI, may offer new opportunities to impact NDD in CHD patients.

LIMS1: A Link Between Innate Immunity and CHD

Besides RANBP2, the other CNV locus at chromosome 2 containing the gene LIMS1 may provide important insight into innate immunity and CHD, based on its significant change of gene expression and previous studies of its role in CHD (Glessner et al., 2014; Saraclar et al., 1996). As shown in this study, the CNV locus at chr2:109113426-109287320 was associated with a significantly decreased expression of LIMS1. In mouse models, a previous study has shown that the LIMS1 protein forms a functional complex with thymosin β4 and integrin-linked kinase, which plays an important role in cardiomyocyte migration, survival, and repair (Bock-Marquette et al., 2004). In humans, LIMS1 polymorphism has been identified in association with monocyte chemotactic protein-1 (MCP1) levels (Ahola-Olli et al., 2017). It is one of the key chemokines that regulate migration and infiltration of monocytes/macrophages (Deshmane et al., 2009), while the critical role of macrophages in neonatal heart regeneration has been demonstrated (Aurora et al., 2014). A previous study on rare CNV in congenital left-sided heart disease also identified a paternally inherited loss of LIMS1 in a case with partially anomalous pulmonary venous return and ligation of patent ductus arteriosus (Hitz et al., 2012). In addition, considering the role of macrophages in resistance to virus infections, LIMS1 might also represent a link of intrauterine virus infections and congenital heart disease.

Conclusions

CNVs have been previously shown to contribute towards CHDs. However, limited sample sizes have prevented these studies from discovering recurrent CNVs in genes contributing towards CHD. To explore the impact of recurrent rare CNVs, we performed the largest family-based CHD structural variation genomic study to date, identifying targets we replicated using different technical platforms based on independent data, and the CNV target regions were further validated through CHD-related knowledge and cardiac tissue gene expression alternations. The new insights into CHD provided by this study offer new opportunities to clarify the development of CHD and its comorbidities.

This study has limitations. For CNVs’ impact on gene expression, we showed unadjusted p values in this paper, which lose significance after correction for multiple testing due to the difficulty in obtaining the statistical power. It is difficult to acquire healthy human heart tissues, and there are limited RNA-Seq data available in public databases such as dbGaP. We could only compare the CNV-impacted gene expression between CHD patients with target CNVs and CHD patients without target CNVs. The tendency of gene expression changes associated with the existence of target CNVs shown in this study warrants validation in future studies. As another limitation of our study, we had only accession to anonymous data of the CHD patient collection, without detailed phenotyping information about the clinical characteristics of CHD. We acknowledge the importance of a clinical investigation on genotype–phenotype analysis for the rare CNVs identified in our study, in the near future, working with the PCGC.

Data Availability

Data can be accessed by the dbGaP accession phs001194.v2.p2 and phs001194.v2.p2.c1.

Ethics Statement

The protocol was approved by the Institutional Review Boards of Boston Children’s Hospital, Brigham and Women’s Hospital, Great Ormond St. Hospital, Children’s Hospital of Los Angeles, Children’s Hospital of Philadelphia, Columbia University Medical Center, Icahn School of Medicine and Mt. Sinai, Rochester School of Medicine and Dentistry, Steven and Alexandra Cohen Children’s Medical Center of New York, and Yale School of Medicine. Written informed consent was obtained from each participating subject or their parent/guardian.

Author Contributions

HH and PS conceived and supervised the project. YL, XC, HQ, and JG designed and implemented the methods. LT, DL, and KN contributed to data acquisition and analysis. YL, HQ, XC, JG, and HH wrote the manuscript. All authors approved the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We thank the CAG staff, including grant support from CHOP and eMERGE as well as dbGaP for data accession phs001194.v2.p2 and phs001194.v2.p2.c1. These data were generated by the Pediatric Cardiac Genomics Consortium (PCGC), under the auspices of the National Heart, Lung, and Blood Institute’s Bench to Bassinet Program (http://www.benchtobassinet.org/). The PCGC program is funded by the National Heart, Lung, and Blood Institute; National Institutes of Health; and U.S. Department of Health and Human Services through grants U01HL098123, U01HL098147, U01HL098153, U01HL098162, U01HL098163, and U01HL098188.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2019.00819/full#supplementary-material

Figure S1 | Corresponding gene expression for target copy number variants (CNVs) that do not have significant alternations. Similar to Figure 2 except the target CNV regions in this plot do not have gene expression alteration significantly while housekeeping genes are stable. For pulmonary valve leaflets tissue, since ACTB is not expressed, B2M was selected as a housekeeping gene.

Table S1 | CNV targets that not occurred in any health controls, and supported by at least one evidence from previous knowledge, biological supports or replication set.

Table S2 | Selected CNV targets based on replication set, previous knowledge and biological supports.

Table S3 | De novo and novel CNV targets identified in the discovery set.

References

Ahola-Olli, A. V., Wurtz, P., Havulinna, A. S., Aalto, K., Pitkanen, N., Lehtimaki, T., et al. (2017). Genome-wide association study identifies 27 loci influencing concentrations of circulating cytokines and growth factors. Am. J. Hum. Genet. 100 (1), 40–50. doi: 10.1016/j.ajhg.2016.11.007

PubMed Abstract | CrossRef Full Text | Google Scholar

Aragam, K. G., Chaffin, M., Levinson, R. T., McDermott, G., Choi, S. H., Shoemaker, M. B., et al. (2018). Phenotypic refinement of heart failure in a national biobank facilitates genetic discovery. Circulation. 139(4), 489–501 doi: 10.1161/CIRCULATIONAHA.119.040940

CrossRef Full Text | Google Scholar

Aslanukov, A., Bhowmick, R., Guruju, M., Oswald, J., Raz, D., Bush, R. A., et al. (2006). RanBP2 Modulates Cox11 and hexokinase I activities and haploinsufficiency of RanBP2 causes deficits in glucose metabolism. PLoS Genet. 2 (10), e177. doi: 10.1371/journal.pgen.0020177

PubMed Abstract | CrossRef Full Text | Google Scholar

Aurora, A. B., Porrello, E. R., Tan, W., Mahmoud, A. I., Hill, J. A., Bassel-Duby, R., et al. (2014). Macrophages are required for neonatal heart regeneration. J. Clin. Investig. 124 (3), 1382–1392. doi: 10.1172/JCI72181

CrossRef Full Text | Google Scholar

Bock-Marquette, I., Saxena, A., White, M. D., Dimaio, J. M., Srivastava, D. (2004). Thymosin beta4 activates integrin-linked kinase and promotes cardiac cell migration, survival and cardiac repair. Nature 432 (7016), 466–472. doi: 10.1038/nature03000

PubMed Abstract | CrossRef Full Text | Google Scholar

Coe, B. P., Witherspoon, K., Rosenfeld, J. A., van Bon, B. W., Vulto-van Silfhout, A. T., Bosco, P., et al. (2014). Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat. Genet. 46 (10), 1063–1071. doi: 10.1038/ng.3092

PubMed Abstract | CrossRef Full Text | Google Scholar

D O’Fallon, B., Durstchi, J., Lewis, T., Close, D. (2018). Algorithmic improvements for discovery of germline copy number variants in next-generation sequencing data. bioRxiv. doi: 10.1101/441378

CrossRef Full Text | Google Scholar

Deshmane, S. L., Kremlev, S., Amini, S., Sawaya, B. E. (2009). Monocyte chemoattractant protein-1 (MCP-1): an overview. J Interferon Cytokine Res. 29 (6), 313–326. doi: 10.1089/jir.2008.0027

PubMed Abstract | CrossRef Full Text | Google Scholar

Fakhro, K. A., Choi, M., Ware, S. M., Belmont, J. W., Towbin, J. A., Lifton, R. P., et al. (2011). Rare copy number variations in congenital heart disease patients identify unique genes in left-right patterning. Proc. Natl. Acad. Sci. U S A 108 (7), 2915–2920. doi: 10.1073/pnas.1019645108

PubMed Abstract | CrossRef Full Text | Google Scholar

Fan, X., Abbott, T. E., Larson, D., Chen, K. (2014). BreakDancer: identification of genomic structural variation from paired-end read mapping. Curr. Protoc. Bioinform. 45 (15.6), 1–11. doi: 10.1002/0471250953.bi1506s45

CrossRef Full Text | Google Scholar

Fromer, M., Moran, J. L., Chambert, K., Banks, E., Bergen, S. E., Ruderfer, D. M., et al. (2012). Discovery and statistical genotyping of copy-number variation from whole-exome sequencing depth. Am. J. Hum Genet. 91 (4), 597–607. doi: 10.1016/j.ajhg.2012.08.005

PubMed Abstract | CrossRef Full Text | Google Scholar

Glessner, J. T., Bick, A. G., Ito, K., Homsy, J. G., Rodriguez-Murillo, L., Fromer, M., et al. (2014). 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. 115 (10), 884–896. doi: 10.1161/CIRCRESAHA.115.304458

PubMed Abstract | CrossRef Full Text | Google Scholar

Greenway, S. C., Pereira, A. C., Lin, J. C., DePalma, S. R., Israel, S. J., Mesquita, S. M., et al. (2009). De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot. Nat. Genet. 41 (8), 931–935. doi: 10.1038/ng.415

PubMed Abstract | CrossRef Full Text | Google Scholar

Hitz, M. P., Lemieux-Perreault, L. P., Marshall, C., Feroz-Zada, Y., Davies, R., Yang, S. W., et al. (2012). Rare copy number variants contribute to congenital left-sided heart disease. PLoS Genet. 8 (9), e1002903. doi: 10.1371/journal.pgen.1002903

PubMed Abstract | CrossRef Full Text | Google Scholar

Hoffman, J. I., Kaplan, S. (2002). The incidence of congenital heart disease. J. Am. Coll. Cardiol. 39 (12), 1890–1900. doi: 10.1016/S0735-1097(02)01886-7

PubMed Abstract | CrossRef Full Text | Google Scholar

Homsy, J., Zaidi, S., Shen, Y., Ware, J. S., Samocha, K. E., Karczewski, K. J., et al. (2015). De novo mutations in congenital heart disease with neurodevelopmental and other congenital anomalies. Science (New York, NY) 350 (6265), 1262–1266. doi: 10.1126/science.aac9396

CrossRef Full Text | Google Scholar

Kim, E. Y., Chen, L., Ma, Y., Yu, W., Chang, J., Moskowitz, I. P., et al. (2012). Enhanced desumoylation in murine hearts by overexpressed SENP2 leads to congenital heart defects and cardiac dysfunction. J. Mol. Cell. Cardiol. 52 (3), 638–649. doi: 10.1016/j.yjmcc.2011.11.011

PubMed Abstract | CrossRef Full Text | Google Scholar

MacDonald, J. R., Ziman, R., Yuen, R. K., Feuk, L., Scherer, S. W. (2014). The database of genomic variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 42 (Database issue), D986–D992. doi: 10.1093/nar/gkt958

PubMed Abstract | CrossRef Full Text | Google Scholar

Marino, B. S., Lipkin, P. H., Newburger, J. W., Peacock, G., Gerdes, M., Gaynor, J. W., et al. (2012). Neurodevelopmental outcomes in children with congenital heart disease: evaluation and management: a scientific statement from the American Heart Association. Circulation 126 (9), 1143–1172. doi: 10.1161/CIR.0b013e318265ee8a

PubMed Abstract | CrossRef Full Text | Google Scholar

Miller, C. L., Pjanic, M., Wang, T., Nguyen, T., Cohain, A., Lee, J. D., et al. (2016). Integrative functional genomics identifies regulatory mechanisms at coronary artery disease loci. Nat. Commun. 7, 12092. doi: 10.1038/ncomms12092

PubMed Abstract | CrossRef Full Text | Google Scholar

Pediatric Cardiac Genomics C, Gelb, B., Brueckner, M., Chung, W., Goldmuntz, E., Kaltman, J., et al. (2013). The congenital heart disease genetic network study: rationale, design, and early results. Circ. Res. 112 (4), 698–706. doi: 10.1161/CIRCRESAHA.111.300297

PubMed Abstract | CrossRef Full Text | Google Scholar

Ripoll, C., Rivals, I., Yahya-Graison, E., Dauphinot, L., Paly, E., Mircher, C., et al. (2012). Molecular signatures of cardiac defects in Down syndrome lymphoblastoid cell lines suggest altered ciliome and hedgehog pathways. PLoS One 7 (8), e41616. doi: 10.1371/journal.pone.0041616

PubMed Abstract | CrossRef Full Text | Google Scholar

Sadedin, S. P., Oshlack, A., Ellis, J. A., Masters, S. L. (2018). Ximmer: a system for improving accuracy and consistency of CNV calling from exome data. GigaScience. 7 (10), 1–11. doi: 10.1093/gigascience/giy112

CrossRef Full Text | Google Scholar

Sakin, V., Richter, S. M., Hsiao, H. H., Urlaub, H., Melchior, F. (2015). Sumoylation of the GTPase ran by the RanBP2 SUMO E3 ligase complex. J. Biol. Chem. 290 (39), 23589–23602. doi: 10.1074/jbc.M115.660118

PubMed Abstract | CrossRef Full Text | Google Scholar

Sanchez-Castro, M., Eldjouzi, H., Charpentier, E., Busson, P. F., Hauet, Q., Lindenbaum, P., et al. (2016). Search for rare copy-number variants in congenital heart defects identifies novel candidate genes and a potential role for FOXC1 in patients with coarctation of the aorta. Circ. Cardiovasc. Genet. 9 (1), 86–94. doi: 10.1161/CIRCGENETICS.115.001213

PubMed Abstract | CrossRef Full Text | Google Scholar

Saraclar, M., Cil, E., Ozkutlu, S. (1996). Echocardiography for the diagnosis of congenital cardiac anomalies with multiple lesions. Pediatr. Cardiol. 17 (5), 308–313. doi: 10.1007/s002469900067

PubMed Abstract | CrossRef Full Text | Google Scholar

Sifrim, A., Hitz, M. P., Wilsdon, A., Breckpot, J., Turki, S. H., Thienpont, B., et al. (2016). Distinct genetic architectures for syndromic and nonsyndromic congenital heart defects identified by exome sequencing. Nat. Genet. 48 (9), 1060–1065. doi: 10.1038/ng.3627

PubMed Abstract | CrossRef Full Text | Google Scholar

Silversides, C. K., Lionel, A. C., Costain, G., Merico, D., Migita, O., Liu, B., et al. (2012). Rare copy number variations in adults with tetralogy of Fallot implicate novel risk gene pathways. PLoS Genet. 8 (8), e1002843 doi: 10.1371/journal.pgen.1002843.

CrossRef Full Text | Google Scholar

Soemedi, R., Wilson, I. J., Bentham, J., Darlay, R., Topf, A., Zelenika, D., et al. (2012). Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease. Am. J. Hum. Genet. 91 (3), 489–501. doi: 10.1016/j.ajhg.2012.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Stankiewicz, P., Lupski, J. R. (2010). Structural variation in the human genome and its role in disease. Annu. Rev Med. 61, 437–455. doi: 10.1146/annurev-med-100708-204735

PubMed Abstract | CrossRef Full Text | Google Scholar

Trapnell, C., Williams, B. A., Pertea, G., Mortazavi, A., Kwan, G., van Baren, M. J., et al. (2010). Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. Nat. Biotechnol. 28 (5), 511–515. doi: 10.1038/nbt.1621

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Bom, T., Zomer, A. C., Zwinderman, A. H., Meijboom, F. J., Bouma, B. J., Mulder, B. J. (2011). The changing epidemiology of congenital heart disease. Nat. Rev. Cardiol. 8 (1), 50–60. doi: 10.1038/nrcardio.2010.166

PubMed Abstract | CrossRef Full Text | Google Scholar

van der Linde, D., Konings, E. E., Slager, M. A., Witsenburg, M., Helbing, W. A., Takkenberg, J. J., et al. (2011). Birth prevalence of congenital heart disease worldwide: a systematic review and meta-analysis. J. Am. Coll. Cardiol. 58 (21), 2241–2247. doi: 10.1016/j.jacc.2011.08.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, K., Chen, Z., Tadesse, M. G., Glessner, J., Grant, S. F., Hakonarson, H., et al. (2008). Modeling genetic inheritance of copy number variations. Nucleic Acids Res. 36 (21), e138. doi: 10.1093/nar/gkn641

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Chen, L., Wen, S., Zhu, H., Yu, W., Moskowitz, I. P., et al. (2011). Defective sumoylation pathway directs congenital heart disease. Birth Defects Res. A. Clin. Mol. Teratol. 91 (6), 468–476. doi: 10.1002/bdra.20816

PubMed Abstract | CrossRef Full Text | Google Scholar

Zaidi, S., Choi, M., Wakimoto, H., Ma, L., Jiang, J., Overton, J. D., et al. (2013). De novo mutations in histone-modifying genes in congenital heart disease. Nature 498 (7453), 220–223. doi: 10.1038/nature12141

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: copy number variants, next-generation sequencing, congenital heart defects, large trios study, genomics

Citation: Liu Y, Chang X, Glessner J, Qu H, Tian L, Li D, Nguyen K, Sleiman PMA and Hakonarson H (2019) Association of Rare Recurrent Copy Number Variants With Congenital Heart Defects Based on Next-Generation Sequencing Data From Family Trios. Front. Genet. 10:819. doi: 10.3389/fgene.2019.00819

Received: 12 June 2019; Accepted: 07 August 2019;
Published: 10 September 2019.

Edited by:

Andrew Landstrom, Duke University, United States

Reviewed by:

Feng Zhang, Fudan University, China
Sheng Wang, University of California, United States

Copyright © 2019 Liu, Chang, Glessner, Qu, Tian, Li, Nguyen, Sleiman and Hakonarson. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Hakon Hakonarson, hakonarson@email.chop.edu

These authors have contributed equally to this work.