Exome sequencing reveals a high genetic heterogeneity on familial Hirschsprung disease

Comprehensive characterization of the genetic landscape of familial Hirschsprung's disease

BACKGROUND

Hirschsprung's disease (HSCR) is one of the most common congenital digestive tract malformations and can cause stubborn constipation or gastrointestinal obstruction after birth, causing great physical and mental pain to patients and their families. Studies have shown that more than 20 genes are involved in HSCR, and most cases of HSCR are sporadic. However, the overall rate of familial recurrence in 4331 cases of HSCR is about 7.6%. Furthermore, familial HSCR patients show incomplete dominance. We still do not know the penetrance and genetic characteristics of these known risk genes due to the rarity of HSCR families.

METHODS

To find published references, we used the title/abstract terms "Hirschsprung" and "familial" in the PubMed database and the MeSH terms "Hirschsprung" and "familial" in Web of Science. Finally, we summarized 129 HSCR families over the last 40 years.

RESULTS

The male-to-female ratio and the percentage of short segment-HSCR in familial HSCR are much lower than in sporadic HSCR. The primary gene factors in the syndromic families are ret proto-oncogene (RET) and endothelin B receptor gene (EDNRB). Most families show incomplete dominance and are relevant to RET, and the RET mutation has 56% penetrance in familial HSCR. When one of the parents is a RET mutation carrier in an HSCR family, the offspring's recurrence risk is 28%, and the incidence of the offspring does not depend on whether the parent suffers from HSCR.

CONCLUSION

Our findings will help HSCR patients obtain better genetic counseling, calculate the risk of recurrence, and provide new insights for future pedigree studies.

Whole genome sequencing reveals epistasis effects within RET for Hirschsprung disease

Common variants in RET and NRG1 have been associated with Hirschsprung disease (HSCR), a congenital disorder characterised by incomplete innervation of distal gut, in East Asian (EA) populations. However, the allelic effects so far identified do not fully explain its heritability, suggesting the presence of epistasis, where effect of one genetic variant differs depending on other (modifier) variants. Few instances of epistasis have been documented in complex diseases due to modelling complexity and data challenges. We proposed four epistasis models to comprehensively capture epistasis for HSCR between and within RET and NRG1 loci using whole genome sequencing (WGS) data in EA samples. 65 variants within the Topologically Associating Domain (TAD) of RET demonstrated significant epistasis with the lead enhancer variant (RET+3; rs2435357). These epistatic variants formed two linkage disequilibrium (LD) clusters represented by rs2506026 and rs2506028 that differed in minor allele frequency and the best-supported epistatic model. Intriguingly, rs2506028 is in high LD with one cis-regulatory variant (rs2506030) highlighted previously, suggesting that detected epistasis might be mediated through synergistic effects on transcription regulation of RET. Our findings demonstrated the advantages of WGS data for detecting epistasis, and support the presence of interactive effects of regulatory variants in RET for HSCR.

Circ-ITCH overexpression promoted cell proliferation and migration in Hirschsprung disease through miR-146b-5p/RET axis

BACKGROUND

Hirschsprung disease (HSCR) is a congenital intestinal disease caused by the abnormal proliferation and migration of enteric nerve cells (ENCC). Research suggested critical roles for circular RNA (circRNA) itchy E3 ubiquitin protein ligase (ITCH) in gastrointestinal malignancies progression. However, the function of circ-ITCH in HSCR remains poorly defined.

METHODS

The related genes expression in 30 HSCR patients and 30 controls without HSCR were detected using qRT-PCR. Cell proliferation was assessed by CCK-8 assay and EdU assay. Cell migration was detected with wound-healing assay and transwell assay. The interactions among circ-ITCH, miR-146b-5p, and RET were confirmed by Dual luciferase reporter assay.

RESULTS

Circ-ITCH and RET expressions were downregulated in HSCR patients and cells, while the miR-146b-5p expression was upregulated. Circ-ITCH overexpression facilitated cell proliferation, migration, and activated MAPK pathway, which were reversed by circRNA-ITCH knockdown. Circ-ITCH negatively regulated miR-146b-5p expression. MiR-146b-5p overexpression abolished the promoting effects of circ-ITCH overexpression on cell proliferation and migration. MiR-146b-5p inhibited RET expression. RET overexpression eliminated the inhibitory effects of miR-146b-5p overexpression on cell proliferation and migration.

