Whole genome sequencing reveals epistasis effects within RET for Hirschsprung disease.

Common variants in RET and NRG1 have been associated with Hirschsprung disease (HSCR), a congenital disorder characterised by incomplete innervation of distal gut, in East Asian (EA) populations. However, the allelic effects so far identified do not fully explain its heritability, suggesting the presence of epistasis, where effect of one genetic variant differs depending on other (modifier) variants. Few instances of epistasis have been documented in complex diseases due to modelling complexity and data challenges. We proposed four epistasis models to comprehensively capture epistasis for HSCR between and within RET and NRG1 loci using whole genome sequencing (WGS) data in EA samples. 65 variants within the Topologically Associating Domain (TAD) of RET demonstrated significant epistasis with the lead enhancer variant (RET+3; rs2435357). These epistatic variants formed two linkage disequilibrium (LD) clusters represented by rs2506026 and rs2506028 that differed in minor allele frequency and the best-supported epistatic model. Intriguingly, rs2506028 is in high LD with one cis-regulatory variant (rs2506030) highlighted previously, suggesting that detected epistasis might be mediated through synergistic effects on transcription regulation of RET. Our findings demonstrated the advantages of WGS data for detecting epistasis, and support the presence of interactive effects of regulatory variants in RET for HSCR.

NGS4THAL, a One-Stop Molecular Diagnosis and Carrier Screening Tool for Thalassemia and Other Hemoglobinopathies by Next-Generation Sequencing.

Thalassemia is one of the most common genetic diseases and a major health threat worldwide. Accurate, efficient, and scalable analysis of next-generation sequencing (NGS) data is much needed for its molecular diagnosis and carrier screening. We developed NGS4THAL, a bioinformatics analysis pipeline analyzing NGS data to detect pathogenic variants for thalassemia and other hemoglobinopathies. NGS4THAL realigns ambiguously mapped NGS reads derived from the homologous Hb gene clusters for accurate detection of point mutations and small insertions/deletions. It uses a combination of complementary structural variant (SV) detection tools and an in-house database of control data containing specific SVs to achieve accurate detection of the complex SV types. Detected variants are matched with those in HbVar (A Database of Human Hemoglobin Variants and Thalassemia Mutations), allowing recognition of known pathogenic variants, including disease modifiers. Tested on simulation data, NGS4THAL achieved high sensitivity and specificity. For targeted NGS sequencing data from samples with laboratory-confirmed pathogenic Hb variants, it achieved 100% detection accuracy. Application of NGS4THAL on whole genome sequencing data from unrelated studies revealed thalassemia mutation carrier rates for Hong Kong Chinese and Northern Vietnamese that were consistent with previous reports. NGS4THAL is a highly accurate and efficient molecular diagnosis tool for thalassemia and other hemoglobinopathies based on tailored analysis of NGS data and may be scaled for population carrier screening.

Sequencing of a Chinese tetralogy of Fallot cohort reveals clustering mutations in myogenic heart progenitors.

Tetralogy of Fallot (TOF) is the most common cyanotic heart defect, yet the underlying genetic mechanisms remain poorly understood. Here, we performed whole-genome sequencing analysis on 146 nonsyndromic TOF parent-offspring trios of Chinese ethnicity. Comparison of de novo variants and recessive genotypes of this data set with data from a European cohort identified both overlapping and potentially novel gene loci and revealed differential functional enrichment between cohorts. To assess the impact of these mutations on early cardiac development, we integrated single-cell and spatial transcriptomics of early human heart development with our genetic findings. We discovered that the candidate gene expression was enriched in the myogenic progenitors of the cardiac outflow tract. Moreover, subsets of the candidate genes were found in specific gene coexpression modules along the cardiomyocyte differentiation trajectory. These integrative functional analyses help dissect the pathogenesis of TOF, revealing cellular hotspots in early heart development resulting in cardiac malformations.

