Integrated bioinformatics analysis of potential pathway biomarkers using abnormal proteins in clubfoot

Background: As one of the most common major congenital distal skeletal abnormalities, congenital talipes equinovarus (clubfoot) affects approximately one in one thousandth newborns. Although several etiologies of clubfoot have been proposed and several genes have been identified as susceptible genes, previous studies did not further explore signaling pathways and potential upstream and downstream regulatory networks. Therefore, the aim of the present investigation is to explore abnormal pathways and their interactions in clubfoot using integrated bioinformatics analyses. Methods: KEGG, gene ontology (GO), Reactome (REAC), WikiPathways (WP) or human phenotype ontology (HP) enrichment analysis were performed using WebGestalt, g:Profiler and NetworkAnalyst. Results: A large number of signaling pathways were enriched e.g. signal transduction, disease, metabolism, gene expression (transcription), immune system, developmental biology, cell cycle, and ECM. Protein-protein interactions (PPIs) and gene regulatory networks (GRNs) analysis results indicated that extensive and complex interactions occur in these proteins, enrichment pathways, and TF-miRNA coregulatory networks. Transcription factors such as SOX9, CTNNB1, GLI3, FHL2, TGFBI and HOXD13, regulated these candidate proteins. Conclusion: The results of the present study supported previously proposed hypotheses, such as ECM, genetic, muscle, neurological, skeletal, and vascular abnormalities. More importantly, the enrichment results also indicated cellular or immune responses to external stimuli, and abnormal molecular transport or metabolism may be new potential etiological mechanisms of clubfoot.