De novo large rare copy-number variations contribute to conotruncal heart disease in Chinese patients

Abstract

Conotruncal heart anomalies (CTDs) are particularly prevalent congenital heart diseases (CHD) in Hong Kong. We surveyed large (>500 kb), rare (<1% frequency in controls) copy-number variations (CNVs) in Chinese patients with CTDs to identify potentially disease-causing variations. Adults who tested negative for 22q11.2 deletions were recruited from the adult CHD clinic in Hong Kong. Using a stringent calling criteria, high-confidence CNV calls were obtained, and a large control set comprising 3,987 Caucasian and 1,945 Singapore Chinese subjects was used to identify rare CNVs. Ten large rare CNVs were identified, and 3 in 108 individuals were confirmed to harbour de novo CNVs. All three patients were syndromic with a more complex phenotype, and each of these CNVs overlapped regions likely to be important in CHD. One was a 611 kb deletion at 17p13.3, telomeric to the Miller-Dieker syndrome (MDS) critical region, overlapping the NXN gene. Another was a 5 Mb deletion at 13q33.3, within a previously described critical region for CHD. A third CNV, previously unreported, was a large duplication at 2q22.3 overlapping the ZEB2 gene. The commonly reported 1q21.1 recurrent duplication was not observed in this Chinese cohort. We provide detailed phenotypic and genotypic descriptions of large rare genic CNVs that may represent CHD loci in the East Asian population. Larger samples of Chinese origin will be required to determine whether the genome-wide distribution differs from that found in predominantly European CHD cohorts.

http://www.nature.com/npjgenmed

Web End =www.nature.com/npjgenmed

ARTICLE OPEN

De novo large rare copy-number variations contribute to conotruncal heart disease in Chinese patients

Christopher CY Mak1, Pak Cheong Chow1, Anthony PY Liu1, Kelvin YK Chan2, Yoyo WY Chu1, Gary TK Mok1, Gordon KC Leung1, Kit San Yeung1, Adolphus KT Chau1, Chelsea Lowther3, Stephen W Scherer4, Christian R Marshall4, Anne S Bassett3and Brian HY Chung1

Conotruncal heart anomalies (CTDs) are particularly prevalent congenital heart diseases (CHD) in Hong Kong. We surveyed large (4500 kb), rare (o1% frequency in controls) copy-number variations (CNVs) in Chinese patients with CTDs to identify potentially disease-causing variations. Adults who tested negative for 22q11.2 deletions were recruited from the adult CHD clinic in Hong Kong. Using a stringent calling criteria, high-condence CNV calls were obtained, and a large control set comprising 3,987 Caucasian and 1,945 Singapore Chinese subjects was used to identify rare CNVs. Ten large rare CNVs were identied, and 3 in 108 individuals were conrmed to harbour de novo CNVs. All three patients were syndromic with a more complex phenotype, and each of these CNVs overlapped regions likely to be important in CHD. One was a 611 kb deletion at 17p13.3, telomeric to the MillerDieker syndrome (MDS) critical region, overlapping the NXN gene. Another was a 5 Mb deletion at 13q33.3, within a previously described critical region for CHD. A third CNV, previously unreported, was a large duplication at 2q22.3 overlapping the ZEB2 gene. The commonly reported 1q21.1 recurrent duplication was not observed in this Chinese cohort. We provide detailed phenotypic and genotypic descriptions of large rare genic CNVs that may represent CHD loci in the East Asian population. Larger samples of Chinese origin will be required to determine whether the genome-wide distribution differs from that found in

predominantly European CHD cohorts.

npj Genomic Medicine (2016) 1, 16033; doi:http://dx.doi.org/10.1038/npjgenmed.2016.33

Web End =10.1038/npjgenmed.2016.33 ; published online 14 September 2016

INTRODUCTIONCongenital heart diseases (CHDs) are the most common birth defect and a major cause of morbidity and mortality.1 Conotruncal heart defects (CTDs) are CHDs affecting the cardiac outow tract. These include tetralogy of Fallot (ToF), pulmonary atresia with VSD (PAVSD), truncus arteriosus, interrupted aortic arch, transposition of the great arteries and double outlet right ventricle. Overall, CTDs have an estimated prevalence of 11.6 per 10,000 live births.2 Figures available for ToF, the most common CTD, indicate a prevalence of 2.7 per 10,000 live births in Europe2 and 4.7 per 10,000 live births in the United States.3 In Taiwan, the prevalence has been reported to be 6.26 per 10,000 live births.4 Although Hong Kong does not have the prevalence data for CTD, in terms of the proportion of all CHD, pulmonary outow tract obstruction alone accounts for 31.1% of all symptomatic Chinese neonates with CHD in Hong Kong, almost double the proportion reported in western studies,5 implying that this group of congenital abnormalities may have greater representation in Asian populations.

Although it is thought that genetic factors have an important role in CHD, only about 11% of patients receive a genetic diagnosis.1 In addition to chromosomal abnormalities and rare single-nucleotide variants, a growing proportion of the molecular diagnosis is attributed to rare copy-number variations (CNVs). Multiple studies show an increased burden of rare CNV

in patients affected with ToF compared with controls, and several recurrent loci have been reported to be pathogenic for ToF, the most notable being the deletion responsible for 22q11.2 deletion syndrome.612 Nevertheless, Chinese studies are lacking.

