Machine Learning-Based Identification of Hub Genes and Temporal Regulation Mechanisms in Zebrafish Fin Regeneration

Background/Objectives: Zebrafish fin regeneration serves as a classic model for investigating vertebrate tissue regeneration, yet the core regulatory networks and their crosstalk with the immune microenvironment remain incompletely characterized. This study aimed to identify hub genes, and elucidate the underlying molecular mechanisms and immune microenvironment dynamics during zebrafish fin regeneration. Methods: We integrated multiple bulk RNA-seq datasets of zebrafish fin regeneration from the GEO database, followed by data standardization with batch effect removal. Hub genes were screened via differential expression analysis, weighted gene co-expression network analysis (WGCNA), and predictive models constructed with 13 classic machine learning algorithms. Functional enrichment, time-ordered gene co-expression network (TO-GCN) method, immune infiltration analyses and RT-qPCR validation were further performed. Results: We identified upregulated differentially expressed genes, regeneration-correlated gene modules and their overlapping genes, including 82 candidate genes and 10 hub genes enriched in cytoskeleton remodeling, extracellular matrix organization, and focal adhesion. Temporal analysis uncovered hierarchical gene regulation and functional switching during regeneration. Hub gene expression was significantly correlated with the infiltration of B cells, M1/M2 macrophages and CD8+ T cells, revealing a stage-specific immune microenvironment. RT-qPCR validation showed high consistency with the multi-omics data. Conclusions: This study provides potential gene targets for understanding zebrafish fin regeneration, and offers a valuable reference for investigating the crosstalk between regulatory networks and the immune microenvironment in vertebrate tissue regeneration.