Fine mapping of the NRG1 Hirschsprung's disease locus

The interplay of common genetic variants NRG1 rs2439302 and RET rs2435357 increases the risk of developing Hirschsprung's disease

Introduction: As a congenital and genetically related disease, many single nucleotide polymorphisms (SNPs) have been reported to be associated with the risk of HSCR. Our previous research showed that SNP rs2439302 (NRG1) interacted with rs2435357 (RET) to increase the risk of HSCR development. However, the underlying molecular mechanism is still not well understood. Methods: SNP rs2439302 (NRG1) and rs2435357 (RET) were genotyped in 470 HSCR cases. The expression of NRG1 and RET was investigated in the colon of HSCR patients. Knockdown of the NRG1 and RET homologs was performed in zebrafish to investigate their synergistic effect on ENS development. The effect of SNP rs2439302 and rs2435357 polymorphism on neuron proliferation, migration, and differentiation were investigated in SHSY-5Y cells and IPSCs. Results: Significant downregulation of NRG1 and RET expression was noticed in the aganglionic segment of HSCR patients and SHSY-5Y cells with rs2439302 GG/rs2435357 TT genotype. NRG1 and RET double mutants caused the most severe reduction in enteric neuron numbers than NRG1 single mutant or RET single mutant in the hindgut of zebrafish. SHSY-5Y cells and IPSCs with rs2439302 GG/rs2435357 TT genotype exhibited a decreased proliferative, migration, and differentiative capacity. CTCF showed a considerably higher binding ability to SNP rs2439302 CC than GG. NRG1 reduction caused a further decrease in SOX10 expression via the PI3K/Akt pathway, which regulates RET expression by directly binding to rs2435357. Discussion: SNP rs2439302 (NRG1) GG increases the risk of developing HSCR by affecting the binding of transcription factor CTCF and interacting with rs2435357 (RET) to regulate RET expression via the PI3K/Akt/SOX10 pathway.

The impact of NRG1 expressions and methylation on multifactorial Hirschsprung disease

Background

Hirschsprung disease (HSCR) is a complex genetic disorder characterized by the lack of ganglion cells in the intestines. A current study showed that the NRG1 rare variant frequency in Indonesian patients with HSCR is only 0.9%. Here, we investigated the impact of NRG1 expressions and methylation patterns on the pathogenesis of HSCR.

Methods

This cross-sectional study determined NRG1 type I (HRGα, HRGβ1, HRGβ2, HRGβ3, HRGγ, and NDF43 isoforms), type II and type III expressions in both ganglionic and aganglionic colons of 20 patients with HSCR and 10 control colons by real-time polymerase chain reaction (qPCR). For methylation studies, we treated the extracted gDNA from 16 HSCR patients’ and 17 control colons with sodium bisulfate and analyzed the methylation pattern of NRG1 exon 1 with methylation-specific PCR. The samples were collected and analyzed at our institution from December 2018 to December 2020.

Results

NRG1 types I, II and III expressions were upregulated (17.2-, 3.2-, and 7.2-fold, respectively) in the ganglionic colons compared with control colons (type I: 13.32 ± 1.65 vs. 17.42 ± 1.51, p < 0.01; type II: 13.73 ± 2.02 vs. 16.29 ± 2.19, p < 0.01; type III: 13.47 ± 3.01 vs. 16.32 ± 2.58, p = 0.03; respectively); while only type I (7.7-fold) and HRGβ1/HRGβ2 (3.3-fold) isoforms were significantly upregulated in the aganglionic colons compared to the controls (type I: 14.47 ± 1.66 vs. 17.42 ± 1.51, p < 0.01; HRGβ1/HRGβ2: 13.62 ± 3.42 vs 14.75 ± 1.26, p = 0.01). Moreover, the frequency of partially methylated NRG1 was higher in the ganglionic (81%) and aganglionic (75%) colons than in the controls (59%).

Conclusions

