A syndromic neurodevelopmental disorder caused by rare variants in PPFIA3.

PPFIA3 encodes the protein-tyrosine phosphatase, receptor-type, F-polypeptide-interacting-protein-alpha-3 (PPFIA3), which is a member of the LAR-protein-tyrosine phosphatase-interacting-protein (liprin) family involved in synapse formation and function, synaptic vesicle transport, and presynaptic active zone assembly. The protein structure and function are evolutionarily well conserved, but human diseases related to PPFIA3 dysfunction are not yet reported in OMIM. Here, we report 20 individuals with rare PPFIA3 variants (19 heterozygous and 1 compound heterozygous) presenting with developmental delay, intellectual disability, hypotonia, dysmorphisms, microcephaly or macrocephaly, autistic features, and epilepsy with reduced penetrance. Seventeen unique PPFIA3 variants were detected in 18 families. To determine the pathogenicity of PPFIA3 variants in vivo, we generated transgenic fruit flies producing either human wild-type (WT) PPFIA3 or five missense variants using GAL4-UAS targeted gene expression systems. In the fly overexpression assays, we found that the PPFIA3 variants in the region encoding the N-terminal coiled-coil domain exhibited stronger phenotypes compared to those affecting the C-terminal region. In the loss-of-function fly assay, we show that the homozygous loss of fly Liprin-α leads to embryonic lethality. This lethality is partially rescued by the expression of human PPFIA3 WT, suggesting human PPFIA3 function is partially conserved in the fly. However, two of the tested variants failed to rescue the lethality at the larval stage and one variant failed to rescue lethality at the adult stage. Altogether, the human and fruit fly data reveal that the rare PPFIA3 variants are dominant-negative loss-of-function alleles that perturb multiple developmental processes and synapse formation.

Msx1 deficiency interacts with hypoxia and induces a morphogenetic regulation during mouse lip development.

Nonsyndromic clefts of the lip and palate are common birth defects resulting from gene-gene and gene-environment interactions. Mutations in human MSX1 have been linked to orofacial clefting and we show here that Msx1 deficiency causes a growth defect of the medial nasal process (Mnp) in mouse embryos. Although this defect alone does not disrupt lip formation, Msx1-deficient embryos develop a cleft lip when the mother is transiently exposed to reduced oxygen levels or to phenytoin, a drug known to cause embryonic hypoxia. In the absence of interacting environmental factors, the Mnp growth defect caused by Msx1 deficiency is modified by a Pax9-dependent 'morphogenetic regulation', which modulates Mnp shape, rescues lip formation and involves a localized abrogation of Bmp4-mediated repression of Pax9 Analyses of GWAS data revealed a genome-wide significant association of a Gene Ontology morphogenesis term (including assigned roles for MSX1, MSX2, PAX9, BMP4 and GREM1) specifically for nonsyndromic cleft lip with cleft palate. Our data indicate that MSX1 mutations could increase the risk for cleft lip formation by interacting with an impaired morphogenetic regulation that adjusts Mnp shape, or through interactions that inhibit Mnp growth.

Next-generation phenotyping contributing to the identification of a 4.7 kb deletion in KANSL1 causing Koolen-de Vries syndrome.

Next-generation phenotyping (NGP) is an application of advanced methods of computer vision on medical imaging data such as portrait photos of individuals with rare disorders. NGP on portraits results in gestalt scores that can be used for the selection of appropriate genetic tests, and for the interpretation of the molecular data. Here, we report on an exceptional case of a young girl that was presented at the age of 8 and 15 and enrolled in NGP diagnostics on the latter occasion. The girl had clinical features associated with Koolen-de Vries syndrome (KdVS) and a suggestive facial gestalt. However, chromosomal microarray (CMA), Sanger sequencing, multiplex ligation-dependent probe analysis (MLPA), and trio exome sequencing remained inconclusive. Based on the highly indicative gestalt score for KdVS, the decision was made to perform genome sequencing to also evaluate noncoding variants. This analysis revealed a 4.7 kb de novo deletion partially affecting intron 6 and exon 7 of the KANSL1 gene. This is the smallest reported structural variant to date for this phenotype. The case illustrates how NGP can be integrated into the iterative diagnostic process of test selection and interpretation of sequencing results.

Truncating Mutations in UBAP1 Cause Hereditary Spastic Paraplegia.