Size matters: Large copy number losses in Hirschsprung disease patients reveal genes involved in enteric nervous system development.

Hirschsprung disease (HSCR) is a complex genetic disease characterized by absence of ganglia in the intestine. HSCR etiology can be explained by a unique combination of genetic alterations: rare coding variants, predisposing haplotypes and Copy Number Variation (CNV). Approximately 18% of patients have additional anatomical malformations or neurological symptoms (HSCR-AAM). Pinpointing the responsible culprits within a CNV is challenging as often many genes are affected. Therefore, we selected candidate genes based on gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics. Next, we used a zebrafish model to investigate whether loss of these genes affects enteric neuron development in vivo. This study included three groups of patients, two groups without coding variants in disease associated genes: HSCR-AAM and HSCR patients without associated anomalies (HSCR-isolated). The third group consisted of all HSCR patients in which a confirmed pathogenic rare coding variant was identified. We compared these patient groups to unaffected controls. Predisposing haplotypes were determined, confirming that every HSCR subgroup had increased contributions of predisposing haplotypes, but their contribution was highest in isolated HSCR patients without RET coding variants. CNV profiling proved that specifically HSCR-AAM patients had larger Copy Number (CN) losses. Gene enrichment strategies using mouse enteric nervous system transcriptomes and constraint metrics were used to determine plausible candidate genes located within CN losses. Validation in zebrafish using CRISPR/Cas9 targeting confirmed the contribution of UFD1L, TBX2, SLC8A1, and MAPK8 to ENS development. In addition, we revealed epistasis between reduced Ret and Gnl1 expression and between reduced Ret and Tubb5 expression in vivo. Rare large CN losses-often de novo-contribute to HSCR in HSCR-AAM patients. We proved the involvement of six genes in enteric nervous system development and Hirschsprung disease.

Identification of a wide spectrum of ciliary gene mutations in nonsyndromic biliary atresia patients implicates ciliary dysfunction as a novel disease mechanism.

BACKGROUND: Biliary atresia (BA) is the most common obstructive cholangiopathy in neonates, often progressing to end-stage cirrhosis. BA pathogenesis is believed to be multifactorial, but the genetic contribution, especially for nonsyndromic BA (common form: > 85%) remains poorly defined. METHODS: We conducted whole exome sequencing on 89 nonsyndromic BA trios to identify rare variants contributing to BA etiology. Functional evaluation using patients' liver biopsies, human cell and zebrafish models were performed. Clinical impact on respiratory system was assessed with clinical evaluation, nasal nitric oxide (nNO), high speed video analysis and transmission electron microscopy. FINDINGS: We detected rare, deleterious de novo or biallelic variants in liver-expressed ciliary genes in 31.5% (28/89) of the BA patients. Burden test revealed 2.6-fold (odds ratio (OR) [95% confidence intervals (CI)]= 2.58 [1.15-6.07], adjusted p = 0.034) over-representation of rare, deleterious mutations in liver-expressed ciliary gene set in patients compared to controls. Functional analyses further demonstrated absence of cilia in the BA livers with KIF3B and TTC17 mutations, and knockdown of PCNT, KIF3B and TTC17 in human control fibroblasts and cholangiocytes resulted in reduced number of cilia. Additionally, CRISPR/Cas9-engineered zebrafish knockouts of KIF3B, PCNT and TTC17 displayed reduced biliary flow. Abnormally low level of nNO was detected in 80% (8/10) of BA patients carrying deleterious ciliary mutations, implicating the intrinsic ciliary defects. INTERPRETATION: Our findings support strong genetic susceptibility for nonsyndromic BA. Ciliary gene mutations leading to cholangiocyte cilia malformation and dysfunction could be a key biological mechanism in BA pathogenesis. FUNDING: The study is supported by General Research Fund, HMRF Commissioned Paediatric Research at HKCH and Li Ka Shing Faculty of Medicine Enhanced New Staff Start-up Fund.

