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.

Introduction

Hirschsprung disease (HSCR: MIM# 142623) is a complex genetic disorder, characterized by the lack of ganglion cells in the bowel, resulting in a functional obstruction during infancy (1). HSCR is categorized into the following types: short-segment, long-segment, and total colonic aganglionosis (1, 2).

HSCR incidence differs among ethnic groups, with 1.5, 2.1, and 2.8 cases per 10,000 live births in European, African, and Asian ancestry cases, respectively (1, 2). There are at least 17 genes responsible for the development of HSCR, with most of them being members of the RET and EDNRB signaling pathways (1, 2). Two genetic risk factors are the RET rs2435357 and rs2506030 variants (3, 4). Our recent studies showed that the RET rs2435357 and rs2506030 risk alleles have higher frequency in Indonesian ancestry cases as compared with European ancestry cases (5, 6), which might relate to the higher incidence of HSCR in Indonesia (3.1 cases per 10,000 live births) than other populations (7).

The third signaling pathway of HSCR pathogenesis includes class 3 semaphorins (SEMA3s), involving SEMA3D (4, 8, 9) SEMA3D has been implicated in the development of HSCR and contributes to risk through both common and rare variants in European ancestries (4, 8, 9), as evidenced by (1) the detection of both common and rare SEMA3D variants in HSCR patients; (2) the expression of SEMA3D in the human, mouse, and zebrafish intestines and, particularly, the enteric nervous system (ENS); and (3) the joint effect of Ret and Sema3d loss of function in an aganglionosis animal model. However, our recent study showed that the effect of SEMA3 rs11766001 common variant on HSCR depends on the ethnic background (10). In addition, the allele frequencies of common variants might differ among Asians, since the North Asians, Han Chinese, Japanese, and Southeast Asians can be distinguished based on their Y chromosome variants (11). Moreover, alterations in the expression of specific genes have been implicated in the development of HSCR (12–15). Therefore, we wished to investigate the role of SEMA3D variants, both rare and common variants, as well as its mRNA expression in Indonesian HSCR patients.

Materials and Methods

Patients for SEMA3D Variant Screening

We identified 54 HSCR patients: 38 males and 16 females (Table 1). We diagnosed HSCR in these patients in Dr. Sardjito Hospital, Yogyakarta, Indonesia, after evaluating clinical findings, contrast enema, and histopathology. For histopathological findings, we used hematoxylin-eosin staining and S100 immunohistochemistry (5–7, 10, 15, 16).

TABLE 1

Table 1. Clinical features of the HSCR patients for SEMA3D sequencing analysis.

All parents signed a written informed consent form before participating in this study. The Institutional Review Board of the Faculty of Medicine, Public Health, and Nursing, Universitas Gadjah Mada/Dr. Sardjito Hospital gave approval for this study (KE/FK/1356/EC/2015). All experiments were performed in accordance with relevant guidelines and regulations.

Polymerase Chain Reaction (PCR) and DNA Sequencing

A QIAamp DNA Extraction Kit (QIAGEN, Hilden, Germany) was used to extract genomic DNA from whole blood from each individual, according to the manufacturer's instructions. We stored the extracted DNA samples at −20°C until analysis. PCR was conducted using a Swift Maxi thermal cycler (Esco Micro Pte. Ltd., Singapore), followed by Sanger sequencing analysis to identify sequence variants in all 17 exons of the SEMA3D gene in HSCR patients using BigDye Terminator V3.1 Cycle Sequencing Kits (Applied Biosystems, Foster City, CA) and a 3730xl Genetic Analyzer (Applied Biosystems), with DNA Sequencing Analysis Software (Applied Biosystems) 0.1 (7). The primer sequences for SEMA3D rare variant analysis were chosen based on a previous study (4).

DNA Genotyping

DNA genotyping was performed using Sanger sequencing analysis. The SEMA3D rs7800072:A>C (chr7: g. 84,628,989A>C) variant was identified during the Sanger sequencing analysis to find a rare variant in Indonesian HSCR patients. The risk allele (C) was determined according to the 1,000 Genomes Project and ExAC population databases (17, 18).

RNA Extraction and Quantitative Real-Time PCR (qPCR)

The ganglionic and aganglionic intestinal specimens were collected at pull-through surgeries from 18 HSCR patients, while control intestinal specimens were obtained at colostomy closure from 13 anorectal malformation (ARM) patients (Table 2).

TABLE 2

Table 2. Clinical characteristics of the HSCR patients for SEMA3D expression study.

Total RNA was isolated from colonic specimens using the total RNA Mini Kit (Tissue) (Geneaid Biotech Ltd., New Taipei City, Taiwan). RNA was quantified by a NanoDrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA). The OD260/280 ratios ranged from 1.8 to 2.0, indicating high RNA purity.

SEMA3D gene expression was quantified using the BioRad CFX Real-Time PCR System (California, USA) and the SensiFAST^TM SYBR^® No-ROX One-Step Kit (Bioline, Meridian Bioscience, Memphis, USA) using the 5′-CAACGCAGCCTGATAAACAA-3′ (forward) and 5′-TCTTTCATCTCTTGTGGGGAGT-3′ (reverse) (19) The primers were designed to bridge SEMA3D exon 8 and 9 junctions (19). Glyceraldehyde-3-phosphate dehydrogenase (GAPDH), a housekeeping gene, was used as an endogenous control. The GAPDH primers were 5′-GCACCGTCAAGGCTGAGAAC-3′ (forward) and 5′-TGGTGAAGACGCCAGTGGA-3′ (reverse). qPCR reactions contained SensiFAST^TM SYBR^® No-ROX One-Step mix (2×) 10 μL, RiboSafe RNase Inhibitor 0.4 μL, reverse transcriptase 0.2 μL, forward primer (10 μM) 0.8 μL, reverse primer (10 μM) 0.8 μL, and total RNA 50 ng, with final volume of 20 μL. qPCR was performed for 10 minutes (min) at 45°C for reverse transcription process, followed by 2 min at 95°C and 39 cycles for 5 s at 95°C, 10 s at 58°C and, 5 s at 72°C, and 1 cycle for 5 s at 65°C, according to the manufacturer's instructions. We performed the gel electrophoresis for the qPCR of SEMA3D and GAPDH (Supplemental Figure 1).

The Livak method was utilized to determine the SEMA3D mRNA expression level (20). This method is designed to calculate a relative gene expression and referred to as the ΔC_T method. The (log) expression is proportional to the negative C_T value (the lower the C_T, the higher the expression) (20).

Statistical Analysis

SEMA3D expression was described as a mean value ± SD. The Kolmogorov–Smirnov test was used to determine the data distribution, and a one-way ANOVA was utilized to assess statistical differences between groups. A chi-square test was used to establish p-value for the case-control association analysis for SEMA3D rs7800072 variant. A p < 0.05 was considered significant.

Results

Most of our HSCR patients were male (70%) and short-segment aganglionosis type (98%), with the mean age at diagnosis being 34.6 ± 44.5 months (Table 1). Among the 54 HSCR patients, 49, 3, and 2 children underwent a definitive surgery, a colostomy, and a full-thickness rectal biopsy waiting a pull-through procedure, respectively, with the most common definitive surgery conducted being transanal endorectal pull-through (43%) (Table 1); and our controls were six males and seven females, with a mean age during stoma closure of 47 ± 45.1 months.

We could not identify any rare variant in all 17 exons of SEMA3D gene in 54 HSCR patients, but encountered one common variant in exon 17: p.Lys701Gln (rs7800072) (Figure 1). The genotype frequencies for rs7800072 variant among HSCR patients were as follows: AA (32/54, 60%), AC (19/54, 35%), and CC (3/54, 5%). Subsequently, we compared the risk allele (C) frequency of rs7800072 in 54 Indonesian HSCR cases and the 1000 Genomes and ExAC East Asian ancestry controls (Table 3). The risk allele (C) frequency was similar in HSCR cases (23%) and the 1000 Genomes (27%) and ExAC (28%) East Asian ancestry controls, with p-values of 0.49 and 0.41, respectively (Table 3). According to the conservation score prediction using PhyloP, the p.Lys701Gln variant did not reach a deleterious threshold of 0.84, while the predicted deleterious effect (Condel) showed the p.Lys701Gln variant as being neutral (4). Furthermore, the SIFT and PolyPhen-2 analysis of p.Lys701Gln showed the variant as being tolerated and benign, respectively (18).

