Aberrant high expression of NRG1 gene in Hirschsprung disease

The roles of non-coding RNAs in Hirschsprung's disease

Hirschsprung's disease (HSCR) is a congenital disorder characterized by the absence of ganglion cells in the colon, leading to various intestinal complications. The etiology of HSCR stems from complex genetic and environmental interactions, of which the intricate roles of non-coding RNAs (ncRNAs) are a key area of research. However, the roles of ncRNAs in the pathogenesis of HSCR have not been fully elucidated. In order to understand the variety of symptoms caused by HSCR and develop new therapeutic approaches, it is essential to understand the underlying biological genetic basis of HSCR. This review presents a comprehensive overview of the current understanding regarding the involvement of ncRNAs in HSCR, including microRNAs (miRNAs), long noncoding RNAs (lncRNAs), and circular RNAs (circRNAs). Additionally, it provides a summary of the molecular mechanisms through which ncRNAs regulate the expression of genes related to the proliferation, migration, and differentiation of intestinal neural crest cells, thereby contributing to the advancement of HSCR research.

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.

The Endocrine Disruptor Bisphenol A (BPA) Affects the Enteric Neurons Immunoreactive to Neuregulin 1 (NRG1) in the Enteric Nervous System of the Porcine Large Intestine

The enteric nervous system (ENS), located in the wall of the gastrointestinal (GI) tract, is characterized by complex organization and a high degree of neurochemical diversity of neurons. One of the less known active neuronal substances found in the enteric neurons is neuregulin 1 (NRG1), a factor known to be involved in the assurance of normal development of the nervous system. During the study, made up using the double immunofluorescence technique, the presence of NRG1 in the ENS of the selected segment of porcine large intestine (caecum, ascending and descending colon) was observed in physiological conditions, as well as under the impact of low and high doses of bisphenol A (BPA) which is commonly used in the production of plastics. In control animals in all types of the enteric plexuses, the percentage of NRG1-positive neurons oscillated around 20% of all neurons. The administration of BPA caused an increase in the number of NRG1-positive neurons in all types of the enteric plexuses and in all segments of the large intestine studied. The most visible changes were noted in the inner submucous plexus of the ascending colon, where in animals treated with high doses of BPA, the percentage of NRG1-positive neurons amounted to above 45% of all neuronal cells. The mechanisms of observed changes are not entirely clear, but probably result from neurotoxic, neurodegenerative and/or proinflammatory activity of BPA and are protective and adaptive in nature.

Aberrant Expressions and Variant Screening of SEMA3D in Indonesian Hirschsprung Patients

Background: The semaphorin 3D (SEMA3D) gene has been implicated in the pathogenesis of Hirschsprung disease (HSCR), a complex genetic disorder characterized by the loss of ganglion cells in varying lengths of gastrointestinal tract. We wished to investigate the role of SEMA3D variants, both rare and common variants, as well as its mRNA expression in Indonesian HSCR patients. Methods: Sanger sequencing was performed in 54 HSCR patients to find a pathogenic variant in SEMA3D. Next, we determined SEMA3D expression in 18 HSCR patients and 13 anorectal malformation colons as controls by quantitative real-time polymerase chain reaction (qPCR). Results: No rare variant was found in the SEMA3D gene, except one common variant in exon 17, p.Lys701Gln (rs7800072). The risk allele (C) frequency at rs7800072 among HSCR patients (23%) was similar to those reported for the 1,000 Genomes (27%) and ExAC (28%) East Asian ancestry controls (p = 0.49 and 0.41, respectively). A significant difference in SEMA3D expression was observed between groups (p = 0.04). Furthermore, qPCR revealed that SEMA3D expression was strongly up-regulated (5.5-fold) in the ganglionic colon of HSCR patients compared to control colon (ΔC_T 10.8 ± 2.1 vs. 13.3 ± 3.9; p = 0.025). Conclusions: We report the first study of aberrant SEMA3D expressions in HSCR patients and suggest further understanding into the contribution of aberrant SEMA3D expression in the development of HSCR. In addition, this study is the first comprehensive analysis of SEMA3D variants in the Asian ancestry.

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.

Aberrant expression of LncRNA-MIR31HG regulates cell migration and proliferation by affecting miR-31 and miR-31* in Hirschsprung's disease