Whole-genome analysis of noncoding genetic variations identifies multiscale regulatory element perturbations associated with Hirschsprung disease.

It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the RET, RASGEF1A, and PIK3C2B loci. In the PIK3C2B regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.

A random forest-based framework for genotyping and accuracy assessment of copy number variations.

Detection of copy number variations (CNVs) is essential for uncovering genetic factors underlying human diseases. However, CNV detection by current methods is prone to error, and precisely identifying CNVs from paired-end whole genome sequencing (WGS) data is still challenging. Here, we present a framework, CNV-JACG, for Judging the Accuracy of CNVs and Genotyping using paired-end WGS data. CNV-JACG is based on a random forest model trained on 21 distinctive features characterizing the CNV region and its breakpoints. Using the data from the 1000 Genomes Project, Genome in a Bottle Consortium, the Human Genome Structural Variation Consortium and in-house technical replicates, we show that CNV-JACG has superior sensitivity over the latest genotyping method, SV2, particularly for the small CNVs (≤1 kb). We also demonstrate that CNV-JACG outperforms SV2 in terms of Mendelian inconsistency in trios and concordance between technical replicates. Our study suggests that CNV-JACG would be a useful tool in assessing the accuracy of CNVs to meet the ever-growing needs for uncovering the missing heritability linked to CNVs.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations.

INTRODUCTION: Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model. METHODS: We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing. RESULTS: When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing. CONCLUSION: Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Cancer gene mutations in congenital pulmonary airway malformation patients.

BACKGROUND: Newborns affected with congenital pulmonary airway malformations (CPAMs) may present with severe respiratory distress or remain asymptomatic. While surgical resection is the definitive treatment for symptomatic CPAMs, prophylactic elective surgery may be recommended for asymptomatic CPAMs owing to the risk of tumour development. However, the implementation of prophylactic surgery is quite controversial on the grounds that more evidence linking CPAMs and cancer is needed. The large gap in knowledge of CPAM pathogenesis results in uncertainties and controversies in disease management. As developmental genes control postnatal cell growth and contribute to cancer development, we hypothesised that CPAMs may be underlain by germline mutations in genes governing airways development. METHODS: Sequencing of the exome of 19 patients and their unaffected parents. RESULTS: A more than expected number of mutations in cancer genes (false discovery rate q-value <5.01×10-5) was observed. The co-occurrence, in the same patient, of damaging variants in genes encoding interacting proteins is intriguing, the most striking being thyroglobulin (TG) and its receptor, megalin (LRP2). Both genes are highly relevant in lung development and cancer. CONCLUSIONS: The overall excess of mutations in cancer genes may account for the reported association of CPAMs with carcinomas and provide some evidence to argue for prophylactic surgery by some surgeons.

Conditional deletion of platelet derived growth factor receptor alpha (Pdgfra) in urorectal mesenchyme causes mesenchyme apoptosis and urorectal developmental anomalies in mice.

In mammals, urorectal development starts at early embryonic stage, defective urorectal development results in anorectal malformations, which are common congenital developmental defects of the anus and the urethra in newborns. The etiology and embryology of the defects are still largely unknown. Platelet-derived growth factor receptor alpha (Pdgfra) is a cell surface receptor tyrosine kinase, upon binding to its ligands (Pdgfa-d), mediates intracellular signaling and regulates embryonic development. The expression of Pdgfra is tightly regulated in the developing urorectal mesenchyme, and its dysregulation is associated with urorectal defects in animals with urorectal defects. Knockout of Pdgfra induces early embryo lethality which precludes investigation of Pdgfra in urorectal development. To address the temporal requirement of Pdgfra in urorectal development, we conditionally deleted Pdgfra in Pdgfra-expressing tissues using a tamoxifen inducible Cre-loxP approach in mice, examined the urorectal development in Pdgfra conditional knockout (Pdgfra-cKO) embryos. We showed that conditional deletion of Pdgfra in Pdgfra-expressing tissues at E10-E11 caused cloaca septation defect, anteriorly displaced anus, defective urogenital folds development and abnormal urethra tubularization in both male and female mice. Furthermore, we showed that Pdgfra was required for the survival of urorectal mesenchyme, deletion of Pdgfra caused apoptosis in the peri-cloacal, the peri-urethra and the urorectal septum mesenchyme of Pdgfra-cKO mutants, associated with an induction of p53, Ndrg1 and activation of caspase-3 in Pdgfra-cKO embryos. In conclusion, Pdgfra is required for the development and survival of the urorectal mesenchyme in embryo, dysregulated Pdgfra signaling induced urorectal defects in mice resembling human congenital diseases of anorectal malformations and hypospadias. Perturbation of PDGFRA signaling may contribute to anorectal malformations and hypospadias in human.

