Gene network analysis of candidate loci for human anorectal malformations

RNA-Seq Profiling of Circular RNAs During Development of Hindgut in Rat Embryos With Ethylenethiourea-Induced Anorectal Malformations

Anorectal malformations (ARMs) are among the most common congenital terminal digestive tract malformations. Circular RNAs (circRNAs), a novel type of endogenous non-coding RNAs, play roles in the development of the digestive system; however, their contributions to the pathogenesis of ARMs are not well-established. In this study, we explored the mechanism underlying ethylenethiourea (ETU)-induced ARMs by profiling circRNA expression via RNA-seq and constructing a regulatory circRNA-miRNA-mRNA network. Nine pregnant rats were gavage-fed a single dose of 125 mg/kg 1% ETU (ARM group) on gestational day 10 (GD10), and another 9 pregnant rats received a similar dose of saline (normal group) as a control. Embryos were obtained by cesarean section on the key time-points of anorectal development (GD14, GD15, and GD16). Hindgut samples isolated from the fetuses were evaluated by high-throughput sequencing and differentially expressed circRNAs were validated by reverse transcription-quantitative polymerase chain reaction, agarose gel electrophoresis, and Sanger cloning and sequencing. A total of 18295 circRNAs were identified in the normal and ARM groups. Based on the 425 differentially expressed circRNAs (|Fc| > 2, p < 0.05), circRNA-miRNA and miRNA-mRNA pairs were predicted using miREAP, miRanda, and TargetScan. A total of 55 circRNAs (14 up- and 41 downregulated in the ARM group compared to the normal group) were predicted to bind to 195 miRNAs and 947 mRNAs. Competing endogenous RNA networks and a Kyoto Encyclopedia of Genes and Genomes analysis revealed that novel_circ_001042 had the greatest connectivity and was closely related to ARM-associated signaling pathways, such as the Wingless Type MMTV integration site family, mitogen-activated protein kinase, and transforming growth factor-β pathways. These results provide original insight into the roles of circRNAs in ARMs and provide a valuable resource for further analyses of molecular mechanisms and signaling networks.

Exome chip association study excluded the involvement of rare coding variants with large effect sizes in the etiology of anorectal malformations

INTRODUCTION

Anorectal malformations (ARM) are rare congenital malformations, resulting from disturbed hindgut development. A genetic etiology has been suggested, but evidence for the involvement of specific genes is scarce. We evaluated the contribution of rare and low-frequency coding variants in ARM etiology, assuming a multifactorial model.

METHODS

We analyzed 568 Caucasian ARM patients and 1,860 population-based controls using the Illumina HumanExome Beadchip array, which contains >240,000 rare and low-frequency coding variants. GenomeStudio clustering and calling was followed by re-calling of 'no-calls' using zCall for patients and controls simultaneously. Single variant and gene-based analyses were performed to identify statistically significant associations, applying Bonferroni correction. Following an extra quality control step, candidate variants were selected for validation using Sanger sequencing.

RESULTS

When we applied a MAF of ≥1.0%, no variants or genes showed statistically significant associations with ARM. Using a MAF cut-off at 0.4%, 13 variants initially reached statistical significance, but had to be discarded upon further inspection: ten variants represented calling errors of the software, while the minor alleles of the remaining three variants were not confirmed by Sanger sequencing.

CONCLUSION

Our results show that rare and low-frequency coding variants with large effect sizes, present on the exome chip do not contribute to ARM etiology.

Integrative analysis of microRNA and mRNA expression profiles in fetal rat model with anorectal malformation

Background

Anorectal malformations (ARMs) are the most common congenital malformation of the gut, and regulated by multiple signal transduction pathways. The microRNA (miRNA) expression profiles and their biologial functions in anorectal malformations (ARMs) remain unclear. The aim of our study was to evaluate miRNA and mRNA expression profiles in the ARM rats.

Methods and Materials

ARM was induced with ethylenethiourea (ETU) on gestational day 10. Cesarean deliveries were performed to harvest the embryos on gestional day 20. For the extraction of total RNA, 1 cm terminal hindgut samples were obtained from three fetal rats that had similar weights. The microarrays and quantitative RT-PCR analysis were conducted to evaluate the miRNA and mRNA expression profiles in normal fetal rats (n = 3) and ARM fetal rats (n = 3).

Results