The diagnostic gap for rare neurodegenerative diseases is still considerable, despite continuous advances in gene identification. Many novel Mendelian genes have only been identified in a few families worldwide. Here we report the identification of an autosomal-dominant gene for hereditary spastic paraplegia (HSP) in 10 families that are of diverse geographic origin and whose affected members all carry unique truncating changes in a circumscript region of UBAP1 (ubiquitin-associated protein 1). HSP is a neurodegenerative disease characterized by progressive lower-limb spasticity and weakness, as well as frequent bladder dysfunction. At least 40% of affected persons are currently undiagnosed after exome sequencing. We identified pathological truncating variants in UBAP1 in affected persons from Iran, USA, Germany, Canada, Spain, and Bulgarian Roma. The genetic support ranges from linkage in the largest family (LOD = 8.3) to three confirmed de novo mutations. We show that mRNA in the fibroblasts of affected individuals escapes nonsense-mediated decay and thus leads to the expression of truncated proteins; in addition, concentrations of the full-length protein are reduced in comparison to those in controls. This suggests either a dominant-negative effect or haploinsufficiency. UBAP1 links endosomal trafficking to the ubiquitination machinery pathways that have been previously implicated in HSPs, and UBAP1 provides a bridge toward a more unified pathophysiology.

Extending the allelic spectrum at noncoding risk loci of orofacial clefting.

Genome-wide association studies (GWAS) have generated unprecedented insights into the genetic etiology of orofacial clefting (OFC). The moderate effect sizes of associated noncoding risk variants and limited access to disease-relevant tissue represent considerable challenges for biological interpretation of genetic findings. As rare variants with stronger effect sizes are likely to also contribute to OFC, an alternative approach to delineate pathogenic mechanisms is to identify private mutations and/or an increased burden of rare variants in associated regions. This report describes a framework for targeted resequencing at selected noncoding risk loci contributing to nonsyndromic cleft lip with/without cleft palate (nsCL/P), the most frequent OFC subtype. Based on GWAS data, we selected three risk loci and identified candidate regulatory regions (CRRs) through the integration of credible SNP information, epigenetic data from relevant cells/tissues, and conservation scores. The CRRs (total 57 kb) were resequenced in a multiethnic study population (1061 patients; 1591 controls), using single-molecule molecular inversion probe technology. Combining evidence from in silico variant annotation, pedigree- and burden analyses, we identified 16 likely deleterious rare variants that represent new candidates for functional studies in nsCL/P. Our framework is scalable and represents a promising approach to the investigation of additional congenital malformations with multifactorial etiology.

Evaluating shared genetic influences on nonsyndromic cleft lip/palate and oropharyngeal neoplasms.

It has been hypothesised that nonsyndromic cleft lip/palate (nsCL/P) and cancer may share aetiological risk factors. Population studies have found inconsistent evidence for increased incidence of cancer in nsCL/P cases, but several genes (e.g., CDH1, AXIN2) have been implicated in the aetiologies of both phenotypes. We aimed to evaluate shared genetic aetiology between nsCL/P and oral cavity/oropharyngeal cancers (OC/OPC), which affect similar anatomical regions. Using a primary sample of 5,048 OC/OPC cases and 5,450 controls of European ancestry and a replication sample of 750 cases and 336,319 controls from UK Biobank, we estimate genetic overlap using nsCL/P polygenic risk scores (PRS) with Mendelian randomization analyses performed to evaluate potential causal mechanisms. In the primary sample, we found strong evidence for an association between a nsCL/P PRS and increased odds of OC/OPC (per standard deviation increase in score, odds ratio [OR]: 1.09; 95% confidence interval [CI]: 1.04, 1.13; p = .000053). Although confidence intervals overlapped with the primary estimate, we did not find confirmatory evidence of an association between the PRS and OC/OPC in UK Biobank (OR 1.02; 95% CI: 0.95, 1.10; p = .55). Mendelian randomization analyses provided evidence that major nsCL/P risk variants are unlikely to influence OC/OPC. Our findings suggest possible shared genetic influences on nsCL/P and OC/OPC.

Investigating the shared genetics of non-syndromic cleft lip/palate and facial morphology.