Our study provides further insights into the aberrant NRG1 expression in the colons of patients with HSCR, both ganglionic and aganglionic bowel, which might contribute to the development of HSCR, particularly in Indonesia. Furthermore, these aberrant NRG1 expressions might be associated with its methylation pattern.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12887-022-03287-1.

Multiple, independent, common variants at RET, SEMA3 and NRG1 gut enhancers specify Hirschsprung disease risk in European ancestry subjects

PURPOSE

Hirschsprung disease (HSCR) is a developmental disorder of the enteric nervous system (ENS) characterized by congenital aganglionosis arising from coding variants in ENS genes causing partial or total loss-of-function. Low-penetrance, common, noncoding variants at RET, SEMA3 and NRG1 loci are also associated with HSCR, with small-to-moderate loss of gene expression mediated through sequence variants in cis-regulatory elements (CRE) as another causal mechanism. Since these latter variants are common, many individuals carry multiple risk variants. However, the extent and combinatorial effects of all putative CRE variants within and across these loci on HSCR is unknown.

METHODS

Using 583 HSCR subjects, one of the largest samples of European ancestry studied, and genotyping 56 tag variants, we evaluated association of all common variants overlapping putative gut CREs and fine-mapped causal variants at RET, SEMA3 and NRG1.

RESULTS

We demonstrate that 28 and 8 tag variants, several of which are genetically independent, overlap putative-enhancers at the RET and SEMA3 loci, respectively, as well as two fine-mapped tag variants at the NRG1 locus, are significantly associated with HSCR. Importantly, disease risk increases with increasing numbers of risk alleles from multiple variants within and across these loci, varying >25-fold across individuals.

CONCLUSION

This increasing allele number-dependent risk, we hypothesize, arises from HSCR-relevant ENS cells sensing the reduced gene expression at multiple ENS genes since their developmental effects are integrated through gene regulatory networks.

Roles of Enteric Neural Stem Cell Niche and Enteric Nervous System Development in Hirschsprung Disease

The development of the enteric nervous system (ENS) is highly modulated by the synchronized interaction between the enteric neural crest cells (ENCCs) and the neural stem cell niche comprising the gut microenvironment. Genetic defects dysregulating the cellular behaviour(s) of the ENCCs result in incomplete innervation and hence ENS dysfunction. Hirschsprung disease (HSCR) is a rare complex neurocristopathy in which the enteric neural crest-derived cells fail to colonize the distal colon. In addition to ENS defects, increasing evidence suggests that HSCR patients may have intrinsic defects in the niche impairing the extracellular matrix (ECM)-cell interaction and/or dysregulating the cellular niche factors necessary for controlling stem cell behaviour. The niche defects in patients may compromise the regenerative capacity of the stem cell-based therapy and advocate for drug- and niche-based therapies as complementary therapeutic strategies to alleviate/enhance niche-cell interaction. Here, we provide a summary of the current understandings of the role of the enteric neural stem cell niche in modulating the development of the ENS and in the pathogenesis of HSCR. Deciphering the contribution of the niche to HSCR may provide important implications to the development of regenerative medicine for HSCR.

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.

Hirschsprung disease is a congenital disorder characterized by the absence of intestinal neurons in the distal part of the intestine. It is a complex genetic disorder in which multiple variations in our genome combined, result in disease. One of these variations are Copy Number Variations (CNVs): large segments of our genome that are duplicated or deleted. Patients often have Hirschsprung disease without other symptoms. However, a proportion of patients has additional associated anatomical malformations and neurological symptoms. We found that CNVs, present in patients with associated anomalies, are more often larger compared to unaffected controls or Hirschsprung patients without other symptoms. Furthermore, Copy Number (CN) losses are enriched for constrained coding regions (CCR; genes usually not impacted by genomic alterations in unaffected controls) of which the expression is higher in the developing intestinal neurons compared to the intestine. We modelled loss of these candidate genes in zebrafish by disrupting the zebrafish orthologues by genome editing. For several genes this resulted in changes in intestinal neuron development, reminiscent of HSCR observed in patients. The results presented here highlight the importance of Copy Number profiling, zebrafish validation and evaluating all CCR expressed in developing intestinal neurons during diagnostic evaluation.

