The major gene for Hirschsprung disease (HSCR) encodes the receptor tyrosine kinase RET. In a study of 690 European- and 192 Chinese-descent probands and their parents or controls, we demonstrate the ubiquity of a >4-fold susceptibility from a C-->T allele (rs2435357: p = 3.9 x 10(-43) in European ancestry; p = 1.1 x 10(-21) in Chinese samples) that probably arose once within the intronic RET enhancer MCS+9.7. With in vitro assays, we now show that the T variant disrupts a SOX10 binding site within MCS+9.7 that compromises RET transactivation. The T allele, with a control frequency of 20%-30%/47% and case frequency of 54%-62%/88% in European/Chinese-ancestry individuals, is involved in all forms of HSCR. It is marginally associated with proband gender (p = 0.13) and significantly so with length of aganglionosis (p = 7.6 x 10(-5)) and familiality (p = 6.2 x 10(-4)). The enhancer variant is more frequent in the common forms of male, short-segment, and simplex families whereas multiple, rare, coding mutations are the norm in the less common and more severe forms of female, long-segment, and multiplex families. The T variant also increases penetrance in patients with rare RET coding mutations. Thus, both rare and common mutations, individually and together, make contributions to the risk of HSCR. The distribution of RET variants in diverse HSCR patients suggests a "cellular-recessive" genetic model where both RET alleles' function is compromised. The RET allelic series, and its genotype-phenotype correlations, shows that success in variant identification in complex disorders may strongly depend on which patients are studied.

Innervation of the gut is segmentally lost in Hirschsprung disease (HSCR), a consequence of cell-autonomous and non-autonomous defects in enteric neuronal cell differentiation, proliferation, migration, or survival. Rare, high-penetrance coding variants and common, low-penetrance non-coding variants in 13 genes are known to underlie HSCR risk, with the most frequent variants in the ret proto-oncogene (RET). We used a genome-wide association (220 trios) and replication (429 trios) study to reveal a second non-coding variant distal to RET and a non-coding allele on chromosome 7 within the class 3 Semaphorin gene cluster. Analysis in Ret wild-type and Ret-null mice demonstrates specific expression of Sema3a, Sema3c, and Sema3d in the enteric nervous system (ENS). In zebrafish embryos, sema3 knockdowns show reduction of migratory ENS precursors with complete ablation under conjoint ret loss of function. Seven candidate receptors of Sema3 proteins are also expressed within the mouse ENS and their expression is also lost in the ENS of Ret-null embryos. Sequencing of SEMA3A, SEMA3C, and SEMA3D in 254 HSCR-affected subjects followed by in silico protein structure modeling and functional analyses identified five disease-associated alleles with loss-of-function defects in semaphorin dimerization and binding to their cognate neuropilin and plexin receptors. Thus, semaphorin 3C/3D signaling is an evolutionarily conserved regulator of ENS development whose dys-regulation is a cause of enteric aganglionosis.

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.

BackgroundHirschsprung 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.ResultsWe 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.ConclusionsOur 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.Electronic supplementary materialThe online version of this article (doi:10.1186/s13059-017-1174-6) contains supplementary material, which is available to authorized users.

No abstract

No abstract

SummaryThe RET proto-oncogene is the major gene involved in the complex genetics of Hirschsprung disease (HSCR), or aganglionic megacolon, showing causative loss-of-function mutations in 15-30% of the sporadic cases. Several RET polymorphisms and haplotypes have been described in association with the disease, suggesting a role for this gene in HSCR predisposition, also in the absence of mutations in the coding region. Finally, the presence of a functional variant in intron 1 has repeatedly been proposed to explain such findings. Here we report a casecontrol study conducted on 97 Italian HSCR sporadic patients and 85 population matched controls, using 13 RET polymorphisms distributed throughout the gene, from the basal promoter to the 3 UTR. Linkage disequilibrium and haplotype analyses have shown increased recombination between the 5 and 3 portions of the gene and an over-representation, in the cases studied, of two haplotypes sharing a common allelic combination that extends from the promoter up to intron 5. We propose that these two disease-associated haplotypes derive from a single founding locus, extending up to intron 19 and successively rearranged in correspondence with a high recombination rate region located between the proximal and distal portions of the gene. Our results suggests the possibility that a common HSCR predisposing variant, in linkage disequilibrium with such haplotypes, is located further downstream than the previously suggested interval encompassing intron 1.

Hirschsprung disease (HSCR) is an inherited disorder characterised by absence of intrinsic ganglion cells in the distal gastrointestinal tract. Different susceptibility genes, involved in either the Ret-tyrosine kinase or the endothelin signalling pathways, contribute to HSCR phenotype. Interestingly, alterations of these genes are detected in only 30-50% of all HSCR patients, suggesting the involvement of modifier genes and/or additional genetic or environmental risk factors. In complex disorders common polymorphic variants can be associated with the disease phenotype, thus modifying the risk of recurrence. To investigate whether sequence variants of the RET proto-oncogene may be associated with the development of the HSCR phenotype, we analysed 92 Italian patients for the 2508C > T synonymous substitution in exon 14 (S836S) finding that the T allele is clearly less frequent than in control individuals (Fisher exact test P = 0.0002). On the other hand, this RET variant allele is overrepresented in patients affected with medullary thyroid carcinoma. Assuming a direct effect of this single-nucleotide polymorphism in predisposing to RET associated pathologies, we have performed functional tests which excluded any possible involvement of the C and T alleles in DNA-protein binding, transcript stability and RNA splicing and editing.