Gene Network Analysis of Candidate Loci for Human Anorectal Malformations

Anorectal malformations (ARMs) are birth defects that require surgery and carry significant chronic morbidity. Our earlier genome-wide copy number variation (CNV) study had provided a wealth of candidate loci. To find out whether these candidate loci are related to important developmental pathways, we have performed an extensive literature search coupled with the currently available bioinformatics tools. This has allowed us to assign both genic and non-genic CNVs to interrelated pathways known to govern the development of the anorectal region. We have linked 11 candidate genes to the WNT signalling pathway and 17 genes to the cytoskeletal network. Interestingly, candidate genes with similar functions are disrupted by the same type of CNV. The gene network we discovered provides evidence that rare mutations in different interrelated genes may lead to similar phenotypes, accounting for genetic heterogeneity in ARMs. Classification of patients according to the affected pathway and lesion type should eventually improve the diagnosis and the identification of common genes/molecules as therapeutic targets.     

Fine mapping of the NRG1 Hirschsprung's disease locus

The primary pathology of Hirschsprung's disease (HSCR, colon aganglionosis) is the absence of ganglia in variable lengths of the hindgut, resulting in functional obstruction. HSCR is attributed to a failure of migration of the enteric ganglion precursors along the developing gut. RET is a key regulator of the development of the enteric nervous system (ENS) and the major HSCR-causing gene. Yet the reduced penetrance of RET DNA HSCR-associated variants together with the phenotypic variability suggest the involvement of additional genes in the disease. Through a genome-wide association study, we uncovered a ～350 kb HSCR-associated region encompassing part of the neuregulin-1 gene (NRG1). To identify the causal NRG1 variants contributing to HSCR, we genotyped 243 SNPs variants on 343 ethnic Chinese HSCR patients and 359 controls. Genotype analysis coupled with imputation narrowed down the HSCR-associated region to 21 kb, with four of the most associated SNPs (rs10088313, rs10094655, rs4624987, and rs3884552) mapping to the NRG1 promoter. We investigated whether there was correlation between the genotype at the rs10088313 locus and the amount of NRG1 expressed in human gut tissues (40 patients and 21 controls) and found differences in expression as a function of genotype. We also found significant differences in NRG1 expression levels between diseased and control individuals bearing the same rs10088313 risk genotype. This indicates that the effects of NRG1 common variants are likely to depend on other alleles or epigenetic factors present in the patients and would account for the variability in the genetic predisposition to HSCR.     

RET mutational spectrum in Hirschsprung disease: evaluation of 601 Chinese patients

Rare (RVs) and common variants of the RET gene contribute to Hirschsprung disease (HSCR; congenital aganglionosis). While RET common variants are strongly associated with the commonest manifestation of the disease (males; short-segment aganglionosis; sporadic), rare coding sequence (CDS) variants are more frequently found in the lesser common and more severe forms of the disease (females; long/total colonic aganglionosis; familial).Here we present the screening for RVs in the RET CDS and intron/exon boundaries of 601 Chinese HSCR patients, the largest number of patients ever reported. We identified 61 different heterozygous RVs (50 novel) distributed among 100 patients (16.64%). Those include 14 silent, 29 missense, 5 nonsense, 4 frame-shifts, and one in-frame amino-acid deletion in the CDS, two splice-site deletions, 4 nucleotide substitutions and a 22-bp deletion in the intron/exon boundaries and 1 single-nucleotide substitution in the 5' untranslated region. Exonic variants were mainly clustered in RET the extracellular domain. RET RVs were more frequent among patients with the most severe phenotype (24% vs. 15% in short-HSCR). Phasing RVs with the RET HSCR-associated haplotype suggests that RVs do not underlie the undisputable association of RET common variants with HSCR. None of the variants were found in 250 Chinese controls.     

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.     

Prokineticin-1 (Prok-1) works coordinately with glial cell line-derived neurotrophic factor (GDNF) to mediate proliferation and differentiation of enteric neural crest cells

