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

Artificial intelligence-based approaches for the detection and prioritization of genomic mutations in congenital surgical diseases

Genetic mutations are critical factors leading to congenital surgical diseases and can be identified through genomic analysis. Early and accurate identification of genetic mutations underlying these conditions is vital for clinical diagnosis and effective treatment. In recent years, artificial intelligence (AI) has been widely applied for analyzing genomic data in various clinical settings, including congenital surgical diseases. This review paper summarizes current state-of-the-art AI-based approaches used in genomic analysis and highlighted some successful applications that deepen our understanding of the etiology of several congenital surgical diseases. We focus on the AI methods designed for the detection of different variant types and the prioritization of deleterious variants located in different genomic regions, aiming to uncover susceptibility genomic mutations contributed to congenital surgical disorders.

miR-938 rs2505901 T>C polymorphism increases Hirschsprung disease risk: a case-control study of Chinese children

Aim: To explore the association between miR-938 rs2505901 T>C polymorphism and Hirschsprung disease (HSCR) risk in Chinese children. Materials & Methods: We conducted a case-control study in a Chinese population with 1381 cases and 1457 controls. The associated correlation strengths were assessed by adjusted odds ratios (AORs) and 95% CIs. Results: The results revealed that the rs2505901 TC and rs2505901 TC/CC genotype were related to an increased HSCR risk compared to the risk contributed by the rs2505901 TT genotype. A stratification analysis showed that the rs2505901 TC/CC genotype promoted the progression of HSCR more significantly in patients with the short-segment HSCR subtype. Conclusion: Our study indicated that miR-938 rs2505901 T>C polymorphism is significantly associated with HSCR risk in Chinese children. This result needs to be confirmed with well-designed studies.

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.

NOX5 is expressed aberrantly but not a critical pathogenetic gene in Hirschsprung disease

Background

Hirschsprung disease (HSCR) is a congenital disorder characterized by the absence of intramural ganglion cells in the distal gastrointestinal tract (GI), which results in tonic contraction of the aganglionic gut segment and functional intestinal obstruction. Recent studies have suggested NADPH oxidase 5 (NOX5) as a candidate risk gene for HSCR. In this study, we examined the function of NOX5 to verify its role in the development of the enteric nervous system (ENS).

Methods

HSCR tissue specimens (n = 10) were collected at the time of pull-through surgery and control specimens (n = 10) were obtained at the time of colostomy closure in patients. The NOX5 expression in aganglionic and ganglionic segments of HSCR colon and normal colon were analyzed by immunohistochemistry (IHC), western blot and real-time quantitative PCR (qPCR). The gene expression levels and spatiotemporal expression spectrum of NOX5 in different development stages of zebrafish embryo were determined using qPCR and in-situ hybridization (ISH). The enteric nervous system in NOX5 Morpholino (MO) knockdown and wild type (WT) zebrafish embryo was analyzed by whole-mount immunofluorescence (IF). Intestinal transit assay was performed to analyze the gastrointestinal motility in NOX5 knockdown and control larvae.

Results

NOX5 is strongly expressed in the ganglion cells in the proximal segment of HSCR colons and all segments of normal colons. Moreover, the expression of NOX5 is markedly decreased in the aganglionic segment of HSCR colon compared to the ganglionic segment. In zebrafish, NOX5 mRNA level is the highest in the one cell stage embryos and it is decreased overtime with the development of the embryos. Interestingly, the expression of NOX5 appears to be enriched in the nervous system. However, the number of neurons in the GI tract and the GI motility were not affected upon NOX5 knockdown.

Conclusions

Our study shows that NOX5 markedly decreased in the aganglionic segment of HSCR but didn’t involve in the ENS development of zebrafish. It implies that absence of intestinal ganglion cells may lead to down-regulation of NOX5.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12887-021-02611-5.

Phenotypic spectrum associated with pathogenic mutation in the NRG1 gene in Acadian family

NRG1 is a gene that encodes for a protein that binds to a receptor of the tyrosine kinase family which is essential for the survival of the central nervous system development during embryogenesis. Mutation of the NRG1 gene causes aganglionosis, which leads to Hirschsprung disease. Two brothers of Acadian descent presented with a history of Hirschsprung disease, in association with other anomalies including congenital heart disease, learning difficulties, developmental issues, and hypopigmented hair patch. Molecular analysis in both siblings revealed a heterozygous pathogenic mutation in the NGR1 gene (c.235C>T [p.Arg79*]), that was inherited from an unaffected father. This family expands our knowledge about the phenotypic spectrum associated with pathogenic mutation in the NRG1 gene with intrafamilial variability and the likely reduced penetrance for the phenotypic expression.

A random forest-based framework for genotyping and accuracy assessment of copy number variations

Detection of copy number variations (CNVs) is essential for uncovering genetic factors underlying human diseases. However, CNV detection by current methods is prone to error, and precisely identifying CNVs from paired-end whole genome sequencing (WGS) data is still challenging. Here, we present a framework, CNV-JACG, for Judging the Accuracy of CNVs and Genotyping using paired-end WGS data. CNV-JACG is based on a random forest model trained on 21 distinctive features characterizing the CNV region and its breakpoints. Using the data from the 1000 Genomes Project, Genome in a Bottle Consortium, the Human Genome Structural Variation Consortium and in-house technical replicates, we show that CNV-JACG has superior sensitivity over the latest genotyping method, SV², particularly for the small CNVs (≤1 kb). We also demonstrate that CNV-JACG outperforms SV² in terms of Mendelian inconsistency in trios and concordance between technical replicates. Our study suggests that CNV-JACG would be a useful tool in assessing the accuracy of CNVs to meet the ever-growing needs for uncovering the missing heritability linked to CNVs.