CONCLUSION

Circ-ITCH overexpression facilitated cell proliferation and migration in HSCR by regulating miR-146b-5p/RET/MAPK axis.

IMPACT

The expressions of Circ-ITCH and RET were markedly reduced in HSCR, while miR-146b-5p expression was increased in HSCR. Circ-ITCH overexpression enhanced the proliferative and migratory abilities of SH-SY5Y and 293T cells. Circ-ITCH negatively regulated miR-146b-5p expression.

Clinical description and genetic analysis of a novel familial skeletal dysplasia characterized by high bone mass and lucent bone lesions

High bone mass (HBM) disorders are a clinically and genetically heterogeneous subgroup of rare skeletal dysplasias. Here we present a case of a previously unreported familial skeletal dysplasia characterized by HBM and lucent bone lesions that we aimed to clinically characterize and genetically investigate. For phenotyping, we reviewed past clinical records and imaging tests, and performed physical examination (PE), bone densitometry, and mineral panels in affected individuals, including a male proband, his son and daughter, in addition to unaffected controls, including the proband's wife and brother. Affected individuals also underwent impact microindentation (IMI). In an effort to elucidate the disorder's molecular etiology, whole exome sequencing (WES) was performed in all individuals to filter for rare variants present only in affected ones. The cases displayed a unique skeletal phenotype with a mix of sclerotic features and lucent bone lesions, and high IMI values. Bone mineral density was very elevated in the proband and his daughter. The proband's daughter also exhibited idiopathic scoliosis (IS), in addition to mild thrombocytopenia and mild structural thyroid abnormalities, which were the only extra-skeletal abnormalities identified. WES analysis yielded 5 rare putative pathogenic variants in affected members in genes that are associated with bone metabolism including: SEM4AD, TBX18, PTCH1, PTK7, and ADGRE5. The PTK7 variant appeared as possibly implicated in the development of IS while the TBX18 and SEMA4D variants stood out as the strongest candidates for the lucent bone lesions and HBM, respectively, given their high predicted pathogenicity and putative role in bone biology. Variant functionality should be addressed in the future to assess their implication in skeletal metabolism as it is the first time that mutations in TBX18 and SEMA4D have been associated to bone developmental lesions and mineral metabolism in a clinical setting.

Identification of FAT3 as a new candidate gene for adolescent idiopathic scoliosis

In an effort to identify rare alleles associated with adolescent idiopathic scoliosis (AIS) whole-exome sequencing was performed on a discovery cohort of 73 unrelated patients and 70 age-and sex matched controls, all of French-Canadian ancestry. A collapsing gene burden test was performed to analyze rare protein-altering variants using case–control statistics. Since no single gene achieved statistical significance, targeted exon sequencing was performed for 24 genes with the smallest p values, in an independent replication cohort of unrelated severely affected females with AIS and sex-matched controls (N = 96 each). An excess of rare, potentially protein-altering variants was noted in one particular gene, FAT3, although it did not achieve statistical significance. Independently, we sequenced the exomes of all members of a rare multiplex family of three affected sisters and unaffected parents. All three sisters were compound heterozygous for two rare protein-altering variants in FAT3. The parents were single heterozygotes for each variant. The two variants in the family were also present in our discovery cohort. A second validation step was done, using another independent replication cohort of 258 unrelated AIS patients having reach their skeletal maturity and 143 healthy controls to genotype nine FAT3 gene variants, including the two variants previously identified in the multiplex family: p.L517S (rs139595720) and p.L4544F (rs187159256). Interestingly, two FAT3 variants, rs139595720 (genotype A/G) and rs80293525 (genotype C/T), were enriched in severe scoliosis cases (4.5% and 2.7% respectively) compared to milder cases (1.4% and 0.7%) and healthy controls (1.6% and 0.8%). Our results implicate FAT3 as a new candidate gene in the etiology of AIS.

Disorders of the enteric nervous system - a holistic view