The Emerging Genetic Landscape of Hirschsprung Disease and Its Potential Clinical Applications

Hirschsprung disease (HSCR) is the leading cause of neonatal functional intestinal obstruction. It is a rare congenital disease with an incidence of one in 3,500-5,000 live births. HSCR is characterized by the absence of enteric ganglia in the distal colon, plausibly due to genetic defects perturbing the normal migration, proliferation, differentiation, and/or survival of the enteric neural crest cells as well as impaired interaction with the enteric progenitor cell niche. Early linkage analyses in Mendelian and syndromic forms of HSCR uncovered variants with large effects in major HSCR genes including RET, EDNRB, and their interacting partners in the same biological pathways. With the advances in genome-wide genotyping and next-generation sequencing technologies, there has been a remarkable progress in understanding of the genetic basis of HSCR in the past few years, with common and rare variants with small to moderate effects being uncovered. The discovery of new HSCR genes such as neuregulin and BACE2 as well as the deeper understanding of the roles and mechanisms of known HSCR genes provided solid evidence that many HSCR cases are in the form of complex polygenic/oligogenic disorder where rare variants act in the sensitized background of HSCR-associated common variants. This review summarizes the roadmap of genetic discoveries of HSCR from the earlier family-based linkage analyses to the recent population-based genome-wide analyses coupled with functional genomics, and how these discoveries facilitated our understanding of the genetic architecture of this complex disease and provide the foundation of clinical translation for precision and stratified medicine.

Common variants of NRG1 and ITGB4 confer risk of Hirschsprung disease in Han Chinese population

BACKGROUND

Hirschsprung disease (HSCR) is a neurodevelopmental disorder with a strong genetic component. Common variants of NRG1 contributed to HSCR risk in Asians, and rare variants of ERBB2 and ITGB4 were found to be associated with HSCR. ERBB2 and ITGB4 are partners of Nrg1/ErbB pathway, which is important in HSCR pathogenesis. We aimed to investigate whether common variants in NRG1, ERBB2 and ITGB4 were associated with HSCR in Chinese Han population.

METHODS

We genotype 17 single nucleotide polymorphisms (SNPs) of NRG1, ERBB2 and ITGB4 in 420 HSCR patients and 1665 controls, and performed association analysis.

RESULTS

We validated associations of two NRG1 SNPs rs7835688 (P_Allelic = 2.2 × 10^-20, OR = 2.21, 95%CI = 1.86-2.62) and rs16879552 (P_Allelic = 5.6 × 10^-9, OR = 1.57, 95%CI = 1.35-1.83) with risk to HSCR. SNP rs3744000 located 5' upstream of ITGB4 showed association with HSCR (P_Allelic = 2.4 × 10^-3, OR = 1.27, 95%CI = 1.09-1.49). Four SNPs of ERBB2 exhibited no association.

CONCLUSIONS

Our results suggested that common variation of ITGB4 and NRG1 conferred risk to HSCR in Chinese Han population, which further highlighted Nrg-1/ErbB pathway involving in the pathogenesis of HSCR.

Advances in understanding the association between Down syndrome and Hirschsprung disease (DS-HSCR)