In total, 33 miRNAs and 772 mRNAs were significantly and differentially expressed in terminal hindgut tissues of ARM group versus control group. Functional annotation was performed to understand the functions and pathways of differentially expressed mRNAs. Also, we constructed a miRNA-target gene regulatory network including 25 differentially expressed miRNAs and 76 mRNAs. Furthermore, the credibility of the microarray-based results were validated by using qRT-PCR.

Conclusion

The miRNA and mRNA expression in terminal hindgut tissue of ARM fetal rats might provide a basis for further research on the pathogenesis of ARMs.

MSX2 and BCL2 expressions in the development of anorectal malformations in ethylenethiourea-induced rat embryos

BACKGROUND

This study aimed to determine Msh homeobox 2 (MSX2) and B cell lymphoma-2 (BCL2) expression patterns during anorectal development in anorectal malformations (ARM) and normal rat embryos, with the goals of determining the role of MSX2 and BCL2 in ARM pathogenesis.

METHODS

ARM was induced in rat embryos with ethylenethiourea administered to dams on gestational day 10 (GD10). Embryos were harvested by cesarean deliveries from GD14 to GD16. MSX2 and BCL2 expression was evaluated via immunohistochemical staining, immunofluorescence, western blotting and quantitative real-time polymerase chain reaction (qRT-PCR).

RESULTS

Immunohistochemical staining of ARM embryos revealed that MSX2 was mainly expressed in the epithelium of the hindgut and urorectal septum (URS) on GD14. On GD15 and GD16, MSX2-immunolabeled cells were noted in the epithelium of the rectum, fistula and URS. However, in normal embryos, faint immunopositivity for MSX2 was demonstrated in the epithelium of the rectum and URS from GD14 to GD16. As for BCL2, in normal embryos, BCL2-immunopositive cells were extensively expressed in the epithelium of the hindgut and URS on GD14 and GD15. In ARM embryos, weak immunopositivity for BCL2 was detected in the epithelium of hindgut and URS on GD14 and GD15. Immunofluorescence revealed that MSX2 and BCL2 colocalized in the hindgut. In ARM embryos, we observed more MSX2-positive than BCL2-positive cells on GD14; the normal embryos had the opposite pattern. Analyses by western blot and qRT-PCR showed that MSX2 protein and mRNA expression was significantly increased in ARM embryos compared with the normal embryos on GD15 and GD16 (p < 0.05). However, BCL2 protein and mRNA expression was significantly decreased in ARM embryos compared with the normal embryos on GD14 (p < 0.05). The MSX2/BCL2 ratio of protein and mRNA expression level in the ARM group was the highest on GD15.

CONCLUSION

These results indicate that upregulation of MSX2 and downregulation of BCL2 during cloacal development into the rectum and urethra might be related to the ARM development, and MSX2 promoted apoptosis through reduction of BCL2 expression during the development of anorectal development in ARM.

Integrating lncRNAs and mRNAs expression profiles in terminal hindgut of fetal rats with anorectal malformations

BACKGROUND

The detailed embryonic etiology and pathogenesis of anorectal malformations (ARMs) remains unclear. Recent studies have shown that gene expression abnormalities were the key factors that result in ARMs. Long non-coding RNAs (lncRNAs) were reported as the 'transcriptional noise' within the genome. The expression profiles of lncRNA and mRNA remain less characterized in the pathogenesis of ARMs. Furthermore, the function of lncRNAs in the regulation of this process has not been investigated so far. Therefore, this current study was aimed to integrate lncRNA and mRNA expression profiles in terminal hindgut of ethylenethiourea (ETU)-induced ARM rats using Agilents lncRNA and mRNA co-expression microarrays.

METHODS

ARM model was induced with ethylenethiourea (ETU) on gestational day 10. Cesarean deliveries were conducted to collect the embryos on gestational day 20. For the extraction of total RNA, 1-cm terminal hindgut tissues were collected from three fetal rats with similair weights. The microarrays and quantitative RT-PCR analysis were conducted to evaluate the lncRNA and mRNA expression profiles in normal fetal rats and ARM fetal rats.

RESULTS

Compared with control group, 164 lncRNAs were observed to be aberrantly expressed (FC ≥ 2; P < 0.05) in ARM group, including 36 upregulated and 128 downregulated, while 772 mRNAs were observed to be aberrantly expressed (FC ≥ 2; P < 0.05) in the terminal hindgut, including 350 up-regulated and 422 down-regulated. The differential expression profiles between the ARM and the control group were used for gene ontology (GO) and pathway analysis. A subset of those RNAs was identified to be closely related to the development process of ARMs. The four RNAs that were differentially expressed between the two groups were selected for qPCR validation, and the results were in line with the microarray data. In addition, the lncRNAs and mRNA co-expression network was established according to the correlation analysis. We predicted the functions of transregulatory lncRNAs by the TFs (transcription factors) which might modulate their expression. In the core network of lncRNA-TF pairs, the lncRNAs can be classified into 5 categories of pathways governed by Jun, c-Myc, Usf1, Alf2, and Stat3.

CONCLUSION

From the above results, it can be suggested that these aberrant lncRNAs might participate in the pathogenesis of ARM, and our present work may provide new research directions for future studies of ARMs.

Microarray analysis of miRNAs during hindgut development in rat embryos with ethylenethiourea‑induced anorectal malformations

Anorectal malformations (ARMs) are one of the most common congenital malformations of the digestive tract; however, the pathogenesis of this disease remains to be fully elucidated. MicroRNAs (miRNAs) are important in gastrointestinal development and may be involved in the pathogenesis of ARMs. The present study aimed to profile miRNAs and examine their potential functions in rats with ethylenethiourea (ETU)-induced ARMs. Pregnant Wistar rats (n=36) were divided randomly into ETU-treated and control groups. The rats in the ETU-treated group were gavage-fed 1% ETU (125 mg/kg) on gestational day 10 (GD10), whereas the control group rats received a corresponding dose of saline. Embryos were harvested by cesarean section on GD14, GD15 and GD16. Hindgut tissue was isolated from the fetuses for RNA extraction and microarray analysis, followed by bioinformatics analysis and reverse transcription-quantitative polymerase chain reaction (RT-qPCR) validation. Overall, 38 miRNAs were differentially expressed (all upregulated) on GD14, 49 (32 upregulated and 17 downregulated) on GD15, and 42 (all upregulated) on GD16 in the ARM group compared with the normal group. The top 18 miRNAs with |log₂(fold change)| >4.25 were selected for further bioinformatics analysis. Among these miRNAs, five were differentially expressed at two time-points and were involved in ARM-associated signaling pathways. The RT-qPCR analysis revealed that three miRNA (miR), miR-125b-2-3p, miR-92a-2-5p and miR-99a-5p, were significantly differentially expressed in rats with ARMs compared with the normal group. In conclusion, the results suggested that the differential expression of miR-125b-2-3p, miR-92a-2-5p and miR-99a-5p during key time-points of anorectal formation in rats may have functions in the pathogenesis of ARM.

A review of genetic factors contributing to the etiopathogenesis of anorectal malformations

BACKGROUND

Anorectal malformation (ARM) is a common congenital anomaly with a wide clinical spectrum. Recently, many genetic and molecular studies have been conducted worldwide highlighting the contribution of genetic factors in its etiology. We summarize the current literature on such genetic factors.

MATERIALS AND METHODS

Literature search was done using different combinations of terms related to genetics in anorectal malformations. From 2012 to June 2017, articles published in the English literature and studies conducted on human population were included.

OBSERVATIONS AND RESULTS

A paradigm shift was observed from the earlier studies concentrating on genetic aberrations in specific pathways to genome wide arrays exploring single nucleotide polymorphisms (SNPs) and copy number variations (CNVs) in ARM patients. Rare CNVs (including 79 genes) and SNPs have been found to genetically contribute to ARM. Out of disrupted 79 genes one such putative gene is DKK4. Down regulation of CDX-1 gene has also been implicated in isolated ARM patients. In syndromic ARM de novo microdeletion at 17q12 and a few others have been identified.

CONCLUSION

Major genetic aberrations proposed in the pathogenesis of ARM affect members of the Wnt, Hox (homebox) genes, Sonic hedgehog (Shh) and Gli2, Bmp4, Fgf and CDX1 signalling pathways; probable targets of future molecular gene therapy.

Genetic and nongenetic etiology of nonsyndromic anorectal malformations: a systematic review

Congenital anorectal malformations (ARMs) are one of the most frequently observed birth defects of the digestive system. However, their etiology remains elusive. Therefore, we aim to summarize and critically appraise all existing literature on the genetic and nongenetic etiology of nonsyndromic ARM and to conclude with unifying hypotheses and directions for future research. A structured literature search on English language human studies was conducted in PubMed and Embase up to October 1, 2013, resulting in 112 included articles. Research on the identification of genes underlying nonsyndromic ARM is remarkably scarce. Most studies were focused on screening of candidate genes for mutations or single-nucleotide polymorphisms, which did not yield any substantial evidence. Nongenetic factors fairly consistently found to be associated with ARM are assisted reproductive techniques, multiple pregnancy, preterm delivery, low birth weight, maternal overweight or obesity, and preexisting diabetes. This review provides indications for the involvement of both genes and nongenetic risk factors in the etiology of ARM. In future studies, large cohorts of patients with ARM from national and international collaborations are needed to acquire new hypotheses and knowledge through hypothesis-generating approaches. Challenges for future studies may also lie in the investigation of gene-gene and gene-environment interactions.

Anorectal malformation: the etiological factors

Anorectal malformation (ARM) is a congenital anomaly commonly encountered in pediatric surgery practice. Although surgical procedures correct the anatomical anomalies, the post-operative bowel function is not universally satisfactory. The etiology of ARM remains unclear. In this review, we summarize the current understanding of the genetic and epigenetic factors contributing to the pathogenesis of ARM, based on published animal models, human genetics and epidemiological researches. Appreciation of these factors may be helpful in the management of ARM in the future.

The Great Divide: Understanding Cloacal Septation, Malformation, and Implications for Surgeons

The anorectal and urogenital systems arise from a common embryonic structure termed cloaca. Subsequent development leads to the division/septation of the cloaca into the urethra, urinary bladder, vagina, anal canal, and rectum. Defective cloacal development and the resulting anorectal and urogenital malformations are some of the most severe congenital anomalies encountered in children. In the most severe form in females, the rectum, vagina, and urethra fail to develop separately and drain via a single common channel known as a cloaca into the perineum. In this review, we summarize our current knowledge of embryonic cloaca development and malformation, and compare them to what has already been described in the literature. We describe the use of mouse models of cloaca malformation to understand which signaling pathways and cellular mechanisms are involved in the process of normal cloaca development. We also discuss the embryological correlation of the epithelial and stromal histology found in step sections of the common channel in 14 human cloaca malformations. Finally, we highlight the significance of these findings, compare them to prior studies, and discuss their implications for the pediatric surgeons. Understanding and identifying the molecular basis for cloaca malformation could provide foundation for tissue engineering efforts that in the future would reflect better surgical reconstruction and improved quality of life for patients.

Defining the molecular pathologies in cloaca malformation: similarities between mouse and human

Anorectal malformations are congenital anomalies that form a spectrum of disorders, from the most benign type with excellent functional prognosis, to very complex, such as cloaca malformation in females in which the rectum, vagina and urethra fail to develop separately and instead drain via a single common channel into the perineum. The severity of this phenotype suggests that the defect occurs in the early stages of embryonic development of the organs derived from the cloaca. Owing to the inability to directly investigate human embryonic cloaca development, current research has relied on the use of mouse models of anorectal malformations. However, even studies of mouse embryos lack analysis of the earliest stages of cloaca patterning and morphogenesis. Here we compared human and mouse cloaca development and retrospectively identified that early mis-patterning of the embryonic cloaca might underlie the most severe forms of anorectal malformation in humans. In mouse, we identified that defective sonic hedgehog (Shh) signaling results in early dorsal-ventral epithelial abnormalities prior to the reported defects in septation. This is manifested by the absence of Sox2 and aberrant expression of keratins in the embryonic cloaca of Shh knockout mice. Shh knockout embryos additionally develop a hypervascular stroma, which is defective in BMP signaling. These epithelial and stromal defects persist later, creating an indeterminate epithelium with molecular alterations in the common channel. We then used these animals to perform a broad comparison with patients with mild-to-severe forms of anorectal malformations including cloaca malformation. We found striking parallels with the Shh mouse model, including nearly identical defective molecular identity of the epithelium and surrounding stroma. Our work strongly suggests that early embryonic cloacal epithelial differentiation defects might be the underlying cause of severe forms of anorectal malformations in humans. Moreover, deranged Shh and BMP signaling is correlated with severe anorectal malformations in both mouse and humans.