Association of Variants in PLD1, 3p24.1, and 10q11.21 Regions With Hirschsprung's Disease in Han Chinese Population

Background and Aims: Hirschsprung's disease (HSCR) is a rare genetically heterogeneous congenital disorder. A recent study based on whole genome sequencing demonstrated that common variants at four novel loci, which contained two intronic variants on CASQ2 and PLD1, and intergenic variants located between SLC4A7 and EOMES at 3p24.1, and between LINC01518 and LOC283028 at 10q11.21, were associated with HSCR susceptibility. To validate these associations with HSCR susceptibility, we performed a case-control study in a Han Chinese sample set. Methods: We selected four previously identified single nucleotide polymorphisms (SNPs) for replication, along with tag SNPs to cover the four associated regions. In total, 61 SNPs were genotyped in 420 HSCR patients and 1,665 healthy controls from the Han Chinese population. Results: None of the 14 tag SNPs in the CASQ2 gene region, including the previously associated rs9428225, showed an association with HSCR. Among the 24 tag SNPs from the SLC4A7-EOMES region at 3p24.1, rs2642925 [odds ratio (OR) = 1.41, 95% confidence interval (95% CI) = 1.10-1.79; P _Additive = 0.007] and the previously associated SNP rs9851320 showed a suggestive association (OR = 1.22, 95% CI = 1.01-1.47; P _Additive = 0.042). A non-synonymous SNP, rs2287579, in PLD1 showed a suggestive association with HSCR susceptibility (OR = 1.71, 95% CI = 1.18-2.46; P _Additive = 0.004). Additionally, the previously associated PLD1 SNP rs12632766 showed a suggestive significance (OR = 1.20, 95% CI = 1.01-1.42, P _Additive = 0.038). In the LINC01518-LOC283028 region at 10q11.21, three SNPs meet the study-wide significance threshold. Rs17153309 was the most associated SNP (OR = 1.60, 95% CI = 1.34-1.90; P _Additive = 1.13 × 10^-7). The previously associated SNP rs1414027 also showed significant association (OR = 1.43, 95% CI = 1.20-1.70, P _Additive = 3.92 × 10^-5). Two associated SNPs at 10q11.21 (rs1414027 and rs624804) were expression quantitative trait loci in digestive tract tissues from GTEx databases. Conclusions: Our results confirmed that variants of the LINC01518-LOC283028 region were associated with HSCR in the Han Chinese population. Additionally, the susceptibility of SNPs in the LINC01518-LOC283028 region were associated with the expression levels of nearby genes. These results provide new insight into the pathogenesis of HSCR.

miR-618 rs2682818 C>A polymorphism decreases Hirschsprung disease risk in Chinese children

MicroRNAs (miRNAs) are endogenous non-coding small RNAs that play an important role in the development of many malignant tumors. In addition, recent studies have reported that single nucleotide polymorphisms (SNPs) located in the miRNA functional region was inextricably linked to tumor susceptibility. In the present study, we investigated the susceptibility between miR-618 rs2682818 C>A and Hirschsprung disease (HSCR) in the Southern Chinese population (1470 patients and 1473 controls). Odds ratios (ORs) and 95% confidence intervals (CIs) were used for estimating the strength of interrelation between them. We found that the CA/AA genotypes of miR-618 rs2682818 were associated with a decreased risk of HSCR when compared with the CC genotype (OR = 0.84, 95% CI = 0.72–0.99, P=0.032). Based on the stratified analysis of HSCR subtypes, the rs2682818 CA/AA genotypes were able to significantly lessen the risk of HSCR compared with CC genotype in patients with long-segment HSCR (adjusted OR = 0.70, 95% CI = 0.52–0.93, P=0.013). In conclusion, our results indicated that the miR-618 rs2682818 C>A polymorphism was associated with a reduced risk of HSCR in Chinese children, especially in patients with long-segment HSCR (L-HSCR) subtype.

Negative Association Between lncRNA HOTTIP rs3807598 C>G and Hirschsprung Disease

Background

Hirschsprung disease (HSCR) is a congenital disease that arises from defective intestinal neural system. LncRNA HOTTIP is a critical gene in various diseases, including HSCR. No epidemiological studies have explored the correlation between lncRNA HOTTIP single nucleotide polymorphisms (SNPs) and HSCR risk. We here lead as a pioneer to explore whether SNPs in lncRNA HOTTIP impact the risk of HSCR and HSCR subtypes in an unrelated Chinese population.

Methods

We used the TaqMan method to genotype rs3807598 C>G of the lncRNA HOTTIP gene using 1470 HSCR cases and 1473 healthy controls. Of them, 1441 cases and 1434 controls were successfully genotyped. We adopted odds ratios (ORs) and 95% confidence intervals (CIs) to quantify the relationship.

