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.

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.

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.

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.

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.

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.

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.

Meta-analysis of GWAS on two Chinese populations followed by replication identifies novel genetic variants on the X chromosome associated with systemic lupus erythematosus.

Systemic lupus erythematosus (SLE) is a prototypic autoimmune disease that affects mainly females. What role the X chromosome plays in the disease has always been an intriguing question. In this study, we examined the genetic variants on the X chromosome through meta-analysis of two genome-wide association studies (GWAS) on SLE on Chinese Han populations. Prominent association signals from the meta-analysis were replicated in 4 additional Asian cohorts, with a total of 5373 cases and 9166 matched controls. We identified a novel variant in PRPS2 on Xp22.3 as associated with SLE with genome-wide significance (rs7062536, OR = 0.84, P = 1.00E-08). Association of the L1CAM-MECP2 region with SLE was reported previously. In this study, we identified independent contributors in this region in NAA10 (rs2071128, OR = 0.81, P = 2.19E-13) and TMEM187 (rs17422, OR = 0.75, P = 1.47E-15), in addition to replicating the association from IRAK1-MECP2 region (rs1059702, OR = 0.71, P = 2.40E-18) in Asian cohorts. The X-linked susceptibility variants showed higher effect size in males than that in females, similar to results from a genome-wide survey of associated SNPs on the autosomes. These results suggest that susceptibility genes identified on the X chromosome, while contributing to disease predisposition, might not contribute significantly to the female predominance of this prototype autoimmune disease.

Inheritance-mode specific pathogenicity prioritization (ISPP) for human protein coding genes.

