Integrated Bioinformatics Analysis Reveals Potential Pathway Biomarkers and Their Interactions for Clubfoot

BACKGROUND Congenital talipes equinovarus (clubfoot), one of the most regular pediatric congenital skeletal anomalies, seriously affects the normal growth and development of about 1 in 1000 newborns. Although it has been investigated widely, the etiology and pathogenesis of clubfoot are still controversial. MATERIAL AND METHODS g: Profiler, NetworkAnalyst and WebGestalt were used to probe the enriched signaling pathways by using the Gene Ontology (GO), Human Phenotype Ontology (HP), Kyoto Encyclopedia of Genes and Genomes (KEGG), Reactome (REAC), and WikiPathways (WP) databases. Large numbers of enriched signaling pathways were identified using the integrated bioinformatics enrichment analyses. RESULTS Apoptosis or programmed cell death (PCD), disease, muscle contraction, metabolism, and immune system were the top functions. Embryo or organ morphogenesis and development, cell or muscle contraction, and apoptosis were the top biological processes, and cell/muscle contraction and apoptosis were the top molecular functions using enriched GO terms analysis. There were a large number of complex interactions in the genes, enriched pathways, and transcription factor (TF)-miRNA co-regulatory networks. Transcription factors such as FOXN3, GLI3, HOX, and NCOR2 family regulated the gene expression of APAF1, BCL2, BID, CASP, MTHFR, and TPM family. CONCLUSIONS The results of bioinformatics enrichment analysis not only supported the previously proposed hypotheses, e.g., extracellular matrix (ECM) abnormality, fetal movement reducing, genetic abnormality, muscle abnormality, neurological abnormality, skeletal abnormality and vascular abnormality, but also indicated that cellular or immune responses to external stimulus, molecular transport and metabolism may be new etiological mechanisms in clubfoot.