Identification of Genes Associated With Hirschsprung Disease, Based on Whole-Genome Sequence Analysis, and Potential Effects on Enteric Nervous System Development.

BACKGROUND & AIMS: Hirschsprung disease, or congenital aganglionosis, is believed to be oligogenic-that is, caused by multiple genetic factors. We performed whole-genome sequence analyses of patients with Hirschsprung disease to identify genetic factors that contribute to disease development and analyzed the functional effects of these variants. METHODS: We performed whole-genome sequence analyses of 443 patients with short-segment disease, recruited from hospitals in China and Vietnam, and 493 ethnically matched individuals without Hirschsprung disease (controls). We performed genome-wide association analyses and gene-based rare-variant burden tests to identify rare and common disease-associated variants and study their interactions. We obtained induced pluripotent stem cell (iPSC) lines from 4 patients with Hirschsprung disease and 2 control individuals, and we used these to generate enteric neural crest cells for transcriptomic analyses. We assessed the neuronal lineage differentiation capability of iPSC-derived enteric neural crest cells using an in vitro differentiation assay. RESULTS: We identified 4 susceptibility loci, including 1 in the phospholipase D1 gene (PLD1) (P = 7.4 × 10-7). The patients had a significant excess of rare protein-altering variants in genes previously associated with Hirschsprung disease and in the ß-secretase 2 gene (BACE2) (P = 2.9 × 10-6). The epistatic effects of common and rare variants across these loci provided a sensitized background that increased risk for the disease. In studies of the iPSCs, we observed common and distinct pathways associated with variants in RET that affect risk. In functional assays, we found variants in BACE2 to protect enteric neurons from apoptosis. We propose that alterations in BACE1 signaling via amyloid ß precursor protein and BACE2 contribute to pathogenesis of Hirschsprung disease. CONCLUSIONS: In whole-genome sequence analyses of patients with Hirschsprung disease, we identified rare and common variants associated with disease risk. Using iPSC cells, we discovered some functional effects of these variants.

Uncovering the genetic lesions underlying the most severe form of Hirschsprung disease by whole-genome sequencing.

Hirschsprung disease (HSCR) is a complex birth defect characterized by the lack of ganglion cells along a variable length of the distal intestine. A large proportion of HSCR patients remain genetically unexplained. We applied whole-genome sequencing (WGS) on 9 trios where the probands are sporadically affected with the most severe form of the disorder and harbor no coding sequence variants affecting the function of known HSCR genes. We found de novo protein-altering variants in three intolerant to change genes-CCT2, VASH1, and CYP26A1-for which a plausible link with the enteric nervous system (ENS) exists. De novo single-nucleotide and indel variants were present in introns and non-coding neighboring regions of ENS-related genes, including NRG1 and ERBB4. Joint analysis with those inherited rare variants found under recessive and/or digenic models revealed both patient-unique and shared genetic features where rare variants were found to be enriched in the extracellular matrix-receptor (ECM-receptor) pathway (p = 3.4 × 10-11). Delineation of the genetic profile of each patient might help finding common grounds that could lead to the discovery of shared molecules that could be used as drug targets for the currently ongoing cell therapy effort which aims at providing an alternative to the surgical treatment.

