Integrative Analysis of Axolotl Gene Expression Data from Regenerative and Wound Healing Limb Tissues

Axolotl (Ambystoma mexicanum) is a urodele amphibian endowed with remarkable regenerative capacities manifested in scarless wound healing and restoration of amputated limbs, which makes it a powerful experimental model for regenerative biology and medicine. Previous studies have utilized microarrays and RNA-Seq technologies for detecting differentially expressed (DE) genes in different phases of the axolotl limb regeneration. However, sufficient consistency may be lacking due to statistical limitations arising from intra-laboratory analyses. This study aims to bridge such gaps by performing an integrative analysis of publicly available microarray and RNA-Seq data from axolotl limb samples having comparable study designs using the “merging” method. A total of 351 genes were found DE in regenerative samples compared to the control in data of both technologies, showing an adjusted p-value < 0.01 and log fold change magnitudes >1. Downstream analyses illustrated consistent correlations of the directionality of DE genes within and between data of both technologies, as well as concordance with the literature on regeneration related biological processes. qRT-PCR analysis validated the observed expression level differences of five of the top DE genes. Future studies may benefit from the utilized concept and approach for enhanced statistical power and robust discovery of biomarkers of regeneration.