Recognising the importance of CTDs in Chinese patients, and the wide spectrum of CNVs contributing to CTDs in other cohorts, our aim was to discover important CNVs in a Chinese cohort, focusing on the large (4500 kb) CNVs, readily detectable on clinical microarray platforms. We hypothesised that large rare structural variants that overlap the coding sequence of one or more genes would be associated with CTDs, over and above the recurrent 22q11.2 deletion we previously showed to be prevalent in this cohort.8 We present case examples characterising the lifetime phenotype of patients harbouring CNVs that are likely to be the most clinically signicant.

RESULTSPatient demographicsThe 116 patients (62 males, 54 females) with CTD recruited from the adult CHD clinic had a mean age s.d. of 31.0 8.4 years that showed a normal distribution (range 1753 years, skewness 0.74, kurtosis 0.031). All of the patients recruited with CTD had either ToF (81.9%, n = 95) or PAVSD (18.1%, n = 21).

1Department of Paediatrics & Adolescent Medicine, LKS Faculty of Medicine, The University of Hong Kong, Hong Kong, China; 2Department of Obstetrics and Gynecology, Queen Mary Hospital, Hong Kong, China; 3The Clinical Genetics Research Program at The Centre for Addiction and Mental Health, The Dalglish Family 22q Clinic at The University Health Network, and The Department of Psychiatry at The University of Toronto, Toronto, ON, Canada and 4The Centre for Applied Genomics and Program in Genetics and Genome Biology, The Hospital for Sick Children, Toronto, Ontario, Canada.

Correspondence: BHY Chung (mailto:[email protected]

Web End [email protected])

Received 14 April 2016; revised 2 August 2016; accepted 5 August 2016

Published in partnership with the Center of Excellence in Genomic Medicine Research

Conotruncal heart anomalies

CCY Mak et al

2

Large rare CNVsUsing stringent criteria and methods,10 rare CNV calls were obtained for 108 (93.2%) of the 116 subjects and therefore these subjects were used as the CNV discovery cohort. Filtering for only CNVs found in o1% of platform-matched controls, a total of 306 rare autosomal CNVs were identied in these 108 subjects. The subsequent analysis focused on CNVs that were large (dened as 4500 kb), rare (o1% in all platform-matched and ethnicity-matched control groups, total = 5,902), and overlapped the exons of genes. Ten large CNVs from ten patients met these criteria (Table 1). Of note, the sex distribution of these 10 patients was skewed, with 7 females and 3 males, although the overall cohort had an excess of male patients. None of the large CNVs identied overlapped the smaller (o500 kb) rare CNVs found in the sample of 108 patients. The clinical features of these 10 patients are summarised in Table 2.

Clinical cases and candidate genes of signicant large rare CNVs Of these 10 CNVs, three were de novo large rare CNVs found in syndromic patients. First is a large duplication at 2q22.3 (Patient 1), the second, a 5 Mb deletion at 13q33.3 (Patient 9) and the third, a 611 kb deletion at 17p13.3 (Patient 10). Below we discuss in detail the phenotypic and molecular features of these three patients.

2q22.3 duplication overlapping the ZEB2 genePatient 1 (female, age 32 years) was born with ToF and patent foramen ovale. She had mild intellectual disability and adjustment disorder. She also has a history of dysfunctional uterine bleeding. On examination, she was dysmorphic, but did not t the facial features of MowatWilson syndrome. The most remarkable observation was that she had a hypernasal voice (Figure 1a).

A de novo 2.1 Mb copy-number gain was found at 2q22.3, which overlaps the ZEB2 gene. Deletions of this gene are a known cause of MowatWilson syndrome (OMIM#235730) in which conotruncal heart disease have been described.13,14 However, there have been no known reports of a duplication leading to the MowatWilson syndrome phenotype. Such a large gain in copy-number may disrupt regulatory elements, such as the expression levels of ZEB2, which would subsequently affect transcription downstream.

The ZEB2 gene (previously known as ZFHX1B, or SIP1) is widely expressed in human embryological development, including in the heart.15 The protein encoded by this gene, is SMAD-interacting protein-1 (ref. 16), which acts to activate or repress transcription by binding to regulatory sequences of E-boxes.17 SMAD proteins are present in the cytoplasm to mediate transforming growth factor (TGF-) signals from receptors of the cell-surface into the nucleus,18 and hence the TGF- group of cytokines are likely to have an important role in embryological development,18 including cardiac outow tract formation. Loss of ZEB2 in the mouse is associated with variable heart defects,19 similar to the cardiac anomalies with deletions of this gene in human MowatWilson syndrome.

A 13q33.3 deletion in a critical region for CHDPatient 9 (female, age 21 years) has a complex cardiac anatomy of dextrocardia, PAVSD, and patent ductus arteriosis. She was born cyanotic and small for dates, with a birth weight of 2.5 kg. Neurological assessment up to 8 months was apparently normal, however, at age 15 months, motor and speech delays were noted. She also had a brain abscess complicated by focal seizures at age 21 months. The head circumference was below the 3rd centile even before the occurrence of the brain abscess, however, and remained at the same percentile after the resolution of this infection. Subsequently, there was overall slow progress in intellectual function and the patient had attended a special

Huangetal.20 ,

McMahonetal.25 ,Yang

etal.21

Silversidesetal.10

PutativecandidategenesReferences

COL4A2,COL4A1,FAM70B,F10,ING1,RASA3,RAB20,

SPACA7,CHAMP1,PROZ,ATP4B,F7,PCID2,MCF2L,IRS2,

CUL4A,UPF3A,TFDP1,TUBGCP3,CARKD,GAS6

CNVeryrarea ValidatedInheritanceNo.of

genes

Table1.LargeraregenicCNVsidentiedinadultChinesepatientswithconotruncaldefects(listedaccordingtocytoband)

2q22.311442234161462336012010185GainYYesDenovo6ARHGAP15,GTDC1,ZEB2

5q33.12151540888152162921622034GainYYesMaternal1NMUR2

6p21.3133442235434943632521279GainYYesMaternal7SPDEF,TAF11,UHRF1BP1,SNRPC,ANKS1A,PACSIN1

7p21.141689361017502331608722GainYYesNA2AHR,AGR3

7q335133760551134367872607322LossYYesNA6LRGUK,AKR1B10,AKR1B15,BPGM,SLC35B4,AKR1B1

16q23.1676070911770871231016213GainYYesNA2CNTNAP4,MIR4719

12q24.327127967297128566111598815GainYYesNA1FLJ37505

13q33.2,13q33.381060409851075812241540240GainYYesNA5EFNB2,ARGLU1,DAOADavyetal.26

13q33.3,13q3491101628111151083974945587LossYYesDenovo41LAMP1,TMCO3,GRTP1,CARS2,ATP11A,SOX1,ANKRD10,

CDC16,ARHGEF7,GRK1AS1,ADPRHL1,DCUN1D2,

17p13.3105007461112604611859Losso0.1%YesDenovo10RNMTL1,FAM57A,GLOD4,GEMIN4,DBIL5P,ABR,VPS53,

NXN,MIR3183,TIMM22

Abbreviations:CN,copy-numberchange;CNV,copy-numbervariants;NA,notavailable;Y,notfoundatallinanyofthe3,957platformcontrol.

Forallavailableparentalsamples,inheritanceofCNVwasdeterminedandshown.

CytobandPatientStartEndCNVSize

(bp)

Candidategeneswithpotentialimportanceincardiacdevelopmentarehighlightedinboldfont.

Allgenomiccoordinatesareprovidedusinggenomebuildhg19.

a RarityforGroup1controls(AllCNVsarefoundino1%ofgroup2Chinesecontrols)indicatingveryrareCNVs.

npj Genomic Medicine (2016) 16033 Published in partnership with the Center of Excellence in Genomic Medicine Research

Conotruncal heart anomalies CCY Mak et al

3

Othermedicalissues

Abbreviations:CHD,congenitalheartdiseases;CNVs,copy-numbervariations;F,female;M,male;PAVSD,PulmonaryatresiawithVSD;TOF,tetralogyofFallot;VSD,ventricularseptaldefect.

school. On examination, she had mild facial dysmorphic features with low set ears and hypernasal speech (no clinical photos).

In this patient, a large de novo 4.9 Mb deletion was found at 13q33.3, detectable by G-banding karyotype performed after the microarray. The CNV overlaps 41 genes, some notable among these were TFDP1, GAS6, COL4A1, COL4A2 and SOX1 (Figure 2). Given the size, and the known critical region, this CNV was clinically classied as pathogenic.

The 13q33.2-33.3 region has previously been proposed to be a critical region for CHD20,21 (Figure 2). Huang et al. summarised deletions with different CHD phenotypes such as ToF, ostium secundum,22 coarctation of the aorta,23 interauricular communication,23 pulmonary valve stenosis,24 patent ductus arteriosus,24 ASD/VSD20 and single atrium.20 Looking at conotruncal defects only, three cases of ToF have so far been reported to have very large deletions that overlap this critical region.22,23

McMahon et al. also reported a case of double outlet right ventricle within this proposed critical region with a 2.5 Mb deletion at 13q33.3-33.4 (Figure 3).25 They speculated that the collagen genes COL4A1 and COL4A2 are candidate genes for CHD as they are likely to have a role in cardiac development. The CNV found in our patient falls within the proposed critical region, and has a complex phenotype of PAVSD, dextrocardia, patent ductus arteriosus and situs solitus. It is interesting that within the small region of overlap with the CNV reported by McMahon et al. are the genes COL4A1 and COL4A2, which agrees with their hypothesis. However, the 13q33.3 deletion of Patient 9 also overlaps a smaller1.1 Mb deletion with single atrium reported by Yang et al. that did not include either of these two collagen genes.21 Of note, in our cohort, there was another large copy-number gain at 13q33.2-33.3 (patient 8), overlapping ve genes (EFNB2, ARGLU1, DAOA-AS1, DAOA and LINC00460), further suggesting the potential importance of this region and ephrin-B2 signalling in CHD.26

Summarising all the evidence so far, we posit that deletions of 13q33.3 are strongly associated with CHD but that narrow critical regions, or candidate genes, remain to be conrmed with more cases.

17p13.3 deletion telomeric to the critical region of MillerDieker Syndrome (MDS)

Patient 10 (Female, age 27 years) is known to have PAVSD. She required multiple surgeries, including a left modied Blalock-Taussig (LmBT) shunt for left pulmonary artery stenosis at age 7 months. Subsequently, closure of VSD and right ventricular outow tract reconstruction was performed at the age of 7 years. She is currently on warfarin due to atrial brillation. She had a history of epilepsy since age 17 years, initially presenting with petit mal seizures with normal EEG, controlled with sodium valproate. She studied in a special school because of mild intellectual disability. Her other medical history includes scoliosis, left hip dysplasia requiring total hip replacement, severe myopia and amblyopia, as well as primary amenorrhoea (diagnosed at age 17 years). On examination, she had dysmorphic facial features (Figure 1b).

The CNV analysis revealed a de novo 611 kb loss at 17p13.3 overlapping 10 genes, including NXN (Figure 3). The NXN gene is expressed in murine heart and there is evidence that loss of the NXN gene is associated with CHD in mice.27 In humans, the gene is expressed in human heart and has a haploinsufciency prediction score of 37.4%,28 as well as RVIS (Residual Variation Intolerance Score)29 percentile of 11.77% (i.e., the top 11.77% most intolerant of human genes to variation). The role of NXN in the canonical Wnt/-catenin signalling pathway through regulation of the dishevelled protein and the genes subsequent role in second heart eld signalling, further suggest this genes potential involvement in cardiac development.

MildmotordelayBrainabscess,VURwithurinary

incontinence,bilateralgenuvalgum

andmildexternaltibialtorsion

1017p13.3loss612F27PAVSDDysmorphicMildintellectualdisability,

epilepsy

LeftHipdysplasia,Severemyopiaand

amblyopia,scoliosis,primary

amenorrhoea

Eczema

AdjustmentdisorderNormal

25q33.1gain622M23TOFNoNormalNormal

36p21.31gain521M22TOFNoNormalNormal

47p21.1gain609F32TOFNoNormalNormal

57q33loss607M25TOFNoNormalNormal

616q23.1gain1,016F25TOFNoMildintellectualdisabilityPrecociouspuberty

712q24.32gain599F22TOFNoNormalNormal

813q33.2-

development

amenorrhoea

developmentaldelay

PatientCNVSize(kb)SexAgeTypeofCHDDysmorphicHypernasalityNeuropsychiatricmanifestationsGrowthand

DysmorphicHypernasal

speech

1,540F27PAVSDNoNormalSecondary

Global

Mildintellectualdisability,

poorperformancein

mathematics

DysmorphicHypernasal

speech

Table2.Clinicalfeaturesofconotruncalanomalypatientswithlarge(4500kb)raregenicCNVs

12q22.3gain1,846F32TOF,patentforamen

ovale

495F21PAVSD,dextrocardia,

situssolitus

13q33.3gain

913q33.3-

13q34loss

Published in partnership with the Center of Excellence in Genomic Medicine Research npj Genomic Medicine (2016) 16033

Conotruncal heart anomalies CCY Mak et al

4

Figure 1. Patient facial phenotypes. (a) Patient 1 with 1.8 Mb gain at 2q22.3 involving the ZEB2 gene. (b) Patient 10 with 611 kb loss at 17p13.3 involving the NXN gene. No consent for clinical photos for Patient 9.

Notably, a similarly sized (544 kb) rare copy-number gain overlapping the NXN gene was reported by another group to be associated with ToF.10 As PAVSD may be considered to be at the more severe end of the spectrum of the cardiac phenotype in relation to ToF, the loss of one copy at this locus may implicate a more profound cardiac developmental effect for CHD than would a gain of copy. Isolated case reports indicate a link between CHD and large deletions in this region at the resolution of the chromosome band, including two from Taiwan and Japan with ToF30 and

PAVSD,31 respectively, that overlap NXN. This suggests that a critical region for CHD may be located in the region of the NXN gene (Figure 3), in addition to the more proximal MDS region (Figure 3).

DISCUSSIONLarge rare CNVs in conotruncal heart diseaseA number of studies have explored the contribution of CNVs to conotruncal heart defects, particularly ToF. Most of these studies involved samples from predominantly Caucasian popula-

tions.6,7,911,32 Although there is evidence to suggest an increased burden of large rare CNVs in ToF patients compared with the normal population,6,10,11 few loci have been found to be

recurrently associated with the disease, implicating substantial genetic heterogeneity.33 The most notable recurrent CNVs responsible for ToF are 22q11.2 deletions, followed by 1q21.1 duplications.7,10,11 In this study of adult Chinese patients with

CTDs, after excluding those with recurrent 22q11.2 deletions(15.3%), we discovered large rare CNVs in 10 of 108 (9.3%) patients; the previously reported 1q21.1 duplication involving GJA5 gene was not found.7,10,11 We were able to characterise the

clinical features of these 10 patients by clinical re-evaluation and determine that three were syndromic. All three of these patients carried a de novo CNV harbouring genes potentially linked to cardiac development that would be considered pathogenic. The other seven large genic CNVs, all very rare and in patients with few or no reported extracardiac features, would represent variants of unknown clinical signicance (VUS). For the two of these seven with parental samples available the CNVs were inherited, as is often found for complex developmental conditions.34

npj Genomic Medicine (2016) 16033 Published in partnership with the Center of Excellence in Genomic Medicine Research

Conotruncal heart anomalies CCY Mak et al

5

Figure 2. CNVs with conotruncal heart disease in the 13q31-34 region. 13three cases of reported ToF with loss CNVs of approximately 2030 Mb (Ballarati et al. and Quelin et al.).22,23 4deletion in patient with double outlet right ventricle, overlapping genes COL4A1 and COL4A2 (McMahon et al.).25 5deletion in Patient 9 from the current study with PAVSD, overlapping genes COL4A1 and COL4A2. 6deletion in patient with single atrium (Yang et al.).21 7duplication in Patient 8 from the current study with PAVSD. 8CHD critical region proposed by Huang et al.20

Figure 3. CNVs associated with congenital heart disease in the 17p13.2-13.3 region and relationship to the MillerDieker Syndrome (MDS) Critical Region. Deletion in Patient 10 in this study overlapping the NXN gene. Deletion in patient with PAVSD, overriding aorta and PDA (Kowase et al.).31 Deletion in patient with Tetralogy of Fallot (Chen et al.).30 Deletion in DECIPHER Patient 4,350 with PDA, PS. Deletion in patient with PDA (Serrenath Nagamani et al. 2009, patient 5). Deletion in patient with PDA, Pulm Hypertension (Shiff et al. 2010, patient B). Duplication in patient with ToF, CNV also overlapping NXN gene (Silversides et al., Case 34).10 Overlapping region of CNV in Patient 10 with other large deletions associated with CHD, a possible critical region for CHD expression. Previously proposed critical region for MDS within the region of genes YWHAE and PAFAH1B1 (Cardoso et al.).50

Possible ethnic specic CNV distribution in ChineseWhist there appears to be some similarities and overlap of CNVs in our cohort compared with the literature, such as the prevalence of recurrent 22q11.2 deletions and overlap of deletion at 5q35.3 (Supplementary Information S1), the distribution of other large rare CNVs in this Chinese CTD cohort seemed to differ to some

extent from the published literature on ToF in the predominantly Caucasian samples studied. Several publications reported recurrent large (over 1 Mb) rare duplications involving the GJA5 gene.7,10,11 Absence of this CNV in our cohort could be due to

insufcient sample size given the ~ 1% prevalence in European ToF cohorts.7 Alternatively, 1q21.1 duplications may be rarer in

Published in partnership with the Center of Excellence in Genomic Medicine Research npj Genomic Medicine (2016) 16033

Conotruncal heart anomalies CCY Mak et al

6

Figure 4. Study protocol for CNVs associated with conotruncal anomalies.

East Asian populations. It is known that there can be considerable stratication of recurrent rearrangements between Asian and European populations,35 and this may be the case for 1q21.1 duplications. Despite our study being the largest Chinese cohort of CTD reported to date, a larger sample sizes will be needed to determine whether 1q21.1 CNVs are truly rarer in individuals of Chinese origin than in those of European descent.

With respect to the rare de novo CNVs found, we speculate that CNVs in the greater 17p13.3 region may be an important cause of developmental disorders in the Chinese population. In general, CNVs in this region are non-recurrent and may lead to various phenotypes other than that of typical MDS. Supporting this possibility, our group has reported cases of 17p13.3 duplication (with aortic stenosis, microcephaly and dysmorphism)36 and

triplication (with split-hand malformation)37 in Hong Kong. Larger studies of the Chinese population will be needed to better classify pathogenic CNVs in the 17p13.3 region, and to determine whether there are any ethnic differences in prevalence.

Advantages and limitationsOur study was carried out using two large groups of European controls,10 and one of Chinese controls,38 to adjudicate rarity of

CNVs. The array platform and stringent criteria used for CNV calling have also been validated before.10 Another advantage of this study is the utilisation of an adult cohort with the availability of long-term medical information. Data regarding both development and later medical complications could be obtained as a result, giving us more information on the potential phenotypic effects of these large rare CNVs over a long period of time. Performing this study on Han Chinese exclusively also meant that we could better begin to characterize the distribution of CNVs related to CHD in this ethnic group, as well as provide documentation of the clinical presentation for reference. Nevertheless, despite being the largest Chinese study so far of CNVs in CHD, our sample size was not sufcient for CNV burden analysis. Furthermore, not all parental samples were available for determination of de novo status of these CNVs.

CONCLUSIONThis study of 116 Chinese patients with conotruncal heart defects identied ten large (4500 kb) rare genic CNVs, three of which were determined to be de novo and pathogenic and the others of less certain clinical signicance. Notably, we did not nd the previously reported recurrent 1q21.1 duplication in our cohort.

npj Genomic Medicine (2016) 16033 Published in partnership with the Center of Excellence in Genomic Medicine Research

Conotruncal heart anomalies CCY Mak et al

7

The small sample size precludes concluding whether or not the distribution of rare CNVs differs in the Chinese population from that observed in European populations. We also provide further evidence for the pathogenicity of the 17p13.3 and 13q33.3 loci, involving the candidate genes COL4A1 and COL4A2, and NXN, respectively.

MATERIALS AND METHODS Patients and samples

A summary of the study design is provided in Figure 4. Adult patients (418 years of age) with CTDs, but no prior genetic diagnosis, were prospectively recruited between February 2012 to July 2013 from the adult congenital heart disease clinic at Queen Mary Hospital, the only tertiary referral centre for CHDs in Hong Kong. All patients were self-reported to be Han Chinese. On a research basis, patients were identied by the attending cardiologist and recruited if a diagnosis of congenital CTD had been recorded. The diagnoses were determined from previous cardiology clinic notes, echocardiograms and operation records based on the cardiologists experience. The standard care of these patients with CHD in Hong Kong does not involve any genetic testing. Using protocols approved by the local institutional review board, written informed consent was obtained for quantitative uorescence PCR (QF-PCR), chromosomal genome-wide microarray analysis and publication of photos (Patient 1 and Patient 10).

Genomic DNA was extracted from whole blood samples. Patients were rst screened for 22q11.2 deletions by QF-PCR followed by conrmatory uorescence in situ hybridisation (FISH) using a standard probe.8 Patients tested 22q11.2 deletion positive8 (n = 21 of 137, 15%) were excluded. We included only the patients with ToF or PAVSD and excluded further two cases with persistent truncus arteriosus (n = 1) and interrupted aortic arch (n = 1) respectively for a more homogenous discovery cohort. Samples were sent for CNV analysis at The Centre for Applied Genomics (The Hospital for Sick Children, Toronto, ON, Canada) using the Affymetrix6.0 SNP array.

Stringent CNV calling criteria and control groupsUsing multiple CNV calling algorithms (Birdsuite,39 iPattern40 and

Affymetrix Genotyping Console), a stringent criteria was used to call the CNVs,10 i.e., at least 10 kb in length, spanning 5 or more consecutive array probes and called by more than one algorithm. We compared CNVs identied to two different groups of population controls and dened rare CNVs as those with a population frequency of o1% using a 50% reciprocal overlap criterion. First, we compared with CNVs analysed in the same way from 3,957 platform-matched controls,10 which included CNVs from the Ottawa Heart Institute Coronary Heart Disease Study, n = 1,234; German PopGen Project, n = 1,123;41 the Ontario Population Genomics Platform, n = 416;10 and Hapmap3 Project,42 n = 1,184. As individuals of Chinese ancestry only formed a small proportion (4.8%) of these controls, as a second stage, we compared CNVs to 1,945 Chinese subjects from the Singapore SgD-CNVdatabase38 in an attempt to eliminate variants with a frequency 41% that were specic to the Chinese population.

Interpretation of CNVsOur analysis prioritised large (4500 kb) CNVs that overlapped genes since these are more likely to be of clinical signicance. CNVs not observed in any of the total 5,902 controls were classied as very rare. Only autosomal CNVs were examined.

The pathogenicity of the CNVs in relation to CTD was analysed systematically using online databases. Each CNV and gene was manually interrogated for the likelihood of causing a cardiac disease phenotype. The databases used for CNV analysis included Database of Genomic Variants (DGV),43 Online Mendelian Inheritance in Man (OMIM), Decipher44 and ISCA.45 Each gene overlapped by the CNV was studied using evidence from PubMed, OMIM, Mouse Genome Informatics (MGI)46 and ZFIN47 to elucidate any potential relevance to congenital heart anomalies and identied as candidate genes. Genomic parameters used were from GRCh37/hg19.

Validation of clinically relevant CNVsCNVs classied as putatively clinically relevant were validated using

NimbleGen CGX-135K oligonucleotide array (Signature Genomics, Roche

NimbleGen, Madison, WI, USA) or 60K CGX v.2 oligonucleotide array (manufactured by Agilent Technologies, Santa Clara, CA, USA; PerkinElmer, Turku, Finland), as described previously48 and by the manufacturers. These microarray platforms were used for clinical diagnostic purposes in our institution. The former NimbleGen array has an average resolution of 140 kb, with a higher resolution (40 kb or less) in regions thought to be of signicance in human development, while the latter Agilent array has an average resolution 190 kb with a higher resolution of 28 kb in the respective regions of focus. The Genoglyphix software (Signature Genomics, Spokane, WA, USA) was used to analyse and annotate the results. Only a limited number of parental samples were obtained when calling back the patients, and trio aCGH testing was offered where possible. Several large rare CNVs detected on the Affy 6.0 platform could not be validated by the NimbleGen CGX-135K platform or quantitative real-time PCR. These CNVs were either in regions of high density of segmental duplications or in close proximity to the centromere. To eliminate these false positives, all large rare CNVs with 470% segmental duplications overlap or near the centromere (e.g., 4p11) were also excluded from this study.

Detailed phenotyping of patients with large, rare CNVsPatients with potentially clinically relevant CNVs were called back for reassessment by a geneticist and given genetic counselling by a genetic counsellor. Lifetime medical history was reviewed and examination was performed for these patients to systematically characterize any syndromic (two out of three features of learning difculties, global dysmorphic facial features, hypernasal voice49) or other extracardiac features. This encompassed birth history, medical history including thyroid function, hearing and immunodeciency, as well as speech, learning and behavioural difculties. A family history of CHD, intellectual disability and congenital anomalies was also obtained.

Ethics statementEthical approval was obtained from the Institutional Review Board of the

University of Hong Kong and Hospital Authority of Hong Kong West Cluster.

ACKNOWLEDGEMENTS

We thank the patients and families for publication of this report. This work was supported by grants from the Childrens Heart Foundation of Hong Kong, S K Yee Medical Research Fund and Small Project Funding of the University of Hong Kong. The funding sources had no involvement in the study design, in the collection, analysis and interpretation of data, in the writing of the report, or in the decision to submit the article for publication.

CONTRIBUTIONS

A.P.L., B.H.C. and C.C.M. designed the project, with contributions from C.R.M., S.W.S. and A.S.B. C.C.M. performed the analysis, reviewed the literature, and drafted the manuscript. P.C.C. and A.K.C. recruited the patients and provided advice regarding the case descriptions and manuscript. K.Y.C., performed the oligonucleotide array analysis. Y.W.C was involved in patient recruitment and genetic counselling. G.T.M. G. K.L. and K.S.Y. provided support for laboratory work and the management of patient data. C.L. provided advice on the literature review and analysis of CNV data. C.R.M. facilitated the SNP array analysis and bioinformatics processing of the CNV data. C.R.M., S.W.S. and A.S.B. revised the manuscript and provided valuable advice on the analysis. B.H.C. oversaw the production of this manuscript and provided guidance and advice for the project. All authors contributed to the discussion and overall data interpretation, provided intellectual input and approved the nal manuscript.

COMPETING INTERESTS

The authors declare no conict of interest.

REFERENCES1. Pediatric Cardiac Genomics C, Gelb, B. et al. The Congenital Heart Disease Genetic Network Study: rationale, design, and early results. Circ. Res. 112, 698706 (2013).

2. Leirgul, E. et al. Birth prevalence of congenital heart defects in Norway 19942009--a nationwide study. Am. Heart J. 168, 956964 (2014).

Published in partnership with the Center of Excellence in Genomic Medicine Research npj Genomic Medicine (2016) 16033

Conotruncal heart anomalies CCY Mak et al

8

3. Reller, M. D., Strickland, M. J., Riehle-Colarusso, T., Mahle, W. T. & Correa, A. Prevalence of congenital heart defects in metropolitan Atlanta, 19982005.J. Pediatr. 153, 807813 (2008).4. Wu, M. H. et al. Prevalence of congenital heart disease at live birth in Taiwan.J. Pediatr. 156, 782785 (2010).5. Jacobs, E. G., Leung, M. P. & Karlberg, J. Distribution of symptomatic congenital heart disease in Hong Kong. Pediatr. Cardiol. 21, 148157 (2000).

6. Bittel, D. C. et al. Ultra high-resolution gene centric genomic structural analysis of a non-syndromic congenital heart defect, Tetralogy of Fallot. PLoS ONE 9, e87472 (2014).

7. Greenway, S. C. et al. De novo copy number variants identify new genes and loci in isolated sporadic tetralogy of Fallot. Nat. Genet. 41, 931935 (2009).

8. Liu, A. P. et al. Under-recognition of 22q11.2 deletion in adult Chinese patients with conotruncal anomalies: implications in transitional care. Eur. J. Med. Genet. 57, 306311 (2014).

9. Peyvandi, S. et al. 22q11.2 deletions in patients with conotruncal defects: data from 1,610 consecutive cases. Pediatr. Cardiol. 34, 16871694 (2013).

10. Silversides, C. K. et al. Rare copy number variations in adults with tetralogy of Fallot implicate novel risk gene pathways. PLoS Genet. 8, e1002843 (2012).

11. Soemedi, R. et al. Contribution of global rare copy-number variants to the risk of sporadic congenital heart disease. Am. J. Hum. Genet. 91, 489501 (2012).

12. Shaheen, R. et al. Positional mapping of PRKD1, NRP1 and PRDM1 as novel candidate disease genes in truncus arteriosus. J. Med. Genet. 52, 322329 (2015).

13. Ghoumid, J. et al. ZEB2 zinc-nger missense mutations lead to hypomorphic alleles and a mild Mowat-Wilson syndrome. Hum. Mol. Genet. 22, 26522661 (2013).

14. Nagase, T. et al. Prediction of the coding sequences of unidentied human genes.XI. The complete sequences of 100 new cDNA clones from brain which code for large proteins in vitro. DNA Res. 5, 277286 (1998).15. Uhlen, M. et al. Proteomics. Tissue-based map of the human proteome. Science 347, 1260419 (2015).

16. Postigo, A. A. Opposing functions of ZEB proteins in the regulation of the TGFbeta/BMP signaling pathway. EMBO J. 22, 24432452 (2003).

17. Lamouille, S., Xu, J. & Derynck, R. Molecular mechanisms of epithelialmesenchymal transition. Nat. Rev. Mol. Cell Biol. 15, 178196 (2014).

18. Moal, F. D. L. et al. ZFHX1B mutations in patients with Mowat-Wilson syndrome. Hum. Mutat. 28, 313321 (2007).

19. Goossens, S. et al. The EMT regulator Zeb2/Sip1 is essential for murine embryonic hematopoietic stem/progenitor cell differentiation and mobilization. Blood 117, 56205630 (2011).

20. Huang, C. et al. Congenital heart defect and mental retardation in a patient with a 13q33.1-34 deletion. Gene 498, 308310 (2012).

21. Yang, Y. F. et al. A 1.1Mb deletion in distal 13q deletion syndrome region with congenital heart defect and postaxial polydactyly: additional support for a CHD locus at distal 13q34 region. Gene 528, 5154 (2013).

22. Ballarati, L. et al. 13q Deletion and central nervous system anomalies: further insights from karyotype-phenotype analyses of 14 patients. J. Med. Genet. 44, e60 (2007).

23. Quelin, C. et al. Twelve new patients with 13q deletion syndrome: genotypephenotype analyses in progress. Eur. J. Med. Genet. 52, 4146 (2009).

24. Kirchhoff, M. et al. Phenotype and 244k array-CGH characterization of chromo-some 13q deletions: an update of the phenotypic map of 13q21.1-qter. Am. J. Med. Genet. A 149A, 894905 (2009).

25. McMahon, C. J., Breathnach, C., Betts, D. R., Sharkey, F. H. & Greally, M. T. De Novo interstitial deletion 13q33.3q34 in a male patient with double outlet right ventricle, microcephaly, dysmorphic craniofacial ndings, and motor and developmental delay. Am. J. Med. Genet. A 167A, 11341141 (2015).

26. Davy, A. & Soriano, P. Ephrin-B2 forward signaling regulates somite patterning and neural crest cell development. Dev. Biol. 304, 182193 (2007).

27. Funato, Y. et al. Nucleoredoxin sustains Wnt/beta-catenin signaling by retaining a pool of inactive dishevelled protein. Curr. Biol. 20, 19451952 (2010).

28. Huang, N., Lee, I., Marcotte, E. M. & Hurles, M. E. Characterising and predicting haploinsufciency in the human genome. PLoS Genet. 6, e1001154 (2010).

29. Petrovski, S., Wang, Q., Heinzen, E. L., Allen, A. S. & Goldstein, D. B. Genic intolerance to functional variation and the interpretation of personal genomes. PLoS Genet. 9, e1003709 (2013).

30. Chen, C. P. et al. Ventriculomegaly, intrauterine growth restriction, and congenital heart defects as salient prenatal sonographic ndings of miller-dieker lissencephaly syndrome associated with monosomy 17p (17p13.2 -4Pter) in a Fetus.

Taiwan. J. Obstet. Gynecol. 49, 8186 (2010).31. Kowase, K. et al. Novel deletion on the short arm of chromosome 17 in a patient with multiple cardiac anomalies. Jpn Circ. J. 61, 882885 (1997).

32. Warburton, D. et al. The contribution of de novo and rare inherited copy number changes to congenital heart disease in an unselected sample of children with conotruncal defects or hypoplastic left heart disease. Hum. Genet. 133, 1127 (2014).

33. Bansal, V. et al. Outlier-based identication of copy number variations using targeted resequencing in a small cohort of patients with Tetralogy of Fallot.PLoS ONE 9, e85375 (2014).

34. Bassett, A. S., Scherer, S. W. & Brzustowicz, L. M. Copy number variations in schizophrenia: critical review and new perspectives on concepts of genetics and disease. Am. J. Psychiatry 167, 899914 (2010).

35. Antonacci, F. et al. Palindromic GOLGA8 core duplicons promote chromosome 15q13.3 microdeletion and evolutionary instability. Nat. Genet. 46, 12931302 (2014).

36. Ho, A. C. et al. A newborn with a 790 kb chromosome 17p13.3 microduplication presenting with aortic stenosis, microcephaly and dysmorphic facial featuresis cardiac assessment necessary for all patients with 17p13.3 microduplication? Eur.J. Med. Genet. 55, 758762 (2012).37. Luk, H. M. et al. A prenatal case of split-hand malformation associated with 17p13.3 triplicationa dilemma in genetic counseling. Eur. J. Med. Genet. 57, 8184 (2014).

38. Xu, H. et al. SgD-CNV, a database for common and rare copy number variants in three Asian populations. Hum. Mutat. 32, 13411349 (2011).

39. Korn, J. M. et al. Integrated genotype calling and association analysis of SNPs, common copy number polymorphisms and rare CNVs. Nat. Genet. 40, 12531260 (2008).

40. Pinto, D. et al. Comprehensive assessment of array-based platforms and calling algorithms for detection of copy number variants. Nat. Biotechnol. 29, 512520 (2011).

41. Lionel A. C. et al. Rare copy number variation discovery and cross-disorder comparisons identify risk genes for ADHD. Sci. Transl. Med. 3: 95ra75. (2011).

42. International HapMap C, Altshuler, D. M. et al. Integrating common and rare genetic variation in diverse human populations. Nature 467, 5258 (2010).

43. MacDonald J. R., Ziman R., Yuen R. K., Feuk L. & Scherer S. W. The Database of Genomic Variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 42(Database issue) 2014; 42 (Database issue): D986D992.

44. Swaminathan, G. J. et al. DECIPHER: web-based, community resource for clinical interpretation of rare variants in developmental disorders. Hum. Mol. Genet. 21 (R1): R37R44 (2012).

45. Riggs, E. R., Jackson, L., Miller, D. T. & Van Vooren, S. Phenotypic information in genomic variant databases enhances clinical care and research: the International Standards for Cytogenomic Arrays Consortium experience. Hum. Mutat. 33, 787796 (2012).

46. Eppig J. T., Blake J. A., Bult C. J., Kadin J. A. & Richardson J. E. Mouse Genome Database G The Mouse Genome Database (MGD): facilitating mouse as a model for human biology and disease. Nucleic Acids Res. 43 (Database issue): D726D736 (2015).

47. Sprague, J. et al. The Zebrash Information Network: the zebrash model organism database. Nucleic Acids Res. 34 (Database issue) D581D585 (2006).

48. Kan, A. S. et al. Whole-genome array CGH evaluation for replacing prenatal karyotyping in Hong Kong. PLoS ONE 9, e87988 (2014).

49. Fung, W. L., Chow, E. W., Webb, G. D., Gatzoulis, M. A. & Bassett, A. S. Extracardiac features predicting 22q11.2 deletion syndrome in adult congenital heart disease.Int. J. Cardiol. 131, 5158 (2008).

50. Cardoso, C. et al. Renement of a 400-kb critical region allows genotypic differentiation between isolated lissencephaly, Miller-Dieker syndrome, and other phenotypes secondary to deletions of 17p13.3. Am. J. Hum. Genet. 72, 918930 (2003).

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the articles Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

Web End =http://creativecommons.org/licenses/ http://creativecommons.org/licenses/by/4.0/

Web End =by/4.0/

The Author(s) 2016

Supplementary Information accompanies the paper on the npj Genomic Medicine website (http://www.nature.com/npjgenmed)

npj Genomic Medicine (2016) 16033 Published in partnership with the Center of Excellence in Genomic Medicine Research