The clinical association between Trisomy 21 (Down syndrome) and aganglionosis (Hirschsprung disease; DS-HSCR) is well-established, being of the order of 5% and remains the most common congenital association with Hirschsprung disease. However, little consensus exists as to the possible etiologic and genetic factors influencing this association. Recent research has identified a number of levels at which development of the enteric nervous system is potentially affected in Trisomy 21. These include a decreased central pool of available neuroblasts for migration into the enteric nervous system, abnormal neuroblast type, poor synaptic nerve function and early germline gene-related influences on the migrating neuroblasts due to genetic mutations of a number of important developmental genes, and possible somatic mutations resulting from alterations in the local tissue microenvironment. In this paper, we review available evidence for this association. In addition, we provide evidence of both germline and somatic gene mutations suggesting causation. Although the picture is complex, recent associations between specific RET proto-oncogene variations have been shown to be significant in Down syndrome patients with Hirschsprung disease, as they probably interfere with vital RET functions in the development of the autonomic and enteric nervous systems, increasing the risk of disturbed normal function. In addition, we explore potential role of other facilitatory influence of other susceptibility genes as well as potential other chromosome 21 gene actions and the microenvironment on the Down syndrome gastro-intestinal tract. The various ways in which trisomy of chromosome influences the enteric nervous system are becoming clearer. The sum of these effects influences the outcome of surgery in Down syndrome patients with Hirschsprung Disease.

Effects of NRG1 Polymorphisms on Hirschsprung's Disease Susceptibility: A Meta-analysis

Substantial resources have been devoted to evaluate the relationship between NRG1 variants rs7835688 and rs16879552 and Hirschsprung’s Disease (HSCR) but no consistency exists. This meta-analysis aimed to assess the association between the two SNPs and HSCR. PubMed, EMBASE, and Chinese Biological Medicine databases were searched for studies potentially eligible up to March, 2017. The summary odds ratios (ORs) with 95% CIs were calculated from different genetic models. Nine case-control studies (8 for both and 1 for rs16879552 only) involving 1984 HSCR patients and 4220 controls were identified. The combined results showed a significant association between HSCR risk and rs7835688 in all genetic models (per-allele model: OR = 1.66, 95% CI = 1.35–2.05; P = 1.940E-06). Rs16879552 was significantly associated with HSCR in per-allele (OR = 1.50, 95% CI = 1.27–1.76; P = 1.087E-06), additive and recessive model, except for dominant model. Stratified analysis by ethnicity showed that rs7835688 and rs16879552 were only causative for Asians, but not risk locus for Caucasians. Furthermore, pooled data based on segment length indicated that individuals with rs7835688 experienced a significantly higher risk for short-segment HSCR in all genotypes; but rs16879552 was only found to be associated with long-segment HSCR/ total colonic aganglionosis at the allele level.

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.

Neuregulin 1 is involved in enteric nervous system development in zebrafish

BACKGROUND

Hirschsprung's disease (HD, also known as congenital colon aganglionosis) is a congenital disorder characterized by the absence of intramural ganglion cells in the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. Recent studies have indicated neuregulin 1 (NRG1) as a new candidate gene involved in the development of the enteric nervous system (ENS) in humans.

METHODS

In our study, we investigated the role of NRG1 in zebrafish ENS development by assessing NRG1 expression patterns during ENS development. Knockdown, overexpression and rescue zebrafish models of NRG1 were created to evaluate differences in phenotype, numbers of enteric neurons, ENS-related factors and nerve fiber arrangements.

RESULTS

NRG1 was expressed in zebrafish intestine at both the larval and adult stage. We also found that decreased expression of NRG1 resulted in reductions in enteric neuron number and decreased expression of ENS development markers. Moreover, NRG1-knockdown zebrafish exhibited a disordered arrangement of nerve fibers.

CONCLUSIONS

Collectively, these results demonstrated that NRG1 expression might play a role in zebrafish ENS development. In addition, by modulating NRG1 expression, we created a model that may be useful for investigating the mechanism underlying HD pathogenesis.

Genetic Background of Hirschsprung Disease: A Bridge Between Basic Science and Clinical Application

Hirschsprung's disease (HSCR) is a congenital disorder, defined by partial or complete loss of the neuronal ganglion cells in the intestinal tract, which is caused by the failure of neural crest cells to migrate completely during intestinal development during fetal life. HSCR has a multifactorial etiology, and genetic factors play a key role in its pathogenesis; these include mutations within several gene loci. These have been identified by screening candidate genes, or by conducting genome wide association (GWAS) studies. However, only a small portion of them have been proposed as major genetic risk factors for the HSCR. In this review, we focus on those genes that have been identified as either low penetrant or high penetrant variants that determine the risk of Hirschsprung's disease. J. Cell. Biochem. 119: 28-33, 2018. © 2017 Wiley Periodicals, Inc.

