Next Article in Journal
APOE Genotypes Modulate Inflammation Independently of Their Effect on Lipid Metabolism
Next Article in Special Issue
Space Omics and Tissue Response in Astronaut Skeletal Muscle after Short and Long Duration Missions
Previous Article in Journal
Photocatalytic Activity of TiO2 Coatings Obtained at Room Temperature on a Polymethyl Methacrylate Substrate
Previous Article in Special Issue
Differential Analysis of Gly211Val and Gly286Val Mutations Affecting Sarco(endo)plasmic Reticulum Ca2+-ATPase (SERCA1) in Congenital Pseudomyotonia Romagnola Cattle
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
IJMS
Volume 23
Issue 21
10.3390/ijms232112933
ijms-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Review

Genes Whose Gain or Loss of Function Changes Type 1, 2A, 2X, or 2B Muscle Fibre Proportions in Mice—A Systematic Review

by
Gabryela Kuhnen
1,
Tiago Guedes Russomanno
1,
Marta Murgia
2,3,
Nicolas J. Pillon
4,
Martin Schönfelder
1 and
Henning Wackerhage
1,*
1
Department of Sports and Health Sciences, Technical University of Munich, 80809 Munich, Germany
2
Max Planck Institute, Martinsried, 82152 Munich, Germany
3
Department of Biomedical Sciences, University of Padova, Via Ugo Bassi, 58/B, 35131 Padua, Italy
4
Department of Physiology and Pharmacology, Karolinska Institutet, 171 77 Stockholm, Sweden
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2022, 23(21), 12933; https://doi.org/10.3390/ijms232112933
Submission received: 20 September 2022 / Revised: 16 October 2022 / Accepted: 18 October 2022 / Published: 26 October 2022
(This article belongs to the Special Issue Novel Molecular Approaches to Skeletal Muscle Disease and Disuse 2.0)
Browse Figures
Versions Notes

Abstract

:
Adult skeletal muscle fibres are classified as type 1, 2A, 2X, and 2B. These classifications are based on the expression of the dominant myosin heavy chain isoform. Muscle fibre-specific gene expression and proportions of muscle fibre types change during development and in response to exercise, chronic electrical stimulation, or inactivity. To identify genes whose gain or loss-of-function alters type 1, 2A, 2X, or 2B muscle fibre proportions in mice, we conducted a systematic review of transgenic mouse studies. The systematic review was conducted in accordance with the 2009 PRISMA guidelines and the PICO framework. We identified 25 “muscle fibre genes” (Akirin1, Bdkrb2, Bdnf, Camk4, Ccnd3, Cpt1a, Epas1, Esrrg, Foxj3, Foxo1, Il15, Mapk12, Mstn, Myod1, Ncor1, Nfatc1, Nol3, Ppargc1a, Ppargc1b, Sirt1, Sirt3, Thra, Thrb, Trib3, and Vgll2) whose gain or loss-of-function significantly changes type 1, 2A, 2X or 2B muscle fibre proportions in mice. The fact that 15 of the 25 muscle fibre genes are transcriptional regulators suggests that muscle fibre-specific gene expression is primarily regulated transcriptionally. A reanalysis of existing datasets revealed that the expression of Ppargc1a and Vgll2 increases and Mstn decreases after exercise, respectively. This suggests that these genes help to regulate the muscle fibre adaptation to exercise. Finally, there are many known DNA sequence variants of muscle fibre genes. It seems likely that such DNA sequence variants contribute to the large variation of muscle fibre type proportions in the human population.

1. Introduction

Human muscle fibres are up to 20 cm-long cells that produce force and heat [1]. Human muscle fibres are multinucleated cells, also termed syncytia, that develop because of myoblast fusion. Per millimetre of length, human muscle fibres have ≈50–250 myonuclei [2] and so we would expect 10,000–50,000 nuclei in a single, 20 cm-long human muscle fibre. Each skeletal muscle contains from hundreds up to many thousands of muscle fibres depending on its size. For example, the human vastus lateralis of ≈20-year-old males can contain between 393,000 and 903,000 muscle fibres [3].
The identification of different adult skeletal muscle fibres has evolved over time. From 1960–1967, researchers distinguished in-between fast white and slow red muscle fibres; from ≈1967–1975 type 1, 2A, and 2B fibres; and from ≈1986–1991 type 1, 2A, 2X, and 2B muscle fibres [4]. Different muscle fibres were first identified using enzyme assays [5,6] as well as histochemical and microscopic visualisation of ATPase activity after pre-incubation with an acidic or alkaline pH [7]. Myosin heavy chain isoforms (MHC), which determine the contraction velocity, were then visualised by immunocytochemistry utilising myosin heavy chain type 1, 2A, 2X, and 2B-specific antibodies [8], or, alternatively, with electrophoretic separation of myosin heavy chain isoforms [9]. Today, the concentrations of thousands of proteins in a single muscle fibre can be measured by first isolating single muscle fibres followed by an unbiased, proteomic analysis via mass spectrometry [10,11]. Currently, muscle fibres are primarily classified by the myosin heavy chain isoform(s) that they express. The major myosin heavy chain isoforms in humans are slow type 1 (gene Myh7), intermediate type 2A (gene Myh2), and fast 2X (gene Myh1) myosin [8]. Additionally, some rodent muscles—but not human muscles—contain the very fast type 2B myosin heavy chain protein (gene Myh4) [12]. In addition to fibres that exclusively (i.e., >90%) express one myosin heavy chain, there are so-called hybrid fibres that express two myosin heavy chain isoforms [13]. Moreover, embryonal (gene Myh3) and perinatal myosin (gene Myh8) heavy chains are expressed in embryonic, foetal, and regenerating muscle fibres [14].
In mammals, the distribution of muscle fibres varies greatly both intra- and inter-individually. For example, the red soleus is predominantly a type 1 slow twitch muscle, whereas the “whiter” rectus femoris contains faster type 2 muscle fibres [15]. In addition to this inter-muscle variability, muscle fibre percentages vary greatly between individuals. For example, an analysis of 418 human vastus lateralis biopsies reported 15–79% type 1 fibres, 13–77% type 2A fibres, and 0–44% type 2X fibres in the vastus lateralis and similarly large variations of metabolic enzyme activities [16]. Extreme muscle fibre compositions occur in athletes, with endurance athletes having typically a high percentage of type 1 and sprint and power athletes having a high percentage of type 2 fibres pin key locomotory muscles such as the vastus lateralis or gastrocnemius [17,18,19,20]. Extreme fibre type percentages are prerequisites for elite performance in both power/speed and endurance sports.
The distribution of skeletal muscle fibres depends both on an individual’s genetics (i.e., variations of the DNA sequence, or heritability) and environmental factors such as exercise training or diet as well as experimental variability, e.g., genetic manipulation on mice. Regarding humans, Bouchard and colleagues estimated that the variation in the proportion of type 1 fibres depends 45% on genetics, 40% on environmental factors, and 15% on experimental variability [21]. In addition to genetics, more innervation and contraction e.g., due to endurance training or fewer contractions e.g., due to denervation, reinnervation, or immobilisation alter fibre type-related gene expression and muscle fibre proportions. Specifically, exercise training promotes gene expression changes and minor fibre type shifts in a fast-to-slow direction. In their classic study, Gollnick, Saltin, and colleagues found that 5 months with four 1 h training sessions per week increased the percentage of slow twitch (i.e., type 1) fibres in the vastus lateralis non-significantly from 32% to 36% [22]. In the Heritage study, twenty weeks of endurance training increased the percentage of type 1 fibres by 3.5% and decreased the percentage of type 2X fibres by 5.4% [23]. Analysing the whole body of evidence suggests that exercise training over months can convert hybrid fibres (e.g., type 2X/2A) to pure (e.g., 2A) fibres and promote some (i.e., <10%) pure fibre conversions, e.g., from pure 2X to pure 2A fibres [24].
Whilst exercise interventions of less than a year result only in limited fibre transitions, it is unclear whether years of exercise training can cause major fibre type transitions where, e.g., >10% of fibres shift their type within a muscle. In contrast to exercise training, near complete fibre type transformations are caused by either denervation where fibres change in a slow-to-fast direction [25], innervation by a different motor neuron [26], or by chronic electrical low-frequency stimulation where many or almost all fibres change in a fast-to-slow direction [27,28].
Given that approximately half of the variation in human muscle fibre proportions is inherited, the question arises as to what genes and DNA sequence variants influence the proportions of muscle fibre types within a given muscle. Transgenic mouse models with germline mutations or injection of genetic constructs have helped to identify genes whose gain or loss of function significantly alters the proportions of type 1, 2A, 2X, or 2B muscle fibres. An important early study used gene expression, reporter assays, and pharmacological inhibition to identify the calcineurin-Nfat (nuclear factor of activated T-cells) pathway as a regulator of muscle fibre proportions [29,30]. Other transgenic mouse studies around the millennium showed that the overexpression of constitutively active Ras promotes a fast-to-slow fibre type change [31], and that the gain of function of the transcriptional co-activator Pgc-1α in muscle not only promotes mitochondrial biogenesis but also increases the number of type 1 fibres [32].
To date, the International Mouse Phenotyping Consortium (IMPC) has created and phenotyped over 5000 transgenic mouse strains [33]. However, the IMPC does not determine the muscle fibre composition of their mice, leaving a knowledge gap over what genes influence fibre type proportions in what muscles and by how much.
To address this question, we systematically searched the literature to identify genes whose gain or loss of function significantly alters the proportion of type 1, 2A, 2X, or 2B muscle fibres in mice. For simplicity, we refer to these genes, and to the proteins that they encode, from here onwards as “muscle fibre genes”. As the second step, we used online databases and re-analysed published datasets to answer questions such as: “In what tissues are muscle fibre genes expressed?” and “Does exercise alter the expression of muscle fibre genes or the phosphorylation of the proteins that they encode?”