MOTIVATION: Exome sequencing studies have facilitated the detection of causal genetic variants in yet-unsolved Mendelian diseases. However, the identification of disease causal genes among a list of candidates in an exome sequencing study is still not fully settled, and it is often difficult to prioritize candidate genes for follow-up studies. The inheritance mode provides crucial information for understanding Mendelian diseases, but none of the existing gene prioritization tools fully utilize this information. RESULTS: We examined the characteristics of Mendelian disease genes under different inheritance modes. The results suggest that Mendelian disease genes with autosomal dominant (AD) inheritance mode are more haploinsufficiency and de novo mutation sensitive, whereas those autosomal recessive (AR) genes have significantly more non-synonymous variants and regulatory transcript isoforms. In addition, the X-linked (XL) Mendelian disease genes have fewer non-synonymous and synonymous variants. As a result, we derived a new scoring system for prioritizing candidate genes for Mendelian diseases according to the inheritance mode. Our scoring system assigned to each annotated protein-coding gene (N = 18 859) three pathogenic scores according to the inheritance mode (AD, AR and XL). This inheritance mode-specific framework achieved higher accuracy (area under curve = 0.84) in XL mode. CONCLUSION: The inheritance-mode specific pathogenicity prioritization (ISPP) outperformed other well-known methods including Haploinsufficiency, Recessive, Network centrality, Genic Intolerance, Gene Damage Index and Gene Constraint scores. This systematic study suggests that genes manifesting disease inheritance modes tend to have unique characteristics. AVAILABILITY AND IMPLEMENTATION: ISPP is included in KGGSeq v1.0 (http://grass.cgs.hku.hk/limx/kggseq/), and source code is available from (https://github.com/jacobhsu35/ISPP.git). CONTACT: mxli@hku.hkSupplementary information: Supplementary data are available at Bioinformatics online.

Three SNPs in chromosome 11q23.3 are independently associated with systemic lupus erythematosus in Asians.

Systemic lupus erythematosus (SLE) has a complex etiology and is affected by both genetic and environmental factors. Although more than 40 loci have shown robust association with SLE, the details of these loci, such as the independent contributors and the genes involved, are still unclear. In this study, we performed meta-analysis of two existing genome-wide association studies (GWASs) on Chinese Han populations from Hong Kong and Anhui, China, and followed the findings by further replication on three additional Chinese and Thailand cohorts with a total of 4254 cases and 6262 controls matched geographically and ethnically. We discovered multiple susceptibility variants for SLE in the 11q23.3 region, including variants in/near PHLDB1 (rs11603023, P(_combined) = 1.25E-08, OR = 1.20), DDX6 (rs638893, P(_combined) = 5.19E-07, OR = 1.22) and CXCR5 (rs10892301, P(_combined) = 2.51E-08, OR = 0.85). Genetic contributions from the newly identified variants were all independent of SNP rs4639966, whose association was reported from the previous GWAS. In addition, the three newly identified variants all showed independent association with the disease through modeling by both stepwise and conditional logistic regression. The presence of multiple independent variants in this region emphasizes its role in SLE susceptibility, and also hints the possibility that distinct biological mechanisms might be involved in the disease involving this genomic region.

Genome-wide search followed by replication reveals genetic interaction of CD80 and ALOX5AP associated with systemic lupus erythematosus in Asian populations.

OBJECTIVES: Genetic interaction has been considered as a hallmark of the genetic architecture of systemic lupus erythematosus (SLE). Based on two independent genome-wide association studies (GWAS) on Chinese populations, we performed a genome-wide search for genetic interactions contributing to SLE susceptibility. METHODS: The study involved a total of 1 659 cases and 3 398 controls in the discovery stage and 2 612 cases and 3 441 controls in three cohorts for replication. Logistic regression and multifactor dimensionality reduction were used to search for genetic interaction. RESULTS: Interaction of CD80 (rs2222631) and ALOX5AP (rs12876893) was found to be significantly associated with SLE (OR_int=1.16, P_int_all=7.7E-04 at false discovery rate<0.05). Single nuclear polymorphism rs2222631 was found associated with SLE with genome-wide significance (P_all=4.5E-08, OR=0.86) and is independent of rs6804441 in CD80, whose association was reported previously. Significant correlation was observed between expression of these two genes in healthy controls and SLE cases, together with differential expression of these genes between cases and controls, observed from individuals from the Hong Kong cohort. Genetic interactions between BLK (rs13277113) and DDX6 (rs4639966), and between TNFSF4 (rs844648) and PXK (rs6445975) were also observed in both GWAS data sets. CONCLUSIONS: Our study represents the first genome-wide evaluation of epistasis interactions on SLE and the findings suggest interactions and independent variants may help partially explain missing heritability for complex diseases.

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.

Genome-wide copy number variation study in anorectal malformations.

Anorectal malformations (ARMs, congenital obstruction of the anal opening) are among the most common birth defects requiring surgical treatment (2-5/10 000 live-births) and carry significant chronic morbidity. ARMs present either as isolated or as part of the phenotypic spectrum of some chromosomal abnormalities or monogenic syndromes. The etiology is unknown. To assess the genetic contribution to ARMs, we investigated single-nucleotide polymorphisms and copy number variations (CNVs) at genome-wide scale. A total of 363 Han Chinese sporadic ARM patients and 4006 Han Chinese controls were included. Overall, we detected a 1.3-fold significant excess of rare CNVs in patients. Stratification of patients by presence/absence of other congenital anomalies showed that while syndromic ARM patients carried significantly longer rare duplications than controls (P = 0.049), non-syndromic patients were enriched with both rare deletions and duplications when compared with controls (P = 0.00031). Twelve chromosomal aberrations and 114 rare CNVs were observed in patients but not in 868 controls nor 11 943 healthy individuals from the Database of Genomic Variants. Importantly, these aberrations were observed in isolated ARM patients. Gene-based analysis revealed 79 genes interfered by CNVs in patients only. In particular, we identified a de novo DKK4 duplication. DKK4 is a member of the WNT signaling pathway which is involved in the development of the anorectal region. In mice, Wnt disruption results in ARMs. Our data suggest a role for rare CNVs not only in syndromic but also in isolated ARM patients and provide a list of plausible candidate genes for the disorder.

Chinese family with diffuse oesophageal leiomyomatosis: a new COL4A5/COL4A6 deletion and a case of gonosomal mosaicism.

BACKGROUND: Diffuse oesophageal leiomyomatosis (DOL) is a rare disorder characterized by tumorous overgrowth of the muscular wall of the oesophagus. DOL is present in 5 % of Alport syndrome (AS) patients. AS is a rare hereditary disease that involves varying degrees of hearing impairment, ocular changes and progressive glomerulonephritis leading to renal failure. In DOL-AS patients, the genetic defect consists of a deletion involving the COL4A5 and COL4A6 genes on the X chromosome. CASE PRESENTATION: We report a two-generation family (4 individuals; parents and two children, one male and one female) with two members (mother and son) affected with oesophageal leiomyomatosis. Signs of potential renal failure, which characterizes AS, were only apparent in the index patient (son) 2 years and three months after the initial diagnosis of DOL. Blood DNA from the four family members were submitted to exome sequencing and array genotyping to perform a genome wide screening for disease causal single nucleotide (SN) and copy number (CN) variations. Analyses revealed a new 40kb deletion encompassing from intron 2 of COL4A5 to intron 1 of COL4A6 at Xq22.3. The breakpoints were also identified. Possible confounding pathogenic exonic variants in genes known to be involved in other extracellular matrices disorders were also shared by the two affected individuals. Meticulous analysis of the maternal DNA revealed a case of gonosomal mosaicism. CONCLUSIONS: This is the first report of gonadosomal mosaicism associated to DOL-AS.

Genome-wide copy number analysis uncovers a new HSCR gene: NRG3.

Hirschsprung disease (HSCR) is a congenital disorder characterized by aganglionosis of the distal intestine. To assess the contribution of copy number variants (CNVs) to HSCR, we analysed the data generated from our previous genome-wide association study on HSCR patients, whereby we identified NRG1 as a new HSCR susceptibility locus. Analysis of 129 Chinese patients and 331 ethnically matched controls showed that HSCR patients have a greater burden of rare CNVs (p = 1.50 × 10(-5)), particularly for those encompassing genes (p = 5.00 × 10(-6)). Our study identified 246 rare-genic CNVs exclusive to patients. Among those, we detected a NRG3 deletion (p = 1.64 × 10(-3)). Subsequent follow-up (96 additional patients and 220 controls) on NRG3 revealed 9 deletions (combined p = 3.36 × 10(-5)) and 2 de novo duplications among patients and two deletions among controls. Importantly, NRG3 is a paralog of NRG1. Stratification of patients by presence/absence of HSCR-associated syndromes showed that while syndromic-HSCR patients carried significantly longer CNVs than the non-syndromic or controls (p = 1.50 × 10(-5)), non-syndromic patients were enriched in CNV number when compared to controls (p = 4.00 × 10(-6)) or the syndromic counterpart. Our results suggest a role for NRG3 in HSCR etiology and provide insights into the relative contribution of structural variants in both syndromic and non-syndromic HSCR. This would be the first genome-wide catalog of copy number variants identified in HSCR.

Contribution of rare and common variants determine complex diseases-Hirschsprung disease as a model.

Finding genes for complex diseases has been the goal of many genetic studies. Most of these studies have been successful by searching for genes and mutations in rare familial cases, by screening candidate genes and by performing genome wide association studies. However, only a small fraction of the total genetic risk for these complex genetic diseases can be explained by the identified mutations and associated genetic loci. In this review we focus on Hirschsprung disease (HSCR) as an example of a complex genetic disorder. We describe the genes identified in this congenital malformation and postulate that both common 'low penetrant' variants in combination with rare or private 'high penetrant' variants determine the risk on HSCR, and likely, on other complex diseases. We also discuss how new technological advances can be used to gain further insights in the genetic background of complex diseases. Finally, we outline a few steps to develop functional assays in order to determine the involvement of these variants in disease development.

Targeted next-generation sequencing on Hirschsprung disease: a pilot study exploits DNA pooling.

To adopt an efficient approach of identifying rare variants possibly related to Hirschsprung disease (HSCR), a pilot study was set up to evaluate the performance of a newly designed protocol for next generation targeted resquencing. In total, 20 Chinese HSCR patients and 20 Chinese sex-matched individuals with no HSCR were included, for which coding sequences (CDS) of 62 genes known to be in signaling pathways relevant to enteric nervous system development were selected for capture and sequencing. Blood DNAs from eight pools of five cases or controls were enriched by PCR-based RainDance technology (RDT) and then sequenced on a 454 FLX platform. As technical validation, five patients from case Pool-3 were also independently enriched by RDT, indexed with barcode and sequenced with sufficient coverage. Assessment for CDS single nucleotide variants showed DNA pooling performed well (specificity/sensitivity at 98.4%/83.7%) at the common variant level; but relatively worse (specificity/sensitivity at 65.5%/61.3%) at the rare variant level. Further Sanger sequencing only validated five out of 12 rare damaging variants likely involved in HSCR. Hence more improvement at variant detection and sequencing technology is needed to realize the potential of DNA pooling for large-scale resequencing projects.

ELF1 is associated with systemic lupus erythematosus in Asian populations.

Systemic lupus erythematosus (SLE) is an autoimmune disease with a strong genetic involvement. The susceptibility genes identified so far can only explain a small proportion of disease heritability. Through a genome-wide association in a Hong Kong Chinese cohort and subsequent replication in two other Asian populations, with a total of 3164 patients and 4482 matched controls, we identified association of ELF1 (E74-like factor 1) with SLE (rs7329174, OR = 1.26, joint P= 1.47 × 10(-8)). ELF1 belongs to the ETS family of transcription factors and is known to be involved in T cell development and function. Database analysis revealed transcripts making use of three alternative exon1s for this gene. Near equivalent expression levels of distinct transcripts initiated from alternative exon1s were detected in peripheral blood mononuclear cells from both SLE patients and healthy controls. Although a direct association of rs7329174 with the three forms of transcripts for this gene was not detected, these findings support an important role of ELF1 in SLE susceptibility and suggest a potentially tight regulation for the expression of this gene.

Common genetic variants regulating ADD3 gene expression alter biliary atresia risk.

BACKGROUND & AIMS: Biliary atresia (BA) is a rare and most severe cholestatic disease in neonates, but the pathogenic mechanisms are unknown. Through a previous genome wide association study (GWAS) on Han Chinese, we discovered association of the 10q24.2 region encompassing ADD3 and XPNPEP1 genes, which was replicated in Chinese and Thai populations. This study aims to fully characterize the genetic architecture at 10q24.2 and to reveal the link between the genetic variants and BA. METHODS: We genotyped 107 single nucleotide polymorphisms (SNPs) in 10q24.2 in 339 Han Chinese patients and 401 matched controls using Sequenom. Exhaustive follow-up studies of the association signals were performed. RESULTS: The combined BA-association p-value of the GWAS SNP (rs17095355) achieved 6.06×10(-10). Further, we revealed the common risk haplotype encompassing 5 tagging-SNPs, capturing the risk-predisposing alleles in 10q24.2 [p=5.32×10(-11); odds ratio, OR: 2.38; confidence interval, CI: (2.14-2.62)]. Through Sanger sequencing, no deleterious rare variants (RVs) residing in the risk haplotype were found, dismissing the theory of "synthetic" association. Moreover, in bioinformatics and in vivo genotype-expression investigations, the BA-associated potentially regulatory SNPs correlated with ADD3 gene expression (n=36; p=0.0030). Remarkably, the risk haplotype frequency coincides with BA incidences in the population, and, positive selection (favoring the derived alleles that arose from mutations) was evident at the ADD3 locus, suggesting a possible role for the BA-associated common variants in shaping the general population diversity. CONCLUSIONS: Common genetic variants in 10q24.2 can alter BA risk by regulating ADD3 expression levels in the liver, and may exert an effect on disease epidemiology and on the general population.

RET and NRG1 interplay in Hirschsprung disease.

Hirschsprung disease (HSCR, aganglionic megacolon) is a complex genetic disorder of the enteric nervous system (ENS) characterized by the absence of enteric neurons along a variable length of the intestine. While rare variants (RVs) in the coding sequence (CDS) of several genes involved in ENS development lead to disease, the association of common variants (CVs) with HSCR has only been reported for RET (the major HSCR gene) and NRG1. Importantly, RVs in the CDS of these two genes are also associated with the disorder. To assess independent and joint effects between the different types of RET and NRG1 variants identified in HSCR patients, we used 254 Chinese sporadic HSCR patients and 143 ethnically matched controls for whom the RET and/or NRG1 variants genotypes (rare and common) were available. Four genetic risk factors were defined and interaction effects were modeled using conditional logistic regression analyses and pair-wise Kendall correlations. Our analysis revealed a joint effect of RET CVs with RET RVs, NRG1 CVs or NRG1 RVs. To assess whether the genetic interaction translated into functional interaction, mouse neural crest cells (NCCs; enteric neuron precursors) isolated from embryonic guts were treated with NRG1 (ErbB2 ligand) or/and GDNF (Ret ligand) and monitored during the subsequent neural differentiation process. Nrg1 inhibited the Gdnf-induced neuronal differentiation and Gdnf negatively regulated Nrg1-signaling by down-regulating the expression of its receptor, ErbB2. This preliminary data suggest that the balance neurogenesis/gliogenesis is critical for ENS development.