Results

We got an unexpected outcome that lncRNA HOTTIP SNP rs3807598 C>G could not modify the risk of HSCR (CG vs. CC: adjusted OR=0.89, 95% CI=0.74–1.07; GG vs. CC: adjusted OR=1.10, 95% CI=0.89–1.37; GG/CG vs CC: adjusted OR=0.95, 95% CI=0.80–1.13; and GG vs. CC/CG: adjusted OR=1.19, 95% CI=0.99–1.43). What’s more, risk effect of lncRNA HOTTIP rs3807598 C>G is still not obvious in stratification analysis by HSCR subtype.

Conclusion

Our studies did not provide statistical evidence of a correlation between lncRNA HOTTIP SNP rs3807598 C>G and susceptibility of HSCR in the Chinese population that is being studied. Further validation study with a larger sample size covering multi-ethnic groups is warranted.

Copy number variations in candidate genomic regions confirm genetic heterogeneity and parental bias in Hirschsprung disease

Background

Hirschsprung Disease (HSCR) is a congenital defect of the intestinal innervations characterized by complex inheritance. Many susceptibility genes including RET, the major HSCR gene, and several linked regions and associated loci have been shown to contribute to disease pathogenesis. Nonetheless, a proportion of patients still remains unexplained. Copy Number Variations (CNVs) have already been involved in HSCR, and for this reason we performed Comparative Genomic Hybridization (CGH), using a custom array with high density probes.

Results

A total of 20 HSCR candidate regions/genes was tested in 55 sporadic patients and four patients with already known chromosomal aberrations. Among 83 calls, 12 variants were experimentally validated, three of which involving the HSCR crucial genes SEMA3A/3D, NRG1, and PHOX2B. Conversely RET involvement in HSCR does not seem to rely on the presence of CNVs while, interestingly, several gains and losses did co-occur with another RET defect, thus confirming that more than one predisposing event is necessary for HSCR to develop. New loci were also shown to be involved, such as ALDH1A2, already found to play a major role in the enteric nervous system. Finally, all the inherited CNVs were of maternal origin.

Conclusions

Our results confirm a wide genetic heterogeneity in HSCR occurrence and support a role of candidate genes in expression regulation and cell signaling, thus contributing to depict further the molecular complexity of the genomic regions involved in the Enteric Nervous System development. The observed maternal transmission bias for HSCR associated CNVs supports the hypothesis that in females these variants might be more tolerated, requiring additional alterations to develop HSCR disease.

Uncovering the genetic lesions underlying the most severe form of Hirschsprung disease by whole-genome sequencing

Hirschsprung disease (HSCR) is a complex birth defect characterized by the lack of ganglion cells along a variable length of the distal intestine. A large proportion of HSCR patients remain genetically unexplained. We applied whole-genome sequencing (WGS) on 9 trios where the probands are sporadically affected with the most severe form of the disorder and harbor no coding sequence variants affecting the function of known HSCR genes. We found de novo protein-altering variants in three intolerant to change genes-CCT2, VASH1, and CYP26A1-for which a plausible link with the enteric nervous system (ENS) exists. De novo single-nucleotide and indel variants were present in introns and non-coding neighboring regions of ENS-related genes, including NRG1 and ERBB4. Joint analysis with those inherited rare variants found under recessive and/or digenic models revealed both patient-unique and shared genetic features where rare variants were found to be enriched in the extracellular matrix-receptor (ECM-receptor) pathway (p = 3.4 × 10-11). Delineation of the genetic profile of each patient might help finding common grounds that could lead to the discovery of shared molecules that could be used as drug targets for the currently ongoing cell therapy effort which aims at providing an alternative to the surgical treatment.

Association of VAMP5 and MCC genetic polymorphisms with increased risk of Hirschsprung disease susceptibility in Southern Chinese children

Hirschsprung disease (HSCR) is a genetic disorder characterized by the absence of neural crest cells in parts of the intestine. This study aims to investigate the association of vesicle-associated membrane protein 5 (VAMP5) and mutated in colorectal cancer (MCC) genetic polymorphisms and their correlated risks with HSCR. We examined the association in four polymorphisms (rs10206961, rs1254900 and rs14242 in VAMP5, rs11241200 in MCC) and HSCR susceptibility in a Southern Chinese population composed of 1473 cases and 1469 controls. Two variants in VAMP5 were replicated as associated with HSCR. Interestingly, we clarified SNPs rs10206961 and rs1254900 in VAMP5 are more essential for patients with long-segment aganglionosis (LHSCR). Relatively high expression correlation was observed between VAMP5 and MCC using data from public database showing there may exist potential genetic interactions. SNP interaction was cross-examined by logistic regression and multifactor dimensionality reduction analysis revealing that VAMP5 rs1254900 and MCC rs11241200 were interacting significantly, thereby contributing to the risk of HSCR. The results suggest that significant associations of the rs10206961 and rs14242 in VAMP5 with an increased risk of HSCR in Southern Chinese, especially in LHSCR patients. This study provided new evidence of epistatic association of VAMP5 and MCC with increased risk of HSCR.

Common genetic variants in GAL, GAP43 and NRSN1 and interaction networks confer susceptibility to Hirschsprung disease