There is increasing evidence that genetic risk variants for non-syndromic cleft lip/palate (nsCL/P) are also associated with normal-range variation in facial morphology. However, previous analyses are mostly limited to candidate SNPs and findings have not been consistently replicated. Here, we used polygenic risk scores (PRS) to test for genetic overlap between nsCL/P and seven biologically relevant facial phenotypes. Where evidence was found of genetic overlap, we used bidirectional Mendelian randomization (MR) to test the hypothesis that genetic liability to nsCL/P is causally related to implicated facial phenotypes. Across 5,804 individuals of European ancestry from two studies, we found strong evidence, using PRS, of genetic overlap between nsCL/P and philtrum width; a 1 S.D. increase in nsCL/P PRS was associated with a 0.10 mm decrease in philtrum width (95% C.I. 0.054, 0.146; P = 2x10-5). Follow-up MR analyses supported a causal relationship; genetic variants for nsCL/P homogeneously cause decreased philtrum width. In addition to the primary analysis, we also identified two novel risk loci for philtrum width at 5q22.2 and 7p15.2 in our Genome-wide Association Study (GWAS) of 6,136 individuals. Our results support a liability threshold model of inheritance for nsCL/P, related to abnormalities in development of the philtrum.

Impact of Maternal Smoking on Nonsyndromic Clefts: Sex-Specific Associations With Side and Laterality.

OBJECTIVE: To compare the incidence of right-sided versus left-sided, and unilateral versus bilateral, nonsyndromic clefting in the affected offspring of smoking and nonsmoking mothers. DESIGN: Self-report data on periconceptual and first trimester smoking behavior were collected from 842 mothers of children with nonsyndromic orofacial clefting. Differences in the incidence of left- versus right-sided clefts, and of unilateral versus bilateral clefts, were analyzed between the children of smoking and nonsmoking mothers. SETTING: Interviews and clinical examinations took place at 8 specialist centers in Germany. PATIENTS AND PARTICIPANTS: Children with nonsyndromic clefts were recruited during the course of surgical or orthodontic treatment, or within the context of the annual control consultation. Patients with cleft palate only or missing data were excluded. The final cohort comprised 842 patients (540 males and 302 females) with unilateral or bilateral clefts. The respective mothers were interviewed. MAIN OUTCOME MEASURE: Side and laterality of nonsyndromic clefts were the main outcome measures. RESULTS: Children of smoking mothers more often had right-sided clefts than children of nonsmoking mothers (42% right-sided clefts in children of smoking mothers vs 31% of nonsmoking mothers). Children of smoking mothers more often had bilateral clefts than children of nonsmoking mothers (35% bilateral clefts in children of smoking mothers vs 29% of nonsmoking mothers). Sex-specific analyses confirmed substantially and statistically significant associations only for girls. CONCLUSIONS: The results suggest that maternal smoking is a sex-specific, exogenous determinant of laterality and side in nonsyndromic clefts.

Nine newly identified individuals refine the phenotype associated with MYT1L mutations.

Both point mutations and deletions of the MYT1L gene as well as microdeletions of chromosome band 2p25.3 including MYT1L are associated with intellectual disability, obesity, and behavioral problems. Thus, MYT1L is assumed to be the-at least mainly-causative gene in the 2p25.3 deletion syndrome. Here, we present comprehensive descriptions of nine novel individuals bearing MYT1L mutations; most of them single nucleotide variants (SNVs). This increases the number of known individuals with causative deletions or SNVs of MYT1L to 51. Since eight of the nine novel patients bear mutations affecting MYT1L only, the total number of such individuals now nearly equals the number of individuals with larger microdeletions affecting additional genes, allowing for a comprehensive phenotypic comparison of these two patient groups. For example, 55% of the individuals with mutations affecting MYT1L only were overweight or obese as compared to 86% of the individuals with larger microdeletions. A similar trend was observed regarding short stature with 5 versus 35%, respectively. However, these differences were nominally significant only after correction for multiple testing, further supporting the hypothesis that MYT1L haploinsufficiency is central to the 2p25.3 deletion phenotype. Most importantly, the large number of individuals with MYT1L mutations presented and reviewed here allowed for the delineation of a more comprehensive clinical picture. Seizures, postnatal short stature, macrocephaly, and microcephaly could be shown to be over-represented among individuals with MYT1L mutations.

Imputation of orofacial clefting data identifies novel risk loci and sheds light on the genetic background of cleft lip ± cleft palate and cleft palate only.