Enteric neural crest cells (NCC) are multipotent progenitors which give rise to neurons and glia of the enteric nervous system (ENS) during fetal development. Glial cell line-derived neurotrophic factor (GDNF)/RET receptor tyrosine kinase (Ret) signaling is indispensable for their survival, migration and differentiation. Using microarray analysis and isolated NCCs, we found that 45 genes were differentially expressed after GDNF treatment (16 h), 29 of them were up-regulated including 8 previously undescribed genes. Prokineticin receptor 1 (PK-R1), a receptor for Prokineticins (Prok), was identified in our screen and shown to be consistently up-regulated by GDNF in enteric NCCs. Further, PK-R1 was persistently expressed at a lower level in the enteric ganglions of the c-Ret deficient mice when compared to that of the wild-type littermates. Subsequent functional analysis showed that GDNF potentiated the proliferative and differentiation effects of Prok-1 by up-regulating PK-R1 expression in enteric NCCs. In addition, expression analysis and gene knock-down experiments indicated that Prok-1 and GDNF signalings shared some common downstream targets. More importantly, Prok-1 could induce both proliferation and expression of differentiation markers of c-Ret deficient NCCs, suggesting that Prok-1 may also provide a complementary pathway to GDNF signaling. Taken together, these findings provide evidence that Prok-1 crosstalks with GDNF/Ret signaling and probably provides an additional layer of signaling refinement to maintain proliferation and differentiation of enteric NCCs.     

Genetic Analyses of a Three Generation Family Segregating Hirschsprung Disease and Iris Heterochromia

We present the genetic analyses conducted on a three-generation family (14 individuals) with three members affected with isolated-Hirschsprung disease (HSCR) and one with HSCR and heterochromia iridum (syndromic-HSCR), a phenotype reminiscent of Waardenburg-Shah syndrome (WS4). WS4 is characterized by pigmentary abnormalities of the skin, eyes and/or hair, sensorineural deafness and HSCR. None of the members had sensorineural deafness. The family was screened for copy number variations (CNVs) using Illumina-HumanOmni2.5-Beadchip and for coding sequence mutations in WS4 genes (EDN3, EDNRB, or SOX10) and in the main HSCR gene (RET). Confocal microscopy and immunoblotting were used to assess the functional impact of the mutations. A heterozygous A/G transition in EDNRB was identified in 4 affected and 3 unaffected individuals. While in EDNRB isoforms 1 and 2 (cellular receptor) the transition results in the abolishment of translation initiation (M1V), in isoform 3 (only in the cytosol) the replacement occurs at Met91 (M91V) and is predicted benign. Another heterozygous transition (c.-248G/A; -predicted to affect translation efficiency-) in the 5′-untranslated region of EDN3 (EDNRB ligand) was detected in all affected individuals but not in healthy carriers of the EDNRB mutation. Also, a de novo CNVs encompassing DACH1 was identified in the patient with heterochromia iridum and HSCRSince the EDNRB and EDN3 variants only coexist in affected individuals, HSCR could be due to the joint effect of mutations in genes of the same pathway. Iris heterochromia could be due to an independent genetic event and would account for the additional phenotype within the family.     

Mutations in the NRG1 gene are associated with Hirschsprung disease

Hirschsprung disease (HSCR, congenital colon aganglionosis) is a relatively common complex genetic condition caused by abnormal development of the enteric nervous system (ENS). Through a recent genome-wide association study conducted on Chinese HSCR patients, we identified a new HSCR contributing locus, neuregulin 1 (NRG1; 8p12), a gene known to be involved in the development of the ENS. As genes in which disease-associated common variants are found are to be considered as candidates for the search of deleterious rare variants (RVs) in the coding sequences, we sequenced the NRG1 exons of 358 sporadic HSCR patients and 333 controls. We identified a total of 13 different heterozygous RVs including 8 non-synonymous (A28G, E134K, V266L, H347Y, P356L, V486M, A511T, P608A) and 3 synonymous amino acid substitutions (P24P, T169T, L483L), a frameshift (E239fsX10), and a c.503-4insT insertion. Functional analysis of the most conserved non-synonymous substitutions, H347Y and P356L, showed uneven intracellular distribution and aberrant expression of the mutant proteins. Except for T169T and V486M, all variants were exclusive to HSCR patients. Overall, there was a statistically significant over-representation of NRG1 RVs in HSCR patients (p?=?0.008). We show here that not only common, but also rare variants of the NRG1 gene contribute to HSCR. This strengthens the role of NRG1.