Hirschsprung disease (HSCR) is a severe multifactorial genetic disorder. Microarray studies indicated GAL,GAP43 and NRSN1 might contribute to the altered risk in HSCR. Thus, we focused on genetic variations in GAL,GAP43 and NRSN1, and the gene‐gene interactions involved in HSCR susceptibility. We recruited a strategy combining case‐control study and MassArray system with interaction network analysis. For GAL,GAP43 and NRSN1, a total of 18 polymorphisms were assessed in 104 subjects with sporadic HSCR and 151 controls of Han Chinese origin. We found statistically significant differences between HSCR and control groups at 5 genetic variants. For each gene, the haplotypes combining all polymorphisms were the most significant. Based on SNPsyn, MDR and GeneMANIA analyses, we observed significant gene‐gene interactions among GAL,GAP43,NRSN1 and our previous identified RELN,GABRG2 and PTCH1. Our study for the first time indicates that genetic variants within GAL,GAP43 and NRSN1 and related gene‐gene interaction networks might be involved in the altered susceptibility to HSCR in the Han Chinese population, which might shed more light on HSCR pathogenesis.

Significance of neurexin and neuroligin polymorphisms in regulating risk of Hirschsprung's disease

By performing a basic case-control study among a Chinese population, the aims of this study were to explore if single nucleotide polymorphisms (SNPs) within neurexin and neuroligin were associated with susceptibility to Hirschsprung's disease (HD). Eleven SNPs within neurexin and neuroligin were selected in this basic case-control study, and this study recruited 210 children with HD and 187 healthy children. The t-test and Χ² test were used to find the difference between case and control in their clinical variables. OR and 95% CI were used to assess the association between HD susceptibility and neurexin/neuroligin polymorphisms/haplotypes. Several SNPs were significantly associated with altered risk of HD in the Chinese Han population, including rs1421589 within NRXN1, rs11795613 and rs4844285 within NLGN3, as well as rs5961397, rs7157669 and rs724373 within NLGX4X (all P<0.05). Further studies presented that the effects of rs1421589 within NRXN1, rs4844285 and rs11795613 within NLGN3, as well as rs5961397 within NLGX4X on HD phenotypes were also statistically significant (all P<0.05). Conclusively, the polymorphisms and haplotypes situated within neurexin and neuroligin were markedly associated with the onset of HD, implying that mutations of neurexin and neuroligin might serve as the treatment target for HD for the Chinese children.

Neuregulin 3 and its roles in schizophrenia risk and presentation

Neuregulins, a four-member family of epidermal growth factor-like signaling molecules, have been studied for over two decades. They were first implicated in schizophrenia in 2002 with the detection of linkage and association at the NRG1 locus followed after a few years by NRG3. However, the associations with disease have not been very consistently observed. In contrast, association of NGR3 variants with disease presentation, specifically the presence of delusions, has been more consistent. This appears to be mediated by quantitative changes in the alternative splicing of the gene, which has also been consistently observed. Additional diseases and phenotypes, psychiatric or not, have also been connected with NRG3. These results demonstrate two important aspects of behavioral genetics research. The first is that if we only consider simple risk and fail to examine the details of each patient's individual phenotype, we will miss important insights on the disease biology. This is an important aspect of the goals of precision medicine. The second is that the functional consequences of variants are often more complex than simple alterations in levels of transcription of a particular gene, including, among others, regulation of alternative splicing. To accurately model and understand the biological consequences of phenotype-associated genetic variants, we need to study the biological consequences of each specific variant. Simply studying the consequences of a null allele of the orthologous gene in a model system, runs the risk of missing the many nuances of hypomorphic and/or gain of function variants in the genome of interest.

Self-domestication in Homo sapiens: Insights from comparative genomics

This study identifies and analyzes statistically significant overlaps between selective sweep screens in anatomically modern humans and several domesticated species. The results obtained suggest that (paleo-)genomic data can be exploited to complement the fossil record and support the idea of self-domestication in Homo sapiens, a process that likely intensified as our species populated its niche. Our analysis lends support to attempts to capture the "domestication syndrome" in terms of alterations to certain signaling pathways and cell lineages, such as the neural crest.

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.

Defining the transcriptomic landscape of the developing enteric nervous system and its cellular environment

BACKGROUND

Motility and the coordination of moving food through the gastrointestinal tract rely on a complex network of neurons known as the enteric nervous system (ENS). Despite its critical function, many of the molecular mechanisms that direct the development of the ENS and the elaboration of neural network connections remain unknown. The goal of this study was to transcriptionally identify molecular pathways and candidate genes that drive specification, differentiation and the neural circuitry of specific neural progenitors, the phox2b expressing ENS cell lineage, during normal enteric nervous system development. Because ENS development is tightly linked to its environment, the transcriptional landscape of the cellular environment of the intestine was also analyzed.

RESULTS