Nonsyndromic cleft lip with or without cleft palate (nsCL/P) is among the most common human birth defects with multifactorial etiology. Here, we present results from a genome-wide imputation study of nsCL/P in which, after adding replication cohort data, four novel risk loci for nsCL/P are identified (at chromosomal regions 2p21, 14q22, 15q24 and 19p13). On a systematic level, we show that the association signals within this high-density dataset are enriched in functionally-relevant genomic regions that are active in both human neural crest cells (hNCC) and mouse embryonic craniofacial tissue. This enrichment is also detectable in hNCC regions primed for later activity. Using GCTA analyses, we suggest that 30% of the estimated variance in risk for nsCL/P in the European population can be attributed to common variants, with 25.5% contributed to by the 24 risk loci known to date. For each of these, we identify credible SNPs using a Bayesian refinement approach, with two loci harbouring only one probable causal variant. Finally, we demonstrate that there is no polygenic component of nsCL/P detectable that is shared with nonsyndromic cleft palate only (nsCPO). Our data suggest that, while common variants are strongly contributing to risk for nsCL/P, they do not seem to be involved in nsCPO which might be more often caused by rare deleterious variants. Our study generates novel insights into both nsCL/P and nsCPO etiology and provides a systematic framework for research into craniofacial development and malformation.

Sequencing the GRHL3 Coding Region Reveals Rare Truncating Mutations and a Common Susceptibility Variant for Nonsyndromic Cleft Palate.

Nonsyndromic cleft lip with/without cleft palate (nsCL/P) and nonsyndromic cleft palate only (nsCPO) are the most frequent subphenotypes of orofacial clefts. A common syndromic form of orofacial clefting is Van der Woude syndrome (VWS) where individuals have CL/P or CPO, often but not always associated with lower lip pits. Recently, ∼5% of VWS-affected individuals were identified with mutations in the grainy head-like 3 gene (GRHL3). To investigate GRHL3 in nonsyndromic clefting, we sequenced its coding region in 576 Europeans with nsCL/P and 96 with nsCPO. Most strikingly, nsCPO-affected individuals had a higher minor allele frequency for rs41268753 (0.099) than control subjects (0.049; p = 1.24 × 10(-2)). This association was replicated in nsCPO/control cohorts from Latvia, Yemen, and the UK (pcombined = 2.63 × 10(-5); ORallelic = 2.46 [95% CI 1.6-3.7]) and reached genome-wide significance in combination with imputed data from a GWAS in nsCPO triads (p = 2.73 × 10(-9)). Notably, rs41268753 is not associated with nsCL/P (p = 0.45). rs41268753 encodes the highly conserved p.Thr454Met (c.1361C>T) (GERP = 5.3), which prediction programs denote as deleterious, has a CADD score of 29.6, and increases protein binding capacity in silico. Sequencing also revealed four novel truncating GRHL3 mutations including two that were de novo in four families, where all nine individuals harboring mutations had nsCPO. This is important for genetic counseling: given that VWS is rare compared to nsCPO, our data suggest that dominant GRHL3 mutations are more likely to cause nonsyndromic than syndromic CPO. Thus, with rare dominant mutations and a common risk variant in the coding region, we have identified an important contribution for GRHL3 in nsCPO.

PEDIA: prioritization of exome data by image analysis.

PURPOSE: Phenotype information is crucial for the interpretation of genomic variants. So far it has only been accessible for bioinformatics workflows after encoding into clinical terms by expert dysmorphologists. METHODS: Here, we introduce an approach driven by artificial intelligence that uses portrait photographs for the interpretation of clinical exome data. We measured the value added by computer-assisted image analysis to the diagnostic yield on a cohort consisting of 679 individuals with 105 different monogenic disorders. For each case in the cohort we compiled frontal photos, clinical features, and the disease-causing variants, and simulated multiple exomes of different ethnic backgrounds. RESULTS: The additional use of similarity scores from computer-assisted analysis of frontal photos improved the top 1 accuracy rate by more than 20-89% and the top 10 accuracy rate by more than 5-99% for the disease-causing gene. CONCLUSION: Image analysis by deep-learning algorithms can be used to quantify the phenotypic similarity (PP4 criterion of the American College of Medical Genetics and Genomics guidelines) and to advance the performance of bioinformatics pipelines for exome analysis.

Male infant with paternal uniparental diploidy mosaicism and a 46,XX/46,XY karyotype.