FIGURE 1

Figure 1. Sanger sequencing of exon 17 SEMA3D gene in a HSCR patient. Arrow indicates a common variant, p.Lys701Gln (rs7800072).

TABLE 3

Table 3. SEMA3D common variant frequency in Indonesian HSCR and population databases (17, 21).

Next, we compared SEMA3D expression in 18 aganglionic and ganglionic colons of HSCR patients and 13 control colons. The distribution of SEMA3D expression data were normal for ganglionic, aganglionic, and control colons (p = 0.2, 0.17, and 0.2, respectively). A significant different of SEMA3D expression was observed between groups (p = 0.04). Furthermore, qPCR showed that the SEMA3D expression was strongly up-regulated (5.5-fold) (Figure 2) 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) (Table 4, Figure 3), while the SEMA3D expression was not significantly different between the aganglionic colon of HSCR patients and the control colon (ΔC_T 13.1 ± 3.0 vs. 13.3 ± 3.9; p = 0.89) (Table 4, Figure 3).

FIGURE 2

Figure 2. The SEMA3D expression was strongly up-regulated (5.5-fold) in the ganglionic colon of HSCR patients compared to control colon (p = 0.025), while the SEMA3D expression was not significantly different between the aganglionic colon of HSCR patients and the control colon (p = 0.89). *p < 0.05.

TABLE 4

Table 4. SEMA3D expressions in colon of HSCR patient and control.

FIGURE 3

Figure 3. Box-plot graph of ΔC_T value of the SEMA3D expressions in HSCR ganglionic colon (ΔC_T 10.8 ± 2.1), HSCR aganglionic colon (ΔC_T 13.4 ± 2.9), and control colon (ΔC_T 13.3 ± 3.9). Box-plot graph of ΔC_T value reveals the median values as lines across the box. Lower and upper boxes are representing the 25th percentile to the 75th percentile, while whiskers indicate the maximum and minimum values. *p < 0.05.

Discussion

In this study, we have performed an in-depth genetic and gene expression study of the SEMA3D in Indonesian HSCR patients. We did not detect any rare variants in this gene, although previous studies in European ancestry cases have identified rare coding variants in SEMA3D associated with HSCR. Our results indicate that the association of such variants to the disease may be restricted to specific ethnic groups. It may still be the case that other semaphorin 3 genes could play a role in HSCR pathogenesis (4, 9). A previous study showed the existence of other genetic factors conferring risk to HSCR in specific ethnic populations (3). For example, there were two different RET haplotypes involving the enhancer mutation that were over-transmitted to the HSCR offspring in the European sample, while in the Chinese sample, only one of those haplotypes was present (3). It might be speculated that the enhancer mutation arose in one haplotype, which after the Asian–Caucasian split, rearranged to also give the other haplotype, but exclusively in the Caucasian (22). Another example is the NRG1 rs7835688 genetic marker, which has been originally discovered in Chinese HSCR patients (23) and has been shown in other Asian ancestry cases (5, 24), but is rare and shows no effect in the Caucasian population (8, 22). Makhmudi et al. (25) also demonstrated that the MTHFR c.677C>T is a genetic risk factor for Indonesian gastroschisis, but not seen in Caucasians (26, 27). A recent study showed that loss of Sema3d in null homozygotes mice had no impact on the intestinal transcriptome (28). The authors were unable to find evidence for Ret and Sema3d interaction affecting survival, presence of myenteric plexus, or intestine transcriptome (28). Furthermore, Tang et al. revealed that the effects of RET and NRG1 variants are universal across Caucasian and Asian ancestries, but the impact of SEMA3 variant was restricted to Caucasian ancestries (29). It should be noted that our screening method did not cover the promoter and enrich regions of CpG of SEMA3D.