Thousands of zebrafish intestines were manually dissected from a transgenic line expressing green fluorescent protein under the phox2b regulatory elements [Tg(phox2b:EGFP) ^w37 ]. Fluorescence-activated cell sorting was used to separate GFP-positive phox2b expressing ENS progenitor and derivatives from GFP-negative intestinal cells. RNA-seq was performed to obtain accurate, reproducible transcriptional profiles and the unbiased detection of low level transcripts. Analysis revealed genes and pathways that may function in ENS cell determination, genes that may be identifiers of different ENS subtypes, and genes that define the non-neural cellular microenvironment of the ENS. Differential expression analysis between the two cell populations revealed the expected neuronal nature of the phox2b expressing lineage including the enrichment for genes required for neurogenesis and synaptogenesis, and identified many novel genes not previously associated with ENS development. Pathway analysis pointed to a high level of G-protein coupled pathway activation, and identified novel roles for candidate pathways such as the Nogo/Reticulon axon guidance pathway in ENS development.

CONCLUSION

We report the comprehensive gene expression profiles of a lineage-specific population of enteric progenitors, their derivatives, and their microenvironment during normal enteric nervous system development. Our results confirm previously implicated genes and pathways required for ENS development, and also identify scores of novel candidate genes and pathways. Thus, our dataset suggests various potential mechanisms that drive ENS development facilitating characterization and discovery of novel therapeutic strategies to improve gastrointestinal disorders.

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.

A genome-wide association analysis of chromosomal aberrations and Hirschsprung disease

Hirschsprung disease (HSCR) is a neurocristopathy characterized by the absence of intramural ganglion cells along variable lengths of the gastrointestinal tract. Although the RET proto-oncogene is considered to be the main risk factor for HSCR, only about 30% of the HSCR cases can be explained by variations in previously known genes including RET. Recently, copy number variation (CNV) and loss of heterozygosity (LOH) have emerged as new ways to understand human genomic variation. The goal of this present study is to identify new HSCR genetic factors related to CNV in Korean patients. In the genome-wide genotyping, using Illumina's HumanOmni1-Quad BeadChip (1,140,419 markers), of 123 HSCR patients and 432 unaffected subjects (total n = 555), a total of 8,188 CNVs (1 kb ∼ 1 mb) were identified by CNVpartition. As a result, 16 CNV regions and 13 LOH regions were identified as associated with HSCR (minimum P = 0.0005). Two top CNV regions (deletions at chr6:32675155-32680480 and chr22:20733495-21607293) were successfully validated by additional real-time quantitative polymerase chain reaction analysis. In addition, 2 CNV regions (6p21.32 and 22q11.21) and 2 LOH regions (3p22.2 and 14q23.3) were discovered to be unique to the HSCR patients group. Regarding the large-scale chromosomal aberrations (>1 mb), 11 large aberrations in the HSCR patients group were identified, which suggests that they may be a risk factor for HSCR. Although further replication in a larger cohort is needed, our findings may contribute to the understanding of the etiology of HSCR.

Exome-Wide Association Study Identified New Risk Loci for Hirschsprung's Disease

Hirschsprung disease (HSCR) is a rare congenital disease caused by impaired proliferation and migration of neural crest cells. In this study, we aimed to investigate the genetic loci involved in the pathogenesis of HSCR. The exome-wide scan was performed to screen the genetic variants with minor allele frequency (MAF) < 0.05 in exonic regions. Candidate mutation type and the wild type were overexpressed to investigate the affection on cell proliferation and migration. We found that ten variants were associated with HSCR at P < 10^-4 in the single-variant analysis while ten genes were also associated with HSCR at P < 10^-4 in the optimized sequence kernel association test (SKAT-O) test analysis. Among these SNPs, the missense variants catechol-O-methyltransferase (COMT) (rs6267) and armadillo repeat gene deleted in velocardiofacial syndrome (ARVCF) (rs80068543) indicated an ectopic expression in colon tissues of HSCR patients. The Ala72Ser variant in COMT induced proliferation suppression through NOTCH signal pathway, while the ARVCF affected cell migration via the downregulating of RHOA and ROC. In conclusion, this exome array study identified the COMT and ARVCF missense coding variants as candidate loci for HSCR. The finding implies the abnormal variant of COMT and ARVCF may account for the pathogenesis of HSCR.

Identification of GLI Mutations in Patients With Hirschsprung Disease That Disrupt Enteric Nervous System Development in Mice

BACKGROUND & AIMS

Hirschsprung disease is characterized by a deficit in enteric neurons, which are derived from neural crest cells (NCCs). Aberrant hedgehog signaling disrupts NCC differentiation and might cause Hirschsprung disease. We performed genetic analyses to determine whether hedgehog signaling is involved in pathogenesis.

METHODS

We performed deep-target sequencing of DNA from 20 patients with Hirschsprung disease (16 men, 4 women), and 20 individuals without (controls), and searched for mutation(s) in GLI1, GLI2, GLI3, SUFU, and SOX10. Biological effects of GLI mutations were tested in luciferase reporter assays using HeLa or neuroblastoma cell lines. Development of the enteric nervous system was studied in Sufu(f/f), Gli3(Δ699), Wnt1-Cre, and Sox10(NGFP) mice using immunohistochemical and whole-mount staining procedures to quantify enteric neurons and glia and analyze axon fasciculation, respectively. NCC migration was studied using time-lapse imaging.

RESULTS

