Time-series transcriptomics reveals distinctive mRNA expression dynamics associated with gene ontology specificity and protein expression in skeletal muscle after electrical stimulation-induced resistance exercise

Resistance exercise upregulates and downregulates the expression of a wide range of genes in skeletal muscle. However, detailed analysis of mRNA dynamics such as response rates and temporal patterns of the transcriptome after resistance exercise has not been performed. We aimed to clarify the dynamics of time-series transcriptomics after resistance exercise. We used electrical stimulation-induced muscle contraction as a resistance exercise model (5 sets x 10 times of 3 s of 100-Hz electrical stimulation) on the tibialis anterior muscle of rats and measured the transcriptome in the muscle before and at 0, 1, 3, 6, and 12 h after muscle contractions by RNA sequencing. We also examined the relationship between the parameters of mRNA dynamics and the increase in protein expression at 12 h after muscle contractions. We found that the function of the upregulated genes differed after muscle contractions depending on their response rate. Genes related to muscle differentiation and response to mechanical stimulus were enriched in the sustainedly upregulated genes. Furthermore, there was a positive correlation between the magnitude of upregulated mRNA expression and the corresponding protein expression level at 12 h after muscle contractions. Although it has been theoretically suggested, this study experimentally demonstrated that the magnitude of the mRNA response after electrical stimulation-induced resistance exercise contributes to skeletal muscle adaptation via increases in protein expression. These findings suggest that mRNA expression dynamics such as response rate, a sustained upregulated expression pattern, and the magnitude of the response contribute to mechanisms underlying adaptation to resistance exercise.