Abstract
Physical exercise has been shown to induce epigenetic modifications with various health implications, directly affect DNA methylation (DNAm), as well as reverse the epigenetic age. Hence, we aimed to identify differential methylation changes and assess the epigenetic age in the saliva of 7-9-year-old school children following a 3-month integrated neuromuscular training (INT), as well as to explore if any of the methylation changes are in core genes. Core genes are defined as genes of high relevance and essential importance within the human genome. Forty children (17 boys and 23 girls) were recruited from schools in Girona, Spain, and allocated into control (N = 20) or INT (N = 20) group. The INT group performed a 3-month INT as a warm-up during the physical education (PE) classes, encompassing strength, coordination, dynamic stabilization, plyometrics, speed, and agility exercises, whereas the control group performed traditional warm-up activities, encompassing aerobic exercises that will prepare the cardiovascular system and increase the joint mobility for the upcoming effort during the class. Genome-wide DNAm analysis was performed with the Illumina 900 K microarray. Core genes were recognized based on the accomplishment of a rigorous and widely accepted 3-point criteria: participation in the enriched pathways, high connectivity (≥ 10), and target genes of key transcription factors. There were 1200 differentially methylated positions (DMPs) in the control group and 414 DMPs in the INT group (FDR < 0.05, p < 0.05, Aβ < |0.1|), suggesting a non-significant trend of epigenetic age acceleration in the control group (1.18 months, p > 0.05) and a non-significant 1-month decrease of the epigenetic age in the INT group (p > 0.05). The genes with DMPs in the control group showed low similarity between enriched pathways and low interconnectivity, encompassing distinct pathways, mostly development and growth-related. Additionally, no core genes were identified in the control group. Interestingly, the genes with DMPs in the INT group showed high similarity between enriched pathways and high interconnectivity, encompassing related pathways involving signaling mechanisms, as well as hormone and protein metabolism pathways. Moreover, 17 DMPs in the children from the INT group were in core genes. The main findings of the present study are suggesting an integrated response to the training stimulus in 7-9-year-old school children that performed a 3-month INT, including epigenetic modifications in 17 genes considered as core genes. Trial Registration: The study protocol was registered in the ISRCTN registry (ISRCTN16744821).