We identified 3 mutations in GLI in 5 patients with Hirschsprung disease but no controls; all lead to increased transcription of SOX10 in cell lines. SUFU, GLI, and SOX10 form a regulatory loop that controls the neuronal vs glial lineages and migration of NCCs. Sufu mutants mice had high Gli activity, due to loss of Sufu, disrupting the regulatory loop and migration of enteric NCCs, leading to defective axonal fasciculation, delayed gut colonization, or intestinal hypoganglionosis. The ratio of enteric neurons to glia correlated inversely with Gli activity.

CONCLUSIONS

We identified mutations that increase GLI activity in patients with Hirschsprung disease. Disruption of the SUFU-GLI-SOX10 regulatory loop disrupts migration of NCCs and development of the enteric nervous system in mice.

A collagen VI-dependent pathogenic mechanism for Hirschsprung's disease

Hirschsprung's disease (HSCR) is a severe congenital anomaly of the enteric nervous system (ENS) characterized by functional intestinal obstruction due to a lack of intrinsic innervation in the distal bowel. Distal innervation deficiency results from incomplete colonization of the bowel by enteric neural crest cells (eNCCs), the ENS precursors. Here, we report the generation of a mouse model for HSCR--named Holstein--that contains an untargeted transgenic insertion upstream of the collagen-6α4 (Col6a4) gene. This insertion induces eNCC-specific upregulation of Col6a4 expression that increases total collagen VI protein levels in the extracellular matrix (ECM) surrounding both the developing and the postnatal ENS. Increased collagen VI levels during development mainly result in slower migration of eNCCs. This appears to be due to the fact that collagen VI is a poor substratum for supporting eNCC migration and can even interfere with the migration-promoting effects of fibronectin. Importantly, for a majority of patients in a HSCR cohort, the myenteric ganglia from the ganglionated region are also specifically surrounded by abundant collagen VI microfibrils, an outcome accentuated by Down syndrome. Collectively, our data thus unveil a clinically relevant pathogenic mechanism for HSCR that involves cell-autonomous changes in ECM composition surrounding eNCCs. Moreover, as COL6A1 and COL6A2 are on human Chr.21q, this mechanism is highly relevant to the predisposition of patients with Down syndrome to HSCR.

Exome sequencing reveals a high genetic heterogeneity on familial Hirschsprung disease

Hirschsprung disease (HSCR; OMIM 142623) is a developmental disorder characterized by aganglionosis along variable lengths of the distal gastrointestinal tract, which results in intestinal obstruction. Interactions among known HSCR genes and/or unknown disease susceptibility loci lead to variable severity of phenotype. Neither linkage nor genome-wide association studies have efficiently contributed to completely dissect the genetic pathways underlying this complex genetic disorder. We have performed whole exome sequencing of 16 HSCR patients from 8 unrelated families with SOLID platform. Variants shared by affected relatives were validated by Sanger sequencing. We searched for genes recurrently mutated across families. Only variations in the FAT3 gene were significantly enriched in five families. Within-family analysis identified compound heterozygotes for AHNAK and several genes (N = 23) with heterozygous variants that co-segregated with the phenotype. Network and pathway analyses facilitated the discovery of polygenic inheritance involving FAT3, HSCR known genes and their gene partners. Altogether, our approach has facilitated the detection of more than one damaging variant in biologically plausible genes that could jointly contribute to the phenotype. Our data may contribute to the understanding of the complex interactions that occur during enteric nervous system development and the etiopathology of familial HSCR.

Genotyping analysis of 3 RET polymorphisms demonstrates low somatic mutation rate in Chinese Hirschsprung disease patients

BACKGROUND

Genetic mosaicism has been reported for both coding and non-coding sequences in the RET gene in Hirschsprung disease (HSCR) patients. This study aimed to investigate somatic mutation rate in Chinese population by comparing both homozygous genotype percentage and risk allele frequency of 3 RET single nucleotide polymorphisms (SNPs) among blood and colon samples.

METHODS

DNA was extracted from 59 HSCR blood samples, 59 control blood samples and 76 fresh frozen colon tissue samples (grouped into ganglionic, transitional and aganglionic level). Genotype status of rs2435357 and rs2506030 was examined by competitive allele specific hydrolysis probes (Taqman) PCR technology, and rs2506004 was examined by Sanger sequencing. Homozygous genotype percentage and risk allele frequency were calculated for each type of sample and compared by chi-square test. P<0.05 was regarded as being statistically significant.

RESULTS

Colon tissue DNA samples showed similar frequency of SNPs as that of the blood DNA samples in HSCR patients, both of which are significantly higher than the control blood group (rs2435357 TT genotype: 71.2%, 74.7% versus 22.0% in HSCR blood, HSCR colon and control blood DNA respectively, P=0.000; rs2506004 AA genotype: 72.4%, 83.1% versus 25.5%, P=0.000; rs2506030 GG genotype: 79.7%, 77.2% versus 54.2%, P=0.000 and 0.004). With respect to DNA extracted from ganglionic, transitional and aganglionic levels, no statistically significant difference was demonstrated in those 3 regions (rs2435357: P=0.897; rs2506004: P=0.740; rs2506030: P=0.901).

CONCLUSION

Our data does not support the notion that high frequency of somatic changes as an underlying etiology of Chinese HSCR population.

A common polymorphism in pre-miR-146a underlies Hirschsprung disease risk in Han Chinese

BACKGROUND