In addition, the observed SEMA3D rs7800072 variant frequency in Indonesian HSCR patients is similar to those reported for the 1,000 Genomes Project and ExAC East Asian ancestry controls (17, 18). Therefore, we might conclude that this common variant does not have a role on the development of HSCR in Indonesia.

ENS development is a complex process, regulated by a large range of molecules and signaling pathways. This strictly controlled process needs the correct regulation of ENS-specific gene expression. The expressions of several genes implicated in HSCR development have been shown to be regulated by epigenetic mechanisms. RET expression was regulated by retinoic acid through DNA methylation on this 5′-CG-3′-rich enhancer region (30), while a significantly lower level of EDNRB methylation was also detected in HSCR patients (31). Interestingly, our study clearly demonstrates that SEMA3D gene expression was strongly up-regulated in the ganglionic intestines of HSCR patients as compared to controls. To the best of our knowledge, our study is the first report of the aberrant expressions of SEMA3D in HSCR patients. Another novelty in our findings is we conducted the study on patients with Indonesian ancestry [vs. European (4) ancestry cases]. It has been shown that reduced sema3d expression by morpholino resulted in severe effects on the intestine and its innervation of zebrafish (4). In addition, sema3d bound to nrp1a to facilitate the axonal guidance and contributed to peripheral axon outgrowth interdependently with dpysl3 (32, 33) Binding of Sema3d to neuropilin and plexin receptors introduced biochemical responses in specific neurons and stimulated the neuron migration (34, 35). Therefore, we might hypothesize that the aberrant expressions of SEMA3D will have an impact in our HSCR patients by affecting the neuronal guidance during ENS development.

It has been shown that some HSCR patients have persistent bowel symptoms, such as constipation, soiling, and enterocolitis, after an appropriate pull-through procedure. Most HSCR patients with persistent bowel symptoms do not have any identifiable etiology for their ongoing bowel dysmotility (36). The current hypothesis is that the aberrant expression of some genes in the ganglionic colon of HSCR patients includes SK3, Cx26, ChAT, and nNOS (37–40). Our study presented the altered SEMA3D expressions in the ganglionic colon of HSCR patients. Therefore, we might hypothesize that the aberrant SEMA3D expressions in the ganglionic colon involve in the pathogenesis of persistent bowel symptoms in HSCR patients following a properly performed pull-through surgery.

Our results should be interpreted with some caution, however, because they are based on overall means without accounting for other factors generating variation in the data such as gender, age, and degree of aganglionosis. Further studies of the methylation pattern of the SEMA3D gene are required to investigate whether the aberrant expression of SEMA3D in HSCR patients is due to abnormal DNA methylation. It is also necessary to compare the SEMA3D protein expression level between HSCR patients and controls and to identify its location in the colon tissue to prove the increased SEMA3D expressions. Unfortunately, we do not have any data on the methylation pattern, protein expression, and immunohistochemistry of SEMA3D due to resource limitations in our institution. Furthermore, our study utilized ARM patients' colon as controls. It has been shown that most ARM patients have abnormal colonic motility (41) and low expression of interstitial cells of Cajal marker (42). Therefore, further research with more proper control colon (e.g., autopsy bowel material from healthy infants or trauma patients) is necessary to better confirm the role of SEMA3D expression in the pathogenesis of HSCR. Noteworthy, however, our small sample size is a limitation of our study and suggests that a multicenter larger sample population needs to be studied to clarify our findings.

In conclusion, 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.

Data Availability Statement

All data generated or analyzed during this study are included in the submission. The raw data are available from the corresponding author on reasonable request.

Ethics Statement

The Institutional Review Board of the Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada/Dr. Sardjito Hospital gave approval for this study (KE/FK/1356/EC/2015). All parents signed a written informed consent form before participating in this study.

Author Contributions