The enteric nervous system (ENS) is the largest division of the peripheral nervous system and closely resembles components and functions of the central nervous system. Although the central role of the ENS in congenital enteric neuropathic disorders, including Hirschsprung disease and inflammatory and functional bowel diseases, is well acknowledged, its role in systemic diseases is less understood. Evidence of a disordered ENS has accumulated in neurodegenerative diseases ranging from amyotrophic lateral sclerosis, Alzheimer disease and multiple sclerosis to Parkinson disease as well as neurodevelopmental disorders such as autism. The ENS is a key modulator of gut barrier function and a regulator of enteric homeostasis. A 'leaky gut' represents the gateway for bacterial and toxin translocation that might initiate downstream processes. Data indicate that changes in the gut microbiome acting in concert with the individual genetic background can modify the ENS, central nervous system and the immune system, impair barrier function, and contribute to various disorders such as irritable bowel syndrome, inflammatory bowel disease or neurodegeneration. Here, we summarize the current knowledge on the role of the ENS in gastrointestinal and systemic diseases, highlighting its interaction with various key players involved in shaping the phenotypes. Finally, current flaws and pitfalls related to ENS research in addition to future perspectives are also addressed.

Planar cell polarity (PCP) proteins support spermatogenesis through cytoskeletal organization in the testis

Few reports are found in the literature regarding the role of planar cell polarity (PCP) in supporting spermatogenesis in the testis. Yet morphological studies reported decades earlier have illustrated the directional alignment of polarized developing spermatids, most notably step 17-19 spermatids in stage V-early VIII tubules in the testis, across the plane of the epithelium in seminiferous tubules of adult rats. Such morphological features have unequivocally demonstrated the presence of PCP in developing spermatids, analogous to the PCP noted in hair cells of the cochlea in mammals. Emerging evidence in recent years has shown that Sertoli and germ cells express numerous PCP proteins, mostly notably, the core PCP proteins, PCP effectors and PCP signaling proteins. In this review, we discuss recent findings in the field regarding the two core PCP protein complexes, namely the Van Gogh-like 2 (Vangl2)/Prickle (Pk) complex and the Frizzled (Fzd)/Dishevelled (Dvl) complex. These findings have illustrated that these PCP proteins exert their regulatory role to support spermatogenesis through changes in the organization of actin and microtubule (MT) cytoskeletons in Sertoli cells. For instance, these PCP proteins confer PCP to developing spermatids. As such, developing haploid spermatids can be aligned and orderly packed within the limited space of the seminiferous tubules in the testes for the production of sperm via spermatogenesis. Thus, each adult male in the mouse, rat or human can produce an upward of 30, 50 or 300 million spermatozoa on a daily basis, respectively, throughout the adulthood. We also highlight critical areas of research that deserve attention in future studies. We also provide a hypothetical model by which PCP proteins support spermatogenesis based on recent studies in the testis. It is conceivable that the hypothetical model shown here will be updated as more data become available in future years, but this information can serve as the framework by investigators to unravel the role of PCP in spermatogenesis.

ChIP-Seq-Based Approach in Mouse Enteric Precursor Cells Reveals New Potential Genes with a Role in Enteric Nervous System Development and Hirschsprung Disease

Hirschsprung disease (HSCR) is a neurocristopathy characterized by intestinal aganglionosis which is attributed to a failure in neural crest cell (NCC) development during the embryonic stage. The colonization of the intestine by NCCs is a process finely controlled by a wide and complex gene regulatory system. Several genes have been associated with HSCR, but many aspects still remain poorly understood. The present study is focused on deciphering the PAX6 interaction network during enteric nervous system (ENS) formation. A combined experimental and computational approach was performed to identify PAX6 direct targets, as well as gene networks shared among such targets as potential susceptibility factors for HSCR. As a result, genes related to PAX6 either directly (RABGGTB and BRD3) or indirectly (TGFB1, HRAS, and GRB2) were identified as putative genes associated with HSCR. Interestingly, GRB2 is involved in the RET/GDNF/GFRA1 signaling pathway, one of the main pathways implicated in the disease. Our findings represent a new contribution to advance in the knowledge of the genetic basis of HSCR. The investigation of the role of these genes could help to elucidate their implication in HSCR onset.