Effects of RET, NRG1 and NRG3 Polymorphisms in a Chinese Population with Hirschsprung Disease

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: clinical dysmorphology, genes, micro-RNAs, and future perspectives

On the occasion of the 100th anniversary of Dr. Harald Hirschsprung's death, there is a worldwide significant research effort toward identifying and understanding the role of genes and biochemical pathways involved in the pathogenesis as well as the use of new therapies for the disease harboring his name (Hirschsprung disease, HSCR). HSCR (aganglionic megacolon) is a frequent diagnostic and clinical challenge in perinatology and pediatric surgery, and a major cause of neonatal intestinal obstruction. HSCR is characterized by the absence of ganglia of the enteric nervous system, mostly in the distal gastrointestinal tract. This review focuses on current understanding of genes and pathways associated with HSCR and summarizes recent knowledge related to micro RNAs (miRNAs) and HSCR pathogenesis. While commonly sporadic, Mendelian patterns of inheritance have been described in syndromic cases with HSCR. Although only half of the patients with HSCR have mutations in specific genes related to early embryonic development, recent pathway-based analysis suggests that gene modules with common functions may be associated with HSCR in different populations. This comprehensive profile of functional gene modules may serve as a useful resource for future developmental, biochemical, and genetic studies providing insights into the complex nature of HSCR.

Association Analysis of SLC6A20 Polymorphisms With Hirschsprung Disease

PURPOSE

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder, which is caused by no neuronal ganglion cells in part or all of distal gastrointestinal tract. Recently, our genome-wide association study has identified solute carrier family 6, proline IMINO transporter, member 20 (SLC6A20) as one of the potential risk factors for HSCR development. This study performed a replication study for the association of SLC6A20 polymorphisms with HSCR and an extended analysis to investigate further associations for subgroups and haplotypes.

METHODS

For the replication study, a total of 40 single nucleotide polymorphisms (SNPs) of SLC6A20 were genotyped in 187 HSCR subjects composed of 121 short-segment HSCR, 45 long-segment HSCR (L-HSCR), 21 total colonic aganglionosis, and 283 unaffected controls. Imputation was performed using genotype data from our genome-wide association study and this replication study.

RESULTS

Imputed meta-analysis revealed that 13 SLC6A20 SNPs (minimum P = 0.0002 at rs6770261) were significantly associated with HSCR even after correction for multiple comparisons using false discovery rate (FDR) (minimum PFDR =  .005). In further subgroup analysis, SLC6A20 polymorphisms appeared to have increased associations with L-HSCR. Moreover, haplotype analysis also showed significant associations between 2 haplotypes (BL3_ht2 and BL4_ht2) and HSCR susceptibility (PFDR <  .05).

CONCLUSIONS

Although further replications and functional evaluations are required, our results suggest that SLC6A20 may have roles in HSCR development and in the extent of aganglionic segment during enteric nervous system development.

Systematic assessment of imputation performance using the 1000 Genomes reference panels

Genotype imputation has been widely adopted in the postgenome-wide association studies (GWAS) era. Owing to its ability to accurately predict the genotypes of untyped variants, imputation greatly boosts variant density, allowing fine-mapping studies of GWAS loci and large-scale meta-analysis across different genotyping arrays. By leveraging genotype data from 90 whole-genome deeply sequenced individuals as the evaluation benchmark and the 1000 Genomes Project data as reference panels, we systematically examined four important issues related to genotype imputation practice. First, in a study of imputation accuracy, we found that IMPUTE2 and minimac have the best imputation performance among the three popular imputing software evaluated and that using a multi-population reference panel is beneficial. Second, the optimal imputation quality cutoff for removing poorly imputed variants varies according to the software used. Third, the major contributing factors to consistently poor imputation are low variant heterozygosity, high sequence similarity to other genomic regions, high GC content, segmental duplication and being far from genotyping markers. Lastly, in an evaluation of the imputability of all known GWAS regions, we found that GWAS loci associated with hematological measurements and immune system diseases are harder to impute, as compared with other human traits. Recommendations made based on the above findings may provide practical guidance for imputation exercise in future genetic studies.

