The RET proto-oncogene was identified as a major locus involved in Hirschsprung disease (HSCR). A genome-wide association study (GWAS) and whole exome sequencing identified NRG1 and NRG3 as additional HSCR susceptibility loci. We investigated the effects of RET (rs2506030 and rs2435357), NRG1 (rs2439302, rs16879552 and rs7835688) and NRG3 (rs10748842, rs10883866 and rs6584400) polymorphisms in a Chinese population with HSCR. We assessed single nucleotide polymorphisms (SNPs) in the RET, NRG1 and NRG3 genes in a cohort of 362 sporadic HSCR patients and 1,448 normal controls using a TaqMan genotyping assay. Significant associations were found between HSCR risk and rs2506030, rs2435357, rs2439302 and rs7835688 (odds ratio [OR] 1.64, P = 1.72E-06; 2.97, P = 5.15E-33; 1.84, P = 9.36E-11; and 1.93, P = 1.88E-12, respectively). Two locus analyses of SNPs indicated increased disease risks of HSCR between NRG1 rs2439302 and RET rs2435357 or rs2506030. RET rs2506030 (GG genotype) and rs2435357 (TT genotype), in combination with NRG1 rs2439302 (GG genotype), were strongly associated with the highest risk of HSCR (OR = 56.53, P = 4.50E-07) compared with the two loci or a single SNP of either RET or NRG1. Our results support the association between genetic variation of RET and NRG1 and susceptibility to HSCR in the Chinese population.

Hirschsprung's disease (HSCR) is a complex developmental defect characterized by the absence of enteric ganglia in the gastrointestinal tract. Although the genetic defect of enteric nervous system (ENS) was identified to play a critical role in the progress of HSCR, the systemic genetic dissection of HSCR still needs to be clarified. In this study, we firstly performed exome sequencing of two HSCR patients from a Han Chinese family, including the affected mother and son. After the initial quality filtering (coverage  ≥ 5X and SNP quality score ≥ 40) of the raw data, we identified 13,948 and 13,856 single nucleotide variants (SNVs), respectively. We subsequently compared the SNVs against public databases (dbSNP130, HapMap, and 1000 Genome Project) and obtained a total of 15 novel nonsynonymous SNVs in 15 genes, which were shared between these two patients. Follow-up Sanger sequencing and bioinformatics analysis highlighted variant c.853G>A (p.E285K) in NRG3, a gene involved in the development of ENS. In the validation phase, we sequenced all nine exons of NRG3 in 96 additional sporadic HSCR cases and 110 healthy individuals and identified another nonsynonymous variant c.1329G>A (p.M443I) and two synonymous variants c.828G>A (p.T276T) and c.1365T>A (p.P455P) only in the cases. Our results indicated that NRG3 may be a susceptibility gene for HSCR in a Chinese population.

No abstract