2. Methods

2.1. Literature Search

We conducted a systematic review to identify genes whose gain or loss of function results in a statistically significant change in the percentage of at least one muscle fibre type or a significant change in myosin heavy chain isoform protein abundance in a mouse muscle. For the systematic literature review, we followed the 2009 preferred reporting items for systematic reviews and meta-analysis guideline (PRISMA [34]). The PICO framework [35] was also used to choose the key search word list, which was formulated as follows: population/problem; intervention; comparison; and outcome. Based on this, we used the following MeSH (medical subjects headings) listed in Table 1.
We used this search strategy to search the PubMed–Medline database for studies published up until November 2020. PubMed–Medline was used as the standard database due to its ability to match the full text of manuscripts with advance researching filtering and specific tools.

2.2. Inclusion and Exclusion Criteria

After the literature search, we screened all titles and abstracts to remove studies that did not report the effects of transgenesis on muscle fibre proportions in mice. After that, we selected studies from peer-reviewed journals, written in English, that reported muscle fibre type distribution and/or myosin heavy chain expression in gene-manipulated mouse models. We included studies that reported a measure of skeletal muscle fibre type numbers, percentages, or quantified myosin heavy chain isoform expression. We excluded studies as follows: (1) rat or in vitro study; (2) no transgenesis or double mutation; (3) miRNA manipulation; (4) mice with disease or pathological changes; (5) no statistically significant effect on muscle fibre distribution or myosin heavy chain isoform expression; (6) no comparison to wildtype or any other valid control; and (7) the gene manipulation resulted in a disease.

2.3. Data Collection, Extraction, and Analyses

From all included studies we extracted the following information: author; gene name; protein name; method of gene manipulation; animal acclimation; output measure; section of paper from where the data were extracted (figure or table); muscle(s) studied; number of fibre type or MHC expression in transgenic and control mice; difference between transgenic mice and control mice in percentage (this analysis was performed in Rstudio version 3.6.3, then exported to an excel sheet for further analysis); age of mice; mouse strain; additional measurements; and remarks. Additional columns for PubMed identificatory (PMID) paper, official gene name by the universal protein resource (UniProt, NCBI) and aliases regarding mice and humans. When the output measured appears only in bar graph or in stained muscle sections, we estimated the relative difference between transgenic and control group using ImageJ Software [36], and indicated this with an (*) on worksheet S2 in the supplementary Table S1.

2.4. Bioinformatic Analyses