Hirschsprung disease (HSCR) is a rare multigenic congenital disorder characterized by the absence of the enteric ganglia. To date, single nucleotide polymorphisms (SNPs) in pre-miRNAs have been confirmed related with some diseases. Thus, we hypothesized that pre-miRNA polymorphisms might contribute to HSCR susceptibility. We investigated whether rs2910164 and rs11614913 of pre-miR-146a and pre-miR-196a2, are associated with HSCR.

METHODS

Polymorphisms were genotyped using the Taqman method. Real-time PCR was used for detecting the expression level of miR-146a and its target gene ROBO1 in CC and GG genotypes.

RESULTS

Significant differences were found in the genotype distribution of rs2910164 and rs11614913 polymorphism between HSCR cases and controls (p = 0.023 and 0.041, respectively). Furthermore, G allele of rs2910164 might increase the risk of HSCR (OR, 1.54; 95% CI, 1.06-2.23). Moreover, the expression level of miR-146a for homozygote GG was also higher than homozygote CC (p = 0.0193). In contrast, the expression level of its target gene ROBO1 predicted in bioinformatics for homozygote GG was much lower than homozygote CC (p = 0.0096).

CONCLUSIONS

Our results showed that the polymorphism rs2910164 in pre-miR-146a might alter the production of mature miR-146a and then down-regulate the target gene ROBO1, which plays an important role in pathogenesis of HSCR.

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.

Involvement of DNMT3B in the pathogenesis of Hirschsprung disease and its possible role as a regulator of neurogenesis in the human enteric nervous system

PURPOSE

Hirschsprung disease (OMIM 142623) is a neurocristopathy attributed to a failure of cell proliferation or migration and/or failure of the enteric precursors along the gut to differentiate during embryonic development. Although some genes involved in this pathology are well characterized, many aspects remain poorly understood. In this study, we aimed to identify novel genes implicated in the pathogenesis of Hirschsprung disease.

METHODS

We compared the expression patterns of genes involved in human stem cell pluripotency between enteric precursors from controls and Hirschsprung disease patients. We further evaluated the role of DNMT3B in the context of Hirschsprung disease by inmunocytochemistry, global DNA methylation assays, and mutational screening.

RESULTS

Seven differentially expressed genes were identified. We focused on DNMT3B, which encodes a DNA methyltransferase that performs de novo DNA methylation during embryonic development. DNMT3B mutational analysis in our Hirschsprung disease series revealed the presence of potentially pathogenic mutations (p.Gly25Arg, p.Arg190Cys, and p.Gly198Trp).

CONCLUSION

DNMT3B may be regulating enteric nervous system development through DNA methylation in the neural crest cells, suggesting that aberrant methylation patterns could have a relevant role in Hirschsprung disease. Moreover, the synergistic effect of mutations in both DNMT3B and other Hirschsprung disease-related genes may be contributing to a more severe phenotype in our Hirschsprung disease patients.

A genomic copy number variant analysis implicates the MBD5 and HNRNPU genes in Chinese children with infantile spasms and expands the clinical spectrum of 2q23.1 deletion

Background

Infantile spasms (IS) is a specific type of epileptic encephalopathy associated with severe developmental disabilities. Genetic factors are strongly implicated in IS, however, the exact genetic defects remain unknown in the majority of cases. Rare mutations in a single gene or in copy number variants (CNVs) have been implicated in IS of children in Western countries. The objective of this study was to dissect the role of copy number variations in Chinese children with infantile spasms.

Methods

We used the Agilent Human Genome CGH microarray 180 K for genome-wide detection of CNVs. Real-time qPCR was used to validate the CNVs. We performed genomic and medical annotations for individual CNVs to determine the pathogenicity of CNVs related to IS.

Results

We report herein the first genome-wide CNV analysis in children with IS, detecting a total of 14 CNVs in a cohort of 47 Chinese children with IS. Four CNVs (4/47 = 8.5%) (1q21.1 gain; 1q44, 2q31.1, and 17p13 loss) are considered to be pathogenic. The CNV loss at 17p13.3 contains PAFAH1B1 (LIS1), a causative gene for lissencephaly. Although the CNVs at 1q21.1, 1q44, and 2q23.1 have been previously implicated in a wide spectrum of clinical features including autism spectrum disorders (ASD) and generalized seizure, our study is the first report identifying them in individuals with a primary diagnosis of IS. The CNV loss in the 1q44 region contains HNRNPU, a strong candidate gene recently suggested in IS by the whole exome sequencing of children with IS. The CNV loss at 2q23.1 includes MBD5, a methyl-DNA binding protein that is a causative gene of ASD and a candidate gene for epileptic encephalopathy. We also report a distinct clinical presentation of IS, microcephaly, intellectual disability, and absent hallux in a case with the 2q23.1 deletion.

Conclusion

Our findings strongly support the role of CNVs in infantile spasms and expand the clinical spectrum associate with 2q23.1 deletion. In particular, our study implicates the HNRNPU and MBD5 genes in Chinese children with IS. Our study also supports that the molecular mechanisms of infantile spasms appear conserved among different ethnic backgrounds.

A novel stop mutation in the EDNRB gene in a family with Hirschsprung's disease associated with multiple sclerosis

PURPOSE