Clinical and genetic correlations of familial Hirschsprung's disease

BACKGROUND

The risk of familial transmission in Hirschsprung's disease (HSCR) currently lacks correlation between the clinical phenotype and the underlying genetic factors. The aim of this study was to clinically evaluate familial HSCR transmission and to correlate with the genetic background.

METHODS

Clinical and gene analysis of familial HSCR patients were explored. DNA from 45 patients (35 kindreds) was screened for genetic variations of the RET, and EDNRB genes were screened for genetic variation by semi-automated bi-directional sequencing analysis and matched to controls.

MAIN RESULTS

Male:female ratio (3:1) had a female proband in 4 families. Aganglionosis was significantly more frequent with total colonic aganglionosis (TCA) in 40% familial cases (viz: 17/43 (43%) vs. 19/342 non-familial patients (5.6%) (p<0.01)). Transmission of S-HSCR was observed in 13 (31%), which was associated with EDNRB variation. RET gene promoter variation correlated with extended aganglionosis in 6/35 kindreds (17%). In 3 kindreds, both significant EDNRB and RET mutations were identified and where present were associated with increased penetrance in succeeding generations. An increased penetrance with succeeding generations occurred in 6 (14%). In a further 3 generation family, extensive variations in exon 6, 13, and 18 affected 3 males with progressive penetration and aganglionic length, including total intestinal aganglionosis in the further offspring. RET and MEN association was noted in 5 kindreds (14.3%) related to RET variations at Cysteine sites.

CONCLUSIONS

Cumulative effects of the RET and EDNRB genes contribute to long-segment and total colonic aganglionosis.

A genome-wide association study identifies potential susceptibility loci for Hirschsprung disease

Hirschsprung disease (HSCR) is a congenital and heterogeneous disorder characterized by the absence of intramural nervous plexuses along variable lengths of the hindgut. Although RET is a well-established risk factor, a recent genome-wide association study (GWAS) of HSCR has identified NRG1 as an additional susceptibility locus. To discover additional risk loci, we performed a GWAS of 123 sporadic HSCR patients and 432 unaffected controls using a large-scale platform with coverage of over 1 million polymorphic markers. The result was that our study replicated the findings of RET-CSGALNACT2-RASGEF1A genomic region (_rawP = 5.69×10⁻¹⁹ before a Bonferroni correction; _corrP = 4.31×10⁻¹³ after a Bonferroni correction) and NRG1 as susceptibility loci. In addition, this study identified SLC6A20 (_adjP = 2.71×10⁻⁶), RORA (_adjP = 1.26×10⁻⁵), and ABCC9 (_adjP = 1.86×10⁻⁵) as new potential susceptibility loci under adjusting the already known loci on the RET-CSGALNACT2-RASGEF1A and NRG1 regions, although none of the SNPs in these genes passed the Bonferroni correction. In further subgroup analysis, the RET-CSGALNACT2-RASGEF1A genomic region was observed to have different significance levels among subgroups: short-segment (S-HSCR, _corrP = 1.71×10⁻⁵), long-segment (L-HSCR, _corrP = 6.66×10⁻⁴), and total colonic aganglionosis (TCA, _corrP>0.05). This differential pattern in the significance level suggests that other genomic loci or mechanisms may affect the length of aganglionosis in HSCR subgroups during enteric nervous system (ENS) development. Although functional evaluations are needed, our findings might facilitate improved understanding of the mechanisms of HSCR pathogenesis.

Where genotype is not predictive of phenotype: towards an understanding of the molecular basis of reduced penetrance in human inherited disease