G and KI conceived the study. G drafted the manuscript. KI and JV critically revised the manuscript for important intellectual content. AK, NB, and FR performed total RNA extraction and qPCR. HH, AM, MF, and DY conducted the experimental PCR-based work for Sanger sequencing. All authors have read and approved the manuscript and agreed to be accountable for all aspects of the work in ensuring that questions related to the accuracy or integrity of any part of the work are appropriately investigated and resolved.

Funding

This work was supported by a grant from the Faculty of Medicine, Public Health and Nursing, Universitas Gadjah Mada, Yogyakarta, Indonesia (Teuku Jacob Fellowship 2015-2017: #UPPM/202/M/05/04/09.15 and #UPPM/256/05/04/05.17 to G).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

We are grateful to Program Bantuan Seminar Luar Negeri, Ditjen Penguatan Riset dan Pengembangan, Kemenristekdikti, Indonesia for providing the conference grant to present the abstract at the 52nd Annual Meeting of PAPS in Christchurch, New Zealand at March 10–14, 2019. We would like to thank the patients, their families, and all those who were involved and provided excellent technical support and assistance during the study. We also thank Professor Aravinda Chakravarti (New York University, School of Medicine) for critically reading the manuscript and his suggestions.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fped.2020.00060/full#supplementary-material

Supplemental Figure S1. Gel electrophoresis of qPCR: (A) GAPDH and (B) SEMA3D expressions for HSCR ganglionic, aganglionic, and control colons. M: 100 bp DNA ladder; G: ganglionic colon; A: aganglionic colon; C: control colon; W: no RNA.

References

1. Tilghman JM, Ling AY, Turner TN, Sosa MX, Krumm N, Chatterjee S, et al. Molecular genetic anatomy and risk profile of Hirschsprung's disease. N Engl J Med. (2019) 380:1421–32. doi: 10.1056/NEJMoa1706594

PubMed Abstract | CrossRef Full Text | Google Scholar

2. Amiel J, Sproat-Emison E, Garcia-Barcelo M, Lantieri F, Burzynski G, Borrego S, et al. Hirschsprung disease, associated syndromes and genetics: a review. J Med Genet. (2008) 45:1–14. doi: 10.1136/jmg.2007.053959

PubMed Abstract | CrossRef Full Text | Google Scholar

3. Emison ES, Garcia-Barcelo M, Grice EA, Lantieri F, Amiel J, Burzynski G, et al. Differential contributions of rare and common, coding and noncoding ret mutations to multifactorial hirschsprung disease liability. Am J Hum Genet. (2010) 87:60–74. doi: 10.1016/j.ajhg.2010.06.007

PubMed Abstract | CrossRef Full Text | Google Scholar

4. Jiang Q, Arnold S, Heanue T, Kilambi KP, Doan B, Kapoor A, et al. Functional loss of semaphorin 3C and/or semaphorin 3D and their epistatic interaction with ret are critical to hirschsprung disease liability. Am J Hum Genet. (2015) 96:581–96. doi: 10.1016/j.ajhg.2015.02.014

PubMed Abstract | CrossRef Full Text | Google Scholar

5. Gunadi, Kapoor A, Ling AY, Rochadi, Makhmudi A, Herini ES, et al. Effects of ret and NRG1 polymorphisms in Indonesian patients with hirschsprung disease. J Pediatr Surg. (2014) 49:1614–8. doi: 10.1016/j.jpedsurg.2014.04.011

PubMed Abstract | CrossRef Full Text | Google Scholar

6. Gunadi, Iskandar K, Makhmudi A, Kapoor A. Combined genetic effect of RET and NRG1 variants on multifactorial hirschsprung disease in Indonesia. J Surg Res. (2019) 233:96–9. doi: 10.1016/j.jss.2018.07.067

CrossRef Full Text | Google Scholar

7. Gunadi, Karina SM, Dwihantoro A. Outcomes in patients with hirschsprung disease following definitive surgery. BMC Res Notes. (2018) 11:644. doi: 10.1186/s13104-018-3751-5

PubMed Abstract | CrossRef Full Text | Google Scholar

8. Kapoor A, Jiang Q, Chatterjee S, Chakraborty P, Sosa MX, Berrios C, et al. Population variation in total genetic risk of Hirschsprung disease from common RET, SEMA3 and NRG1 susceptibility polymorphisms. Hum Mol Genet. (2015) 24:2997–3003. doi: 10.1093/hmg/ddv051

PubMed Abstract | CrossRef Full Text | Google Scholar

9. Luzón-Toro B, Fernández RM, Torroglosa A, de Agustín JC, Méndez-Vidal C, Segura DI, et al. Mutational spectrum of semaphorin 3A and semaphorin 3D genes in Spanish Hirschsprung patients. PLoS ONE. (2013) 8:e54800. doi: 10.1371/journal.pone.0054800

PubMed Abstract | CrossRef Full Text | Google Scholar

10. Gunadi, Makhmudi A, Agustriani N, Rochadi. Effects of SEMA3 polymorphisms in Hirschsprung disease patients. Pediatr Surg Int. (2016) 32:1025–8. doi: 10.1007/s00383-016-3953-7

PubMed Abstract | CrossRef Full Text | Google Scholar

11. Tajima A, Pan IH, Fucharoen G, Fucharoen S, Matsuo M, Tokunaga K, et al. Three major lineages of Asian Y chromosomes: implications for the peopling of east and Southeast Asia. Hum Genet. (2002) 110:80–8. doi: 10.1007/s00439-001-0651-9

PubMed Abstract | CrossRef Full Text | Google Scholar

12. Torroglosa A, Alves MM, Fernández RM, Antiñolo G, Hofstra RM, Borrego S. Epigenetics in ENS development and hirschsprung disease. Dev Biol. (2016) 417:209–16. doi: 10.1016/j.ydbio.2016.06.017

PubMed Abstract | CrossRef Full Text | Google Scholar

13. O'Donnell AM, Coyle D, Puri P. Decreased expression of NEDL2 in Hirschsprung's disease. J Pediatr Surg. (2016) 51:1839–42. doi: 10.1016/j.jpedsurg.2016.06.016

PubMed Abstract | CrossRef Full Text | Google Scholar

14. Tang W, Li B, Xu X, Zhou Z, Wu W, Tang J, et al. Aberrant high expression of NRG1 gene in Hirschsprung disease. J Pediatr Surg. (2012) 47:1694–8. doi: 10.1016/j.jpedsurg.2012.03.061

PubMed Abstract | CrossRef Full Text | Google Scholar

15. Gunadi, Budi NYP, Kalim AS, Santiko W, Musthofa FD, Iskandar K, et al. Aberrant expressions of miRNA-206 target, FN1, in multifactorial hirschsprung disease. Orphanet J Rare Dis. (2019) 14:5. doi: 10.1186/s13023-018-0973-5

PubMed Abstract | CrossRef Full Text | Google Scholar

16. Widyasari A, Pravitasari WA, Dwihantoro A, Gunadi. Functional outcomes in hirschsprung disease patients after transabdominal soave and duhamel procedures. BMC Gastroenterol. (2018) 18:56. doi: 10.1186/s12876-018-0783-1

CrossRef Full Text | Google Scholar

17. The (1000) Genomes Project Consortium, Abecasis GR, Altshuler D, Auton A, Brooks LD, Durbin RM, Gibbs RA, et al. A map of human genome variation from population-scale sequencing. Nature. (2010) 467:1061–73. doi: 10.1038/nature09534

PubMed Abstract | CrossRef Full Text | Google Scholar

18. Lek M, Karczewski KJ, Minikel EV, Samocha KE, Banks E, Fennell T, et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature. (2016) 536:285–91. doi: 10.1038/nature19057

PubMed Abstract | CrossRef Full Text | Google Scholar

19. Sanchez-Castro M, Pichon O, Briand A, Poulain D, Gournay V, David A, et al. Disruption of the SEMA3D gene in a patient with congenital heart defects. Hum Mutat. (2015) 36:30–3. doi: 10.1002/humu.22702

PubMed Abstract | CrossRef Full Text | Google Scholar

20. Livak KJ, Schmittgen TD. Analysis of relative gene expression data using real-time quantitative PCR and the 2(-Delta Delta C(T)) method. Methods. (2001) 25:402–8. doi: 10.1006/meth.2001.1262

PubMed Abstract | CrossRef Full Text | Google Scholar

21. Gunadi, Miura K, Ohta M, Sugano A, Lee MJ, Sato Y, et al. Two novel mutations in the ED1 gene in Japanese families with X-linked hypohidrotic ectodermal dysplasia. Pediatr Res. (2009) 65:453–7. doi: 10.1203/PDR.0b013e3181991229

PubMed Abstract | CrossRef Full Text | Google Scholar

22. Luzón-Toro B, Torroglosa A, Núñez-Torres R, Enguix-Riego MV, Fernández RM, de Agustín JC, et al. Comprehensive analysis of NRG1 common and rare variants in hirschsprung patients. PLoS ONE. (2012) 7:e36524. doi: 10.1371/journal.pone.0036524

PubMed Abstract | CrossRef Full Text | Google Scholar

23. Garcia-Barcelo MM, Tang CS, Ngan ES, Lui VC, Chen Y, So MT, et al. Genome-wide association study identifies NRG1 as a susceptibility locus for hirschsprung's disease. Proc Natl Acad Sci USA. (2009) 106:2694–9. doi: 10.1073/pnas.0809630105

PubMed Abstract | CrossRef Full Text | Google Scholar

24. Phusantisampan T, Sangkhathat S, Phongdara A, Chiengkriwate P, Patrapinyokul S, Mahasirimongkol S. Association of genetic polymorphisms in the RET-protooncogene and NRG1 with hirschsprung disease in Thai patients. J Hum Genet. (2012) 57:286–93. doi: 10.1038/jhg.2012.18

PubMed Abstract | CrossRef Full Text | Google Scholar

25. Makhmudi A, Sadewa AH, Aryandono T, Chatterjee S, Heij HA, Gunadi. Effects of MTHFR c.677C>T, F2 c.20210G>A and F5 Leiden polymorphisms in gastroschisis. J Invest Surg. (2016) 29:88–92. doi: 10.3109/08941939.2015.1077908

PubMed Abstract | CrossRef Full Text | Google Scholar

26. Cardonick E, Broth R, Kaufmann M, Seaton J, Henning D, Roberts N, et al. Genetic predispositions for thromboembolism as a possible etiology for gastroschisis. Am J Obstet Gynecol. (2005) 193:426–8. doi: 10.1016/j.ajog.2004.12.010

PubMed Abstract | CrossRef Full Text | Google Scholar

27. Torfs CP, Christianson RE, Iovannisci DM, Shaw GM, Lammer EJ. Selected gene polymorphisms and their interaction with maternal smoking, as risk factors for gastroschisis. Birth Defects Res A Clin Mol Teratol. (2006) 76:723–30. doi: 10.1002/bdra.20310

PubMed Abstract | CrossRef Full Text | Google Scholar

28. Kapoor A, Auer DR, Lee D, Chatterjee S, Chakravarti A. Testing the Ret and SEMA3D genetic interaction in mouse enteric nervous system development. Hum Mol Genet. (2017) 26:1811–20. doi: 10.1093/hmg/ddx084

PubMed Abstract | CrossRef Full Text | Google Scholar

29. Tang CS, Gui H, Kapoor A, Kim JH, Luzón-Toro B, Pelet A, et al. Trans-ethnic meta-analysis of genome-wide association studies for Hirschsprung disease. Hum Mol Genet. (2016) 25:5265–75. doi: 10.1093/hmg/ddw333

PubMed Abstract | CrossRef Full Text | Google Scholar

30. Angrisano T, Sacchetti S, Natale F, Cerrato A, Pero R, Keller S, et al. Chromatin and DNA methylation dynamics during retinoic acid-induced RET gene transcriptional activation in neuroblastoma cells. Nucleic Acids Res. (2011) 39:1993–2006. doi: 10.1093/nar/gkq864

PubMed Abstract | CrossRef Full Text | Google Scholar

31. Tang W, Li B, Tang J, Liu K, Qin J, Wu W, et al. Methylation analysis of EDNRB in human colon tissues of Hirschsprung's disease. Pediatr Surg Int. (2013) 29:683–8. doi: 10.1007/s00383-013-3308-6

PubMed Abstract | CrossRef Full Text | Google Scholar

32. Wolman MA, Liu Y, Tawarayama H, Shoji W, Halloran MC. Repulsion and attraction of axons by semaphorin3D are mediated by different neuropilins in vivo. J Neurosci. (2004) 24:8428–35. doi: 10.1523/JNEUROSCI.2349-04.2004

PubMed Abstract | CrossRef Full Text | Google Scholar

33. Tanaka H, Nojima Y, Shoji W, Sato M, Nakayama R, Ohshima T, et al. Islet1 selectively promotes peripheral axon outgrowth in Rohon-Beard primary sensory neurons. Dev Dyn. (2011) 240:9–22. doi: 10.1002/dvdy.22499

PubMed Abstract | CrossRef Full Text | Google Scholar

34. Tran TS, Kolodkin AL, Bharadwaj R. Semaphorin regulation of cellular morphology. Annu Rev Cell Dev Biol. (2007) 23:263–92. doi: 10.1146/annurev.cellbio.22.010605.093554

PubMed Abstract | CrossRef Full Text | Google Scholar

35. Osborne NJ, Begbie J, Chilton JK, Schmidt H, Eickholt BJ. Semaphorin/neuropilin signaling influences the positioning of migratory neural crest cells within the hindbrain region of the chick. Dev Dyn. (2005) 232:939–49. doi: 10.1002/dvdy.20258

PubMed Abstract | CrossRef Full Text | Google Scholar

36. Zimmer J, Tomuschat C, Puri P. Long-term results of transanal pull-through for Hirschsprung's disease: a meta-analysis. Pediatr Surg Int. (2016) 32:743–9. doi: 10.1007/s00383-016-3908-z

PubMed Abstract | CrossRef Full Text | Google Scholar

37. Coyle D, O'Donnell AM, Puri P. Altered distribution of small-conductance calcium-activated potassium channel SK3 in Hirschsprung's disease. J Pediatr Surg. (2015) 50:1659–64. doi: 10.1016/j.jpedsurg.2015.01.013

PubMed Abstract | CrossRef Full Text | Google Scholar

38. Gunadi, Sunardi M, Budi NYP, Kalim AS, Iskandar K, Dwihantoro A. The impact of down-regulated SK3 expressions on Hirschsprung disease. BMC Med Genet. (2018) 19:24. doi: 10.1186/s12881-018-0539-3

PubMed Abstract | CrossRef Full Text | Google Scholar

39. Coyle D, Doyle B, Murphy JM, O'Donnell AM, Gillick J, Puri P. Expression of connexin 26 and connexin 43 is reduced in Hirschsprung's disease. J Surg Res. (2016) 206:242–51. doi: 10.1016/j.jss.2016.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

40. Coyle D, O'Donnell AM, Gillick J, Puri P. Altered neurotransmitter expression profile in the ganglionic bowel in Hirschsprung's disease. J Pediatr Surg. (2016) 51:762–9. doi: 10.1016/j.jpedsurg.2016.02.018

PubMed Abstract | CrossRef Full Text | Google Scholar

41. Levitt MA, Peña A. Anorectal malformations. Orphanet J Rare Dis. (2007) 2:33. doi: 10.1186/1750-1172-2-33

PubMed Abstract | CrossRef Full Text | Google Scholar

42. Xiao H, Huang R, Cui DX, Xiao P, Diao M, Li L. Histopathologic and immunohistochemical findings in congenital anorectal malformations. Medicine. (2018) 97:e11675. doi: 10.1097/MD.0000000000011675

PubMed Abstract | CrossRef Full Text | Google Scholar