Actionable secondary findings from whole-genome sequencing of 954 East Asians.

Recently, the American College of Medical Genetics (ACMG) recommended the return of actionable secondary findings detected from clinical sequencing. The reported frequency of secondary findings in Asian populations were highly variable and it is unclear whether the uniformity in coverage offered by whole-genome sequencing (WGS) may impact the estimate. In this analysis, we aimed to refine the rate of secondary findings on East Asians through a large-scale WGS study. We classified 1256 protein-altering or splicing variants of the 59 actionable genes detected from WGS of 954 East Asians in strict accordance with the ACMG and the Association for Molecular Pathology guidelines. A total of 21 pathogenic or likely pathogenic variants were detected in 24 of the 954 East Asian genomes with an estimate of 2.5% of East Asians carrying actionable variants. Although the overall estimate of secondary findings was consistent with those reported for non-East Asian ethnicities, genetic and allelic heterogeneity was observed. WGS offers a wider breadth of coverage over WES, which highlights the need to further investigate the variable sensitivity of WES and WGS in the detection of secondary findings. Identifying secondary findings in populations underrepresented in previous genetic literature might improve variant interpretation and has a profound impact on local decision-making with regard to the cost-effectiveness of returning the secondary findings from clinical sequencing.

De novo mutations in Caudal Type Homeo Box transcription Factor 2 (CDX2) in patients with persistent cloaca.

The cloaca is an embryonic cavity that is divided into the urogenital sinus and rectum upon differentiation of the cloacal epithelium triggered by tissue-specific transcription factors including CDX2. Defective differentiation leads to persistent cloaca in humans (PC), a phenotype recapitulated in Cdx2 mutant mice. PC is linked to hypo/hyper-vitaminosis A. Although no gene has ever been identified, there is a strong evidence for a genetic contribution to PC. We applied whole-exome sequencing and copy-number-variants analyses to 21 PC patients and their unaffected parents. The damaging p.Cys132* and p.Arg237His de novo CDX2 variants were identified in two patients. These variants altered the expression of CYP26A1, a direct CDX2 target encoding the major retinoic acid (RA)-degrading enzyme. Other RA genes, including the RA-receptor alpha, were also mutated. Genes governing the development of cloaca-derived structures were recurrently mutated and over-represented in the basement-membrane components set (q-value < 1.65 × 10-6). Joint analysis of the patients' profile highlighted the extracellular matrix-receptor interaction pathway (MsigDBID: M7098, FDR: q-value < 7.16 × 10-9). This is the first evidence that PC is genetic, with genes involved in the RA metabolism at the lead. Given the CDX2 de novo variants and the role of RA, our observations could potentiate preventive measures. For the first time, a gene recapitulating PC in mouse models is found mutated in humans.

Patient complexity and genotype-phenotype correlations in biliary atresia: a cross-sectional analysis.