A male patient with mosaic paternal uniparental diploidy (PUD) is presented. After birth, the patient presented with hypoglycemia, hemihypertrophy, umbilical hernia, and hepatomegaly. Afterward pancreatic hypertrophy, liver hemangiomas, and cysts were detected sonographically. At the age of 3.5 months, hepatoblastoma was diagnosed. To investigate suspected Beckwith-Wiedemann syndrome (BWS), extensive genetic analyses were performed using DNA from chorionic villus sampling, amniocentesis, and peripheral blood lymphocytes (chromosome analysis, methylation-specific multiplex ligation-dependent probe amplification assays, microsatellite analyses, and single nucleotide polymorphism array analysis). These analyses led to the detection of mosaic PUD. In peripheral blood lymphocytes, a male cell line (46,XY[27]/46,XX[5]) predominated, suggesting a mixture of uniparental isodisomy and heterodisomy. The genetic analyses suggest that the mosaic PUD status was attributable to fertilization of an oocyte by two sperms, with subsequent triploidy rescue giving rise to haploidy, which in turn was rescued. Notably, in the majority of the 28 mosaic PUD patients reported to date, BWS was initially suspected. Mosaic PUD status is associated with a higher risk for a broad range of malignant and benign tumors than in BWS. As tumors can also occur after childhood surveillance into adolescence is indicated. Mosaic PUD must therefore be considered in patients with suspected BWS.

Cleft Palate as Distinguishing Feature in a Patient with GABRB3 Epileptic Encephalopathy.

Mutations in GABAA-receptor subunit genes are associated with a heterogeneous spectrum of epilepsies. Patients with epilepsy caused by mutations in a specific GABAA-receptor (GABRA3) occasionally present with orofacial dysmorphism (e.g., cleft palates). While cleft palates have been described in Gabrb3 knockout mice and in humans with GABRB3 variants without epilepsy, the specific combination of epilepsy and cleft palate in humans with GABRB3 mutations has not yet been reported.We describe a patient with epileptic encephalopathy (EE) who presented with therapy-refractory neonatal-onset myoclonic seizures and severe developmental delay. Electroencephalogram showed burst suppression pattern at neonatal age and hypsarrhythmia at infantile age. Initial magnetic resonance imaging was unremarkable. As he additionally presented with a cleft palate, we were curious whether cleft palate and EE had the same genetic origin. Whole exome sequencing of the index patient revealed a novel pathogenic heterozygous de novo mutation in GABRB3 (c.899T > C; p.I300T). In consistency with Gabrb3 knockout mice data, this is the first report of cleft palate in a patient with GABRB3 associated EE.We suggest to add cleft palate to the phenotypic GABRB3 spectrum and to screen for mutations in GABAA-receptors in patients with EE and orofacial dysmorphism.

Meta-analysis Reveals Genome-Wide Significance at 15q13 for Nonsyndromic Clefting of Both the Lip and the Palate, and Functional Analyses Implicate GREM1 As a Plausible Causative Gene.

Nonsyndromic orofacial clefts are common birth defects with multifactorial etiology. The most common type is cleft lip, which occurs with or without cleft palate (nsCLP and nsCLO, respectively). Although genetic components play an important role in nsCLP, the genetic factors that predispose to palate involvement are largely unknown. In this study, we carried out a meta-analysis on genetic and clinical data from three large cohorts and identified strong association between a region on chromosome 15q13 and nsCLP (P = 8.13 × 10(-14) for rs1258763; relative risk (RR): 1.46, 95% confidence interval (CI): 1.32-1.61)) but not nsCLO (P = 0.27; RR: 1.09 (0.94-1.27)). The 5 kb region of strongest association maps downstream of Gremlin-1 (GREM1), which encodes a secreted antagonist of the BMP4 pathway. We show during mouse embryogenesis, Grem1 is expressed in the developing lip and soft palate but not in the hard palate. This is consistent with genotype-phenotype correlations between rs1258763 and a specific nsCLP subphenotype, since a more than two-fold increase in risk was observed in patients displaying clefts of both the lip and soft palate but who had an intact hard palate (RR: 3.76, CI: 1.47-9.61, Pdiff<0.05). While we did not find lip or palate defects in Grem1-deficient mice, wild type embryonic palatal shelves developed divergent shapes when cultured in the presence of ectopic Grem1 protein (P = 0.0014). The present study identified a non-coding region at 15q13 as the second, genome-wide significant locus specific for nsCLP, after 13q31. Moreover, our data suggest that the closely located GREM1 gene contributes to a rare clinical nsCLP entity. This entity specifically involves abnormalities of the lip and soft palate, which develop at different time-points and in separate anatomical regions.