A complementary study approach unravels novel players in the pathoetiology of Hirschsprung disease

Hirschsprung disease (HSCR, OMIM 142623) involves congenital intestinal obstruction caused by dysfunction of neural crest cells and their progeny during enteric nervous system (ENS) development. HSCR is a multifactorial disorder; pathogenetic variants accounting for disease phenotype are identified only in a minority of cases, and the identification of novel disease-relevant genes remains challenging. In order to identify and to validate a potential disease-causing relevance of novel HSCR candidate genes, we established a complementary study approach, combining whole exome sequencing (WES) with transcriptome analysis of murine embryonic ENS-related tissues, literature and database searches, in silico network analyses, and functional readouts using candidate gene-specific genome-edited cell clones. WES datasets of two patients with HSCR and their non-affected parents were analysed, and four novel HSCR candidate genes could be identified: ATP7A, SREBF1, ABCD1 and PIAS2. Further rare variants in these genes were identified in additional HSCR patients, suggesting disease relevance. Transcriptomics revealed that these genes are expressed in embryonic and fetal gastrointestinal tissues. Knockout of these genes in neuronal cells demonstrated impaired cell differentiation, proliferation and/or survival. Our approach identified and validated candidate HSCR genes and provided further insight into the underlying pathomechanisms of HSCR.

Whole-genome analysis of noncoding genetic variations identifies multiscale regulatory element perturbations associated with Hirschsprung disease

It is widely recognized that noncoding genetic variants play important roles in many human diseases, but there are multiple challenges that hinder the identification of functional disease-associated noncoding variants. The number of noncoding variants can be many times that of coding variants; many of them are not functional but in linkage disequilibrium with the functional ones; different variants can have epistatic effects; different variants can affect the same genes or pathways in different individuals; and some variants are related to each other not by affecting the same gene but by affecting the binding of the same upstream regulator. To overcome these difficulties, we propose a novel analysis framework that considers convergent impacts of different genetic variants on protein binding, which provides multiscale information about disease-associated perturbations of regulatory elements, genes, and pathways. Applying it to our whole-genome sequencing data of 918 short-segment Hirschsprung disease patients and matched controls, we identify various novel genes not detected by standard single-variant and region-based tests, functionally centering on neural crest migration and development. Our framework also identifies upstream regulators whose binding is influenced by the noncoding variants. Using human neural crest cells, we confirm cell stage-specific regulatory roles of three top novel regulatory elements on our list, respectively in the RET, RASGEF1A, and PIK3C2B loci. In the PIK3C2B regulatory element, we further show that a noncoding variant found only in the patients affects the binding of the gliogenesis regulator NFIA, with a corresponding up-regulation of multiple genes in the same topologically associating domain.

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.

Identification of New Potential LncRNA Biomarkers in Hirschsprung Disease

Hirschsprung disease (HSCR) is a neurocristopathy defined by intestinal aganglionosis due to alterations during the development of the Enteric Nervous System (ENS). A wide spectrum of molecules involved in different signaling pathways and mechanisms have been described in HSCR onset. Among them, epigenetic mechanisms are gaining increasing relevance. In an effort to better understand the epigenetic basis of HSCR, we have performed an analysis for the identification of long non-coding RNAs (lncRNAs) by qRT-PCR in enteric precursor cells (EPCs) from controls and HSCR patients. We aimed to test the presence of a set lncRNAs among 84 lncRNAs in human EPCs, which were previously related with crucial cellular processes for ENS development, as well as to identify the possible differences between HSCR patients and controls. As a result, we have determined a set of lncRNAs with positive expression in human EPCs that were screened for mutations using the exome data from our cohort of HSCR patients to identify possible variants related to this pathology. Interestingly, we identified three lncRNAs with different levels of their transcripts (SOCS2-AS, MEG3 and NEAT1) between HSCR patients and controls. We propose such lncRNAs as possible regulatory elements implicated in the onset of HSCR as well as potential biomarkers of this pathology.

Family-based exome sequencing combined with linkage analyses identifies rare susceptibility variants of MUC4 for gastric cancer