Hirschsprung's disease (HSCR) is a birth defect that causes a failure of the enteric nervous system to cover the distal gut during early embryonic development. Evidence shows that long non-coding RNAs (lncRNA) play important roles in HSCR. The MIR31 host gene (MIR31HG), also known as Loc554202, is a long non-coding RNA (lncRNA), which acts as the host gene of (microRNA) miR-31 and miR-31*. There have been no studies regarding its function in early developmental defects during pregnancy, and its downstream genetic receptors. We report that downregulation of MIR31HG inhibited migration and proliferation in 293T and SH-SY5Y cell lines, by suppressing miR-31 and miR-31*. Moreover, the downregulation of miR-31 and miR-31* enhanced inter-α-trypsin inhibitor heavy chain 5 (ITIH5) and the phosphatidylinositol-4, 5-bisphosphate 3-kinase catalytic gamma subunit (PIK3CG), respectively with reductions of cell migration and proliferation in 293T and SH-SY5Y cell lines. In addition, synergistic actions were observed between miR-31 and miR-31* in cell migration and proliferation. Our results demonstrated that the MIR31HG-miR-31/31*-ITIH5/PIK3CG pathway plays a role in the pathogenesis of HSCR.

Downregulation of lncRNA MEG3 and miR-770-5p inhibit cell migration and proliferation in Hirschsprung's disease

The long noncoding RNA (lncRNA) MEG3 is involved in various biological processes including cell migration and cell proliferation. In present study, it was found that MEG3 and the intronic miR-770-5p were decreased in samples from HSCR patients. Besides, knockdown of MEG3 and miR-770-5p suppressed cell migration and proliferation, while cell cycle and apoptosis were not affected in human 293T and SH-SY5Y cells. SRGAP1 mRNA and protein upregulation was inversely correlated with miR-770-5p expression in tissue samples and cell lines, which was confirmed to be a target gene of miR-770-5p by dual-luciferase reporter assay. Moreover, silencing of SRGAP1 rescued the inhibition of cell migration and proliferation induced by MEG3 siRNA and miR-770-5p inhibition. The present study elucidates a novel mechanism of the development of HSCR and shows that the MEG3/miR-770-5p/SRGAP1 pathway plays a vital role in the pathogenesis of HSCR.

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.

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.

Expression and function of Neuregulin 1 and its signaling system ERBB2/3 in the enteric nervous system

Neuregulin 1 (NRG1) is suggested to promote the survival and maintenance of the enteric nervous system (ENS). As deficiency in its corresponding receptor signaling complex ERBB2/ERBB3 leads to postnatal colonic hypo/aganglionosis we assessed the distributional and expressional pattern of the NRG1-ERBB2/ERBB3 system in the human colon and explored the neurotrophic capacity of NRG1 on cultured enteric neurons. Site-specific mRNA expression of the NRG1-ERBB2/3 system was determined in microdissected samples harvested from enteric musculature and ganglia. Localization of NRG1, ERBB2 and ERBB3 was determined by dual-label-immunohistochemistry using pan-neuronal and pan-glial markers. Morphometric analysis was performed on NRG1-stimulated rat enteric nerve cultures to evaluate neurotrophic effects. mRNA expression of the NRG1-ERBB2/3 system was determined by qPCR. Co-localization of NRG1 with neuronal or synaptic markers was analyzed in enteric nerve cultures stimulated with glial cell line-derived neurotrophic factor (GDNF). The NRG1 system was expressed in both neurons and glial cells of enteric ganglia and in nerve fibers. NRG1 significantly enhanced growth parameters in enteric nerve cell cultures and ErB3 mRNA expression was down-regulated upon NRG1 stimulation. GDNF negatively regulates ErbB2 and ErbB3 mRNA expression. The NRG1-ERBB2/3 system is physiologically present in the human ENS and NRG1 acts as a neurotrophic factor for the ENS. The down-regulation of ErbB3/ErbB2 in GDNF stimulated nerve cell cultures points to an interaction of both neurotrophic factors. Thus, the data may provide a basis to assess disturbed signaling components of the NRG1 system in enteric neuropathies.

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.

MiR-195 affects cell migration and cell proliferation by down-regulating DIEXF in Hirschsprung's disease

Background

Hirschsprung’s disease (HSCR) is the most common congenital gut motility disorder. We aimed to investigate the roles of miR-195 in the pathogenesis of HSCR.

Methods

In this study, we measured the expression levels of miRNA, mRNA, and protein in colon tissues from 78 patients with HSCR and 66 controls without HSCR. Transwell, Cell Counting Kit-8 (CCK-8) and flow cytometry assay were employed to detect the function role of miR-195 in vitro.

Results

Our results showed that expression levels of miR-195 from patients with HSCR were significantly higher than control group; along with aberrant lower expression levels of digestive-organ expansion factor (DIEXF) were tested. Increased level of miR-195 could suppress the level of DIEXF in cell, which induced the impairment of cell migration and proliferation.

Conclusions

Aberrant expression of miR-195 may involved in the pathogenesis of HSCR by down-regulated the level of DIEXF.