Some individuals with a particular disease-causing mutation or genotype fail to express most if not all features of the disease in question, a phenomenon that is known as 'reduced (or incomplete) penetrance'. Reduced penetrance is not uncommon; indeed, there are many known examples of 'disease-causing mutations' that fail to cause disease in at least a proportion of the individuals who carry them. Reduced penetrance may therefore explain not only why genetic diseases are occasionally transmitted through unaffected parents, but also why healthy individuals can harbour quite large numbers of potentially disadvantageous variants in their genomes without suffering any obvious ill effects. Reduced penetrance can be a function of the specific mutation(s) involved or of allele dosage. It may also result from differential allelic expression, copy number variation or the modulating influence of additional genetic variants in cis or in trans. The penetrance of some pathogenic genotypes is known to be age- and/or sex-dependent. Variable penetrance may also reflect the action of unlinked modifier genes, epigenetic changes or environmental factors. At least in some cases, complete penetrance appears to require the presence of one or more genetic variants at other loci. In this review, we summarize the evidence for reduced penetrance being a widespread phenomenon in human genetics and explore some of the molecular mechanisms that may help to explain this enigmatic characteristic of human inherited disease.

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.

Tcof1 acts as a modifier of Pax3 during enteric nervous system development and in the pathogenesis of colonic aganglionosis

Hirschsprung disease (HSCR) is a human congenital disorder, defined by the absence of ganglia from variable lengths of the colon. These ganglia comprise the enteric nervous system (ENS) and are derived from migratory neural crest cells (NCCs). The inheritance of HSCR is complex, often non-Mendelian and characterized by variable penetrance. Although extensive research has identified many key players in the pathogenesis of Hirschsprung disease, a large number of cases remain genetically undefined. Therefore, additional unidentified genes or modifiers must contribute to the etiology and pathogenesis of Hirschsprung disease. We have discovered that Tcof1 may be one such modifier. Haploinsufficiency of Tcof1 in mice results in a reduction of vagal NCCs and their delayed migration along the length of the gut during early development. This alone, however, is not sufficient to cause colonic aganglionosis as alterations in the balance of NCC proliferation and differentiation ensures NCC colonize the entire length of the gut of Tcof1(+/-) mice by E18.5. In contrast, Tcof1 haploinsufficiency is able to sensitize Pax3(+/-) mice to colonic aganglionosis. Although, Pax3 heterozygous mice do not show ENS defects, compound Pax3;Tcof1 heterozygous mice exhibit cumulative apoptosis which severely reduces the NCC population that migrates into the foregut. In addition, the proliferative capacity of these NCC is also diminished. Taken together with the opposing effects of Pax3 and Tcof1 on NCC differentiation, the synergistic haploinsufficiency of Tcof1 and Pax3 results in colonic aganglionosis in mice and may contribute to the pathogenesis of Hirschsprung disease.

Comprehensive analysis of NRG1 common and rare variants in Hirschsprung patients

Hirschsprung disease (HSCR, OMIM 142623) is a developmental disorder characterized by the absence of ganglion cells along variable lengths of the distal gastrointestinal tract, which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. The RET proto-oncogene is the major gene for HSCR with differential contributions of its rare and common, coding and noncoding mutations to the multifactorial nature of this pathology. Many other genes have been described to be associated with the pathology, as NRG1 gene (8p12), encoding neuregulin 1, which is implicated in the development of the enteric nervous system (ENS), and seems to contribute by both common and rare variants. Here we present the results of a comprehensive analysis of the NRG1 gene in the context of the disease in a series of 207 Spanish HSCR patients, by both mutational screening of its coding sequence and evaluation of 3 common tag SNPs as low penetrance susceptibility factors, finding some potentially damaging variants which we have functionally characterized. All of them were found to be associated with a significant reduction of the normal NRG1 protein levels. The fact that those mutations analyzed alter NRG1 protein would suggest that they would be related with HSCR disease not only in Chinese but also in a Caucasian population, which reinforces the implication of NRG1 gene in this pathology.