Genome-wide association studies of gastric cancer (GC) cases have revealed common gastric cancer susceptibility loci with low effect size. We investigated rare variants with high effect size via whole-exome sequencing (WES) of subjects with familial clustering of gastric cancer. WES of DNAs from the blood of 19 gastric cancer patients and 36 unaffected family members from 14 families with two or more gastric cancer patients were tested. Linkage analysis combined with association tests were performed using Pedigree Variant Annotation, Analysis, and Search Tool (pVAAST) software. Based on the logarithm of odds (LOD) and permutation-based composite likelihood ratio test (CLRT) from pVAAST, MUC4 was identified as a predisposing gene (LOD P-value = 1.9×10-5; permutation-based P-value of CLRT ≤ 9.9×10-9). In a larger cohort consisting of 597 GC patients and 9,759 healthy controls genotyped with SNP array, we discovered common variants in MUC4 regions (rs148735556, rs11717039, and rs547775645) significantly associated with GC supporting the association of MUC4 with gastric cancer. And the MUC4 variants were found in higher frequency in The Cancer Genome Atlas Study (TCGA) germline samples of patients with multiple cancer types. Immunohistochemistry indicated that MUC4 was downregulated in the noncancerous gastric mucosa of subjects with MUC4 germline missense variants, suggesting that loss of the protective function of MUC4 predisposes an individual to gastric cancer. Rare variants in MUC4 can be novel gastric cancer susceptibility loci in Koreans possessing the familial clustering of gastric cancer.

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.

Loss of Tbx3 in murine neural crest reduces enteric glia and causes cleft palate, but does not influence heart development or bowel transit

Transcription factors that coordinate migration, differentiation or proliferation of enteric nervous system (ENS) precursors are not well defined. To identify novel transcriptional regulators of ENS development, we performed microarray analysis at embryonic day (E) 17.5 and identified many genes that were enriched in the ENS compared to other bowel cells. We decided to investigate the T-box transcription factor Tbx3, which is prominently expressed in developing and mature ENS. Haploinsufficiency for TBX3 causes ulnar-mammary syndrome (UMS) in humans, a multi-organ system disorder. TBX3 also regulates several genes known to be important for ENS development. To test the hypothesis that Tbx3 is important for ENS development or function, we inactivated Tbx3 in all neural crest derivatives, including ENS progenitors using Wnt1-Cre and a floxed Tbx3 allele. Tbx3 fl/fl; Wnt1-Cre conditional mutant mice die shortly after birth with cleft palate and difficulty feeding. The ENS of mutants was well-organized with a normal density of enteric neurons and nerve fiber bundles, but small bowel glial cell density was reduced. Despite this, bowel motility appeared normal. Furthermore, although Tbx3 is expressed in cardiac neural crest, Tbx3 fl/fl; Wnt1-Cre mice had structurally normal hearts. Thus, loss of Tbx3 within neural crest has selective effects on Tbx3-expressing neural crest derivatives.

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.

Hirschsprung disease - integrating basic science and clinical medicine to improve outcomes

Hirschsprung disease is defined by the absence of enteric neurons at the end of the bowel. The enteric nervous system (ENS) is the intrinsic nervous system of the bowel and regulates most aspects of bowel function. When the ENS is missing, there are no neurally mediated propulsive motility patterns, and the bowel remains contracted, causing functional obstruction. Symptoms of Hirschsprung disease include constipation, vomiting, abdominal distension and growth failure. Untreated disease usually causes death in childhood because bloodstream bacterial infections occur in the context of bowel inflammation (enterocolitis) or bowel perforation. Current treatment is surgical resection of the bowel to remove or bypass regions where the ENS is missing, but many children have problems after surgery. Although the anatomy of Hirschsprung disease is simple, many clinical features remain enigmatic, and diagnosis and management remain challenging. For example, the age of presentation and the type of symptoms that occur vary dramatically among patients, even though every affected child has missing neurons in the distal bowel at birth. In this Review, basic science discoveries are linked to clinical manifestations of Hirschsprung disease, including partial penetrance, enterocolitis and genetics. Insights into disease mechanisms that might lead to new prevention, diagnostic and treatment strategies are described.