BACKGROUND: Biliary Atresia (BA) is rare and genetically complex, and the pathogenesis is elusive. The disease course is variable and can represent heterogeneity, which hinders effective disease management. Deciphering the BA phenotypic variance is a priority in clinics and can be achieved by the integrative analysis of genotype and phenotype. We aim to explore the BA phenotypic features and to delineate the source of its variance. METHODS: The study is a cross-sectional observational study collating with case/control association analysis. One-hundred-and-eighty-one type III non-syndromic BA patients and 431 controls were included for case-control association tests, including 89 patients (47.19% males, born June 15th, 1981 to September 17th, 2007) have detailed clinical records with follow-up of the disease course (median ~17.2 years). BA-association genes from the genome-wide gene-based association test on common genetic variants (CV) and rare copy-number-variants (CNVs) from the genome-wide survey, the later comprise only CNVs > 100 kb and found in the BA patients but not in the local population (N = 1,381) or the database (N = 11,943). Hereby comorbidity is defined as a chronic disease that affects the BA patients but has no known relationship with BA or with the BA treatment. We examined genotype-phenotype correlations of CNVs, connectivity of these novel variants with BA-associated CVs, and their role in the BA candidate gene network. RESULTS: Of the 89 patients, 41.57% have comorbidities, including autoimmune-allergic disorders (22.47%). They carried 29 BA-private CNVs, including 3 CNVs underpinning the carriers' immunity comorbidity and one JAG1 micro-deletion. The BA-CNV-intersected genes (N = 102) and the CV-tagged genes (N = 103) were both enriched with immune-inflammatory pathway genes (FDR q < 0.20), and the two gene sets were interconnected (permutation p = 0.039). The molecular network representing CVs and rare-CNV association genes fit into a core/periphery structure, the immune genes and their related modules are found at the coherence core of all connections, suggesting its dominant role in the BA pathogenesis pathway. CONCLUSIONS: The study highlights a patient-complexity phenomenon as a novel BA phenotypic feature, which is underpinned by rare-CNVs that biologically converge with CVs into the immune-inflammatory pathway and drives the BA occurrence and the likely BA association with immune diseases in clinics.

Whole exome sequencing coupled with unbiased functional analysis reveals new Hirschsprung disease genes.

BACKGROUND: Hirschsprung disease (HSCR), which is congenital obstruction of the bowel, results from a failure of enteric nervous system (ENS) progenitors to migrate, proliferate, differentiate, or survive within the distal intestine. Previous studies that have searched for genes underlying HSCR have focused on ENS-related pathways and genes not fitting the current knowledge have thus often been ignored. We identify and validate novel HSCR genes using whole exome sequencing (WES), burden tests, in silico prediction, unbiased in vivo analyses of the mutated genes in zebrafish, and expression analyses in zebrafish, mouse, and human. RESULTS: We performed de novo mutation (DNM) screening on 24 HSCR trios. We identify 28 DNMs in 21 different genes. Eight of the DNMs we identified occur in RET, the main HSCR gene, and the remaining 20 DNMs reside in genes not reported in the ENS. Knockdown of all 12 genes with missense or loss-of-function DNMs showed that the orthologs of four genes (DENND3, NCLN, NUP98, and TBATA) are indispensable for ENS development in zebrafish, and these results were confirmed by CRISPR knockout. These genes are also expressed in human and mouse gut and/or ENS progenitors. Importantly, the encoded proteins are linked to neuronal processes shared by the central nervous system and the ENS. CONCLUSIONS: Our data open new fields of investigation into HSCR pathology and provide novel insights into the development of the ENS. Moreover, the study demonstrates that functional analyses of genes carrying DNMs are warranted to delineate the full genetic architecture of rare complex diseases.

Genetic study of congenital bile-duct dilatation identifies de novo and inherited variants in functionally related genes.

BACKGROUND: Congenital dilatation of the bile-duct (CDD) is a rare, mostly sporadic, disorder that results in bile retention with severe associated complications. CDD affects mainly Asians. To our knowledge, no genetic study has ever been conducted. METHODS: We aim to identify genetic risk factors by a "trio-based" exome-sequencing approach, whereby 31 CDD probands and their unaffected parents were exome-sequenced. Seven-hundred controls from the local population were used to detect gene-sets significantly enriched with rare variants in CDD patients. RESULTS: Twenty-one predicted damaging de novo variants (DNVs; 4 protein truncating and 17 missense) were identified in several evolutionarily constrained genes (p < 0.01). Six genes carrying DNVs were associated with human developmental disorders involving epithelial, connective or bone morphologies (PXDN, RTEL1, ANKRD11, MAP2K1, CYLD, ACAN) and four linked with cholangio- and hepatocellular carcinomas (PIK3CA, TLN1 CYLD, MAP2K1). Importantly, CDD patients have an excess of DNVs in cancer-related genes (p < 0.025). Thirteen genes were recurrently mutated at different sites, forming compound heterozygotes or functionally related complexes within patients. CONCLUSIONS: Our data supports a strong genetic basis for CDD and show that CDD is not only genetically heterogeneous but also non-monogenic, requiring mutations in more than one genes for the disease to develop. The data is consistent with the rarity and sporadic presentation of CDD.