We identified a girl with Hirschsprung's disease (HSCR) whose mother and grandmother had HSCR associated with multiple sclerosis (MS). The aim of this study was to outline mutations in HSCR-related genes and MS susceptibility alleles in these three individuals.

METHODS

The phenotypes were reviewed based on medical records. The three subjects had rectosigmoid HSCR verified with histopathology. The mother and grandmother fulfilled the McDonald criteria for MS. DNA was isolated from EDTA-preserved blood according to standard procedures. Exome sequencing aiming mainly at analyzing HSCR associated genes as well as Sanger sequencing for confirmation was performed.

RESULTS

All affected individuals carry a novel heterozygous nonsense mutation in the EDNRB gene (c.C397T,p.R133X,refNM_000115), changing an arginine at position 133 into a premature stop codon. None of the subjects were homozygous for the HLA risk alleles for MS.

CONCLUSION

We report a novel non-sense EDNRB gene mutation in a girl with HSCR and her mother and grandmother with HSCR and MS. We propose that this EDNRB gene mutation plays a role in the etiology of HSCR and also makes the subjects susceptible to MS.

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.

Common genetic variations in Patched1 (PTCH1) gene and risk of hirschsprung disease in the Han Chinese population

Hirschsprung disease (HSCR) is the most frequent genetic cause of congenital intestinal obstruction with an incidence of 1:5000 live births. In a pathway-based epistasis analysis of data generated by genome-wide association study on HSCR, specific genotype of Patched 1 (PTCH1) has been linked to an increased risk for HSCR. The aim of the present study is to examine the contribution of genetic variants in PTCH1 to the susceptibility to HSCR in Han Chinese. Accordingly, we assessed 8 single nucleotide polymorphisms (SNPs) within PTCH1 gene in 104 subjects with sporadic HSCR and 151 normal controls of Han Chinese origin by the Sequenom MassArray technology (iPLEX GOLD). Two of the eight genetic markers were found to be significantly associated with Hirschsprung disease (rs357565, allele P = 0.005; rs2236405, allele P = 0.002, genotype P = 0.003). Both the C allele of rs357565 and the A allele of rs2236405 served as risk factors for HSCR. During haplotype analysis, one seven-SNP-based haplotype was the most significant, giving a global P = 0.0036. Our results firstly suggest common variations of PTCH1 may be involved in the altered risk for HSCR in the Han Chinese population, providing potential molecular markers for early diagnosis of Hirschsprung disease.

GStream: improving SNP and CNV coverage on genome-wide association studies

We present GStream, a method that combines genome-wide SNP and CNV genotyping in the Illumina microarray platform with unprecedented accuracy. This new method outperforms previous well-established SNP genotyping software. More importantly, the CNV calling algorithm of GStream dramatically improves the results obtained by previous state-of-the-art methods and yields an accuracy that is close to that obtained by purely CNV-oriented technologies like Comparative Genomic Hybridization (CGH). We demonstrate the superior performance of GStream using microarray data generated from HapMap samples. Using the reference CNV calls generated by the 1000 Genomes Project (1KGP) and well-known studies on whole genome CNV characterization based either on CGH or genotyping microarray technologies, we show that GStream can increase the number of reliably detected variants up to 25% compared to previously developed methods. Furthermore, the increased genome coverage provided by GStream allows the discovery of CNVs in close linkage disequilibrium with SNPs, previously associated with disease risk in published Genome-Wide Association Studies (GWAS). These results could provide important insights into the biological mechanism underlying the detected disease risk association. With GStream, large-scale GWAS will not only benefit from the combined genotyping of SNPs and CNVs at an unprecedented accuracy, but will also take advantage of the computational efficiency of the method.

Enteric nervous system development: migration, differentiation, and disease

The enteric nervous system (ENS) provides the intrinsic innervation of the bowel and is the most neurochemically diverse branch of the peripheral nervous system, consisting of two layers of ganglia and fibers encircling the gastrointestinal tract. The ENS is vital for life and is capable of autonomous regulation of motility and secretion. Developmental studies in model organisms and genetic studies of the most common congenital disease of the ENS, Hirschsprung disease, have provided a detailed understanding of ENS development. The ENS originates in the neural crest, mostly from the vagal levels of the neuraxis, which invades, proliferates, and migrates within the intestinal wall until the entire bowel is colonized with enteric neural crest-derived cells (ENCDCs). After initial migration, the ENS develops further by responding to guidance factors and morphogens that pattern the bowel concentrically, differentiating into glia and neuronal subtypes and wiring together to form a functional nervous system. Molecules controlling this process, including glial cell line-derived neurotrophic factor and its receptor RET, endothelin (ET)-3 and its receptor endothelin receptor type B, and transcription factors such as SOX10 and PHOX2B, are required for ENS development in humans. Important areas of active investigation include mechanisms that guide ENCDC migration, the role and signals downstream of endothelin receptor type B, and control of differentiation, neurochemical coding, and axonal targeting. Recent work also focuses on disease treatment by exploring the natural role of ENS stem cells and investigating potential therapeutic uses. Disease prevention may also be possible by modifying the fetal microenvironment to reduce the penetrance of Hirschsprung disease-causing mutations.

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.

Exome sequencing identified NRG3 as a novel susceptible gene of Hirschsprung's disease in a 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.