Abstract
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 Χ2 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.
Significance of this study
What is already known about this subject?
The morbidity of Hirschsprung’s disease (HD) was ranked high at second of all disorders relevant to newborn congenital digestive tract malformation.
About 20% of patients with HD were accompanied with family heredity.
The development of HD may be related with migration of neural crest cells in the embryonic period.
What are the new findings?
The C allele of NRXN1 rs1421589 acted as a protective element for HD risk.
The mutant alleles of rs11795613 (A>G) and rs4844285 (G>A) within NLGN3 were both associated with significantly lower susceptibility to HD than their wide alleles.
The C allele of NLGX4X rs5961397 could reduce the risk of HD in comparison to the T allele.
The subjects carrying mutant alleles of rs7157669 and rs724373 were more prone to suffer from HD than carriers of the corresponding wide alleles.
How might these results change the focus of research or clinical practice?
The results strongly evidenced the hypothesis that the aberrant neuroligin and neurexin expressions caused by genetic mutations participated in the connections between intestinal neurons, leading to developmental immaturity of enteric nervous system and thus the presence of such diseases as HD.
Introduction
Hirschsprung’s disease (HD), a developmental disorder, was featured by congenital malformations in the gastrointestinal tract, and it could be commonly found among cases of pediatric surgery. The morbidity of HD (ie, 1/2000–1/5000) was ranked high at second of all disorders relevant to newborn congenital digestive tract malformation, and males possessed threefold more susceptibility to HD than females.1 The pathological changes of HD lied in the dysfunctional enteric nervous system (ENS), which was particularly manifested as shortage of intestinal ganglion cells within myenteric nerve plexus or intestinal submucosa of distal spastic colon and disorderly arranged nerve fibers. The above symptoms would render intestinal canal to take the spastic shape, making the fecal sediment difficult to empty the dung and facilitating compensatory dilation of proximal intestine.2 3
Up to now, the familial incidence of HD was demonstrated to account for around 3.6%–7.8% of the overall incidence of HD, and about 20% of patients with HD were accompanied with family heredity.4 Thus, the consensus has been reached that one dominating etiologic factor of HD was genetic variation of ENS. For instance, single nucleotide polymorphisms (SNPs) within RET (ie, rs2435357 and rs2506030) and SEMA3 (ie, rs12707682) could impart tremendous risk of HD development and disorders of phenotypic conditions among a Chinese population.5 A Chinese genome-wide association study disclosed that NRG3 mattered a lot to the onset of HD,6 yet later a case–control study did not uncover the association of three NRG3 SNPs with HD within a Chinese crowd.7 Nonetheless, it was documented that rs78356888 of NRG1, as well as rs2506030 and rs2435357 of RET, was significantly correlated with susceptibility to HD.7 Since the study results about the association of SNPs with HD risk were quite limited and were sometimes in dispute, finding novel target SNPs or conduction of replicated studies were critical to hindering the HD onset.
Of note, the achievements related with the development of central nervous system were vital to the exploration with respect to development of ENS, for that ENS and central nervous system were both originating from embryonic neural crest.8 Since production of synapse symbolized complete development of central nervous system, the sensory and motor functions of ENS should also be guaranteed by the appropriate synaptic actions. It has been verified that neurexin family and neuroligin family were located within presynaptic membranes or across the subsynaptic membranes.9 10 Furthermore, neurexin might interact with neuroligin to mediate synaptic junctions within the central nervous system via affecting scaffolding protein.11–13 Hence, it was hypothesized that mutants of neurexin and neuroligin SNPs were highly potential factors accounting for ENS-relevant disorders, such as HD.
Accordingly, the current study was aimed to explore whether SNPs within neurexin and neuroligin were associated with susceptibility to HD.
Methods
Subjects
Totally 210 children with HD were retrospectively recruited from Zhoukou Central Hospital from March 2011 to February 2016. They all experienced the transanal pull-through operation for the first time. According to the diagnostic criteria for HD declared by the fourth international conference on HD and related neural crest diseases,14 the children with HD were preoperatively managed with barium enema examination, rectum mucosa biopsy, H&E staining and acetylcholinesterase staining. Consecutive sections were performed for each sample, and HD was confirmed when hardly any neurons were observed and obvious hyperplasia of nerve plexus were found within submucosa. Since there existed a cell area featured by physiological hyporesponsiveness and aganglionosis beyond the normal rectal line, biopsies for HD diagnosis were obtained based on age: (1) 2 cm beyond the line for newborns; (2) 2.5 cm beyond the line for the children aged <1 year; (3) 3 cm beyond the line for the children aged 1–3 years; (4) 3.5 cm beyond the line for the children aged >4 years. In addition, HD was divided into the short-segment HD, the long-segment HD and the total-segment HD. In particular, the lesions of short-segment HD were limited to be within distal sigmoid, rectal, as well as the border between distal sigmoid and rectal; the lesions of long-segment HD covered proximal sigmoid, ascending colon and even a wider range; the lesions of total-segment HD ran through rectum, total colon and terminal ileum (<30 cm far from ileocecal valve). Besides, 187 normal children receiving the physical examination in our hospital were included in the control group, and they all had no malformations relevant to digestive tract and neural crest. Parents of all subjects have signed informed consents.
Selection of SNPs
The investigated SNPs were partly selected with aid of the bioinformation database (http://www.ncbi.nlm.nih.gov/SNP/). Specifically, with the SNP browser software, tag SNPs were chosen if their r2 was more than 0.8 and their minor allelic frequency was at least 5%, in accordance with the pairwise r2 method. Other SNPs were deemed as eligible based on related documentations.15 16 Finally, SNPs within NRXN1 (ie, rs1363032 and rs1424589), NRXN3 (ie, rs11624704, rs7157669 and rs724373), NLGN1 (rs13074723 and rs1488547), NLGN3 (rs11795613 and rs4844285) and NLGX4X (rs6529901 and rs5961379) were chosen and arranged for the following studies.
Genotyping
Peripheral blood (volume: 2 mL) was extracted from HD and healthy children and 3.8% sodium citrate was used for anticoagulation. The Qiagen DNA extraction kit was utilized to extract the DNA, and specific procedures were carried out strictly in accordance with the manufacturers’ protocols. The primers were designed based on the primer 5 software and were synthesized with assistance of Shanghai Invitrogen Biotechnology Co (Supplementary table 1). All the primers were all purified via polyacrylamide gel electrophoresis and were then recycled using C18 column.
Supplementary file 1
Then PCR-restriction fragment length polymorphism (RFLP) was implemented to genotype the SNPs. The PCR reaction system (50 µL) mainly included double distilled water (32.5 µL), 10× buffer (5.0 µL), 25 mmol/L magnesium chloride (4.0 µL), 10 mmol/L deoxyribonucleotide triphosphate (4.0 µL), 20 μmol/L sense primer (1.0 µL), 20 μmol/L antisense primer (1.0 µL), 5 U/µL Taq DNA polymerase and 100 ng DNA template (2.0 µL). Moreover, the PCR reaction was proceeded according to the following process: (1) predenaturation at 95°C for 5 min; (2) 35 cycles of 95°C for 40 s, 60°C for 30 s and 72°C for 30 s and (3) extension at 72°C for 10 min and (4) conservation at 4°C.
Operation
The operation methods were selected in line with the age of children and the length of pathological intestinal tract. For the children aged between 6 months and 1 year, they would be treated with the transanal endorectal pull-through operation if they were diagnosed as short-segment HD. The children diagnosed as long-segment HD were alternatively treated with transabdominal heart anastomosis. As for the children with total colonic aganglionosis, they were first managed with colostomy and were then treated with Martin surgery or ileoanal anastomosis when they were around 1 year old. The children aged more than 1 year experienced the improved Duhamel surgery or transabdominal heart anastomosis.
Evaluation of HD phenotypes
The biopsies of all HD children were retrieved. In line with the range of the affected bowels, the samples were classified as short segment type, long segment type and total type (ie, the most serious type).17
Statistical analysis
The measurement data were exhibited in the form of mean±SD, and the enumeration data were presented as percentages (%). Besides, Χ2 test was conducted to compare the enumeration data, and t-test or Mann-Whitney test was utilized to analyze the measurement data. Furthermore, the Hardy-Weinberg equilibrium law was applied to examine whether the gene and genotype frequencies were balanced among the study samples. The linkage disequilibrium (LD) among neurexin/neuroligin SNPs was fitted based on Shesis software, and the haplotypes of the SNPs were also analyzed with the Shesis software.18 In addition, the correlations between potential SNPs/haplotypes and risk/phenotype of HD were shown in the form of ORs and 95% CIs. It would be deemed as statistically significant when P was less than 0.05. All the statistical analyses were handled with SPSS V.13.0 software.
Results
Baseline characteristics of the subjects
There showed hardly any significant distinctions between the HD group and healthy controls in the aspects of mean age, mean weight and mean height (P>0.05) (table 1). Moreover, the HD children altogether possessed 69 (32.86%) short segment type, 101 (40.09%) long segment type and 40 (19.05%) total type.
Baseline characteristics of the children with Hirschsprung’s disease (case group) and healthy children (control group)
Correlation between SNPs/haplotypes within neurexin/neuroligin and susceptibility to HD
Regarding the studied population, the LD of the selected neurexin/neuroligin SNPs was shown as figure 1. The C allele of NRXN1 rs1421589 acted as a protective element for HD risk when compared with T allele (OR 0.59, 95% CI 0.44 to 0.81, P=0.001) (table 2). Similarly, the mutant alleles of rs11795613 (A>G) and rs4844285 (G>A) within NLGN3 were both associated with significantly lower susceptibility to HD than their wide alleles (OR 0.39, 95% CI 0.29 to 0.51, P<0.001; OR 0.55, 95% CI 0.41 to 0.73, P<0.001). Moreover, the C allele of NLGX4X rs5961397 could reduce the risk of HD in comparison to the T allele (OR 0.51, 95% CI 0.38 to 0.68, P<0.001). In contrast, the subjects carrying mutant alleles of rs7157669 and rs724373 were more prone to suffer from HD than carriers of the corresponding wide alleles (A vs C: OR 1.44, 95% CI 1.03 to 2.02, P<0.034; C vs T: OR 2.66, 95% CI: 1.88 to 3.78, P<0.001).
- Download figure
- Open in new tab
- Download powerpoint
The linkage disequilibrium of single nucleotide polymorphism within (A) NRXN1 (ie, rs1363032 (chromosome 2: 49993858) and rs1421589 (chromosome 2: 50063015)), (B) NRXN3 ([ie, rs11624704 (chromosome 14: 78319734), rs7157669 (chromosome 14: 78471334) and rs724373 (chromosome 14: 78476555)), (C) NLGN1 (ie rs13074723 (chromosome 3: 173804307) and rs1488547 (chromosome 3: 173807978)) and (D) NLGN3 (ie rs11795613 (chromosome X: 71147478) and rs4844285 (chromosome X: 71150394) and NLGX4X (ie rs6529901 (chromosome X: 5937979) and rs5961397 (chromosome X: 6004425)).
Correlation between SNPs within neurexin/neuroligin and susceptibility to HD
After investigating the combined role of the significant SNPs mentioned above in facilitating HD risk (table 3), it was observed that haplotypes CT of NRXN1, haplotypes AAC, ACT and CAC of NRXN3, haplotypes AG and GA of NLGN1, as well as haplotypes AGCT, AGTT and GACC of NLGN3 and NLGX4X, respectively, served as the hazard parameter for susceptibility to HD (OR 1.74, 95% CI 1.25 to 2.42, P=0.001; OR 1.97, 95% CI 1.47 to 2.63, P<0.001; OR 2.05, 95% CI 1.26 to 3.31, P=0.003; OR 2.22, 95% CI 1.41 to 3.50, P<0.001; OR 2.97, 95% CI 2.20 to 3.99, P<0.001; OR 2.37, 95% CI 1.74 to 3.22, P<0.001; OR 6.34, 95% CI 3.53 to 11.40, P<0.001; OR 10.02, 95% CI 6.00 to 16.75, P<0.001; OR 3.40, 95% CI 2.31 to 5.01, P<0.001). On the contrary, carriers with haplotypes ACC of NRXN3, haplotype GG of NLGN1 and haplotypes AATT or GACT of NLGN3 and NLGX4X possessed lower possibility of attacking HD (OR 0.28, 95% CI 0.17 to 0.46, P<0.001; OR 0.14, 95% CI 0.08 to 0.24, P<0.001; OR 0.00, 95% CI 0.00 to 0.01, P<0.001; OR 0.24, 95% CI 0.13 to 0.46, P<0.001).
Correlation between haplotypes within neurexin/neuroligin and susceptibility to HD
Correlation between SNPs within neurexin/neuroligin and HD phenotypes
With respect to NRXN3 rs724373 (table 4), the distribution of its genotypes (ie, TT, TC and CC) within the long-segment HD population and the total-segment HD population, respectively, differed from the short-segment HD ones (P=0.024 and P=0.045). The rs11795613 was also significantly linked with altered phenotypes of HD (long-segment HD vs short-segment HD: P=0.003; total-segment HD vs short-segment HD: P=0.004). Ultimately, the genotyping frequencies of NLGN3 rs4844285 (G>A) was prominently discrepant between the total-segment HD versus short-segment HD groups (P=0.014), and NLGX4X rs5961397 (T>C) also displayed different genotyping frequencies among the long-segment HD group, when compared with the short-segment HD group (P=0.001).
Correlation between SNPs within neurexin/neuroligin and HD phenotypes
Discussion
ENS was composed of myenteric plexus and submucous plexus, which interacted with each other to regulate the motor and secretion function of the intestinal tract.19 In-depth studies also indicated that ENS was crucial to the development of the diseases relevant to the intestinal dynamics and gastrointestinal functions. The normal functioning of enteric nerves was dependent on the integral connection of intestinal nerve cells, intestinal muscle cells and epithelial cells. It was demonstrated that genetic mutations could induce disordered gastrointestinal dynamics, by way of rendering loss of interneuron functions and disconnecting between neurons and their effectors.20 It was the synaptic connection that prompted intestinal tubes to perform segmental motion, reflection secretion and neural immune responses.21 Hence, it was hypothesized that the genetic mutations that modified the synaptic function were involved in the development and maturation of ENS.
Neurexin and neuroligin have been verified to mediate the connection and differentiation of synapses through glycosylation, oligomerization and calcium binding (figure 2).10 22 It has been confirmed in a glycosylation experiment that the O-bound and N-bound saccharides could regulate neuroligin that was rich in Ser–Thr domains,23 and the ChE-like domain (CLD) structural domain residues of neuroligin-1 might also affect the interaction between neuroligin and neurexin.22 24 Of note, the mutations of EDNRB could contribute to less inflow of calcium ions, and the attenuated binding of Ca2+-dependent neuroligin and neurexin, which finally hindered the regular shaping of synapses and also ENS. In the meantime, mutations of neuroligin and neurexin would also aggravate this condition via damaging the aggregation of postsynaptic proteins and neuroligins.22
- Download figure
- Open in new tab
- Download powerpoint
The mechanism underlying the effects of neurexin and neuroligin on synapse-related information transfer. As an action potential is initiated (①), the Ca2+-channel of the presynaptic terminal would be open (②), which triggers inflowing of Ca2+. Subsequently, the synaptic vesicles within the presynaptic cell would be fused with the presynaptic plasma membrane (③), and neurotransmitters are released from vesicles into the synaptic cleft (④). Ultimately, the neurotransmitters would contact with the postsynaptic receptors, through which the process of information transfer is accomplished (⑤).
Virtually, certain SNPs appeared to decisively modify the functions of neuroligin and neurexin. For instance, a haplotype composed of six NLGN4X SNPs, namely, rs663857 (G), rs3810687 (G), rs3810688 (T), rs1882260 (T), rs3810686 (C) and rs5916269 (G), was regarded as a hazard for autism risk among an Italian population.25 Additional SNPs have also been affirmed with their impacts on the central nervous system-induced disorders, including rs3747333 or rs3747334 of NLGN4X and rs2303298 of NRXN-1.26–28. Accordingly, in our study, it was well founded that the mutated SNPs of NRXN1 (ie, rs1421859), NRXN3 (ie, rs7157669 and rs724373), NLGN3 (ie, rs11795613 and rs4844285) and NLGX4X (ie, rs5961397) were involved in the etiology of ENS immaturation-caused HD. Nevertheless, the SNPs considered as crucial to HD were, in some degree, distinct from ones that exhibited significance in encouraging the development of other diseases. This phenomenon could be attributed to that SNPs that interfered with the occurrence of different diseases were also disparate, for that they might participate in different molecular mechanisms when diverse diseases were considered.
Similarly, diverse SNPs of the genes that impacted the functioning of nervous system have been confirmed to be associated with incremental susceptibility to HD. For instance, recent studies documented the association of SNPs within catechol-O-methyltransferase (COMT) (ie, rs6267) and armadillo repeat gene deleted in velocardiofacial syndrome (ARVCF) (ie, rs80068543) with HD.29 The COMT was involved with canalization of brain catecholamine neurotransmitters and COMT interacted with ARVCF to modify the development of neuronal system.30 Furthermore, as high as 50% of familial HD cases were discovered with missing of RET functional mutations, and it was demonstrated that the intron 1 of RET gene shared rs2435357 with a critical transcription factor of ENS.31 32 In addition, γ-aminobutyricacid A receptor gamma 2 (GABRG2) mediated the neuron-inhibition process of the mammalian central nervous system,33 and its SNPs (ie, rs209350 and rs169793) exhibited significantly distinct frequencies between HD subjects and controls.34 Finally, RELN was correlated with migration and positioning of neuronal cells, and its tag SNPs also displayed marked association with HD risk, including rs802788, rs6977616 and rs56345626.34 35 It was further deeply stressed that whether the nervous system could work properly was vital to regulating the risk of HD.
To sum up, this hypothesis-generating study was the first one to investigate the correlation between SNPs relevant to NLGN or NRXN and risk of HD. Nonetheless, the selected SNPs might not have the largest representation of NLGN and NRXN functioning among all the genetic mutations, and more meaningful SNPs could be missed. Moreover, the research subjects were merely a Chinese Han population recruited from the hospital, so a large-scale prospective study targeting the entire population were demanded to replicate the study results. In addition, the limited size and confined ethnicities of the included children less able to be generalized to other populations, and the smaller population managed for the haplotype analysis made the results less convincing and less recommended for clinical application. Furthermore, thorough analyses regarding the interaction of exposure factors and SNPs on HD risk should also be completed in future. Finally, how the relationship between the SNPs and pathogenesis of HD was established remained unknown. Hence, in vivo and in vitro studies should be carried out to explore the potential mechanisms.
In conclusion, certain SNPs/haplotypes within NLGN or NRXN were believed as the biomarkers for risk of HD, though more diverse studies were in need to confirm the study result.
Footnotes
Contributors YL analyzed the data. HL and YD drafted the manuscript. All the authors contributed to the conception and design of the experiments and performed the experiments. All authors read and approved the final manuscript.
Funding This research received no specific grant from any funding agency in the public, commercial or not-for-profit sectors.
Competing interests None declared.
Patient consent Obtained.
Ethics approval The ethics committee of Zhoukou Central Hospital.
Provenance and peer review Not commissioned; externally peer reviewed.
This is an Open Access article distributed in accordance with the Creative Commons Attribution Non Commercial (CC BY-NC 4.0) license, which permits others to distribute, remix, adapt, build upon this work non-commercially, and license their derivative works on different terms, provided the original work is properly cited and the use is non-commercial. See: http://creativecommons.org/licenses/by-nc/4.0/