Sacral agenesis: a pilot whole exome sequencing and copy number study.

BACKGROUND: Caudal regression syndrome (CRS) or sacral agenesis is a rare congenital disorder characterized by a constellation of congenital caudal anomalies affecting the caudal spine and spinal cord, the hindgut, the urogenital system, and the lower limbs. CRS is a complex condition, attributed to an abnormal development of the caudal mesoderm, likely caused by the effect of interacting genetic and environmental factors. A well-known risk factor is maternal type 1 diabetes. METHOD: Whole exome sequencing and copy number variation (CNV) analyses were conducted on 4 Caucasian trios to identify de novo and inherited rare mutations. RESULTS: In this pilot study, exome sequencing and copy number variation (CNV) analyses implicate a number of candidate genes, including SPTBN5, MORN1, ZNF330, CLTCL1 and PDZD2. De novo mutations were found in SPTBN5, MORN1 and ZNF330 and inherited predicted damaging mutations in PDZD2 (homozygous) and CLTCL1 (compound heterozygous). Importantly, predicted damaging mutations in PTEN (heterozygous), in its direct regulator GLTSCR2 (compound heterozygous) and in VANGL1 (heterozygous) were identified. These genes had previously been linked with the CRS phenotype. Two CNV deletions, one de novo (chr3q13.13) and one homozygous (chr8p23.2), were detected in one of our CRS patients. These deletions overlapped with CNVs previously reported in patients with similar phenotype. CONCLUSION: Despite the genetic diversity and the complexity of the phenotype, this pilot study identified genetic features common across CRS patients.

Trans-ethnic meta-analysis of genome-wide association studies for Hirschsprung disease.

Hirschsprung disease (HSCR) is the most common cause of neonatal intestinal obstruction. It is characterized by the absence of ganglia in the nerve plexuses of the lower gastrointestinal tract. So far, three common disease-susceptibility variants at the RET, SEMA3 and NRG1 loci have been detected through genome-wide association studies (GWAS) in Europeans and Asians to understand its genetic etiologies. Here we present a trans-ethnic meta-analysis of 507 HSCR cases and 1191 controls, combining all published GWAS results on HSCR to fine-map these loci and narrow down the putatively causal variants to 99% credible sets. We also demonstrate that the effects of RET and NRG1 are universal across European and Asian ancestries. In contrast, we detected a European-specific association of a low-frequency variant, rs80227144, in SEMA3 [odds ratio (OR) = 5.2, P = 4.7 × 10-10]. Conditional analyses on the lead SNPs revealed a secondary association signal, corresponding to an Asian-specific, low-frequency missense variant encoding RET p.Asp489Asn (rs9282834, conditional OR = 20.3, conditional P = 4.1 × 10-14). When in trans with the RET intron 1 enhancer risk allele, rs9282834 increases the risk of HSCR from 1.1 to 26.7. Overall, our study provides further insights into the genetic architecture of HSCR and has profound implications for future study designs.

Genetics of enteric neuropathies.

Abnormal development or disturbed functioning of the enteric nervous system (ENS), the intrinsic innervation of the gastrointestinal tract, is associated with the development of neuropathic gastrointestinal motility disorders. Here, we review the underlying molecular basis of these disorders and hypothesize that many of them have a common defective biological mechanism. Genetic burden and environmental components affecting this common mechanism are ultimately responsible for disease severity and symptom heterogeneity. We believe that they act together as the fulcrum in a seesaw balanced with harmful and protective factors, and are responsible for a continuum of symptoms ranging from neuronal hyperplasia to absence of neurons.