Sporadic Hirschsprung Disease: Mutational Spectrum and Novel Candidate Genes Revealed by Next-generation Sequencing

Hirschsprung disease (HSCR) is a common cause of functional colonic obstruction in children. The currently available genetic testing is often inadequate as it mainly focuses on RET and several other genes, accounting for only 15–20% of cases. To identify novel, potentially pathogenic variants, we isolated a panel of genes from a whole-exome sequencing study and from the published mouse aganglionosis phenotypes, enteric nervous system development, and a literature review. The coding exons of 172 genes were analyzed in 83 sporadic patients using next-generation sequencing. Rare stop-gain, splice-site variants, frameshift and in-frame insertions/deletions and non-synonymous variants (conserved and predicted to be deleterious) were prioritized as the most promising variants to have an effect on HSCR and subjected to burden analysis. GeneMANIA interaction database was used to identify protein–protein interaction-based networks. In addition, 6 genes (PTPN13, PHKB, AGL, ZFHX3, LAMA1, and AP3B2) were prioritized for follow-up studies: both their time-space expression patterns in mouse and human colon showed that they are good candidates for predicting pathogenicity. The results of this study broaden the mutational spectrum of HSCR candidate genes, and they provide an insight into the relative contributions of individual genes to this highly heterogeneous disorder.

Upregulation of MiR-369-3p suppresses cell migration and proliferation by targeting SOX4 in Hirschsprung's disease

BACKGROUND

Hirschsprung disease (HSCR) is a congenital digestive disease in the new born. miR-369-3p has been reported to be involved in many human diseases. However, the relationship between miR-369-3p and HSCR remains largely unknown.

METHODS

In this study, qRT-PCR was used to detect the relative expression of miR-369-3p in 60 HSCR bowel tissue samples and 47 matched controls. Bioinformatic analysis and dual-luciferase reporter assay were performed to evaluate the target for miR-369-3p. Cell Counting Kit-8 (CCK-8) assay, Transwell assay, wound healing assay and flow cytometry were employed to investigate the biological function of miR-369-3p in human SH-SY5Y and 293T cell lines.

RESULTS

We found that ganglion cell numbers were remarkably reduced while miR-369-3p was significantly upregulated in HSCR tissues compared to that in adjacent normal tissues (P<0.01). Dual-luciferase reporter assay showed that the 3'-UTR of SOX4 was a direct target to miR-369-3p. Moreover, an increased level of miR-369-3p was inversely correlated with decreased levels of SOX4 mRNA and protein (P<0.05, respectively). Dysregulation of miR-369-3p and SOX4 significantly suppressed cell proliferation and migration in SH-SY5Y and 293T cell lines in vitro (P<0.05, respectively).

CONCLUSION

Our study demonstrates that aberrant expression of miR-369-3p might play a crucial role in the development HSCR by regulating SOX4 expression, which may infer that it is an effective diagnostic target in the pathogenesis of HSCR, but investigation is still needed to explore the underlying mechanism.

Overexpression of DNMT3b target genes during Enteric Nervous System development contribute to the onset of Hirschsprung disease

Hirschsprung disease (HSCR) is attributed to a failure of neural crest cells (NCCs) to migrate, proliferate, differentiate and/or survive in the bowel wall during embryonic Enteric Nervous System (ENS) development. ENS formation is the result from a specific gene expression pattern regulated by epigenetic events, such DNA methylation by the DNA methyltransferases (DNMTs), among other mechanisms. Specifically, DNMT3b de novo methyltransferase is associated with NCCs development and has been shown to be implicated in ENS formation and in HSCR. Aiming to elucidate the specific mechanism underlying the DNMT3b role in such processes, we have performed a chromatin immunoprecipitation coupled with massively parallel sequencing analysis to identify the DNMT3B target genes in enteric precursor cells (EPCs) from mice. Moreover, the expression patterns of those target genes have been analyzed in human EPCs from HSCR patients in comparison with controls. Additionally, we have carried out a search of rare variants in those genes in a HSCR series. Through this approach we found 9 genes showing a significantly different expression level in both groups. Therefore, those genes may have a role in the proper human ENS formation and a failure in their expression pattern might contribute to this pathology.

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.