- Research article
- Open Access
- Open Peer Review
CNV analysis in 169 patients with bladder exstrophy-epispadias complex
BMC Medical Geneticsvolume 17, Article number: 35 (2016)
The bladder exstrophy-epispadias complex (BEEC) represents the severe end of the congenital uro-rectal malformation spectrum. Initial studies have implicated rare copy number variations (CNVs), including recurrent duplications of chromosomal region 22q11.21, in BEEC etiology.
To detect further CNVs, array analysis was performed in 169 BEEC patients. Prior to inclusion, 22q11.21 duplications were excluded using multiplex ligation-dependent probe amplification.
Following the application of stringent filter criteria, seven rare CNVs were identified: n = 4, not present in 1307 in-house controls; n = 3, frequency of <0.002 in controls. These CNVs ranged from 1 to 6.08 Mb in size. To identify smaller CNVs, relaxed filter criteria used in the detection of previously reported BEEC associated chromosomal regions were applied. This resulted in the identification of six additional rare CNVs: n = 4, not present in 1307 in-house controls; n = 2, frequency <0.0008 in controls. These CNVs ranged from 0.03–0.08 Mb in size. For 10 of these 13 CNVs, confirmation and segregation analyses were performed (5 of maternal origin; 5 of paternal origin). Interestingly, one female with classic bladder extrophy carried a 1.18 Mb duplication of 22q11.1, a chromosomal region that is associated with cat eye syndrome.
A number of rare CNVs were identified in BEEC patients, and these represent candidates for further evaluation. Rare inherited CNVs may constitute modifiers of, or contributors to, multifactorial BEEC phenotypes.
The bladder-exstrophy-epispadias complex (BEEC; MIM %600057) represents the severe end of the uro-rectal malformation spectrum, and has a profound impact on continence and sexual and renal functions. The BEEC is an anterior wall midline defect with variable phenotypic expression. The phenotype ranges from epispadias (E) and classic bladder exstrophy (CBE), to the most severe form, cloacal exstrophy (CE). The latter is often referred to as the OEIS complex (omphalocele, exstrophy, imperforate anus, and spinal defects) [1–3]. Around one third of BEEC patients present with associated urological malformations, e.g. ectopic kidney, renal agenesis, and hydronephrosis. BEEC has an overall prevalence of 2.07 in 100,000 live births, and is more common in males . For the specific subtypes, estimated birth prevalences after the inclusion of terminated pregnancies are 1 in 117,000 in males and 1 in 484,000 in females for E; 1 in 37,000 for CBE; and 1 in 200,000 to 1 in 400,000 for CE [2–7]. Although BEEC can occur as part of a complex malformation syndrome, approximately 98.5 % of cases are isolated [4, 8, 9].
Extensive recent research has implicated both inherited and de novo genetic factors in BEEC etiology. These factors include common single nucleotide polymorphisms (SNPs) [10, 11] and rare larger genomic aberrations, such as chromosomal aberrations and copy number variations (CNVs) (Table 1 [12–28]). Genetic risk factors involving larger genomic regions typically show stronger individual effects on disease causation, and are more likely to have a de novo occurrence. The largest systematic array-based genome-wide CNV study of BEEC to date investigated 110 patients, and identified a de novo 0.9 Mb microduplication on chromosome 19p13.12 in a single CBE patient . Two earlier array-based genome-wide CNV studies, which included a total of 102 CBE patients, identified a duplication of 22q11.21 in four individuals [23, 25]. An additional case report described an array-based CNV analysis in a single CBE patient with a duplication of 22q11.21 . Following array-based genome-wide CNV studies, Draaken et al.  used a multiplex ligation-dependent probe amplification (MLPA) based approach to perform a regional screen for 22q11.21 duplications in 244 independent BEEC patients. The authors identified four novel duplications of variable size in four unrelated CBE patients.
The aim of the present study was to detect further BEEC-associated CNVs by performing a state-of-the-art genome-wide single nucleotide polymorphism (SNP)-array based analysis in 169 BEEC patients. Standardized filter criteria for a genome-wide approach were applied. To detect smaller CNVs, we then conducted a high-resolution analysis of genomic regions previously implicated in BEEC phenotypes using relaxed filter criteria.
Patients, controls, and DNA isolation
The present study was part of an ongoing multicenter investigation of the molecular genetic causes of BEEC. In an earlier study, our group performed a regional screen for 22q11.21 duplications in 244 previously unreported BEEC patients using MLPA . For 169 of these 244 patients, the DNA samples were suitable for genome-wide array-based CNV analysis. None of these 169 unrelated patients carried a 22q11.21 duplication. These 169 patients were therefore used as the cohort for the present analyses. Of these 169 patients (E, n = 17; CBE, n = 126; CE, n = 26), 109 were male and 60 were female. Patients were of Central European (n = 128); Spanish (n = 24); Italian (n = 7); Bosnian (n = 2); Croatian (n = 1); Portuguese (n = 1); and Turkish origin (n = 6). The patients were recruited by one of four experienced physicians. For 125 patients, DNA from both parents was available. All patients had a negative family history of BEEC. A total of 1,307 population-based controls were drawn from the Heinz Nixdorf Recall Study (HNR) . The study was approved by the ethics committee of the Medical Faculty of the University of Bonn. Written informed consent was obtained from all participants prior to inclusion.
Blood or saliva samples were obtained from patients, the population-based control group, and (when possible) from the parents of the present BEEC patients. Isolation of genomic DNA from blood was carried out using a Chemagic Magnetic Separation Module I (Chemagen, Baesweiler, Germany). Isolation of genomic DNA from saliva samples was carried out using the Oragene DNA Kit (DNA Genotek Inc., Kanata, Canada).
Array-based molecular karyotyping
For CNV detection, Illumina’s HumanOmniExpress-12 v1.1 microarray (San Diego, California, USA) was used. This comprises 719,665 markers, and has a median marker spacing of 2.1 Kb. The controls were genotyped using Illumina’s HumanOmniExpress-12 v1.0 microarray. The v1.1 and v1.0 microarrays have an overall marker overlap of 99.98 %. A DNA sample was considered to have failed if less than 99 % of the markers were called on the respective microarray.
CNVs were predicted using the program QuantiSNP (v2.2, www.well.ox.ac.uk/QuantiSNP/). This program applies the Objective-Bayes Hidden-Markov model . The following quality control (QC) criteria were used to exclude CNVs: (I) log Bayes factor <30; and (II) regions with <5 consecutive aberrant markers. In a subsequent step, samples which still presented with >10 CNVs (twice the standard deviation), or with a standard deviation of the log R ratio of >0.3, were excluded.
In the subsequent analysis, CNVs were excluded if they: (I) covered both equivocal telomeric regions and HLA-loci; (II) presented without gene content; (III) affected segmental duplications only; (IV) had a frequency in the present control cohort of >1 %; or (V) had >10 entries in the Database of Genomic Variants (DGV; http://dgv.tcag.ca/dgv/app/home).
Of the remaining CNVs, only those with a length of >1 Mb were considered. Filtering was performed using the package ‘intervals’, as implemented in R (R: A Language and Environment for Statistical Computing; http://www.R-project.org), and the UCSC Human Genome Browser assembly hg19 (http://genome.ucsc.edu/) .
Irrespective of the above filter criteria, a separate analysis was performed for regions/CNVs previously associated with BEEC phenotypes (see Table 1). This analysis included CNVs with a length of <1 Mb.
Following the application of the above mentioned criteria, all remaining CNVs were visually inspected using GenomeStudio genotyping module (v2011.1, www.illumina.com/). Possible candidate genes within these regions were evaluated for their expression in BEEC-relevant tissues during the respective critical embryonic time frame in mice (E9.5–14.5), as indicated in the Mouse Genome Informatics Database (MGI; http://www.informatics.jax.org/). Further information on the function of these candidate genes was obtained from the Uniprot Database (http://www.uniprot.org/), and via an NCBI literature research of PubMed and OMIM (http://www.ncbi.nlm.nih.gov).
Quantitative polymerase chain reaction (qPCR)
Confirmation of the remaining visually inspected CNVs was carried out using qPCR and SYBR Green or TaqMan. To detect the origin of each CNV, the respective parents were screened. The qPCR was performed on an ABI Prism 7900HT Fast Real-Time PCR System with SYBR Green (Applied Biosystems, Foster City, USA), as described elsewhere . All primer sequences are available upon request. Two CNVs were confirmed by qPCR using TaqMan Copy Number Assay Hs07478160_cn (Applied Biosystem, Foster City, CA, USA), and LightCycler® 480 Instrument II (384 well version; Roche Diagnostics GmbH, Mannheim, Germany). This probe locates at nucleotides chr1:1,398,345-1,398,369 (hg19). Copy numbers were calculated using the ΔΔCt method, as implemented in the CopyCaller Software (v2.0, http://www.appliedbiosystems.com/support/software/copycaller/).
EFNB1 sequence analysis
In all 25 female patients with CE, sequence analysis of the EFNB1 gene was performed. Of these patients, 23 had undergone previous microarray analysis, without detection of any disease-associated CNVs [22, 23]. All five exons with their adjacent splice sites were amplified by PCR (oligonucleotide sequences obtainable on request). For mutational analysis, PCR-amplified DNA products were subjected to direct automated sequencing (3130XL Genetic Analyzer, Applied Biosystems, Foster City, USA) and sequencing was performed for both strands of each amplicon. Nucleotides were numbered according to GenBank entry NM_004429.4.
Conventional cytogenetic analysis was performed using standard procedures (data not shown).
QuantiSNP array analysis in the initial 169 samples detected 13,767 putative CNVs. The samples of 18 patients failed to meet initial QC criteria, and were excluded from further analysis.
Using the primary filter criteria, six rare CNVs were identified (Table 2). All six reside in regions not yet implicated in BEEC. These six CNVs comprised five duplications and one deletion, and were identified in a total of seven patients. Examination of CNVs of <1 Mb in regions previously associated with BEEC (Table 1) revealed six additional CNVs in a further six patients (Table 2), and comprised deletions only. For three of these 13 CNVs, confirmation of their presence was impossible due to their partial overlap with segmental duplications.
The six larger CNVs included a 1.7 Mb duplication comprising seven RefSeq genes, which was detected at Xp22.31 in two CBE females (Pat. 15 and 16). These two identical duplications were of maternal origin (Pat. 15), and paternal origin (Pat. 16) respectively. Since the mother of Patient 15, and the father of Patient 16, were unaffected, it is unlikely that this CNV was a highly penetrant genetic causal factor. However, it may nonetheless contribute to disease development.
Furthermore, we detected a 6.08 Mb duplication at Xq11.1-q13.1 in a CBE female, which had been transmitted from the non-affected father. This female patient also presented with persistent foramen ovale and bilateral inguinal hernia. The duplication represents the largest CNV detected in the present study, and encompasses 43 RefSeq genes. This region contains the Ephrin B1 (EFNB1) gene. EFNB1 has previously been associated with craniofrontonasal syndrome (MIM #304110), a severe craniofacial midline defect that is only expressed in female carriers. Interestingly, two reports in the literature describe the co-occurrence of CE—the most severe form of the BEEC—and craniofrontonasal syndrome in two unrelated female patients . Therefore, although the female patient with the EFNB1 comprising duplication displayed CBE and not CE, the subsequent sequence analysis focused on all female CE patients in our cohort (n = 25 CE females). Sequence analysis of all five EFNB1 exons and their adjacent splice sites revealed no mutation in any of these 25 CE females. In one patient, an extremely rare but silent variant was detected in exon 5 (rs143341175, p. Ser281=). No minor allele frequency (MAF) for this variant is given in dbSNP. In four patients, a common polymorphism was detected in the 3′-UTR (rs2230423, C/T, MAF 0.1 in the European population).
Interestingly, one female CBE patient who additionally showed coxa valga (Pat. 14), carried a 1.18 Mb duplication on chromosome 22q11.1 (Fig. 1), which involves a region typically amplified in cat eye syndrome (CES; #115470). Karyotype analysis detected no supernumerary marker chromosome. Due to the partial overlap of this CNV with segmental duplications, qPCR could not be performed in the mother. As we did not had a paternal sample, it was not further investigated, whether this CNV had been inherited. The breakpoints did not coincide with the known low copy repeat (LCR) regions, as this CNV is proximal to LCR-A. CES conventional cytogenetic analysis from peripheral blood revealed a normal female karyotype (46,XX) in 30 metaphases. No supernumerary marker chromosome 22 was detected. The region affected by this duplication harbors six pseudogenes, and four genes encoding the transcripts for POTE ankyrin domain family member H (POTEH); olfactory receptor 11H1 (OR11H1); putative T-complex protein 1 subunit theta-like 2 (CCT8L2); and XK-related protein 3 (XKR3) (Fig. 1).
In a male with epispadias and penoscrotal transposition (Pat. 5), a 1 Mb duplication was detected at 4q26. This was of paternal origin. This duplication affects the translocation associated membrane protein 1-like 1 (TRAM1L1) gene and one microRNA (MIR1973). A CBE female (Pat. 6) was found to carry a 1.06 Mb duplication on chromosome region 5q22.2, which involved seven genes. The duplicated genes were adenomatous polyposis coli (APC); signal recognition particle 19 kDa (SRP19); U2 small nuclear ribonucleoprotein auxiliary factor 35 kDa subunit-related protein 1 (ZRSR1); receptor accessory protein 5 (REEP5); decapping mRNA 2 (DCP2); mutated in colorectal cancers (MCC); and testis-specific serine kinase 1B (TSSK1B).
Finally, a 1.68 Mb deletion at 13q33.1-q33.2 was detected in a male CBE patient (Pat. 11). This affected the D-amino acid oxidase activator (DAOA) gene and its antisense RNA (DAOA-AS1), as well as a long intergenic non-coding RNA (LINC00343).
In chromosomal regions previously associated with BEEC, we identified six deletions with an unknown effect. A deletion of chromosomal region 1p36.33 was observed in three patients. Two of these patients had CBE (Pat. 19 and 20, Table 2) and these individuals carried the same 40 Kb deletion. Patient 19 was male, and had inherited the CNV from his healthy father. Patient 20 was female, and had inherited the CNV from her healthy mother. A smaller, overlapping 30 Kb 1p36.33 microdeletion was detected in a male patient (Pat. 21) with epispadias. However, confirmation of this microdeletion was impossible due to the presence of a partially overlapping segmental duplication. The larger CNVs of this 1q36.33 region encompass two members of the family of mitochondrial AAA + −ATPase ATAD3 genes, i.e. ATAD3B and ATAD3C. The smaller CNV affects ATAD3B only.
In addition, we detected a small, maternally inherited deletion of chromosomal region 1q41 in a male CBE patient (Pat. 2). This deletion affects the Usherin (USH2A) gene, which is mutated in mild autosomal recessive usher syndrome 2A.
A maternally transmitted 10 Kb deletion at 9q34.2, which encompasses the ABO blood group gene, was detected in a CBE female (Pat. 9). A further CBE female patient (Pat. 22) carried an 80 kb deletion at 19q13.42. The genes affected by this 19q13.42 deletion encode two zinc finger proteins (ZNF761 and ZNF813), and the tropomyosin 3 pseudogene 9 (TPM3P9) non-coding RNA. However, confirmation of this CNV was impossible due to the presence of a partially overlapping segmental duplication.
The largest CNV detected in the present study was a paternally inherited 6.08 Mb duplication. This contains 43 RefSeq genes, and was found in a female CBE patient (Pat. 17) (Table 2). Previous authors have reported that an increased dosage of X-linked genes impacts normal neurocognitive development . The present 6.08 Mb duplication is listed in Decipher and the Database of Genomic Variants (DGV; http://dgv.tcag.ca/) as having shown association in several patients with neurocognitive impairment. In contrast to patients from these previous reports, the present female CBE patient was otherwise healthy and showed no neurocognitive impairment. Of the 43 RefSeq genes within the duplicated region, we considered EFNB1, encoding the ephrin receptor ligand ephrin-B1, to be a promising candidate gene. Although heterozygous mutations in this gene cause craniofrontonasal syndrome, the genetic defect causes no—or only mild—abnormalities in male carriers . If a duplication of EFNB1 exerts a similar effect, this might explain our observations in a healthy carrier father. This hypothesis is supported by the familial hypertelorism study of Babbs et al. , which identified a duplication of EFNB1 in three affected females. A duplication model led to an imbalance in murine Ephrin-B1 expression and abnormal cell sorting. Interestingly, around 10 % of mice—whether heterozygous, homozygous, or hemizygous for the conditional EfnB1 Lox allele—died within 24 h due to severe cleft palate . The literature also includes at least two reports of female patients with craniofrontonasal syndrome and CE, thus suggesting a common etiology [37, 38]. Moreover, from embryonic day 10.5, Efnb1 expression has been detected in the renal, urinary, and reproductive systems of the mouse . Research has also shown, that in humans another member of the family of ephrin receptor ligands, ephrin-B2, acts as a signaling molecule in uro-rectal development . However, we detected no potential causal EFNB1 variant for CE in the present cohort of 25 female patients, although the sample size may have been too small to detect rare causal mutational events. Furthermore, we cannot exclude the possibility that the method applied in the present study overlooked mutations in the promoter region, as-yet-unknown regulatory sequences, or non-coding regions.
In that context of ephrin receptor ligands, Walczak-Sztulpa et al.  also reported genital malformations in patients with deletions of 13q33-34, where EFNB2 is located. The authors suggested, that this chromosomal region harbors a gene for male genital development. Of note, the EFNB2 gene is directly adjacent to the 13q33.1-q33.2 deletion found in our male patient 11 (Table 2). EFNB2 has also been analyzed as a candidate gene in 13 patients with OEIS complex in the study by Vlangos et al.  however, no mutations were identified. Hence, further studies are warranted to investigate a potential dose effect of EFNB1 and EFNB2 in the etiology of BEEC, and to determine whether EFNB1 shares functions with EFNB2.
In previous studies, our group and others have generated strong evidence for the involvement of 22q11.21 duplications in the etiology of BEEC [23–26] and thus, the present finding of a 1.18 Mb duplication of the neighboring chromosomal region 22q11.1 (Fig. 1) in a female CBE patient (Pat. 14) is interesting, particularly since this CNV resides within the region typically amplified in cat eye syndrome (CES; MIM #115470). This duplication has been detected in numerous (apparently healthy) controls [42–44], and may thus represent a benign variant. However, it remains possible that it is causally related to the phenotype but with incomplete penetrance, as has been observed for the duplication 22q11.21 [23–25]. Interestingly, the CNVs on 22q11.1 and Xp22.31 were previously reported in a patient with OEIS complex . Although these CNVs are not identical, the duplicated 22q11.1 region is partially encompassed within the duplicated region in our patient 14, and the region on Xp22.31 is a much smaller deletion that is completely encompassed within the region, that is duplicated in our patients 15 and 16.
Of the six deletions identified in previously implicated (Table 1) chromosomal regions, the deletion of chromosomal region at 1p36.33 was detected in three patients. This region is involved in one of the most common terminal subtelomeric microdeletion syndromes, i.e. the 1p36 contiguous gene deletion syndrome, which typically presents with central nervous system involvement, cardiac defects, and dysmorphic craniofacial features . Two deletions of chromosomal region 1p36 in two different CBE patients involved ATAD3B and ATAD3C, while in the patient with epispadias only, ATAD3B was deleted. Interestingly, de novo deletions affecting this chromosomal band and deleting all three ATAD3 genes were detected in two previously reported patients with CE, i.e. the severest form of BEEC [12, 13]. However, the patients presented in these separate reports may actually represent one (the same) single patient as the description of the genital phenotype is strikingly similar. Nevertheless, this finding suggests an additive effect of ATAD3 genes in BEEC etiology. While Atad genes are expressed in early embryonic development [46, 47], Atad3 deficient mice usually die at E7.5, and heterozygotes display no urogenital anomalies . Moreover, around 70 heterozygous deletions, which involve all human ATAD3 genes, have been deposited in Decipher and the Database of Genomic Variants . With the exception of one individual with hypospadias, none of these patients presented with BEEC, thus rendering a contribution of these genes to disease formation unlikely.
Available data suggest that disease causing CNVs other than duplications of chromosomal region 22q11.21 are a rare cause of BEEC. Around 98.5 % of cases with BEEC are isolated, and yet many of the described CNVs in this study and by others are inherited from a supposedly healthy parent. This argues either, that non-penetrance is extremely common, or that the CNVs detected are unrelated. Further research is warranted to determine the role of the presently identified CNVs in BEEC etiology. Some of these rare inherited CNVs might at least constitute modifiers or contributors in a multifactorial mode of inheritance.
bladder exstrophy-epispadias complex
classic bladder exstrophy
copy number variation
minor allele frequency
polymerase chain reaction
single nucleotide polymorphism
Husmann DA, Vandersteen DR. Anatomy of the cloacal exstrophy. In: Gearhart JP, Mathews R, editors. The exstrophy-epispadias complex. New York: Kluwer Academic/Plenum Publishers; 1999. p. 199–206.
Hurst JA. Anterior abdominal wall defects. In: Firth HV, Hall JG, editors. Oxford desk reference. Clinical genetics. New York: Oxford University Press; 2012. p. 566.
Mahfuz I, Darling T, Wilkins S, White S, Cheng W. New insights into the pathogenesis of bladder exstrophy-epispadias complex. J Pediatr Urol. 2013;9(6):996–1005.
Sifel C, Correa A, Amar E, Bakker MK, Bermejo-Sánchez E, Bianca S, et al. Bladder exstrophy: An epidemiologic study from the International Clearinghouse for Birth Defects Surveillance and Research, and an overview of the literature. Am J Med Genet C Semin Med Genet. 2011;157C(4):321–32.
Hurwitz RS, Manzon GA, Ransley PG, Stephens FD. Cloacal exstrophy: a report of 34 cases. J Urol. 1987;138(4):1060–4.
Wiesel A, Queisser-Luft A, Clementi M, Bianca S, Stoll C, EUROSCAN Study Group. Prenatal detection of congenital renal malformations by fetal ultrasonographic examination: an analysis of 709,030 births in 12 European countries. Eur J Med Genet. 2005;48(2):131–44.
Ebert AK, Reutter H, Ludwig M, Rösch WH. The exstrophy-epispadias complex. Orphanet J Rare Dis. 2009;4(23):1–17.
Boyadjiev SA, Dodson JL, Radford CL, Ashrafi GH, Beaty TH, Mathews RI, et al. Clinical and molecular characterization of the bladder exstrophy-epispadias complex: analysis of 232 families. BJU Int. 2004;94(9):1337–43.
Gambhir L, Holler T, Muller M, Schott G, Vogt H, Detlefsen B, et al. Epidemiological survey of 214 families with bladder exstrophy-epispadias complex. J Urol. 2008;179(4):1539–43.
Reutter H, Draaken M, Pennimpede T, Wittler L, Brockschmidt FF, Ebert AK, et al. Genome-wide association study and mouse expression data identify a highly conserved 32 kb intergenic region between WNT3 and WNT9b as possible susceptibiliy locus for isolated classic exstrophy of the bladder. Hum Mol Genet. 2014;23(20):5536–44.
Draaken M, Knapp M, Pennimpede T, Schmidt JM, Ebert AK, Rösch W, et al. Genome-wide association study and meta-analysis identify ISL1 as genome-wide significant susceptibility gene for bladder exstrophy. PLoS Genet. 2015;11(3):e1005024.
Tannour-Louet M, Han S, Corbett ST, et al. Identification of de novo copy number variants associated with human disorders of sexual development. PLoS ONE. 2010;5(10):e15392.
El-Hattab AW, Skorupski JC, Hsieh MH, Breman AM, Patel A, Cheung SW, et al. OEIS complex associated with chromosome 1p36 deletion: A case report and review. Am J Med Genet A. 2010;152A(2):504–11.
Rotmensch S, Liberati M, Luo JS, Tallini G, Mahoney MJ, Hobbins JC. Prenatal diagnosis of a fetus with terminal deletion of chromosome 1 (q41). Prenat Diagn. 1991;11(11):867–73.
Zaki MS, Gillesseen-Kaesbach G, Vater I, Caliebe A, Siebert R, Kamel AK, et al. Bladder exstrophy and extreme genital anomaly in a patient with pure terminal 1q deletion: Expansion of phenotypic spectrum. Eur J Med Genet. 2012;55(1):43–8.
Jorgez CJ, Rosenfeld JA, Wilken NR, Vangapandu HV, Sahin A, Pham D, et al. Genitourinary defects associated with genomic deletions in 2p15 encompassing OTX1. PLoS ONE. 2014;9(9):e107028.
Kosaki R, Fukuhara Y, Kosuga M, Okuyama T, Kawashima N, Honna T, et al. OEIS complex with del (3) (q12.2q13.2). Am J Med Genet A. 2005;135A(2):224–6.
Nicholls G, Duffy PG. Anatomical correction of the exstrophy-epispadias complex: analysis of 34 patients. Br J Urol. 1998;82(6):865–9.
Battaglia A, Carey JC, Cederholm P, Viskochil DH, Brothman AR, Galasso C. Natural history of Wolf-Hirschhorn syndrome: experience with 15 cases. Pediatrics. 1999;103(4 Pt 1):830–6.
Chipail A, Constantinescu V, Covic M, Angheloni T. Phenotypic and cytogenetic analysis of an unusual malformative syndrome (trisomy 9 p+). Rev Pediatr Obstet Ginecol Pediatr. 1976;25(3):201–10.
Thauvin-Robinet C, Faivre L, Cusin V, Khau Van Kien P, Callier P, Parker KL, et al. Cloacal exstrophy in an infant with 9q34.1-qter deletion resulting from a de novo unbalanced translocation between chromosome 9q and Yq. Am J Med Genet A. 2004;126A(3):303–7.
Draaken M, Mughal SS, Pennimpede T, et al. Isolated bladder exstrophy associated with a de novo 0.9 Mb microduplication on chromosome 19p13.12. Birth Defects Res A Clin Mol Teratol. 2013;97(3):133–9.
Draaken M, Reutter H, Schramm C, Bartels E, Boemers TM, Ebert AK, et al. Microduplications at 22q11.21 are associated with non-syndromic classic bladder exstrophy. Eur J Med Genet. 2010;53(2):55–60.
Draaken M, Baudisch F, Timmermann B, Kuhl H, Kerick M, Proske J, et al. Classic bladder exstrophy: Frequent 22q11.21 duplications and definition of a 414 kb phenocritical region. Birth Defects Res A Clin Mol Teratol. 2014;100(6):512–7.
Lundin J, Soderhäll C, Lundén L, Hammarsjö A, White I, Schoumans J, et al. 22q11.2 microduplication in two patients with bladder exstrophy and hearing impairment. Eur J Med Genet. 2010;53(2):61–5.
Pierquin G, Uwineza A. 22q11.2 microduplication in a patient with bladder exstrophy and delayed psychomotor development. Eur J Hum Genet. 2012;20 suppl 1:89. abstr.
Soderhäll C, Lundin J, Lagerstedt-Robinson K, Grigelioniene G, Lackgren G, Clementson Kockum C, et al. A case with bladder exstrophy and unbalanced X chromosome rearrangement. Eur J Pediatr Surg. 2014;24(4):353–9.
Vlangos CN, Siuniak A, Ackley T, van Bokhoven H, Veltman J, Iyer R, et al. Comprehensive genetic analysis of OEIS complex reveals no evidence for a recurrent microdeletion or duplication. Am J Med Genet A. 2011;155A(1):38–49.
Schmermund A, Möhlenkamp S, Stang A, Grönemeyer D, Seibel R, Hirche H, et al. Assessment of clinically silent atherosclerotic disease and established and novel risk factors for predicting myocardial infarction and cardiac death in healthy middleaged subjects: rationale and design of the Heinz Nixdorf RECALL Study. Risk factors, evaluation of coronary calcium and lifestyle. Am Heart J. 2002;144(2):212–8.
Colella S, Yau C, Taylor JM, Mirza G, Butler H, Clouston P, et al. An Objective Bayes Hidden-Markov Model to detect and accurately map copy number variation using SNP genotyping data. Nucleic Acids Res. 2007;35(6):2013–25.
Meyer LR, Zweig AS, Hinrichs AS, Karolchik D, Kuhn RM, Wong M, et al. The UCSC Genome Browser database: extensions and updates 2013. Nucleic Acids Res. 2013;41(D1):D64–9.
Draaken M, Giesen CA, Kesselheim AL, Jabs R, Aretz S, Kugaudo M, et al. Maternal de novo triple mosaicism for two single OCRL nucleotide substitutions (c.1736A > T, c.1736A > G) in a Lowe syndrome family. Hum Genet. 2011;129(5):513–9.
Ludwig M, Ching B, Reutter H, Boyadjiev SA. Bladder exstrophy-epispadias complex. Birth Defects Res Part A Clin Mol Teratol. 2009;85(6):509–22.
Froyen G, van Esch H, Bauters M, Hollanders K, Frints SG, Vermeesch JR, et al. Detection of genomic copy number changes in patients with idiopathic mental retardation by high-resolution X-array-CGH: Important role for increased gene dosage of XLMR genes. Hum Mutat. 2007;28(10):1034–42.
Wieland I, Jakkubiczka S, Muschke P, Cohen M, Thiele H, Gerlach KL, et al. Mutations of the ephrin-B1 gene cause craniofrontonasal syndrome. Am J Hum Genet. 2004;74(6):1209–15.
Babbs C, Stewart HS, Williams LJ, Connell L, Goriely A, Twigg SRF, et al. Duplication of the EFNB1 gene in familial hypertelorism: Imbalance in Ephrin-B1 expression and abnormal phenotypes in humans and mice. Hum Mutat. 2011;32(8):930–8.
Neidlich JA, Whitaker LA, Natowicz M, McDonald DM, Schnur R, Zackai EH. Aglossia with congenital absence of the mandibular rami and other craniofacial abnormalities. Am J Med Genet. 1988;4 suppl 4:161–6.
Robin NH, Neidich JA, Bason LD, Whitaker LA, McDonald-McGinn D, Hunter J, et al. Frontonasal malformation and cloacal exstrophy: A previously unreported asociation. Am J Med Genet. 1996;61(1):75–8.
Visel A, Thaller C, Eichele G. GenePaint.org: an atlas of gene expression patterns in the mouse embryo. Nucleic Acids Res. 2004;32(Database issue):D552–6.
Dravis C, Yokoyama N, Chumley MJ, Cowan CA, Silvany RE, Shay J, et al. Bidirectional signaling mediated by ephrin-B2 and EphB2 controls urorectal development. Dev Biol. 2004;271(2):272–90.
Walczak-Sztulpa J, Wisniewska M, Latos-Bielenska A, Linné M, Kelbova C, Belitz B, et al. Chromosome deletions in 13q33-34: Report of four patients and review of the literature. Am J Med Genet A. 2008;146A(3):337–42.
Pinto D, Marshall C, Feuk L, Scherer SW. Copy-number variation in control population cohorts. Hum Mol Genet. 2007;16(2):R168–73.
Redon R, Ishikawa S, Fitch KR, Feuk L, Perry GH, Andrews TD, et al. Global variation in copy number in the human genome. Nature. 2006;444(7118):444–54.
Zogopoulos G, Ha KC, Naqib F, Moore S, Kim H, Montpetit A, et al. Germ-line DNA copy number variation frequencies in a large North American population. Hum Genet. 2007;122(3–4):345–53.
Giannikou K, Fryssira H, Oikonomakis V, Syrmou A, Kosma K, Tzetis M, et al. Further delineation of novel 1p36 rearrangements by array-CGH analysis: Narrowing the breakpoints and clarifying the “extended” phenotype. Gene. 2012;506(2):360–8.
Goller T, Seibold UK, Kremmer E, Voos W, Kolanus W. Atad3 function is essential for early post-implantation development in the mouse. PLoS Genet. 2013;8(1):e54799.
Li S, Lamarche F, Charton R, Delphin C, Gires O, Hubstenberger A, et al. Expression analysis of ATAD3 isoforms in rodent and human cell lines and tissues. Gene. 2014;535(1):60–9.
MacDonald JR, Ziman R, Yuen RK, Feuk L, Scherer SW. The database of genomic variants: a curated collection of structural variation in the human genome. Nucleic Acids Res. 2013;42(Database issue):D986–92.
HR is supported by grant RE 1723/1-1 from the German Research Foundation (Deutsche Forschungsgemeinschaft, DFG). MMN received support from the Alfried Krupp von Bohlen und Halbach-Stiftung, and is a member of the DFG-funded Excellence Cluster Immunosensation. We thank all patients and their parents for their cooperation. We acknowledge the assistance of the German, Spanish, and Italian self-help organizations in the recruitment of individuals with BEEC, and thank Enrika Bartels and Sara Trbojevic for their work on the patient registry. We thank Peter Teßmann and Christina Ergang for their excellent technical assistance. The study used data generated by the Decipher Consortium (http://decipher.sanger.ac.uk/).
The authors have no commercial associations, such as consultancies, stock ownership, or other equity interests or patent licensing arrangements, to declare.
The authors declare that they have no competing interests.
HR, MMN, and ML designed the study. AKE, WR, RS, TMB, KH, CM, WFJF, AB, NM, MDG, and SM obtained informed consent and collected clinical data and blood samples from patients. CVL, AH, RZ, FM, and MD performed the genetic studies. MD, HR, and ML coordinated the research and drafted the manuscript. MMN, HR, and ML revised the manuscript for important intellectual content. All authors have read and approved the final manuscript and its submission for publication.