Genetic testing in patients with retinitis pigmentosa: Features of unsolved cases.

IMPORTANCE: Uncommon characteristics in genetically unsolved retinitis pigmentosa (RP) patients may indicate an incorrect clinical diagnosis or as yet unknown genetic causes resulting in specific retinal phenotypes. The diagnostic yield of targeted next-generation sequencing may be increased by a reasonable preselection of RP-patients. BACKGROUND: To systematically evaluate and compare features of genetically solved and unsolved RP-patients. DESIGN: Retrospective, observational study. PARTICIPANTS: One-hundred and twelve consecutive RP-patients who underwent extensive molecular genetic analysis. METHODS: Characterization of patients based on multimodal imaging and medical history. MAIN OUTCOME MEASURES: Differences between genetically solved and unsolved RP-patients. RESULTS: Compared to genetically solved patients (n = 77), genetically unsolved patients (n = 35) more frequently had an age of disease-onset above 30 years (60% vs 8%; P < 0.0001), showed atypical fundus features (49% vs 8%; P < 0. 0001) and indicators for phenocopies (eg, autoimmune diseases) (17% vs 0%; P < 0. 001). Evidence for a particular inheritance pattern was less common (20% vs 49%; P < 0. 01). The diagnostic yield was 84% (71/85) in patients with first symptoms below 30 years-of-age, compared to 69% (77/112) in the overall cohort. The other selection criteria alone or in combination resulted in limited further increase of the diagnostic yield (up to 89%) while excluding considerably more patients (up to 56%) from genetic testing. CONCLUSIONS AND RELEVANCE: The medical history and retinal phenotype differ between genetically solved and a subgroup of unsolved RP-patients, which may reflect undetected genotypes or retinal conditions mimicking RP. Patient stratification may inform on the individual likelihood of identifying disease-causing mutations and may impact patient counselling.

Whole exome association of rare deletions in multiplex oral cleft families.

By sequencing the exomes of distantly related individuals in multiplex families, rare mutational and structural changes to coding DNA can be characterized and their relationship to disease risk can be assessed. Recently, several rare single nucleotide variants (SNVs) were associated with an increased risk of nonsyndromic oral cleft, highlighting the importance of rare sequence variants in oral clefts and illustrating the strength of family-based study designs. However, the extent to which rare deletions in coding regions of the genome occur and contribute to risk of nonsyndromic clefts is not well understood. To identify putative structural variants underlying risk, we developed a pipeline for rare hemizygous deletions in families from whole exome sequencing and statistical inference based on rare variant sharing. Among 56 multiplex families with 115 individuals, we identified 53 regions with one or more rare hemizygous deletions. We found 45 of the 53 regions contained rare deletions occurring in only one family member. Members of the same family shared a rare deletion in only eight regions. We also devised a scalable global test for enrichment of shared rare deletions.

Genome-wide association study and meta-analysis identify ISL1 as genome-wide significant susceptibility gene for bladder exstrophy.

The bladder exstrophy-epispadias complex (BEEC) represents the severe end of the uro-rectal malformation spectrum, and is thought to result from aberrant embryonic morphogenesis of the cloacal membrane and the urorectal septum. The most common form of BEEC is isolated classic bladder exstrophy (CBE). To identify susceptibility loci for CBE, we performed a genome-wide association study (GWAS) of 110 CBE patients and 1,177 controls of European origin. Here, an association was found with a region of approximately 220kb on chromosome 5q11.1. This region harbors the ISL1 (ISL LIM homeobox 1) gene. Multiple markers in this region showed evidence for association with CBE, including 84 markers with genome-wide significance. We then performed a meta-analysis using data from a previous GWAS by our group of 98 CBE patients and 526 controls of European origin. This meta-analysis also implicated the 5q11.1 locus in CBE risk. A total of 138 markers at this locus reached genome-wide significance in the meta-analysis, and the most significant marker (rs9291768) achieved a P value of 2.13 × 10-12. No other locus in the meta-analysis achieved genome-wide significance. We then performed murine expression analyses to follow up this finding. Here, Isl1 expression was detected in the genital region within the critical time frame for human CBE development. Genital regions with Isl1 expression included the peri-cloacal mesenchyme and the urorectal septum. The present study identified the first genome-wide significant locus for CBE at chromosomal region 5q11.1, and provides strong evidence for the hypothesis that ISL1 is the responsible candidate gene in this region.