After identifying muscle fibre genes, we conducted five bioinformatic analyses to add more information on the tissue of expression, common functions, and properties and their regulation by exercise, etc. More specifically the analyses investigated the following:
(1)
To find out whether muscle fibre proteins are encoded we performed a String analysis (https://string-db.org/; RRID:SCR_005223 [37]).
(2)
To identify common functions and properties of muscle fibre genes, we performed a ToppGene enrichment analysis (https://toppgene.cchmc.org/ [38]).
(3)
To determine whether the skeletal muscle distribution and MHC expression genes are expressed specifically in skeletal muscle or elsewhere, we retrieved expression figures from the Genotype-Tissue Expression (GTEx; RRID:SCR_001618 [39]).
(4)
To identify associations between muscle fibre genes and human phenotypes we searched the Genome-wide association studies (GWAS) catalogue (https://www.ebi.ac.uk/gwas/; RRID:SCR_012745).
(5)
To find out whether muscle fibre genes change their expression after acute endurance exercise, resistance exercise, or in response to activity, we used the Meta-analysis of skeletal muscle response to exercise (MetaMEx) gene expression database to determine expression changes in muscle biopsies after acute endurance exercise, acute resistance exercise, and activity in health subjects (https://www.metamex.eu [40]). We also investigated whether muscle fibre proteins become phosphorylated or dephosphorylated after exercise. For this, we retrieved supplementary data from two phospho-proteome studies. Study 1 investigated protein phosphorylation changes after a single bout of high intensity training in human muscles [41]. Study 2 investigated protein phosphorylation in mouse skeletal muscle after electrically evoked maximal-intensity contractions [42].
More complete data of the above analyses are in the supplementary data.

3. Results and Discussion

3.1. Results

The initial PubMed–Medline search yielded 371 studies published before November 2020. We identified another three studies through other sources. After reading titles and abstracts, we excluded 261 manuscripts and 113 articles remained. We then read the full text of these articles and excluded more articles based on our inclusion and exclusion criteria. In the end, 24 publications were analysed quantitatively. The PRISMA flowchart in Figure 1 summarises the search and subsequent selection of the publications analysed in this study.

Genes Whose Gain or Loss of Function Significantly Changes Muscle Fibre Distribution in Mice

Overall, we identified 25 genes whose gain or loss of function significantly changed the percentages of type 1, 2A, 2X, or 2B muscle fibres or myosin heavy chain abundance in at least one muscle in mice. A summary of the studies included, and some of their main findings, is given in Table 2.
Genes whose gain or loss of function significantly changed the percentage of type 1, 2A, or 2X fibres are presented in Figure 2 and the genes that affect myosin heavy chain isoform expression are shown in Figure 3.
Specifically, we identified 13 genes (Bdkrb2, Bdnf, Camk4, Ccnd3, Cpt1a, Foxj3, Mapk12, Mstn, Myod1, Nfatc1, Nol3, Thra, and Thrb) whose loss of and two genes whose gain of function (Foxo1, Ppargc1a) significantly changed the proportions of at least one muscle fibre type. Additionally, we identified one gene whose knockout (Ncor1) and five genes whose gain of function (Esrrg, Il15, Ppargc1b, Sirt3, and Trib3) significantly altered the expression of at least one myosin heavy chain isoform in mice. Moreover, loss of function of two genes (Epas1 and Vgll2) and the gain of function of two other genes (Akirin1 and Sirt1) significantly altered both the proportions of at least one muscle fibre type and significantly changed the expression of at least one myosin heavy chain isoform. The effect sizes of the gene manipulations range from a reduction in a fibre type by 37% (Ppargc1a, Plantaris type 2B fibres) to a gain of 28% of a single fibre type (Epas1, soleus type 2B fibres). Genes whose knockout affected more than one fibre type in the soleus are Bdkrb2, Camk4, Mpak12, Nol3, Thra, and Thrb. The muscle fibre genes Bdnf, Ccnd3, and Mstn affected muscle fibre distribution in the tibialis anterior, Bdnf in the extensor digitorum longus, and Nol3 in the plantaris. In addition, Foxo1 affected more than one fibre type in the soleus, and, similarly, Ppargc1a in the plantaris. Ncor1 knockout affected the myosin heavy chain expression in the gastrocnemius and quadriceps. On the other hand, overexpression of Esrrg, Il15, Ppargc1b, and Sirt3 changed muscle fibre proportions in the gastrocnemius; Il15 and Trib3 in the soleus; Il15 and Ppargc1b in the extensor digitorum longus; and Ppargc1b and Trib3 in the tibialis anterior.
Next, we used the list of 25 muscle fibre genes to answer direct research questions through further bioinformatic analyses. Information on the searches and detailed results can be found in the supplemental data.
(1)
Do muscle fibre proteins interact, and do muscle fibre genes share common functional features?
To answer this question, we completed a String protein interaction analysis [37] and a ToppGene enrichment analysis [38]. The String analysis suggests interactions in and between muscle fibre genes. Clusters of muscle fibre genes included a cluster of genes that encoded thyroid (i.e., the expression of all MYH genes that respond to the thyroid hormone [66]) and oestrogen hormone receptors (Thra, Thrb, and Essrg). Further, there were also a cluster with the transcriptional co-factors and transcription factors Ppargc1a, Ppargc1b, Vgll2, Foxo1, Myod1, and Nfatc1; the sirtuins Sirt1 and Sirt3; and a cluster of the circulating factors Mstn and Bdnf, as well as the kinases Mapk12 and Camk4 (Figure 4).
We also used a ToppGene functional enrichment analysis to identify common features and functions among the muscle fibre genes identified. Specifically, we found 15 muscle fibre genes that regulate gene transcription (Foxo1, Ppargc1a, Ncor1, Ppargc1b, Thra, Thrb, Akirin1, Trib3, Nfatc1, Foxj3, Vgll2, Myod1, Epas1, Sirt1, and Esrrg), 5 genes that regulate muscle adaptation to contractile activity (see below), loading conditions, substrate supply, and environmental factors. Of the 25 muscle fibre genes, 13 genes regulate cellular responses to hormones (Ccnd3, Foxo1, Ppargc1a, Ncor1, Ppargc1b, Thra, Thrb, Mstn, Myod1, Bdnf, Sirt1, and Esrrg) as well as 6 genes that are linked to energy metabolism, i.e., in the form of mitochondrial biogenesis (Foxo1, Ppargc1a, Ppargc1b, Camk4, Sirt3, and Sirt1) (see supplemental data worksheet S7 for more detail).
(2)
In what human tissues are muscle fibre genes expressed?
To find out whether muscle fibre genes are primarily expressed in skeletal muscle or elsewhere, we retrieved human gene expression data from the Genotype-Tissue Expression (GTEx) Project website (https://gtexportal.org/home [39]). This analysis revealed that two of the muscle fibre genes, Vgll2 and MyoD1, are exclusively expressed in skeletal muscle. Moreover, Mapk12, Foxo1, and Nol3 are most expressed in the human skeletal muscle, but they are also expressed elsewhere in the body (Figure 5).
(3)
Are muscle fibre genes regulated in response to exercise or inactivity?
To systematically investigate whether muscle fibre genes are regulated by exercise, we retrieved from the MetaMEx website (https://www.metamex.eu/ [40]) meta-analysed human muscle fibre gene expression data, which compared pre and post exercise as well as inactivity data (Figure 5. Worksheet S4 in Supplementary Table S1 [40]). To find out whether muscle fibre proteins are phosphorylated, and whether they become phosphorylated or dephosphorylated after a bout of human exercise, we also retrieved published phosphoproteomics data [41] to analyse (see supplemental data worksheet S8 for more detail). The gene expression analysis identified PPARGC1A, which encodes the mitochondrial biogenesis regulator Pgc-1α; further, it also identified VGLL2 as the gene that roughly doubles its expression after acute bouts of endurance or resistance exercise in the human vastus lateralis muscle, and that which decreases its expression in inactive muscles. EPAS1, which encodes a hypoxia-induced transcription factor, also increases its expression after a bout of endurance and resistance exercise but decreases in response to inactivity. Conversely, MSTN expression decreases after a bout of endurance and resistance exercise but increases in response to inactivity. The expression changes in all muscle fibre genes in response to acute endurance exercise, acute resistance exercise, and inactivity are shown in Figure 6.
This reveals that muscle fibre genes such as PPARGC1A and VGLL2 increase their expression, whereas MSTN decreases its expression especially after a bout of endurance exercise. However, the opposite is true for inactivity.
When analysing muscle protein phosphorylation, we found that Vgll2 Ser261 phosphorylation increased by 30% after maximal muscle contractions in mice (p = 0.07) [42]. In addition, FOXO1, MAPK12, NOL3, NCOR1, and SIRT1 were detected as phosphorylated proteins in human muscle after a single high-intensity exercise bout. However, of these muscle fibre proteins, only MAPK12 Ser362 phosphorylation increased by more than 1.5-fold (p < 0.05) [41]. Collectively, this suggests that several muscle fibre genes are regulated in response to acute endurance or resistance exercise or inactivity.
(4)
What is known about sequence variability of the muscle fibre genes in human exome?
Human fibre type distribution in muscles vary in the human population and this is partially, ≈45%, explained by genetics [21]. To determine the frequency of human DNA sequence variants of muscle fibre genes, we retrieved exome sequence data for 60,706 humans [67]. The analysis of this data revealed that each muscle fibre gene had on average 160 missense, 3 loss-of-function, and 87 synonymous DNA variants (see supplemental data worksheet S9 for more detail). For BDKRB2, CCND3, NOL3, THRA, NCOR1, and EPAS2 homozygous loss-of-function DNA variants are reported [67]. Together, this suggests that exome sequence variability of human muscle fibre genes could at least partially explain the currently poorly explained variability of muscle fibre genes. Moreover, we also searched the term “muscle fibre genes” in the GWAS catalogue (https://www.ebi.ac.uk/gwas/) to find out whether these genes are associated with similar physiological or pathological phenotypes that could be linked to muscle fibre alterations. However, we did not find any systematic pattern (see supplemental data worksheet S3 for more detail).

3.2. Discussion

By conducting a systematic review, we have identified 25 genes whose gain or loss of function significantly changes muscle fibre percentages, or the abundance of myosin heavy chain isoforms in mouse muscles. This confirms that muscle fibre type proportions are a polygenic trait. The effect sizes of the mutated, individual genes range from a 37% reduction in plantaris type 2B fibres (Ppargc1a) to a 28% increase in type 2B fibres in the soleus muscle (Epas1). There is no muscle fibre gene whose manipulation results in a transformation of all muscle fibres to another muscle fibre type. There are known human DNA sequence variants for all muscle fibre genes, and this could be useful to identify in future DNA sequence variants that have a measurable effect size on human muscle fibre distribution.
The fact that 15 of the 25 muscle fibre genes (Foxo1, Ppargc1a, Ncor1, Ppargc1b, Thra, Thrb, Akirin1, Trib3, Nfatc1, Foxj3, Vgll2, Myod1, Epas1, Sirt1, and Esrrg), are transcriptional regulators suggests that muscle fibre-specific gene expression is primarily regulated transcriptionally. How these transcription factors and co-factors combine to regulate the coordinated expression of hundreds of muscle fibre-specific genes—in either a binary (i.e., an on/off regulation such as the case for Myh genes in pure fibres) or graded fashion (i.e., a higher/lower expression, e.g., mitochondrial genes)—is incompletely understood. Generally, the identity of cells, such as in muscle fibre identity, is often regulated by super-enhancers. Such super-enhancers are groups of enhancers bound by master transcription factors such as Brd4 and Med1 that compartmentalise DNA and its genes, allowing a coordinated regulation of the expression of target genes [68,69]. Consistent with this, the muscle-specific knockout of the super-enhancer-associated factor Med1 resulted in increased expression of Myh7 (type 1 myosin heavy chain), Myh2 (type 2A), as well as myoglobin and metabolic gene mRNAs. This suggests a coordinated regulation of muscle fibre identity-related genes [70]. This is consistent but does not proof that Med1 targets a fibre identity-regulating super enhancer.
Super enhancers are often associated with long intergenic non-coding RNA (lincRNAs) [71]. In relation to this, linc-MYH was identified as a lincRNA in “fast fibre” nuclei that inhibits slow-type enhancers and promotes fast/glycolytic muscle fibre gene expression [72]. The regulation of multiple fibre identify associated genes is consistent with a fast/glycolytic fibre super enhancer. Finally, Gunderson and colleagues have investigated the epigenetics of purified myonuclei from the predominantly fast twitch extensor digitorum longus and the slow twitch soleus mouse muscles [73]. They found an overrepresentation of binding sites for Mef2C, Nfatc2, and Ppara in the soleus and of MyoD1 and Sox1 in the EDL. There is knowledge of multiple, muscle fibre type identify regulating transcription factors, the likely existence of a fibre type-regulating super enhancer [72], and methods allowing to analyse transcriptional regulation of nuclei with a fast/glycolytic or slow/oxidative identity [73]. These should enable future researchers to find out why more than ten transcription factors can regulate sets of genes associated with a fast/glycolytic or slow/oxidative muscle fibre identity.
Motor neuron activity and contractions are the main stimuli for muscle fibre-specific gene expression and for muscle fibre transitions. Consistent with this, our analysis reveals several muscle fibre genes that are regulated by exercise or immobilisation. A key mode of regulation for adaptation to muscle contractions and exercise is transcriptional regulation. For example, acute human endurance or resistance exercise increases the expression of muscle fibre genes such as PPARGC1A and VGLL2 but decreases the expression of MSTN [40]. In contrast, inactivity decreases the expression of muscle fibre genes such as VGLL2 but increases the expression of MSTN [40]. A second mode of regulation is post-translational modification of muscle fibre proteins especially by phosphorylation. One example for this is the 1.5-fold increased phosphorylation of MAPK12 at Ser362 [41]. Finally, muscle fibre proteins can also be regulated by translocation. Prime examples of this are Nfat transcription factors. For example, Nfatc1 is more nuclear in slower type 1 than faster type 2 fibres and, therefore, translocates to the nucleus in response electrical stimulation [74]. Whether and how this drives gene expression that is related to muscle fibre identity, however, is poorly understood.
Finally, all muscle fibre genes are affected by human DNA sequence variation. This ranges from rare, homozygous knockouts of the whole gene (such as in the case of MSTN [75]) to single nucleotide polymorphisms that occur in all muscle fibre genes [67], such as Bdnf, Camk4, and Ccnd3. Given that the distribution of muscle fibres is ≈50% inherited [21], and given that we know few muscle fibre-proportion varying DNA sequence variants, the list of muscle fibre genes could be used to find out whether DNA sequence variants of muscle fibre genes occur at a higher or lower frequency in, e.g., elite sprinters or endurance athletes with often extreme fibre type proportions when compared to the sedentary population. Such a targeted, genetic analysis might help to uncover the genetic causes of the variability of muscle fibre proportions in the human population.
The first limitation of this study is that the choice of the targeted genes is subjective, that not all genes in the genome have been manipulated and that, e.g., the IMPC phenotyping pipeline does not include muscle fibre typing. A second limitation is that it is impossible to determine whether gene manipulation affects development, e.g., by blocking the slow-to-fast switch that occurs after birth in the fast muscles or is instead related to adult muscle fibre plasticity by, e.g., promoting a fast-to-slow twitch [4]. A third limitation is the strict inclusion and exclusion criteria used that included the use of “statistical significance” as an inclusion criterion. For example, the manipulation of genes such as Actn3 results in nearly significant fibre type changes (e.g., p = 0.051 in the spinalis) [76] and this is an issue because of the limited p-values [77] used. Further, as a systematic review requires clearcut inclusion and exclusion criteria, and due to the fact that we have chosen statistical significance as an inclusion criterion, this study cannot, therefore, be considered to be a systematic review. Fourth and finally, some publications suggest fibre type shifts based on mRNA expression alone and these were not included. Examples are transgenic mice with transgenes of p43 [78] and Med1 [70].

4. Conclusions

In summary, we detected 25 genes whose gain or loss of function alters muscle fibre proportions in mice. Given that the 25 genes causally affect fibre proportions, functionally relevant DNA sequence variants that are related to these genes should affect human muscle fibre type proportions, provided that the function of these genes is conserved between mice and men. Together with a list of muscle hypertrophy genes [79], endurance genes [80], and glucose uptake genes [81] this additional list of causative genes should help to identify causative DNA variants that influence human sport and exercise-related traits [82] as well as human health—as many exercise-related traits are predictors of health.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms232112933/s1.

Author Contributions

G.K. helped to conceptualise the research, performed and led the systematic literature research and bioinformatics analyses, as well as produced figures and helped to write the original draft. T.G.R. was the second systematic reviewer. M.M. helped to write the original draft, reviewed, and edited the manuscript. N.J.P. performed the bioinformatics analysis in Figure 5 and reviewed and edited the manuscript. M.S. helped to conceptualise the research, reviewed, and edited the manuscript. H.W. conceptualised the research, supervised the work, led the writing of the first draft, and edited later drafts. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are available in the supplementary data file.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

Actn3Actinin alpha 3
Akirin1Akirin 1
Bdkrb2Bradykinin receptor 2
BdnfBrain derived neurotrophic factor
Brd4Bromodomain containing 4
Camk4Calcium-/calmodulin-dependent protein kinase IV
Ccnd3Cyclin D3
Cpt1aCarnitine palmitoyltransferase 1A
EDLExtensor digitorum longus
Epas1Endothelial PAS domain protein 1
EsrrgOestrogen related receptor gamma
Foxj3Forkhead box J3
Foxo1Forkhead box O1
GTExGenotype-Tissue Expression
GWASGenome Wide Association Studies
Il15Interleukin 15
IMPCInternational Mouse Phenotyping Consortium
lincRNAslong intergenic non-coding
Mapk12Mitogen-activated protein kinase 12
Med1Mediator complex subunit 1
Mef2CMyocyte enhancer factor 2C
MeSHMedical Subjects Heading
MetaMExMeta-analysis of Skeletal Muscle Response to Exercise
MHCMyosin Heavy Chain
MstnMyostatin
Myh1Myosin heavy chain 1; encodes myosin 2X protein
Myh2Myosin heavy chain 2; encodes myosin 2A protein
Myh4Myosin heavy chain 4; encodes myosin 2AB protein
Myh7Myosin heavy chain 7; encodes myosin type 1 protein
Myod1Myogenic differentiation 1
NCBINational Center for Biotechnology Information
Ncor1Nuclear receptor corepressor 1
Nfatc1Nuclear factor of activated T cells 1
Nfatc2Nuclear factor of activated T cells 2
Nol3Nuclear protein 3
PICOPopulation/Problem, Intervention, Comparison, Outcome
PMIDPubMed identificatory
Ppargc1aPeroxisome proliferative activated receptor, gamma, coactivator 1 alpha
Ppargc1bPeroxisome proliferative activated receptor, gamma, coactivator 1beta
PparaPeroxisome proliferator activated receptor alpha
PRISMAPreferred Reporting Items for Systematic Review and Meta-Analysis Guideline
Sirt1Sirtuin 1
Sirt3Sirtuin 3
Sox1SRY (sex determining region Y)-box 1
ThraThyroid hormone receptor alpha
ThrbThyroid hormone receptor beta
TPMTranscript per million
Trib3Tribbles pseudokinase 3
UniprotUniversal Protein Resources
Vgll2Vestigial like family member 2

References

  1. Heron, M.I.; Richmond, F.J. In-series fiber architecture in long human muscles. J. Morphol. 1993, 216, 35–45. [Google Scholar] [CrossRef] [PubMed]
  2. Cristea, A.; Qaisar, R.; Edlund, P.K.; Lindblad, J.; Bengtsson, E.; Larsson, L. Effects of aging and gender on the spatial organization of nuclei in single human skeletal muscle cells. Aging Cell 2010, 9, 685–697. [Google Scholar] [CrossRef] [PubMed]
  3. Lexell, J.; Henriksson-Larsen, K.; Winblad, B.; Sjostrom, M. Distribution of different fiber types in human skeletal muscles: Effects of aging studied in whole muscle cross sections. Muscle Nerve 1983, 6, 588–595. [Google Scholar] [CrossRef] [PubMed]
  4. Schiaffino, S. Fibre types in skeletal muscle: A personal account. Acta Physiol. 2010, 199, 451–463. [Google Scholar] [CrossRef] [PubMed]
  5. Barany, M. ATPase activity of myosin correlated with speed of muscle shortening. J. Gen. Physiol. 1967, 50, 197–218. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Lowry, C.V.; Kimmey, J.S.; Felder, S.; Chi, M.M.; Kaiser, K.K.; Passonneau, P.N.; Kirk, K.A.; Lowry, O.H. Enzyme patterns in single human muscle fibers. J. Biol. Chem. 1978, 253, 8269–8277. [Google Scholar] [CrossRef]
  7. Brooke, M.H.; Kaiser, K.K. Three “myosin adenosine triphosphatase” systems: The nature of their pH lability and sulfhydryl dependence. J. Histochem. Cytochem. 1970, 18, 670–672. [Google Scholar] [CrossRef]
  8. Schiaffino, S. Muscle fiber type diversity revealed by anti-myosin heavy chain antibodies. FEBS J. 2018, 285, 3688–3694. [Google Scholar] [CrossRef]
  9. LaFramboise, W.A.; Daood, M.J.; Guthrie, R.D.; Moretti, P.; Schiaffino, S.; Ontell, M. Electrophoretic separation and immunological identification of type 2X myosin heavy chain in rat skeletal muscle. Biochim. Et Biophys. Acta (BBA) Gen. Subj. 1990, 1035, 109–112. [Google Scholar] [CrossRef]
  10. Lang, F.; Khaghani, S.; Türk, C.; Wiederstein, J.L.; Hölper, S.; Piller, T.; Nogara, L.; Blaauw, B.; Günther, S.; Müller, S.; et al. Single Muscle Fiber Proteomics Reveals Distinct Protein Changes in Slow and Fast Fibers during Muscle Atrophy. J. Proteome Res. 2018, 17, 3333–3347. [Google Scholar] [CrossRef]
  11. Murgia, M.; Nogara, L.; Baraldo, M.; Reggiani, C.; Mann, M.; Schiaffino, S. Protein profile of fiber types in human skeletal muscle: A single-fiber proteomics study. Skelet. Muscle 2021, 11, 24. [Google Scholar] [CrossRef] [PubMed]
  12. Smerdu, V.; Karsch-Mizrachi, I.; Campione, M.; Leinwand, L.; Schiaffino, S. Type IIx myosin heavy chain transcripts are expressed in type IIb fibers of human skeletal muscle. Am. J. Physiol. 1994, 267 Pt 1, C1723–C1728. [Google Scholar] [CrossRef]
  13. Pette, D.; Staron, R.S. Myosin isoforms, muscle fiber types, and transitions. Microsc. Res. Tech. 2000, 50, 500–509. [Google Scholar] [CrossRef]
  14. Schiaffino, S.; Rossi, A.C.; Smerdu, V.; Leinwand, L.A.; Reggiani, C. Developmental myosins: Expression patterns and functional significance. Skelet. Muscle 2015, 5, 22. [Google Scholar] [CrossRef] [Green Version]
  15. Johnson, M.A.; Polgar, J.; Weightman, D.; Appleton, D. Data on the distribution of fibre types in thirty-six human muscles. An autopsy study. J. Neurol. Sci. 1973, 18, 111–129. [Google Scholar] [CrossRef]
  16. Simoneau, J.A.; Bouchard, C. Human variation in skeletal muscle fiber-type proportion and enzyme activities. Am. J. Physiol. 1989, 257 Pt 1, E567–E572. [Google Scholar] [CrossRef]
  17. Costill, D.L.; Daniels, J.; Evans, W.; Fink, W.; Krahenbuhl, G.; Saltin, B. Skeletal muscle enzymes and fiber composition in male and female track athletes. J. Appl. Physiol. 1976, 40, 149–154. [Google Scholar] [CrossRef]
  18. Costill, D.L.; Fink, W.J.; Pollock, M.L. Muscle fiber composition and enzyme activities of elite distance runners. Med. Sci. Sports 1976, 8, 96–100. [Google Scholar] [CrossRef] [PubMed]
  19. Coyle, E.F.; Bell, S.; Costill, D.L.; Fink, W.J. Skeletal muscle fiber characteristics of world class shot-putters. Res. Q. 1978, 49, 278–284. [Google Scholar] [CrossRef]
  20. Foster, C.; Costill, D.L.; Daniels, J.T.; Fink, W.J. Skeletal muscle enzyme activity, fiber composition and VO2 max in relation to distance running performance. Eur. J. Appl. Physiol. Occup. Physiol. 1978, 39, 73–80. [Google Scholar] [CrossRef] [PubMed]
  21. Simoneau, J.A.; Bouchard, C. Genetic determinism of fiber type proportion in human skeletal muscle. FASEB J. 1995, 9, 1091–1095. [Google Scholar] [CrossRef] [PubMed]
  22. Gollnick, P.D.; Armstrong, R.B.; Saltin, B.; Saubert, C.W.; Sembrowich, W.L.; Shepherd, R.E. Effect of training on enzyme activity and fiber composition of human skeletal muscle. J. Appl. Physiol. 1973, 34, 107–111. [Google Scholar] [CrossRef] [Green Version]
  23. Rico-Sanz, J.; Rankinen, T.; Joanisse, D.R.; Leon, A.S.; Skinner, J.S.; Wilmore, J.H.; Rao, D.C.; Bouchard, C. Familial resemblance for muscle phenotypes in the HERITAGE Family Study. Med. Sci. Sports Exerc. 2003, 35, 1360–1366. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Plotkin, D.L.; Roberts, M.D.; Haun, C.T.; Schoenfeld, B.J. Muscle Fiber Type Transitions with Exercise Training: Shifting Perspectives. Sports 2021, 9, 127. [Google Scholar] [CrossRef] [PubMed]
  25. Vrbova, G. The effect of motoneurone activity on the speed of contraction of striated muscle. J. Physiol. 1963, 169, 513–526. [Google Scholar] [CrossRef] [Green Version]
  26. Buller, A.J.; Eccles, J.C.; Eccles, R.M. Interactions between motoneurones and muscles in respect of the characteristic speeds of their responses. J. Physiol. 1960, 150, 417–439. [Google Scholar] [CrossRef] [PubMed]
  27. Pette, D.; Vrbova, G. What does chronic electrical stimulation teach us about muscle plasticity? Muscle Nerve 1999, 22, 666–677. [Google Scholar] [CrossRef]
  28. Windisch, A.; Gundersen, K.; Szabolcs, M.J.; Gruber, H.; Lomo, T. Fast to slow transformation of denervated and electrically stimulated rat muscle. J. Physiol. 1998, 510 Pt 1, 623–632. [Google Scholar] [CrossRef]
  29. Chin, E.R.; Olson, E.N.; Richardson, J.A.; Yang, Q.; Humphries, C.; Shelton, J.M.; Wu, H.; Zhu, W.; Bassel-Duby, R.; Williams, R.S. A calcineurin-dependent transcriptional pathway controls skeletal muscle fiber type. Genes Dev. 1998, 12, 2499–2509. [Google Scholar] [CrossRef] [Green Version]
  30. Serrano, A.L.; Murgia, M.; Pallafacchina, G.; Calabria, E.; Coniglio, P.; Lomo, T.; Schiaffino, S. Calcineurin controls nerve activity-dependent specification of slow skeletal muscle fibers but not muscle growth. Proc. Natl. Acad. Sci. USA 2001, 98, 13108–13113. [Google Scholar] [CrossRef]
  31. Murgia, M.; Serrano, A.L.; Calabria, E.; Pallafacchina, G.; Lomo, T.; Schiaffino, S. Ras is involved in nerve-activity-dependent regulation of muscle genes. Nat. Cell Biol. 2000, 2, 142–147. [Google Scholar] [CrossRef] [PubMed]
  32. Lin, J.; Wu, H.; Tarr, P.T.; Zhang, C.Y.; Wu, Z.; Boss, O.; Michael, L.F.; Puigserver, P.; Isotani, E.; Olson, E.N.; et al. Transcriptional co-activator PGC-1 alpha drives the formation of slow-twitch muscle fibres. Nature 2002, 418, 797–801. [Google Scholar] [CrossRef] [PubMed]
  33. Birling, M.C.; Yoshiki, A.; Adams, D.J.; Ayabe, S.; Beaudet, A.L.; Bottomley, J.; Bradley, A.; Brown, S.D.M.; Bürger, A.; Bushell, W.; et al. A resource of targeted mutant mouse lines for 5061 genes. Nat. Genet. 2021, 53, 416–419. [Google Scholar] [CrossRef] [PubMed]
  34. Page, M.J.; McKenzie, J.E.; Bossuyt, P.M.; Boutron, I.; Hoffmann, T.C.; Mulrow, C.D.; Shamseer, L.; Tetzlaff, J.M.; Akl, E.A.; Brennan, S.E.; et al. The PRISMA 2020 statement: An updated guideline for reporting systematic reviews. BMJ 2021, 372, n71. [Google Scholar] [CrossRef] [PubMed]
  35. Schardt, C.; Adams, M.B.; Owens, T.; Keitz, S.; Fontelo, P. Utilization of the PICO framework to improve searching PubMed for clinical questions. BMC Med. Inform. Decis. Mak. 2007, 7, 16. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Schneider, C.A.; Rasband, W.S.; Eliceiri, K.W. NIH Image to ImageJ: 25 years of image analysis. Nat. Methods 2012, 9, 671–675. [Google Scholar] [CrossRef]
  37. Szklarczyk, D.; Gable, A.L.; Lyon, D.; Junge, A.; Wyder, S.; Huerta-Cepas, J.; Simonovic, M.; Doncheva, N.T.; Morris, J.H.; Bork, P.; et al. STRING v11: Protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets. Nucleic Acids Res. 2019, 47, D607–D613. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Chen, J.; Bardes, E.E.; Aronow, B.J.; Jegga, A.G. ToppGene Suite for gene list enrichment analysis and candidate gene prioritization. Nucleic Acids Res. 2009, 37, W305–W311. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. GTExConsortium; Ardlie, K.G.; Deluca, D.S.; Segre, A.V.; Sullivan, T.J.; Young, T.R.; Gelfand, E.T.; Trowbridge, C.A.; Maller, J.B.; Tukiainen, T.; et al. Human genomics. The Genotype-Tissue Expression (GTEx) pilot analysis: Multitissue gene regulation in humans. Science 2015, 348, 648–660. [Google Scholar]
  40. Pillon, N.J.; Gabriel, B.M.; Dollet, L.; Smith, J.A.B.; Sardón Puig, L.; Botella, J.; Bishop, D.J.; Krook, A.; Zierath, J.R. Transcriptomic profiling of skeletal muscle adaptations to exercise and inactivity. Nat. Commun. 2020, 11, 470. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  41. Hoffman, N.J.; Parker, B.L.; Chaudhuri, R.; Fisher-Wellman, K.H.; Kleinert, M.; Humphrey, S.J.; Yang, P.; Holliday, M.; Trefely, S.; Fazakerley, D.J.; et al. Global Phosphoproteomic Analysis of Human Skeletal Muscle Reveals a Network of Exercise-Regulated Kinases and AMPK Substrates. Cell Metab. 2015, 22, 922–935. [Google Scholar] [CrossRef] [PubMed]
  42. Potts, G.K.; McNally, R.M.; Blanco, R.; You, J.-S.; Hebert, A.S.; Westphall, M.S.; Coon, J.J.; Hornberger, T.A. A map of the phosphoproteomic alterations that occur after a bout of maximal-intensity contractions. J. Physiol. 2017, 595, 5209–5226. [Google Scholar] [CrossRef] [Green Version]
  43. Rao, V.V.; Sangiah, U.; Mary, K.A.; Akira, S.; Mohanty, A. Role of Akirin1 in the regulation of skeletal muscle fiber-type switch. J. Cell. Biochem. 2019, 120, 11284–11304. [Google Scholar] [CrossRef]
  44. Reis, F.C.G.; Haro, A.S.; Bacurau, A.V.N.; Hirabara, S.M.; Wasinski, F.; Ormanji, M.S.; Moreira, J.B.N.; Kiyomoto, B.H.; Bertoncini, C.R.A.; Brum, P.C.; et al. Deletion of Kinin B2 Receptor Alters Muscle Metabolism and Exercise Performance. PLoS ONE 2015, 10, e0134844. [Google Scholar] [CrossRef] [Green Version]
  45. Delezie, J.; Weihrauch, M.; Maier, G.; Tejero, R.; Ham, D.J.; Gill, J.F.; Karrer-Cardel, B.; Rüegg, M.A.; Tabares, L.; Handschin, C. BDNF is a mediator of glycolytic fiber-type specification in mouse skeletal muscle. Proc. Natl. Acad. Sci. USA 2019, 116, 16111–16120. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Akimoto, T.; Ribar, T.J.; Williams, R.S.; Yan, Z. Skeletal muscle adaptation in response to voluntary running in Ca2+/calmodulin-dependent protein kinase IV-deficient mice. Am. J. Physiology. Cell Physiol. 2004, 287, C1311–C1319. [Google Scholar] [CrossRef] [Green Version]
  47. Giannattasio, S.; Giacovazzo, G.; Bonato, A.; Caruso, C.; Luvisetto, S.; Coccurello, R.; Caruso, M. Lack of cyclin D3 induces skeletal muscle fiber-type shifting, increased endurance performance and hypermetabolism. Sci. Rep. 2018, 8, 12792. [Google Scholar] [CrossRef] [Green Version]
  48. Hénique, C.; Mansouri, A.; Vavrova, E.; Lenoir, V.; Ferry, A.; Esnous, C.; Ramond, E.; Girard, J.; Bouillaud, F.; Prip-Buus, C.; et al. Increasing mitochondrial muscle fatty acid oxidation induces skeletal muscle remodeling toward an oxidative phenotype. FASEBJ. Off. Publ. Fed. Am. Soc. Exp. Biol. 2015, 29, 2473–2483. [Google Scholar] [CrossRef] [PubMed]
  49. Alexander, M.S.; Shi, X.; Voelker, K.A.; Grange, R.W.; Garcia, J.A.; Hammer, R.E.; Garry, D.J. Foxj3 transcriptionally activates Mef2c and regulates adult skeletal muscle fiber type identity. Dev. Biol. 2010, 337, 396–404. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Kamei, Y.; Miura, S.; Suzuki, M.; Kai, Y.; Mizukami, J.; Taniguchi, T.; Mochida, K.; Hata, T.; Matsuda, J.; Aburatani, H.; et al. Skeletal muscle FOXO1 (FKHR) transgenic mice have less skeletal muscle mass, down-regulated Type I (slow twitch/red muscle) fiber genes, and impaired glycemic control. J. Biol. Chem. 2004, 279, 41114–41123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  51. Foster, W.H.; Tidball, J.G.; Wang, Y. p38γ activity is required for maintenance of slow skeletal muscle size. Muscle Nerve 2012, 45, 266–273. [Google Scholar] [CrossRef]
  52. Hennebry, A.; Berry, C.; Siriett, V.; O’Callaghan, P.; Chau, L.; Watson, T.; Sharma, M.; Kambadur, R. Myostatin regulates fiber-type composition of skeletal muscle by regulating MEF2 and MyoD gene expression. Am. J. Physiol. Cell Physiol. 2009, 296, C525–C534. [Google Scholar] [CrossRef]
  53. Macharia, R.; Otto, A.; Valasek, P.; Patel, K. Neuromuscular junction morphology, fiber-type proportions, and satellite-cell proliferation rates are altered in MyoD(-/-) mice. Muscle Nerve 2010, 42, 38–52. [Google Scholar] [CrossRef] [PubMed]
  54. Ehlers, M.L.; Celona, B.; Black, B.L. NFATc1 controls skeletal muscle fiber type and is a negative regulator of MyoD activity. Cell Rep. 2014, 8, 1639–1648. [Google Scholar] [CrossRef] [Green Version]
  55. Mitchell, A.S.; Smith, I.C.; Gamu, D.; Donath, S.; Tupling, A.R.; Quadrilatero, J. Functional, morphological, and apoptotic alterations in skeletal muscle of ARC deficient mice. Apoptosis 2015, 20, 310–326. [Google Scholar] [CrossRef] [PubMed]
  56. Yu, F.; Göthe, S.; Wikström, L.; Forrest, D.; Vennström, B.; Larsson, L. Effects of thyroid hormone receptor gene disruption on myosin isoform expression in mouse skeletal muscles. Am. J. Physiol.-Regul. Integr. Comp. Physiol. 2000, 278, R1545–R1554. [Google Scholar] [CrossRef] [Green Version]
  57. Rasbach, K.A.; Gupta, R.K.; Ruas, J.L.; Wu, J.; Naseri, E.; Estall, J.L.; Spiegelman, B.M. PGC-1alpha regulates a HIF2alpha-dependent switch in skeletal muscle fiber types. Proc. Natl. Acad. Sci. USA 2010, 107, 21866–21871. [Google Scholar] [CrossRef] [Green Version]
  58. Rangwala, S.M.; Wang, X.; Calvo, J.A.; Lindsley, L.; Zhang, Y.; Deyneko, G.; Beaulieu, V.; Gao, J.; Turner, G.; Markovits, J. Estrogen-related receptor gamma is a key regulator of muscle mitochondrial activity and oxidative capacity. J. Biol. Chem. 2010, 285, 22619–22629. [Google Scholar] [CrossRef] [Green Version]
  59. Quinn, L.S.; Anderson, B.G.; Conner, J.D.; Wolden-Hanson, T. IL-15 overexpression promotes endurance, oxidative energy metabolism, and muscle PPARdelta, SIRT1, PGC-1alpha, and PGC-1beta expression in male mice. Endocrinology 2013, 154, 232–245. [Google Scholar] [CrossRef] [Green Version]
  60. Yamamoto, H.; Williams, E.G.; Mouchiroud, L.; Canto, C.; Fan, W.; Downes, M.; Heligon, C.; Barish, G.D.; Desvergne, B.; Evans, R.M.; et al. NCoR1 is a conserved physiological modulator of muscle mass and oxidative function. Cell 2011, 147, 827–839. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Arany, Z.; Lebrasseur, N.; Morris, C.; Smith, E.; Yang, W.; Ma, Y.; Chin, S.; Spiegelman, B.M. The transcriptional coactivator PGC-1beta drives the formation of oxidative type IIX fibers in skeletal muscle. Cell Metab. 2007, 5, 35–46. [Google Scholar] [CrossRef] [PubMed]
  62. Lin, L.; Chen, K.; Abdel Khalek, W.; Ward, J.L., 3rd; Yang, H.; Chabi, B.; Wrutniak-Cabello, C.; Tong, Q. Regulation of skeletal muscle oxidative capacity and muscle mass by SIRT3. PLoS ONE 2014, 9, e85636. [Google Scholar] [CrossRef]
  63. An, D.; Lessard, S.J.; Toyoda, T.; Lee, M.-Y.; Koh, H.-J.; Qi, L.; Hirshman, M.F.; Goodyear, L.J. Overexpression of TRB3 in muscle alters muscle fiber type and improves exercise capacity in mice. Am. J. Physiol.-Regul. Integr. Comp. Physiol. 2014, 306, R925–R933. [Google Scholar] [CrossRef] [Green Version]
  64. Honda, M.; Hidaka, K.; Fukada, S.-I.; Sugawa, R.; Shirai, M.; Ikawa, M.; Morisaki, T. Vestigial-like 2 contributes to normal muscle fiber type distribution in mice. Sci. Rep. 2017, 7, 7168. [Google Scholar] [CrossRef] [Green Version]
  65. Chalkiadaki, A.; Igarashi, M.; Nasamu, A.S.; Knezevic, J.; Guarente, L. Muscle-specific SIRT1 gain-of-function increases slow-twitch fibers and ameliorates pathophysiology in a mouse model of duchenne muscular dystrophy. PLoS Genet. 2014, 10, e1004490. [Google Scholar] [CrossRef]
  66. Izumo, S.; Nadal-Ginard, B.; Mahdavi, V. All members of the MHC multigene family respond to thyroid hormone in a highly tissue-specific manner. Science 1986, 231, 597–600. [Google Scholar] [CrossRef]
  67. Lek, M.; Karczewski, K.J.; Minikel, E.V.; Samocha, K.E.; Banks, E.; Fennell, T.; O’Donnell-Luria, A.H.; Ware, J.S.; Hill, A.J.; Cummings, B.B.; et al. Analysis of protein-coding genetic variation in 60,706 humans. Nature 2016, 536, 285–291. [Google Scholar] [CrossRef] [Green Version]
  68. Hnisz, D.; Abraham, B.J.; Lee, T.I.; Lau, A.; Saint-André, V.; Sigova, A.A.; Hoke, H.A.; Young, R.A. Super-enhancers in the control of cell identity and disease. Cell 2013, 155, 934–947. [Google Scholar] [CrossRef] [Green Version]
  69. Sabari, B.R.; Dall’Agnese, A.; Boija, A.; Klein, I.A.; Coffey, E.L.; Shrinivas, K.; Abraham, B.J.; Hannett, N.M.; Zamudio, A.V.; Manteiga, J.C.; et al. Coactivator condensation at super-enhancers links phase separation and gene control. Science 2018, 361, eaar3958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  70. Chen, W.; Zhang, X.; Birsoy, K.; Roeder, R.G. A muscle-specific knockout implicates nuclear receptor coactivator MED1 in the regulation of glucose and energy metabolism. Proc. Natl. Acad. Sci. USA 2010, 107, 10196–10201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  71. Soibam, B. Super-lncRNAs: Identification of lncRNAs that target super-enhancers via RNA:DNA:DNA triplex formation. RNA 2017, 23, 1729–1742. [Google Scholar] [CrossRef] [PubMed]
  72. Sakakibara, I.; Santolini, M.; Ferry, A.; Hakim, V.; Maire, P. Six homeoproteins and a Iinc-RNA at the fast MYH locus lock fast myofiber terminal phenotype. PLoS Genet. 2014, 10, e1004386. [Google Scholar] [CrossRef]
  73. Bengtsen, M.; Winje, I.M.; Eftestøl, E.; Landskron, J.; Sun, C.; Nygård, K.; Domanska, D.; Millay, D.P.; Meza-Zepeda, L.A.; Gundersen, K. Comparing the epigenetic landscape in myonuclei purified with a PCM1 antibody from a fast/glycolytic and a slow/oxidative muscle. PLoS Genet. 2021, 17, e1009907. [Google Scholar] [CrossRef] [PubMed]
  74. Tothova, J.; Blaauw, B.; Pallafacchina, G.; Rudolf, R.; Argentini, C.; Reggiani, C.; Schiaffino, S. NFATc1 nucleocytoplasmic shuttling is controlled by nerve activity in skeletal muscle. J. Cell Sci. 2006, 119 Pt 8, 1604–1611. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Schuelke, M.; Wagner, K.R.; Stolz, L.E.; Hubner, C.; Riebel, T.; Komen, W.; Braun, T.; Tobin, J.F.; Lee, S.J. Myostatin mutation associated with gross muscle hypertrophy in a child. N. Engl. J. Med. 2004, 350, 2682–2688. [Google Scholar] [CrossRef] [Green Version]
  76. MacArthur, D.G.; Seto, J.T.; Chan, S.; Quinlan, K.G.; Raftery, J.M.; Turner, N.; Nicholson, M.D.; Kee, A.J.; Hardeman, E.C.; Gunning, P.W.; et al. An Actn3 knockout mouse provides mechanistic insights into the association between {alpha}-actinin-3 deficiency and human athletic performance. Hum. Mol. Genet. 2008, 17, 1076–1086. [Google Scholar] [CrossRef] [Green Version]
  77. Amrhein, V.; Greenland, S.; McShane, B. Scientists rise up against statistical significance. Nature 2019, 567, 305–307. [Google Scholar] [CrossRef] [Green Version]
  78. Casas, F.; Pessemesse, L.; Grandemange, S.; Seyer, P.; Gueguen, N.; Baris, O.; Lepourry, L.; Cabello, G.; Wrutniak-Cabello, C. Overexpression of the mitochondrial T3 receptor p43 induces a shift in skeletal muscle fiber types. PLoS ONE 2008, 3, e2501. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Verbrugge, S.A.J.; Schönfelder, M.; Becker, L.; Yaghoob Nezhad, F.; Hrabě de Angelis, M.; Wackerhage, H. Genes Whose Gain or Loss-Of-Function Increases Skeletal Muscle Mass in Mice: A Systematic Literature Review. Front. Physiol. 2018, 9, 553. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  80. Nezhad, Y.F.; Verbrugge, S.A.J.; Schonfelder, M.; Becker, L.; Hrabe de Angelis, M.; Wackerhage, H. Genes Whose Gain or Loss-of-Function Increases Endurance Performance in Mice: A Systematic Literature Review. Front. Physiol. 2019, 10, 262. [Google Scholar] [CrossRef] [Green Version]
  81. Verbrugge, S.A.J.; Alhusen, J.A.; Kempin, S.; Pillon, N.J.; Rozman, J.; Wackerhage, H.; Kleinert, M. Genes controlling skeletal muscle glucose uptake and their regulation by endurance and resistance exercise. J. Cell. Biochem. 2022, 123, 202–214. [Google Scholar] [CrossRef] [PubMed]
  82. Sarzynski, M.A.; Bouchard, C. World-class athletic performance and genetic endowment. Nat. Metab. 2020, 2, 796–798. [Google Scholar] [CrossRef] [PubMed]
Ijms 23 12933 g001 550
Figure 1. PRISMA flow chart of the study research methodology used on the databases. See also supplemental data worksheet S1.
Figure 1. PRISMA flow chart of the study research methodology used on the databases. See also supplemental data worksheet S1.
Ijms 23 12933 g001
Ijms 23 12933 g002 550
Figure 2. Genes whose gain or loss of function significantly changes the percentages of type 1 (A), type 2A (B), or type 2X/B (C) fibres when compared to wildtype or control mice. Red symbols: loss-of-function; blue symbols: gain-of-function. ● Soleus; ⬣ extensor digitorum longus (EDL); ◆ plantaris; ◼ tibialis anterior; ▲ biceps femoris; ▼ quadriceps; and ◯ gastrocnemius. In (C), filled symbols refer to 2X muscle fibres and symbols without filling indicated 2B fibres. Note that the “◯”-marked data in (A,B) were reported as statistically significant even though they are close to “0” (i.e., no change).
Figure 2. Genes whose gain or loss of function significantly changes the percentages of type 1 (A), type 2A (B), or type 2X/B (C) fibres when compared to wildtype or control mice. Red symbols: loss-of-function; blue symbols: gain-of-function. ● Soleus; ⬣ extensor digitorum longus (EDL); ◆ plantaris; ◼ tibialis anterior; ▲ biceps femoris; ▼ quadriceps; and ◯ gastrocnemius. In (C), filled symbols refer to 2X muscle fibres and symbols without filling indicated 2B fibres. Note that the “◯”-marked data in (A,B) were reported as statistically significant even though they are close to “0” (i.e., no change).
Ijms 23 12933 g002
Ijms 23 12933 g003 550
Figure 3. Genes whose gain or loss of function significantly changes the relative expression of type 1 (A), type 2A (B), type 2X (C), or type 2B (D) myosin heavy chain (MHC) isoforms when compared to wildtype or control mice. Red symbols: loss-of-function; blue symbols: gain-of-function. ● Soleus; ⬣ extensor digitorum longus (EDL); ◼ tibial anterior; ▼ quadriceps; and ◯ gastrocnemius.
Figure 3. Genes whose gain or loss of function significantly changes the relative expression of type 1 (A), type 2A (B), type 2X (C), or type 2B (D) myosin heavy chain (MHC) isoforms when compared to wildtype or control mice. Red symbols: loss-of-function; blue symbols: gain-of-function. ● Soleus; ⬣ extensor digitorum longus (EDL); ◼ tibial anterior; ▼ quadriceps; and ◯ gastrocnemius.
Ijms 23 12933 g003
Ijms 23 12933 g004 550
Figure 4. String analysis for interactions between muscle fibre proteins. Lines between proteins indicate evidence for interaction. See supplemental data S5 for more detail.
Figure 4. String analysis for interactions between muscle fibre proteins. Lines between proteins indicate evidence for interaction. See supplemental data S5 for more detail.
Ijms 23 12933 g004
Ijms 23 12933 g005 550
Figure 5. Heatmap illustrating the expression of muscle fibre genes in different human tissues in transcript per million (TPM) and in the order of highest expression in skeletal muscle. See also supplemental data worksheet S6.
Figure 5. Heatmap illustrating the expression of muscle fibre genes in different human tissues in transcript per million (TPM) and in the order of highest expression in skeletal muscle. See also supplemental data worksheet S6.
Ijms 23 12933 g005
Ijms 23 12933 g006 550
Figure 6. Meta−analysed expression change data of muscle fibre genes after acute bouts of endurance exercise, resistance exercise, and inactivity (“*” indicates FDR < 0.05). See supplemental data worksheet S4 for more detail.
Figure 6. Meta−analysed expression change data of muscle fibre genes after acute bouts of endurance exercise, resistance exercise, and inactivity (“*” indicates FDR < 0.05). See supplemental data worksheet S4 for more detail.
Ijms 23 12933 g006
Table
Table 1. Search Strategy used.
Table 1. Search Strategy used.
DatabaseSearch Formula
Medline (via PubMed)(mice OR “mouse” OR “mouse transgenic” OR “mice transgenic” OR “mouse knockout” OR “mice knockout” OR “mouse model” OR “mice model” OR “mice overexpressed” OR “mouse overexpressed”) AND (“gene expression” OR “gene knockout” OR “gene overexpression” OR “gene knock in” OR “gene transfer techniques” OR “gene deletion”) AND (“muscle fiber distribution” OR “muscle fiber fast twitch” OR “muscle fiber slow twitch” OR “muscle fiber type I” OR “muscle fiber type II” OR “oxidative muscle” OR “oxidative fiber” OR “glycolytic muscle” OR “glycolytic fiber”)
Table
Table 2. Summary of included studies. Main characteristics and outcome.
Table 2. Summary of included studies. Main characteristics and outcome.
StudyGeneKnockout or OverexpressionFibre Type Analysis ProcedureOutput MeasureOutcome (Relative to Wildtype)
[43]Akirin1KnockoutMHC immunofluorescence
MHC analysis		Quadriceps Type 1 ↓
Quadriceps Type 2A ↑
Quadriceps Myh7 1 ↓
Quadriceps Myh2 2A↑
[44]Bdkrb2KnockoutATPase stainingSoleus % type 1 ↑
Soleus % type 2A ↓
Soleus % type Intermediary ↑
[45]BdnfKnockoutMHC immunohistochemistryTA % type 2X↑
TA % type 2B ↓
EDL % type 2X ↑
EDL % type 2B ↓
[46]CamK4KnockoutMHC immunofluorescenceSoleus % type 1 ↑
Soleus % type 2A↓
Soleus % type other↓
[47]Ccnd3KnockoutMHC immunostainingTA % myofiber type 2A↑
TA % myofiber type 2B↓
[48]Cpt1aConditional knockoutMHC immunostaining TA % type 2A↑
[49]FoxJ3KnockoutATPase stainingEDL % type 1 ↓
Soleus % type 1 ↓
[50]Foxo1OverexpressionATPase stainingSoleus number type 1 ↓
Soleus number type 2 ↑
[51]Mapk12KnockoutMHC immunostaining Soleus % type 1↑
Soleus % type 2A ↓
[52]MstnKnockout MHC immunostaining Biceps femoris % type 1 ↓
Biceps femoris % type 2A ↓
Biceps femoris % type 2B/X ↑
TA % type 2A↓
TA % type 2B/X↑
[53]MyoD1KnockoutMHC immunohistochemistry EDL % Type 1↑
[54]Nfatc1Conditional knockoutATPase stainingSoleus number Type 1↓
[55]Nol3KnockoutImmunofluorescence MHCSoleus % type 1↓
Soleus % Type 2A↑
Plantaris % Type 2A↓
Plantaris % Type 2B↑
[32]Ppargc1aOverexpressionATPase stainingPlantaris number Type 1↑
Plantaris number Type 2A↑
Plantaris number Type 2B ↓
[56]ThraKnockoutATPase stainingSoleus % Type 1↑
Soleus % Type 2A↓
ThrbSoleus % Type 1↓
Soleus % Type 2A↑
[57]Epas1Conditional KnockoutMHC analysis Soleus number type 1↓
Soleus number type 2B ↑
Soleus MHC 1 ↓
Soleus MHC 2B ↑
[58]EsrrgConditional overexpressionMHC gene expression analysisGastrocnemius MHC 2A ↑
Gastrocnemius MHC 2B ↓
[59]Il15Conditional overexpressionMHC gene expression analysisSoleus MHC 1 ↑
Soleus MHC 2A ↓
Soleus a MHC 2X ↓
EDL MHC 1↑
EDL MHC 2X ↑
Gastrocnemius MHC 1 ↑
[60]Ncor1Conditional knockoutMHC gene expression analysisGastrocnemius MHC 2A ↑
Gastrocnemius MHC 2X ↑
Quadriceps MHC 1 ↑
Quadriceps MHC 2A ↑
Quadriceps MHC 2X ↑
[61]Ppargc1bConditional overexpressionMHC gene expression analysisGastrocnemius MHC 1 ↓
Gastrocnemius MHC 2X ↑
Gastrocnemius MHC 2B ↓
EDL MHC 1 ↓
EDL MHC 2A ↓
EDL MHC 2X ↑
EDL MHC 2B ↓
TA MHC 2A ↓
TA MHC 2X ↑
TA MHC 2B ↓
[62]Sirt3Conditional overexpressionWestern blotQuadriceps MHC 1↑
Quadriceps MHC 2B↓
[63]Trib3Conditional overexpressionElectrophorese of MHCSoleus MHC 1↑
Soleus MHC 2B↓
TA MHC 2B↓
[64]Vgll2KnockoutqPCR MHC analysisSoleus number type 2A↓
Soleus number type 2B ↑
Soleus Myh7 1↓
Soleus Myh2 2A ↓
Soleus Myh1 2X ↑
Soleus Myh4 2B↑
[65]Sirt1Conditional overexpressionATPase staining and PCR MHC analysisGastrocnemius % Type 1↑
Gastrocnemius % Type 2↓
Gastrocnemius MHC 2A↑
Gastrocnemius MHC 2X↑
Gastrocnemius MHC 2B↓
Legend ↑ Significantly increased proportion, concentration, or expression when compared to control. ↓ Significantly decreased proportion, concentration, or expression when compared to control. For more detail, see also Table S2 in the supplemental data.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Kuhnen, G.; Guedes Russomanno, T.; Murgia, M.; Pillon, N.J.; Schönfelder, M.; Wackerhage, H. Genes Whose Gain or Loss of Function Changes Type 1, 2A, 2X, or 2B Muscle Fibre Proportions in Mice—A Systematic Review. Int. J. Mol. Sci. 2022, 23, 12933. https://doi.org/10.3390/ijms232112933

AMA Style

Kuhnen G, Guedes Russomanno T, Murgia M, Pillon NJ, Schönfelder M, Wackerhage H. Genes Whose Gain or Loss of Function Changes Type 1, 2A, 2X, or 2B Muscle Fibre Proportions in Mice—A Systematic Review. International Journal of Molecular Sciences. 2022; 23(21):12933. https://doi.org/10.3390/ijms232112933

Chicago/Turabian Style

Kuhnen, Gabryela, Tiago Guedes Russomanno, Marta Murgia, Nicolas J. Pillon, Martin Schönfelder, and Henning Wackerhage. 2022. "Genes Whose Gain or Loss of Function Changes Type 1, 2A, 2X, or 2B Muscle Fibre Proportions in Mice—A Systematic Review" International Journal of Molecular Sciences 23, no. 21: 12933. https://doi.org/10.3390/ijms232112933

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop