Copy number variations in candidate genomic regions confirm genetic heterogeneity and parental bias in Hirschsprung disease
Orphanet Journal of Rare Diseases volume 14, Article number: 270 (2019)
Hirschsprung Disease (HSCR) is a congenital defect of the intestinal innervations characterized by complex inheritance. Many susceptibility genes including RET, the major HSCR gene, and several linked regions and associated loci have been shown to contribute to disease pathogenesis. Nonetheless, a proportion of patients still remains unexplained. Copy Number Variations (CNVs) have already been involved in HSCR, and for this reason we performed Comparative Genomic Hybridization (CGH), using a custom array with high density probes.
A total of 20 HSCR candidate regions/genes was tested in 55 sporadic patients and four patients with already known chromosomal aberrations. Among 83 calls, 12 variants were experimentally validated, three of which involving the HSCR crucial genes SEMA3A/3D, NRG1, and PHOX2B. Conversely RET involvement in HSCR does not seem to rely on the presence of CNVs while, interestingly, several gains and losses did co-occur with another RET defect, thus confirming that more than one predisposing event is necessary for HSCR to develop. New loci were also shown to be involved, such as ALDH1A2, already found to play a major role in the enteric nervous system. Finally, all the inherited CNVs were of maternal origin.
Our results confirm a wide genetic heterogeneity in HSCR occurrence and support a role of candidate genes in expression regulation and cell signaling, thus contributing to depict further the molecular complexity of the genomic regions involved in the Enteric Nervous System development. The observed maternal transmission bias for HSCR associated CNVs supports the hypothesis that in females these variants might be more tolerated, requiring additional alterations to develop HSCR disease.
Hirschsprung Disease (HSCR) (OMIM# 142623) is a congenital intestinal aganglionosis caused by a premature arrest of the cranio-caudal migration of neural crest cells during embryogenesis, showing an incidence of around 1/5000 live births . The phenotype is highly variable, with severity classified according to the length of the affected gastrointestinal tract as short-segment HSCR (aganglionosis does not extend beyond the upper sigmoid, S-HSCR: 80% of cases) and long-segment HSCR (L-HSCR: 20% of cases). About 70% of cases presents as isolated, while 30% shows additional anomalies, including chromosomal aberrations . HSCR is characterized by increased sibling recurrence risk varying by gender, length of aganglionosis and familial occurrence (up to 80% of cases are sporadic). These observations, together with a distorted sex ratio (M:F = 4:1), make HSCR a model for complex genetic disease. Genetic heterogeneity in HSCR is demonstrated by involvement of several genes and loci [2,3,4]. The major gene involved in isolated HSCR is the RET proto-oncogene , located on 10q11.2 and linked to HSCR in 90% of the familial forms . In addition, HSCR has been associated with several RET polymorphisms, most of which part of a common risk haplotype encompassing the RET gene from the promoter to exon 2 [5,6,7,8]. An association with NRG1 (OMIM# 142445) and SEMA3A (OMIM# 603961) / SEMA3D (OMIM# 609907) has also been described [9,10,11]. However, a proportion of HSCR patients still remains unexplained as only 50% of familial and 7–35% of sporadic HSCR cases can be explained by RET coding variants . Deletions in genes already known to be involved in HSCR might account for these latter cases. Indeed, the RET gene was discovered following observation of de novo interstitial deletions of 10q11.2 [12, 13] and about 12% of HSCR patients have structural abnormalities .
Copy number variations (CNVs), already proven to be genetic risk factors in disease pathogenesis [14, 15], might thus account for part of the missing heritability in HSCR. Jiang et al. (2011) performed a custom-designed array CGH to examine 67 candidate HSCR genes in 18 HSCR patients, identifying seven CNVs at three loci, all likely hosting regulatory genes in syndromic HSCR patients . On the other hand, Tang et al. (2012) assessed the CNV contribution to HSCR from genome-wide SNP data finding a greater burden for rare CNVs in HSCR cases over controls and larger CNVs in syndromic HSCR than in isolated cases. Only six CNVs overlapped with known HSCR loci, none involving the RET gene . Another study investigated 123 HSCR patients and 432 unaffected subjects, with Illumina’s HumanOmni1-Quad BeadChip, finding 16 CNV regions associated with HSCR . Finally, very recently, Tilghman et al. (2019) have dissected, through both karyotyping and exome sequencing, the differential contribution to HSCR development of three different molecular classes of risk alleles, namely rare coding variants, common non coding variants and large CNVs and chromosomal anomalies. In this latter class, at least 9 loci have been reported, whose aberrations showed a very high odd ratio (63.07) and are involved in 11.4% of the patients .
To further explore genetic heterogeneity in HSCR, we have performed CGH, using a custom array with high density probes and focusing on a total of 20 candidate regions/genes already known to be involved in HSCR, on a selected panel of 55 sporadic HSCR previously genotyped at the RET locus  and four HSCR patients with already known chromosomal aberrations.
HSCR patients and regions analyzed
A total of 55 Italian sporadic HSCR patients fully genotyped at the RET locus  were retrospectively included in the study. To investigate the possible presence of interstitial deletions of the RET region, we selected 52 cases homozygous for the risk haplotype [5,6,7,8] as well as 3 patients carrying very uncommon haplotypes, consistent with possible hemizygosity of the same region. Four additional HSCR patients were included as positive controls: two with a de novo deletion at the centromeric region of chromosome 10 [12, 13], one with an inverted duplication at chromosome 22, and another one with a trisomy 21 (in addition to two HSCR patients already included in the sample set and presenting with Down syndrome, OMIM# 190685). The whole sample analyzed is thus constituted by 59 HSCR patients.
Besides the major RET gene, other candidate genes and loci were selected for the analysis based on i) linkage with HSCR, ii) association with the disease, iii) mutation in syndromic and isolated HSCR patients, iv) involvement in the transcriptional regulation of RET, and v) preliminary evidence, not confirmed later. Finally, other loci were included because altered in disorders presenting HSCR with a higher prevalence than the general population (Table 1).
The selected HSCR patients and positive controls underwent custom aCGH. As reported in Additional file 1, a total of 75 calls were estimated from the raw data using the Agilent Aberration Detection method as described under Methods. Two of these calls corresponded to already known trisomies of chromosome 21, and four calls identified the alterations included as controls.
In addition, we evaluated the profiles of all the samples by visual inspection. This allowed us to add 6 calls to the list of aberrations, undetected by the software. Finally, although we did not expect any new variant in replicates, in a triplicate sample we found two aberrations that had not been detected in the array firstly investigated, but were present in both the two replicates and thus assumed as reliable, for a total of 83 aberrations detected in 64 different chromosomal locations in 44 samples (Additional file 1, Figure S1). Twenty-six of these aberrations had previously been reported on DGV; all of them were common CNVs (frequency > 5%), with the exception of the deletion at SEMA3A/3D, that is compatible with two CNVs detected by sequencing [26, 27] with an overall frequency of less than 1% (Table 2, Additional file 1: Figure S1). Four of the CNVs common in DGV were recurrent in the samples. The variants at 15q11 and 5q13, and the gain at 9p11, common on DGV, showed frequencies roughly similar to those reported for controls in the high resolution gnomAD (https://gnomad.broadinstitute.org/) and deciphering developmental disorders (DDD) (https://decipher.sanger.ac.uk/) databases [28, 29]. Losses at 9p11 and 9q31 had frequencies in between gnomAD and DDD. Interestingly, variants on 22q11.2 were all but one more frequent in our sample (Additional file 2). Moreover, three regions, found to carry anomalies in 5 samples, are compatible with CNVs reported on the Decipher database. The region 22:25672585–25,892,401 was found duplicated in two patients and deleted in a third patient, with anomalies also reported on DGV, with frequencies similar to controls in gnomAD and DDD and regarded as likely benign common CNVs. Also the deletion at 9:113025039–113,029,430 is common on DGV and likely benign, but interestingly it is reported in a patients affected by aganglionic megacolon, intellectual disability and short stature. Finally, the deletion at 1:146638075–149,224,043 is compatible with several deletions reported on decipher including the 1q21.1 recurrent microdeletion (OMIM# 612474).
In conclusion, excluding the control regions and chromosome 21 for the two Down syndrome HSCR patients, we have detected 51 novel aberrations, plus the one reported on DGV with very low frequency (Additional file 1, Figure S1). As two of these variants were recurrent in three patients each, we had a total of 48 distinct variants detected in 25 patients. Several of them looked unlikely at the visual inspection and, as a matter of fact, mainly not confirmed on the replicate, when available. Conversely, most of the variants classified as likely or possible at the visual inspection were also replicated  (Additional file 1).
CNVs already reported in HSCR
Despite 41 aberrations were called at the RET locus in 40 samples, they were barely overlapping and rarely located on the risk haplotype. Only two of them were considered after applying the MALR> 0.30 criteria (see Methods), but no one looked as likely at visual inspection and could be confirmed at validation, proving that those CNVs were false positive. Therefore, neither deletions hypothesized based on homozigosity nor very rare haplotypes compatible with hemizygous conditions were sustained by the present data.
Most of the regions reported as duplicated or deleted in other studies [16, 17] were not included in our selected regions, neither we could detect any aberration comparable with those already reported when covered by probes in our design.
Variant validation and parental origin
We elected to focus on gains/losses that seemed more promising at the visual inspection (classified as likely) and not reported on DGV, in addition to one deletion of the SEMA3A/SEMA3D region, compatible with a CNV reported on DGV but with a low frequency. We included also three aberrations unlikely at the visual inspection but particularly interesting as located in the SEMA3A/SEMA3D region and the RET locus, two master loci in HSCR development [1, 10, 11]. We have thus selected for validation 18 aberrations found in 14 patients (Table 3, Fig. 1, Fig. 2, Additional file 1).
Eleven of these 18 aberrations were confirmed in 11 patients (Table 3), including one variant that gave no conclusive results by validation with a different method but was confirmed on a replicate. Seven were instead not confirmed, mostly of which on RET and/or already considered as unlikely after visual inspection. In addition, one variant initially not selected for validation, because considered unlikely at the visual inspection, was confirmed on a replicate and resulted recurrent in two other samples that have not been further analyzed. Interestingly, some of these 12 true novel aberrations involved the NRG1, SEMA3A/SEMA3D and PHOX2B loci, three of the strongest candidates among our target regions. In particular, the short deletion detected between exon 6 and 7 of NRG1 affected a male isolated patient with L-form HSCR, the ~ 9 kb deletion in SEMA3A/3D was intergenic (between SEMA3A and SEMA3D), and was inherited from the mother by a male S-form HSCR isolated patient. We could also detect a deletion that involved almost the entire PHOX2B gene in a female isolated S-form HSCR patient, shown to be inherited by the mother. Other CNVs were found at known HSCR candidate loci 9q31  (in four patients, although not overlapping between each other and covering a gene in only one patient), 16q23.3 (a short inter-exonic gain) , and 19q12 (not involving any known gene) . In addition, a CNV recurrent in three samples has been detected at 22q11.2, a locus affected in DiGeorge and VeloCardioFacial syndromes, and in der (22) and Cat-Eye syndrome, disorders presenting HSCR with a higher prevalence than the general population . Finally two CNVs, at 1q21 and 15q21, were found outside any of the HSCR candidate loci represented at high-density probes in the array (Table 3). Among these 12 losses and gains found in 11 patients, the M/F rate was 7/4, with an enrichment in females with respect to the whole sample analyzed (from 28.8 to 36.4%), seven patients had an S-form (70.0% with respect to 58.8% for the whole sample analyzed), while three were L and one unknown (Table 4). Nine were isolated cases, while one patient had Down syndrome and another one presented with cardiac and facial malformations. Finally, 5 patients carried either RET variants or RET locus deletions, with increased frequency of RET anomalies with respect to the original sample (45.5% vs 23.7%). Interestingly, the four patients with aberrations at 9q31 were all defective for RET.
No difference in the distribution of patients’ characteristics was detected compared to the whole sample analyzed, when considering both the newly detected 77 aberrations and the “true” 37 CNVs detected. Indeed, in both cases the median size of CNVs was larger among syndromic patients than among isolated cases: 1,5 Mb vs 16.5 Kb in 9 and 68 newly detected CNVs in syndromic and isolated patients respectively, and 2.8 Mb vs 514 Kb in 5 and 32 “true” CNVs, with a borderline p-value (p = 0,0866), as already reported by others . No difference was detected instead with regards to patients gender, while RET negative patients carried CNVs larger on average than RET mutated patients.
While parents were not available for three patients, and only one gain showed to be de novo, all the other seven validated CNVs resulted inherited by unaffected mothers and none by fathers (100%, exact confidence interval from 59 to 100%, p = 0.0078).
We have performed a high density custom array-CGH to search for DNA copy imbalances at selected candidate genes and loci in a total of 59 HSCR patients. Despite our interest on the RET gene, we could not detect any novel variant at this locus, in line with what reported by others [16,17,18, 25, 32]. Moreover, the RET locus presented with false positive calls, confirming the difficulties raised by this subcentromeric region. Also other genes, previously investigated for deletions and amplifications (ZEB2, EDN3 and GDNF), did not show any alteration [16, 32]. At opposite, we could detect CNVs at two loci, 1q21 and 15q21 (Fig. 2), never implicated in HSCR before, by probes randomly selected to cover the whole genome at low density.
Unfortunately, given the uneven probes distribution of the present design, we could not assess the possible CNVs enrichment in HSCR candidate genes with respect to the other chromosomal regions. Nonetheless, data from controls were looked up in the high resolution databases gnomAD and DDD [28, 29] and analyzed to investigate the CNVs detected in our panel of HSCR cases. As shown in the Additional file 2, we were able to demonstrate the presence of several novel deletions/duplications in candidate genes and loci and to suggest an enrichment of common CNVs in 22q11.2 over controls. 9q31 and 9p11 losses have resulted with a frequency significantly different compared to both the control databases. However, as these latter control frequencies are very different from each other, a degree of discordance can be hypothesized between the control sets of these two databases.
Linkage of HSCR to 9q31 was shown in families with no or hypomorphic RET gene mutations, suggesting that these latters would require the action of other defects . Novel chromosomal variants at locus 9q31 were confirmed in our dataset in 4 out of 59 patients, who also carried either heterozygous missense variants of the RET gene or large deletions at the RET locus, associations consistent with the digenic inheritance of HSCR already suggested [3, 19]. Other studies have pointed to 9q31 as a region involved in HSCR, but only a few suggestive causative genes have been identified so far. Among these, IKBKAP (OMIM# 603722) was found associated with HSCR in Chinese samples, especially in patients carrying RET coding variants , suggesting population specificity and implying that, in agreement with our observations, RET variants are found to co-occur with additional chromosomal anomalies. Interestingly, despite lack of concordance about HSCR gene(s) on 9q31, a quantitative linkage analysis carried out on genes likely involved in the enteric nervous system development identified a “master regulator” locus in 9q31 .
We also detected novel CNVs in SEMA3A/3D, NRG1 and PHOX2B (Fig. 2). Class 3 Semaphorins, known to be involved in neuronal migration, proliferation, survival, and axonal guidance , have been demonstrated to be HSCR susceptibility factors [10, 11]. The importance of SEMA3D signaling in the ENS is further supported by gene expression comparison between wild-type and Ret k−/k- mouse gastrointestinal tracts . Neuregulin 1 (NRG1) is essential for the development of the nervous system and the heart and its deregulation has been linked to cancer, schizophrenia and bipolar disorder (BPD) (OMIM# 181500) . NRG1 has also been identified as an additional HSCR susceptibility locus in Asian populations [9, 38]. Such an association, found initially to be below genome-wide significance in Caucasians [22, 39], has been demonstrated also in an European population . Moreover, NRG1 expression has been found to be significantly higher in HSCR than in controls tissues . An interplay between RET and NRG1 has been suggested . We could also detect a deletion of the PHOX2B gene, a gene that codes for a homeodomain transcription factor involved in the development of several noradrenergic neuronal populations in the autonomous nervous system . Different heterozygous mutations of PHOX2B are known to cause Congenital Central Hypoventilation Syndrome (CCHS) (OMIM# 209880), a rare disease characterized by impaired ventilator response to hypercapnia and hypoxia, often associated with HSCR and neuroblastomas . A PHOX2B interstitial deletion, as well as PHOX2B mutations, have been reported in HSCR patients [20, 21]. Moreover, in-frame deleted and common polyA contracted alleles of the PHOX2B gene have shown to either abolish or reduce the transactivation activity of the mutant proteins, respectively [21, 42]. Therefore our results support PHOX2B loss-of-function as a rare cause of the HSCR phenotype.
We also confirmed one aberration affecting 19q12, a locus found to be in linkage with HSCR , and detected a de novo gain at locus 16q23.3 (Fig. 2), previously identified by a genome-wide association study in 43 Mennonite family trios . The only known gene in this region is MPHOSPH6 (M-phase PHOSPHoprotein6) (OMIM# 605500), an exosome-associated protein that is phosphorylated during mitosis .
In addition, we found quite large aberrations in two regions not included among the selected candidate loci. The 15q21.3 locus contains several genes, among which ALDH1A2 (Aldehyde Dehydrogenase 1 Family, Member A2) (OMIM# 603687) is particularly interesting, coding for an enzyme that catalyzes the synthesis of retinoic acid (RA) from retinaldehyde. RA is a hormonal signaling molecule critical during embryonic development, that has already been documented as a regulator of RET expression in cardiac and renal development, to delay the colonization of the hindgut by RET-positive enteric neuroblasts, and to result in ectopic RET expression during embryogenesis. RA has also been proposed to maintain migratory signals and deficiency of its precursor, Vitamin A, and could therefore increase HSCR penetrance and expressivity . Finally, targeted inactivation of the mouse aldh1a2 has been shown to lead to agenesis of the enteric ganglia, a condition reminiscent of human Hirschprung’s disease . Another interesting gene in the same region is ADAM10 (OMIM# 602192), a member of the ADAM family, cell surface proteins with both adhesion and protease domains, which cleaves TNF-alpha, E-cadherin, L1cam and other proteins, besides regulating Notch signaling, a process required for progenitor cell lineage specification and maintenance .
The deletion in 1q21.1-q21.2 spans the 1q microdeletion syndrome region, increasing the risk of delayed development, intellectual disability, physical abnormalities, and neurological and psychiatric problems. Recently, a deletion and two duplications at locus 1q21.1 have been detected in HSCR patients and reported as significantly overrepresented compared to controls, thus confirming the consistency of our result . The patient carrying this maternally inherited deletion is a female with S-form HSCR without any RET coding variant, reported to present ventricular septal defect (VSD), mandibular hypoplasia, and low-set ears. As far as we know, this is the first case reported of 1q21 microdeletion syndrome associated with HSCR.
In our study, we have elected to focus on CNVs never reported on DGV database. Surprisingly, seven out of 8 of such novel and confirmed CNVs were inherited by the unaffected mothers, while the eighth one had a de novo occurrence. Similar observations have already been made in Type 2 Diabetes (T2D) (OMIM# 125853) and long QT syndrome (OMIM# 192500) [46, 47] and it has also been demonstrated that inherited CNVs might be pathogenic . Interestingly, an increase in CNV burden in the mothers of children affected by mental retardation (MR) have been recently reported in the Chinese population, suggesting that females might be more tolerant than males to deleterious variations and that MR manifestation for females might require a higher burden of deleterious variants . Consistent with our observation of an excess of inheritance from unaffected mothers, a parental mutation transmission asymmetry has already been reported at the RET locus . Apparently, this bias in the transmission of RET single base mutations was not due to different expression of the disease depending on the gender of the transmitting parent, but rather to differential reproductive rate between male and female carriers, with mothers carrying a severe mutation that would be more likely than fathers to reproduce and transmit . Considering that HSCR penetrance is less reduced in male than females, with recurrent risk for male sibs higher when the HSCR patient is female [1, 5], the differential paternal and maternal fitness and the supposed greater mutation burden toleration in females seem to be the two faces of the same coin. Therefore, asymptomatic females carrying these HSCR structural variants would be more likely to transmit the causative CNV alleles, as already demonstrated in other diseases, especially in the presence of additional alterations.
Three fundamental genes already involved in HSCR pathogenesis, namely SEMA3A/3D, NRG1 and PHOX2B, have been shown in our study to play a role also through the presence of CNVs, thus suggesting their haploinsufficiency is responsible of damaging effects. Indeed, to our knowledge the presence of structural variants in HSCR patients has never or rarely been reported for NRG1, SEMA3A/3D, and PHOX2B. Conversely RET involvement in HSCR does not seem to rely on the presence of CNVs but, interestingly, several gains and losses did co-occur with another RET defect in our sample, thus sustaining the hypothesis that more than one predisposing event is necessary for HSCR to develop. Our results, not surprisingly for a complex genetic disease like HSCR, support a role of candidate genes in transcription and expression regulation and in ENS development, confirming the known genetic heterogeneity and showing the possible involvement of new loci. Finally, all the CNVs shown to be inherited in our samples were of maternal origin, including the four novel CNVs detected on 9q31 affecting patients who were also carrying variants of the RET proto-oncogene.
HSCR patients and microarray design
A total of 55 Italian sporadic HSCR patients were retrospectively included in the study, as described under the “Results” section.
We performed high-resolution oligonucleotide array-CGH analyses on 20 candidate genes/loci known to be involved in HSCR, using a customized 8x15K array (Agilent Technologies, Santa Clara, CA, USA), in accordance with manufacturers’ instructions, whose details are reported in Table 1. Additional probes were also distributed along the whole genome, including those useful for sample replication and normalization. Genomic positions are based on the Human Genome GRCh37 (hg19) assembly (http://genome-euro.ucsc.edu).
To assess genomic imbalances, we applied the Aberration Detection Methods ADM-2 with a threshold of 6, as recommended by Agilent. We also applied the centralization and the GC correction algorithms and considered as aberrant only those regions with a minimum of 2 consecutive probes and exceeding a Mean Absolute Log2 Ratio of 0.30 (referred as MALR> 0.30). In addition, sample profiles were evaluated at the whole genome level by visual inspection. Such a manual search was performed paying attention to log2 ratios values above 0.5 and below − 1.0 and taking into consideration not only the overall profile of the single individual sample (and its quality) but also the profiles of the entire cohort. The above parameters were kept, with the exception of the MALR> 0.30 filter, and observations were confirmed by a second operator. Based on such inspection, we classified the CNVs detected as “likely”, “possible”, or “unlikely”. Loci with nearby gain or loss intervals and an intervening region of more than 2 probes were considered as separate CNVs, as well as those differing for 2 probes with opposite log2 ratios (log2ratio < − 0.3 for gain and > 0.3 for deletions). The quality of the experiments was evaluated on the basis of the QC metrics generated by the Genomic Workbench 5.0.14 software (Agilent Technologies), such as the DLRSpread (derivative log ratio spread), a measure of the log ratio noise for each sample. DLRSs and the other sample metrics are detailed elsewhere .
Ten samples with bad profiles (DRLS ≥0.3) together with other 16 arbitrarily selected samples, were replicated at least once on another array, for a total of 26 samples with at least one replicate. The successive variants search was performed in the replicates of overall better quality.
Aberrations were compared with CNVs observed in the normal population, as reported in the Database of Genomic Variants (DGV, http://dgv.tcag.ca/dgv/app/home), and with the CNVs reported in the DECIPHER database of phenotypes, v8.7 release (https://decipher.sanger.ac.uk/). Data comparing is challenging as exact boundaries of the aberrations are not known, but only assumed to be between the last “normal” probe (outer) and the first “aberrant” probe (inner), depending on the mean coverage. However, we considered aberrations as consistent with those already reported if they showed an overlap ≥80%, did not differ for more than two probes with compatible log2 ratios (that is ≥ |0.3|), and were of the same kind (gain or loss).
To get further insight into the aberrations thus identified, we also compared their frequencies in our sample with frequencies of corresponding CNVs publically available at a resolution similar to that of our design, considering comparable those CNVs with identical boundaries (Additional file 2). To this end we used the control data from the European population in the GnomAD website (https://gnomad.broadinstitute.org/) and from the DDD database, browsed through the UCSC genome browser (https://genome-euro.ucsc.edu/index.html), and assessed statistical differences through the Fisher’s test or the Chisquare test with Yate’s correction for continuity when more feasible.
Results obtained with the custom aCGH, together with the concordance degree among the replicates on the same design array, showed that the replication rate was not very high, and that the visual inspection outperformed the mere software call . However, a high false positive rate is not surprising as a few studies have shown a not infrequent presence of false positive and false negative results from aCGH [51,52,53,54]. As a matter of the fact, Conrad et al. (2010), using quantitative PCR (qPCR) for initial validation of aCGH experiments on 450 HapMap samples, suggested a preliminary false-discovery rate of ~ 20%, then assessed to 15% when comparing the results with another CNV-typing array .
For this reason, the most promising regions, based on rare presence/absence on DGV and visual inspection, were selected to undergo validation with different approaches. Parents were also investigated, when available, to verify whether the aberrations were inherited or de novo. When the deleted region was sufficiently small, the DNA was amplified with PCR and checked for anomalous bands by electrophoresis. For the majority of the aberrations, validation was carried out applying quantitative PCR. Primer pairs were designed with Primer-blast (http://www.ncbi.nlm.nih.gov/tools/primer-blast/) according to stringent parameters to ensure successful assay. qPCR analysis was performed with the LightCycler 480 Instrument and relative software using the SYBR Green I Master following manufacturer’s recommendations (Roche, Manheim, Germany). In case of inconclusive results, we spotted the sample on a CGH 8x60K array (Agilent).
The probability that 7 out of 7 inherited CNVs were transmitted by mothers was checked by the binomial test.
Finally, we defined variant as “true” when either reported on DGV, confirmed on a second array replicate, or validated with a different approach as described above.
Availability of data and materials
All CNV data generated and analyzed during this study are included in this published article and in particular in the Additional file 1.
Amiel J, Sproat-Emison E, Garcia-Barcelo M, Lantieri F, Burzynski G, Borrego S, et al. Hirschsprung disease, associated syndromes and genetics: a review. J Med Genet. 2008;45:1–14.
Gui H, Schriemer D, Cheng WW, Chauhan RK, Antiňolo G, Berrios C, et al. Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes. Genome Biol. 2017;18:48.
Bolk S, Pelet A, Hofstra RM, Angrist M, Salomon R, Croaker D, et al. A human model for multigenic inheritance: phenotypic expression in Hirschsprung disease requires both the RET gene and a new 9q31 locus. Proc Natl Acad Sci U S A. 2000;97:268–73.
Gabriel SB, Salomon R, Pelet A, Angrist M, Amiel J, Fornage M, et al. Segregation at three loci explains familial and population risk in Hirschsprung disease. Nat Genet. 2002;31:89–93.
Emison ES, McCallion AS, Kashuk CS, Bush RT, Grice E, Lin S, et al. A common sex-dependent mutation in a RET enhancer underlies Hirschsprung disease risk. Nature. 2005;434:857–63.
Lantieri F, Griseri P, Puppo F, Campus R, Martucciello G, Ravazzolo R, et al. Haplotypes of the human RET proto-oncogene associated with Hirschsprung disease in the Italian population derive from a single ancestral combination of alleles. Ann Hum Genet. 2006;70:12–26.
Emison ES, Garcia-Barcelo M, Grice EA, Lantieri F, Amiel J, Burzynski G, et al. Differential contributions of rare and common, coding and noncoding ret mutations to multifactorial Hirschsprung disease liability. Am J Hum Genet. 2010;87:60–74.
Chatterjee S, Kapoor A, Akiyama JA, Auer DR, Lee D, Gabriel S, et al. Enhancer variants synergistically drive dysfunction of a gene regulatory network in Hirschsprung disease. Cell. 2016;167:355–68.
Garcia-Barcelo MM, Tang CS, Ngan ES, Lui VC, Chen Y, So MT, et al. Genome-wide association study identifies NRG1 as a susceptibility locus for Hirschsprung's disease. Proc Natl Acad Sci U S A. 2009;106:2694–9.
Wang LL, Zhang Y, Fan Y, Li H, Zhou FH, Miao JN, et al. SEMA3A rs7804122 polymorphism is associated with Hirschsprung disease in the northeastern region of China. Birth Defects Res A Clin Mol Teratol. 2012;94:91–5.
Jiang Q, Arnold S, Heanue T, Kilambi KP, Doan B, Kapoor A, et al. Functional loss of semaphorin 3C and/or semaphorin 3D and their epistatic interaction with ret are critical to Hirschsprung disease liability. Am J Hum Genet. 2015;96:581–96.
Martucciello G, Bicocchi M, Dodero P, Lerone M, Cirillo MS, Puliti A, et al. Total colonic aganglionosis with interstitial deletion of the long arm of chromosome 10. Pediat Surg Int. 1992;7:308–10.
Fewtrell MS, Tam PK, Thomson AH, Fitchett M, Currie J, Huson SM, et al. Hirschsprung's disease associated with a deletion of chromosome 10 (q11.2q21.2): a further link with the neurocristopathies? J Med Genet. 1994;31:325–7.
Yin CL, Chen HI, Li LH, Chien YL, Liao HM, Chou MC, et al. Genome-wide analysis of copy number variations identifies PARK2 as a candidate gene for autism spectrum disorder. Mol Autism. 2016;7:23.
Kushima I, Aleksic B, Nakatochi M, Shimamura T, Shiino T, Yoshimi A, et al. High-resolution copy number variation analysis of schizophrenia in Japan. Mol Psychiatry. 2017;22:430–40.
Jiang Q, Ho YY, Hao L, Nichols Berrios C, Chakravarti A. Copy number variants in candidate genes are genetic modifiers of Hirschsprung disease. PLoS One. 2011;6:e21219.
Tang CS, Cheng G, So MT, Yip BH, Miao XP, Wong EH, et al. Genome-wide copy number analysis uncovers a new HSCR gene: NRG3. PLoS Genet. 2012;8:e1002687.
Bae JS, Koh I, Cheong HS, Seo JM, Kim DY, Oh JT, et al. A genome-wide association analysis of chromosomal aberrations and Hirschsprung disease. Transl Res. 2016;177:31–40.
Tilghman JM, Ling AY, Turner TN, Sosa MX, Krumm N, Chatterjee S, et al. Molecular genetic anatomy and risk profile of Hirschsprung's disease. N Engl J Med. 2019;380:1421–32.
Benailly HK, Lapierre JM, Laudier B, Amiel J, Attié T, De Blois MC, et al. PMX2B, a new candidate gene for Hirschsprung's disease. Clin Genet. 2003;64:204–9.
Fernández RM, Mathieu Y, Luzón-Toro B, Núñez-Torres R, González-Meneses A, Antiñolo G, et al. Contributions of PHOX2B in the pathogenesis of Hirschsprung disease. PLoS One. 2013;8:e54043.
Kapoor A, Jiang Q, Chatterjee S, Chakraborty P, Sosa MX, Berrios C, et al. Population variation in total genetic risk of Hirschsprung disease from common RET, SEMA3 and NRG1 susceptibility polymorphisms. Hum Mol Genet. 2015;24:2997–3003.
Carrasquillo MM, McCallion AS, Puffenberger EG, Kashuk CS, Nouri N, Chakravarti A. Genome-wide association study and mouse model identify interaction between RET and EDNRB pathways in Hirschsprung disease. Nat Genet. 2002;32:237–44.
Leon TY, Ngan ES, Poon HC, So MT, Lui VC, Tam PK, et al. Transcriptional regulation of RET by Nkx2-1, Phox2b, Sox10, and Pax3. J Pediatr Surg. 2009;44:1904–12.
Núñez-Torres R, Fernández RM, López-Alonso M, Antiñolo G, Borrego S. A novel study of copy number variations in Hirschsprung disease using the multiple ligation-dependent probe amplification (MLPA) technique. BMC Med Genet. 2009;10:119.
Abecasis GR, Auton A, Brooks LD, DePristo MA, Durbin RM, Handsaker RE, et al. An integrated map of genetic variation from 1,092 human genomes. Nature. 2012;491:56–65.
Wong LP, Ong RT, Poh WT, Liu X, Chen P, Li R, et al. Deep whole-genome sequencing of 100 southeast Asian Malays. Am J Hum Genet. 2013;92:52–66.
Coe BP, Witherspoon K, Rosenfeld JA, van Bon BW, Vulto-van Silfhout AT, Bosco P, et al. Refining analyses of copy number variation identifies specific genes associated with developmental delay. Nat Genet. 2014;46(10):1063–71.
Cooper GM, Coe BP, Girirajan S, Rosenfeld JA, Vu TH, Baker C, et al. A copy number variation morbidity map of developmental delay. Nat Genet. 2011;43(9):838–46.
Lantieri F, Malacarne M, Gimelli S, Santamaria G, Coviello D, Ceccherini I. Custom Array comparative genomic hybridization: the importance of DNA quality, an expert eye, and variant validation. Int J Mol Sci. 2017;18:E609.
Kerstjens-Frederikse WS, Hofstra RM, van Essen AJ, Meijers JH, Buys CH. A Hirschsprung disease locus at 22q11? J Med Genet. 1999;36:221–4.
Serra A, Görgens H, Alhadad K, Ziegler A, Fitze G, Schackert HK. Analysis of RET, ZEB2, EDN3 and GDNF genomic rearrangements in 80 patients with Hirschsprung disease (using multiplex ligation-dependent probe amplification). Ann Hum Genet. 2009;73:147–51.
Tang CS, Sribudiani Y, Miao XP, de Vries AR, Burzynski G, So MT, et al. Fine mapping of the 9q31 Hirschsprung's disease locus. Hum Genet. 2010;127:675–83.
Lantieri F, Rydbeck H, Griseri P, Ceccherini I, Devoto M. Incorporating prior biological information in linkage studies increases power and limits multiple testing. BMC Proc. 2007;1(Suppl 1):S89.
Kruger RP, Aurandt J, Guan KL. Semaphorins command cells to move. Nat Rev Mol Cell Biol. 2005;6:789–800.
Heanue TA, Pachnis V. Expression profiling the developing mammalian enteric nervous system identifies marker and candidate Hirschsprung disease genes. Proc Natl Acad Sci U S A. 2006;103:6919–24.
Georgieva L, Dimitrova A, Ivanov D, Nikolov I, Williams NM, Grozeva D, et al. Support for neuregulin 1 as a susceptibility gene for bipolar disorder and schizophrenia. Biol Psychiatry. 2008;64:419–27.
Kim JH, Cheong HS, Sul JH, Seo JM, Kim DY, Oh JT, et al. A genome-wide association study identifies potential susceptibility loci for Hirschsprung disease. PLoS One. 2014;9:e110292.
Tang CS, Ngan ES, Tang WK, So MT, Cheng G, Miao XP, et al. Mutations in the NRG1 gene are associated with Hirschsprung disease. Hum Genet. 2012;131:67–76.
Virtanen VB, Salo PP, Cao J, Löf-Granström A, Milani L, Metspalu A, et al. Noncoding RET variants explain the strong association with Hirschsprung disease in patients without rare coding sequence variant. Eur J Med Genet. 2018;S1769–7212(18):30076–4.
Tang W, Li B, Xu X, Zhou Z, Wu W, Tang J, et al. Aberrant high expression of NRG1 gene in Hirschsprung disease. J Pediatr Surg. 2012;47:1694–8.
Di Zanni E, Adamo A, Belligni E, Lerone M, Martucciello G, Mattioli G, et al. Common PHOX2B poly-alanine contractions impair RET gene transcription, predisposing to Hirschsprung disease. Biochim Biophys Acta. 1863;2017:1770–7.
Chen CY, Gherzi R, Ong SE, Chan EL, Raijmakers R, Pruijn GJ, et al. AU binding proteins recruit the exosome to degrade ARE-containing mRNAs. Cell. 2001;107:451–64.
Fu M, Sato Y, Lyons-Warren A, Zhang B, Kane MA, Napoli JL, et al. Vitamin a facilitates enteric nervous system precursor migration by reducing Pten accumulation. Development. 2010;137:631–40.
Tsai YH, VanDussen KL, Sawey ET, Wade AW, Kasper C, Rakshit S, et al. ADAM10 regulates notch function in intestinal stem cells of mice. Gastroenterology. 2014;147:822–34.
Prasad RB, Lessmark A, Almgren P, Kovacs G, Hansson O, Oskolkov N, et al. Excess maternal transmission of variants in the THADA gene to offspring with type 2 diabetes. Diabetologia. 2016;59:1702–13.
Itoh H, Berthet M, Fressart V, Denjoy I, Maugenre S, Klug D, et al. Asymmetry of parental origin in long QT syndrome: preferential maternal transmission of KCNQ1 variants linked to channel dysfunction. Eur J Hum Genet. 2016;24:1160–6.
Vulto-van Silfhout AT, Hehir-Kwa JY, van Bon BW, Schuurs-Hoeijmakers JH, Meader S, Hellebrekers CJ, et al. Clinical significance of de novo and inherited copy-number variation. Hum Mutat. 2013;34:1679–87.
Wang B, Ji T, Zhou X, Wang J, Wang X, Wang J, et al. CNV analysis in Chinese children of mental retardation highlights a sex differentiation in parental contribution to de novo and inherited mutational burdens. Sci Rep. 2016;6:25954.
Jannot AS, Amiel J, Pelet A, Lantieri F, Fernandez RM, Verheij JB, et al. Male and female differential reproductive rate could explain parental transmission asymmetry of mutation origin in Hirschsprung disease. Eur J Hum Genet. 2012;20:917–20.
Wiszniewska J, Bi W, Shaw C, Stankiewicz P, Kang SH, Pursley AN, et al. Combined array CGH plus SNP genome analyses in a single assay for optimized clinical testing. Eur J Hum Genet. 2014;22:79–87.
Zhang X, Du R, Li S, Zhang F, Jin L, Wang H. Evaluation of copy number variation detection for a SNP array platform. BMC Bioinformatics. 2014;15:50.
La Cognata V, Morello G, Gentile G, D'Agata V, Criscuolo C, Cavalcanti F, et al. A customized high-resolution array-comparative genomic hybridization to explore copy number variations in Parkinson's disease. Neurogenetics. 2016;17:233–44.
Yoshikawa Y, Emi M, Hashimoto-Tamaoki T, Ohmuraya M, Sato A, Tsujimura T, et al. High-density array-CGH with targeted NGS unmask multiple noncontiguous minute deletions on chromosome 3p21 in mesothelioma. Proc Natl Acad Sci U S A. 2016;113:13432–7.
Conrad DF, Pinto D, Redon R, Feuk L, Gokcumen O, Zhang Y, et al. Origins and functional impact of copy number variation in the human genome. Nature. 2010;464:704–12.
We wish to thank the numerous patients, their families and referring physicians that have participated in these studies in the Gaslini Institute. We gratefully acknowledge the help of Mirko Travaglia in creating the Figures.
The authors would like to thank the Genome Aggregation Database (gnomAD) and the groups that provided exome and genome variant data to this resource. A full list of contributing groups can be found at https://gnomad.broadinstitute.org/about.
This study makes also 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 email@example.com. Funding for the project was provided by the Wellcome Trust. Those centres bear no responsibility for the further analysis or interpretation that we carried out in this study.
The present work was carried out with the financial support of the Italian Ministry of Health: grant GR-2011-02347381 (to A.P.P. and F.L.),“Cinque per mille”, and Ricerca Corrente (to the Gaslini Institute).
Ethics approval and consent to participate
Research ethics approval was obtained from the Gaslini Ethics Committee (provvedimento IGG 197/2012). Informed consents were obtained from all subjects, with the exception of those samples dating back to ‘90s.
Consent for publication
Not applicable as there is no identifiable individual data.
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
“Variants detected by software and visual inspection, according to selection criteria”. Whole raw data obtained by applying the a-CGH approach.
“Frequencies of CNVs detected in our sample compared to those for controls from gnomAd and DDD”. Comparison of all the aberrations detected by aCGH in our study and the frequencies reported for controls in two databases.
About this article
Cite this article
Lantieri, F., Gimelli, S., Viaggi, C. et al. Copy number variations in candidate genomic regions confirm genetic heterogeneity and parental bias in Hirschsprung disease. Orphanet J Rare Dis 14, 270 (2019). https://doi.org/10.1186/s13023-019-1205-3