Next Article in Journal
Analysis of SI-Related BoGAPDH Family Genes and Response of BoGAPC to SI Signal in Brassica oleracea L.
Previous Article in Journal
The Expression of RAAS Key Receptors, Agtr2 and Bdkrb1, Is Downregulated at an Early Stage in a Rat Model of Wolfram Syndrome
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Genes
Volume 12
Issue 11
10.3390/genes12111718
genes-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:
Article

Gene Expression Profiling of Skeletal Muscles

by
Sarah I. Alto
1,2,
Chih-Ning Chang
1,2,
Kevin Brown
1,3,
Chrissa Kioussi
1,2,* and
Theresa M. Filtz
1,2
1
Department of Pharmaceutical Sciences, College of Pharmacy, Oregon State University, Corvallis, OR 97331, USA
2
Molecular and Cellular Biology Graduate Program, Graduate School, Oregon State University, Corvallis, OR 97331, USA
3
School of Chemical, Biological, and Environmental Engineering, College of Engineering, Oregon State University, Corvallis, OR 97331, USA
*
Author to whom correspondence should be addressed.
Genes 2021, 12(11), 1718; https://doi.org/10.3390/genes12111718
Submission received: 23 September 2021 / Revised: 21 October 2021 / Accepted: 27 October 2021 / Published: 28 October 2021
(This article belongs to the Section Human Genomics and Genetic Diseases)
Browse Figures
Review Reports Versions Notes

Abstract

:
Next-generation sequencing provides an opportunity for an in-depth biocomputational analysis to identify gene expression patterns between soleus and tibialis anterior, two well-characterized skeletal muscles, and analyze their gene expression profiling. RNA read counts were analyzed for differential gene expression using the R package edgeR. Differentially expressed genes were filtered using a false discovery rate of less than 0.05 c, a fold-change value of more than twenty, and an association with overrepresented pathways based on the Reactome pathway over-representation analysis tool. Most of the differentially expressed genes associated with soleus are coded for components of lipid metabolism and unique contractile elements. Differentially expressed genes associated with tibialis anterior encoded mostly for glucose and glycogen metabolic pathway regulatory enzymes and calcium-sensitive contractile components. These gene expression distinctions partly explain the genetic basis for skeletal muscle specialization, and they may help to explain skeletal muscle susceptibility to disease and drugs and further refine tissue engineering approaches.

1. Introduction

Skeletal muscles (SM) provide voluntary body movement and locomotion, posture and body position, energy production, fatty acid (FA) β-oxidation, carbohydrate metabolism, and soft tissue support. Skeletal muscles are composed of myofibers, which are formed by fused myoblasts. Myofibers have defining metabolic and contractile properties.
SM contribute to glucose, lipid, and protein metabolism. SM of healthy and well-fed mammals rely on carbohydrate or glucose metabolism as their primary sources of energy from aerobic respiration, anaerobic respiration, or both. SM can also use FAs as a fuel source via β-oxidation when mammals consume a high-fat diet or are in a starving state. As lipid resources become depleted, SM proteins can be broken down for energy use. Oxidative myofibers utilize mostly aerobic respiration. Glycolytic myofibers use primarily anaerobic respiration. Myofibers that utilize both are called oxidative-glycolytic.
Based on contractile properties, myofibers are defined by which myosin heavy chain is present. Myosin heavy chains are a significant component of the thick filament in a sarcomere. Myosin heavy chains (MyHCs) differentiate myofibers into types I (MYH7), II(A) (MYH2), II(B) (MYH4), and II(D or X) (MYH1) [1]. The four contractile properties relate to the three metabolic distinctions. Type I myofibers are oxidative [2]. Type II(A) myofibers are oxidative-glycolytic [2]. Types II(B) and II(D) myofibers are glycolytic [2]. Type II(D) myofibers combine with type II(A) or (B) to create hybrids called II(AD) and II(DB).
The hindlimb SM Soleus (So) and Tibialis Anterior (Ta) were used to determine the gene expression profile of oxidative and glycolytic myofibers by RNA-Seq in one-month-old mice (postpartum day 30, P30), the age at which isolation of single myofibers was feasible. So and Ta act in opposition, wherein one muscle contracts while the other relaxes. So is one of the superficial muscles that shape the lower hindlimb’s posterior component and is essential for ankle joint movement to plantarflex the foot and maintain the standing posture [3]. The murine So is composed of approximately 50% type II(A), 40% type I, and 10% type II(D) myofibers [4]. Ta is located along the anterior lateral side of the tibia; the muscle acts in the foot’s dorsiflexion and inversion, stabilizes the ankle when the foot hits the ground, and pulls the foot off the ground [3]. The murine Ta is composed of approximately 50–60% type II(B), 30–35% type II(D), and 5–10% II(A) myofibers [4].
Pathway analysis identified molecular players involved in glucose and glycogen metabolism, ion transport, and contraction as being overrepresented among the differentially expressed (DE) genes between the So and Ta biopsies. Additionally, DE transcripts with a high fold change (FC) were identified as being involved in neuromuscular junctions and the immune system. Data from this study broadens our understanding of molecular diversity in SM and potentially differential susceptibility to myopathies.

2. Materials and Methods

2.1. Mice

All animal experiments were performed following the Institutional and National Health and Medical Research Council guidelines. The experimental protocol was approved by the Institutional Animal Care and Use Committee at Oregon State University. Biopsies from Pitx2FL/Z mice were utilized as previously described for this study [5].

2.2. RNA-Sequencing (RNA-Seq) Data Analysis

Five So and two Ta myofiber biopsies were collected at P30. RNA was extracted, sequenced, and analyzed as previously described [5]. The raw RNA-Seq data files for these samples are publicly available through the National Center for Biotechnology Information (NCBI) Sequence Read Archive (SRA) under the reference numbers SRP127367 (So) and SRP145066 (Ta). The reads were aligned using the mouse genome (mm10) and TopHat2. The genomic feature counts were then obtained using HTSeq [5], resulting in a total of 24,421 elements identified. The read counts were normalized to reduce transcript length bias [6]. Longer transcripts appear more prevalent and overshadow smaller transcripts, if not corrected before performing differential expression analysis. Quantitative changes in gene expression levels between the So and Ta groups were discovered using the quasi-likelihood F test from the R package edgeR [7], which resulted in 6123 differentially expressed (DE) genes (p < 0.05 c, False Discovery Rate (FDR)-corrected) between So and Ta. The “c” in 0.05 c refers to the p-value being FDR-corrected in all instances. A smaller set of DE genes based on a log-fold change of ±1.0 or greater was used, resulting in 2481 DE genes that passed our fold change cutoff.

3. Results

3.1. RNA-Seq Data Quality Check

So and Ta were processed and subject to RNA-Seq as described above. Principal component analysis of standardized samples showed a clear separation of So and Ta based on RNA expression (Figure 1a). We use red font to denote increased expression in So relative to Ta and blue for decreased expression in So relative to Ta. A volcano plot (negative log of the gene’s FDR-corrected p-value vs. log2 fold change) of all 24,421 annotated genes showed that the subset of 2481 genes included many highly significant DE genes with increased expression differences in So and Ta (Figure 1b). A heatmap of sample-standardized RNA counts of the final subset of 2481 genes shows clear clustering of myofibers from So and Ta based on gene expression (Figure 1c), as we expect from the biplot in Figure 1a.
Molecular pathways enriched for DE genes from the So and Ta samples were identified using the Reactome database [8]. The 2481 DE genes were matched to molecular pathways, and those pathways were statistically analyzed for over- or under-representation using the Reactome database and its corresponding ReactomePA R package [8,9]. The number of genes analyzed by the Reactome R package decreased to 2198 after removing genes with no Entrez ID. Hypergeometric testing, provided by the ReactomePA R package, identified 19 Reactome pathways as overrepresented with an adjusted p-value of less than 0.05 c. Three pathways had fewer than ten associated genes and were removed from further analysis. The DE genes related to the 16 Reactome pathways were ranked from highest to lowest log2(FC) values (Figure S1). The ranking revealed that the 16 Reactome pathways could be consolidated into broader categories, such as contraction, ion and amino acid transport, and glucose and glycogen metabolism (Figure S1). We identified 119 DE genes with FC ≥ 20, of which 14 were on the Reactome-curated list (Figure 1d). Examination of known gene functions of the remaining 105 high FC genes revealed that 36 could be manually categorized into lipid metabolism, glycogen metabolism, glucose metabolism, or contraction.

3.2. Lipid Metabolism

SM primarily metabolize glucose. Oxidative muscles can also use FAs along with glucose and glycogen as fuel sources. Free FA molecules are transported from the liver via cholesterol transport, which shuttles cholesterol, triglycerides, and protein-bound FAs. Three cholesterol-related genes, Pon1, Tspo2, and Apoa2, were expressed in So by 68, 51, and 27-fold, respectively (Figure 2, Table 1). Pon1 protects against lipid oxidation for high-density lipoprotein (HDL) cholesterol and low-density lipoprotein (LDL) cholesterol [10]. Tspo2 traffics free cholesterol in erythroid cells [11]. Apoa2 regulates steroid concentrations by modulating cholesterol transport [12], the precursor of steroid hormones.
Once free FA molecules pass through the plasma membrane, they enter the mitochondrial matrix to be broken down through β-oxidation. Specific protein channels handle different FA aliphatic tail lengths ranging from greater than 22 to fewer than 5 carbons. Acsm5 was expressed 30-fold higher in So (Figure 2, Table 1). Acsm5 encodes for an enzyme that catalyzes the activation of FAs with aliphatic tails of 6 to 12 carbons by CoA to produce acyl-CoA, the first step in FA metabolism [13]. Dio1, expressed 69-fold higher in So (Figure 2, Table 1), is involved in FA oxidation and oxidative phosphorylation uncoupling primarily in the liver, kidney, and thyroid [14], in addition to altering the thyroid hormone balance of triiodothyronine (T3) and thyroxine (T4) [15]. Thyroid hormone signaling regulates numerous genes involved in skeletal muscle homeostasis, function, and metabolism [16].
Several of the highest FC DE genes in So were involved in lipid, sphingolipid, and derived lipid biosynthesis. FA molecules can be elongated and modified into complex and derived lipid molecules. Tecrl, involved in FA elongation in polyunsaturated FA biosynthesis [17] and myoblast differentiation upon TNF activation [18], was expressed 55-fold higher in So (Figure 2, Table 1). The increased expression of Tecrl could indicate that the So muscle is still developing in P30 mice.
One of the highest FC genes in this study was Cyp4f39 with a 269-fold higher expression in Ta, which encodes for a FA ω-hydroxylase involved in acyl-ceramide synthesis [19] (Figure 2, Table 1). Cyp4f39 is involved in cell surface protection, cell recognition, signaling, membrane transportation, increased rigidity, hydrolyzation of very long-chain FAs, and steroid hormone synthesis [19]. Ceramides have been linked to insulin resistance in type II diabetes [20].

3.3. Glycogen Metabolism

Glycogen metabolism combines two inverse processes: glycogenesis and glycogenolysis. During glycogenesis, glucose or glucose 6-phosphate molecules are converted into glycogen, using 1 adenosine triphosphate (ATP) and 1 ATP equivalent to add one glucose molecule to glycogen (Figure 3a). Only the glycolytic fibers store glycogen as a source of energy. Glycogen is tapped during times of nutritional insufficiency or rapid demand for muscle activity. This study cannot determine the direction of glycogen metabolism alone, only inferring potentially emphasized steps in either So or Ta. The glycogenesis-involved enzymes phosphoglucomutase (Pgm2, Pgm2l1) and glycogen synthase (Gys2) had increased expression in Ta (Figure 3a; Table 2). Phosphoglucomutase reversibly converts glucose 6-phosphate to glucose 1-phosphate. Increased expression of phosphoglucomutase could contribute to either glycogen metabolic processes.
Glycogen synthase is a critical enzyme in glycogenesis that converts glucose-1-phosphate into glycogen with the assistance of pyro-phosphorylase and a branching enzyme [25]. Increased expression of Gys2 in Ta may imply that glycogenesis is the priority for this muscle (Figure 3a; Table 2). Gys2 activity is regulated by dephosphorylation and phosphorylation by protein phosphatase 1 (PP1) and cyclic-adenosine monophosphate (AMP)-dependent protein kinase. PP1 is an enzyme complex with regulatory and catalytic subunits, of which two regulatory subunits, Ppp1r1c and Ppp1r3g, were expressed higher in So (Figure 3b; Table 2). Ppp1r3g regulates glycogen accumulation [26] and, along with the predicted inhibitory regulatory subunit Ppp1r1c, could be controlling glycogenesis in the So myofibers. The increased expression of the PP1 regulatory subunits in So could indicate a tighter control mechanism that prevents glycogen build-up. Epm2a, a factor that regulates PP1, was increased in Ta (Figure 3b; Table 2). Epm2a encodes for laforin, which interacts with PP1 and regulates the dephosphorylation of glycogen to promote glycogen branching [27]. Laforin is also involved in preventing cytotoxicity by protein ubiquitination [28], and it may control glycogen synthesis by ubiquitinating PP1, which would prevent the activation of glycogen synthase.
Glycogenesis must be inhibited to trigger glycogenolysis or glycogen catabolism. Glycogenolysis neither requires nor produces energy when breaking glycogen into glucose molecules. The process is triggered by the accumulation of gluconeogenic precursors, such as lactate or alanine. Glycogen synthase is inactivated to inhibit glycogenesis. The inactivation triggers the activation of phosphorylase kinase, which phosphorylates glycogen phosphorylase. Phosphorylase kinase is an enzyme complex of four subunits [29], three of which (Phka1, Phkg1, Phkb) also had increased expression in Ta (Figure 3c; Table 2). Phosphorylase kinase could be regulating the glycogenolysis in Ta without inhibiting glycogen synthase by ubiquitination. Glycogen phosphorylase pairs with a debranching enzyme to break off glucose 1-phosphate molecules from glycogen; a gene encoding for each (Pygm and Agl, respectively) was expressed higher in Ta (Figure 3a; Table 2).

3.4. Glucose Metabolism

Glucose metabolism consolidates glycolysis and gluconeogenesis enzymatic processes. The former process breaks glucose into an intermediate pyruvate, and the latter is the reverse process. Glycolysis converts glucose into pyruvate in 10 steps, consuming one glucose, two NAD+, two Pi, four ADP, and two ATP molecules to generate two pyruvate, two NADH, two H+, four ADP, and four ATP molecules. The net products are two pyruvate, two ATP, two NADH, and two H+ molecules. Gluconeogenesis converts two pyruvate, four ATP, two GTP, two NADH, two H+, and two H2O molecules into fructose 6-phosphate or glucose 6-phosphate. Overall, gluconeogenesis consumes six ATP, and glycolysis generates a net two ATP. Therefore, the whole glycolysis/gluconeogenesis cycle costs four ATP.
Three critical enzymes regulate the complete glycolysis process, and they are hexokinase, phosphofructokinase-1 (PFK1), and pyruvate kinase. Genes encoding for PFK1 (Pfkm) and pyruvate kinase (Pkm) had increased expression in Ta (Figure 4a; blue font, Table 2). PFK1 converts fructose 6-phosphate to F1,6-BP and is activated by high AMP concentrations and fructose 2,6-bisphosphate (F2,6-BP) [31]. Pyruvate kinase converts phosphoenolpyruvate to pyruvate when activated by high F1,6-BP concentration and is deactivated by an elevated ATP concentration. Four other isoenzymes involved in glycolysis and gluconeogenesis, aldolase (Aldoa, Aldoc), triosephosphate isomerase (Tpi1), phosphoglycerate kinase (Pgk1), and enolase (Eno3 in Ta and Eno2 in So), were increased in Ta (Figure 4a, Table 2).
The indirect regulator of glycolysis is the 6-phosphofructo-2-kinase/fructose 2,6-bisphosphatase (PFK2/FBPase2) complex. The PFK2/FBPase2 complex pulls fructose 6-phosphate out of the glycolytic process and converts it to F2,6-BP [30]. In Ta, three isoenzymes of PFK2/FBPase2, Pfkfb3, Pfkfb1, and Pfkfb4 were increased (Figure 4b; blue font, Table 2).
When blood glucose concentrations are low, glucose can be generated through gluconeogenesis, an 11-step process for converting pyruvate to glucose. The enzymes pyruvate carboxylase, PEP carboxykinase, and fructose 1,6-bisphosphatase (FBPase1) are critical regulatory points in gluconeogenesis. FBPase1 (Fbp2) and PEP carboxykinase (Pck1) were increased in Ta and So, respectively (Figure 4a, Table 2). FBPase1 converts F1,6-BP to fructose 6-phosphate in the presence of a high concentration of AMP and F2,6-BP, while increased citrate levels deactivate it. PEP carboxykinase turns oxaloacetate into phosphoenolpyruvate and is enabled by a high ADP concentration.
All myofibers utilize glycolysis, whether or not they use aerobic and/or anaerobic respiration. Glycolysis occurs under both aerobic and anaerobic conditions because it does not require oxygen. Most of the DE genes associated with glucose metabolism had prominent expression in Ta compared to So. As a muscle that uses anaerobic metabolism to gain fuel, Ta may better utilize the anaerobic pathway due to the increased expression levels of glycolytic and regulatory enzymes.

3.5. Contraction

The Reactome pathway overrepresentation analysis identified contraction as another significant difference between So and Ta. The sliding filament theory explains the molecular basis of muscle contraction, where myosin and actin filaments interact in a particular manner to produce contractile force. Most contraction-related DE genes exposed were associated with the thin filament, thick filament, and Z-line of the sarcomere (Figure 5a, Table 3).
The thin filament is composed of actin subunits. Actc1 had increased expression in Ta (Figure 5b, Table 3). Tropomyosin wraps outside the actin helix to stabilize the strand and rotate the thin filament to expose the binding sites. Two isoforms of tropomyosin were increased in So (Tpm3) [32] or Ta (Tpm1) [33] (Figure 5b, Table 3). The troponin complex is a regulatory element of the actin filament and is associated with calcium binding, inhibitory regulation, and tropomyosin binding. Tropomyosin-binding element encoding genes Tnnt1 [34] and Tnnt2 [35] were increased in So, while Ta had increased expression of Tnnt3 [36] (Figure 5b, Table 3). The calcium and inhibitory troponin units, Tnnc1, Tnni1, and Tnni3, had an increased expression in So (Figure 5b, Table 3). The actin filaments are capped by tropomodulin Tmod1, which was expressed higher in Ta (Figure 5b, Table 3).
The thick filament consists of several strands of heavy myosin chains, which are two myosin heavy chains wrapped around one another with myosin heads at the ends that attach and flex to move along the thin filament. Myosin heavy chain isoforms unique to type I (Myh7) [1] and type II(B) (Myh4) [37] myofibers were increased in So and Ta, respectively (Figure 5c, Table 3).
Myosin light chains and myosin-binding proteins regulate myosin head movement. Myosin light chains are classified into essential (ELCs) and regulatory (RLCs). ELCs stabilize the myosin head, while RLCs stiffen the myosin neck domain. ELCs encoded by Myl6b, Myl3, and Myl4 [2] and the RLCs encoded by Myl2, Myl10 [38] were increased in So (Figure 5c, Table 3) while only Mylpf [39] was increased in Ta. Mybpc2 [40] and Mybph that encode for myosin-binding proteins that assist in movement along the thin filament were increased in Ta (Figure 5c, Table 3).
The thin and thick filaments are held in an antiparallel configuration by a structure called the Z-line at the sarcomere ends (Figure 5d). The titin cap (Tcap) [41] is involved in mechano-electrical links between Z-lines and T-tubules and was expressed higher in So (Figure 5d, Table 3). Actinin (Actn3) [42] is involved in stabilizing the antiparallel formation and was increased in Ta (Figure 5d, Table 3). These structural differences could be altering the electrical activity and, in turn, the contractile properties of the different myofibers in the So and Ta muscles.
Several of the highest FC DE genes were part of the neuromuscular junction (NMJ) apparatus. Genes involved in the presynaptic membrane side included Grp, Slc30a3, Asphd1, Nrxn3 [43], and Nrsn2, which were increased in So, and Fam19a4 [44], Kcnc4 [45], and Fgfbp1 [46], which were increased in Ta (Figure 5e, Table 3). On the postsynaptic membrane side, Tom1 [47], Cntnap4 [48], Wdr72 [49], and Lypd1 were increased in So (Figure 5e; Table 3). Several of these genes are associated with endocytosis and neurotransmitter receptor signaling.

4. Discussion

Two broad-scale patterns emerged along with the three categories of overrepresented molecular pathways when examining the gene expression profiles of So and Ta muscles (Figure S1). First, DE genes involved in lipid, glucose, and glycogen metabolism pathways associated more with one muscle than the other (Table 1 and Table 2). Second, each muscle expressed unique genes encoding for similarly functioning isoenzymes, particularly for genes involved in the contraction and ion transport categories (Table 3 and Table S1). We put these differences into context below, starting with the pathways that differentiate the muscle types. DE genes with the highest FC difference were not prevalent among genes associated with the overrepresented pathways. Therefore, DE genes with a greater than 20-fold expression difference were included for further investigation. Similarly seen in a previous study [54], we found that approximately 10% of DE transcripts were mainly related to fatty acid metabolism, structural components, and neuromuscular junction assembly.

4.1. Lipid Metabolism

FA lipids are an alternative energy source during fasting, starvation, and endurance exercise in oxidative muscles. Both So and Ta muscles contain oxidative myofibers, and both should express genes associated with FA metabolism (Figure 2). However, DE genes explicitly associated with FA transport and catabolism were increased in expression in So, which contains a higher percentage of oxidative fibers in its overall myofiber composition. In contrast, the 269-fold increased expression of Cyp4f39 in Ta suggests that lipid catabolism to generate ceramides is more active in Ta. Ceramide accumulation has been linked to insulin resistance in type II diabetes [55]. The increased Cyp4f39 expression may allow Ta to utilize ceramides to shift from being insulin-sensitive and from facilitating efficient glucose uptake [20] to becoming insulin-resistant. Alternatively, increased expression of Cyp4f39 and presumed increased ceramide levels in Ta might be associated with differential apoptosis, cell cycling, or autophagy.

4.2. Glycogen and Glycosaminoglycan Metabolism

Glycogen acts as an energy reserve for myofibers by storing excess glucose molecules that are then utilized when blood glucose levels are low. Glycolytic or oxidative-glycolytic myofibers retain glycogen because they rely heavily on glucose as energy, and Ta has a high percentage of those myofibers; therefore, the increased expressions of the glycogen metabolic enzyme genes Gys2, Pygm, Phka1, Phkg1, and Phkb in Ta support the known metabolic composition of the muscle relative to So (Figure 3a; Table 2).
Hydrolysis of muscle glycogen to glucose occurs in lysosomes that engulf glycogen granules. The lysosome-associated DE genes, Slc35d3, Atp6v0d2, and Galns, exhibited increased expression in So, suggesting that lysosomal-related organelles in So and Ta myofibers perform slightly different functions (Table S1). Slc35d3 is associated with the biosynthesis of platelet-dense granules [56]. Platelet-dense granules contain high concentrations of calcium, adenine nucleotides, pyrophosphate, and polyphosphate molecules that enhance autophagy in lysosome-related organelles. Atp6v0d2 is a part of vacuolar ATPases involved in proton translocation into vacuoles, lysosomes, or the Golgi apparatus to lower the pH [57]. The presence of Atp6v0d2 suggested that lysosome-related organelles in So need to reduce their pH levels routinely. If lysosome-related organelles in So myofibers have an enhanced autophagy, the pH levels within those organelles would fluctuate as the cell engulfs and metabolizes more molecules from the extracellular space and may require proton pumps.
Galns is involved in glycosaminoglycan biosynthesis. Glycosaminoglycans are cell surface proteins with branches of sugars projecting into the extracellular matrix to support cell identity, adhesion, and growth. Galns encodes for a protein that breaks keratan sulfate off glycosaminoglycans [58] in lysosomes. Stab2, B3gat1, and Cspg5 are also involved in glycosaminoglycan synthesis, indicating that So and Ta synthesize distinct glycosaminoglycans associated with cell adhesion and proliferation (Tables S1 and S2). Stab2 is involved in the endocytosis of metabolic waste products, including circulating hyaluronic acid (HA) that promotes cell proliferation once the plasma concentration decreases. B3gat1 is involved in generating HNK-1 carbohydrate (CD57) cell surface epitope associated with cell adhesion [59]. Cspg5 is involved in chondroitin sulfate synthesis, another molecule added to proteoglycans and connected to cell adhesion, growth, migration, and receptor binding in the central nervous system (Table S1). The high Cspg5 expression in So suggests that neurons associated with the So muscle are marked differently.

4.3. Glucose Metabolism

Genes associated with glucose metabolism were increased in Ta more than So despite both muscles utilizing glycolysis. Yet, Eno2 and Pck1 were increased in expression in So compared to the 11 isoenzymes in Ta, perhaps due to the higher percentage of type II(B) glycolytic myofibers in Ta (Figure 4, Table 2). Glycolysis-related enzymes encoded by Pfkm, Pkm Pfkfb3, Pfkfb1, and Pfkfb4 were increased in Ta (Figure 4b, Table 2). In contrast, two gluconeogenesis regulatory enzymes that regulate the reversion of pyruvate to glucose in skeletal muscle encoded by Fbp2 and Pck1 were expressed higher in Ta and So, respectively. Gluconeogenesis in Ta might not be as tightly controlled across several steps as glycolysis seems to be based on the 3-to-1 regulatory enzyme ratio. Gluconeogenesis may be controlled by the PFK2/FBPase2 complex instead of relying on the energy-costing enzyme PEP carboxylase seemingly utilized by So. Pfkfb4 had a higher DE than Pfkfb3. Pfkfb4 opposes Pfkfb3 as it redirects glucose to the pentose phosphate pathway to promote the detoxification of reactive oxygen species and lipid and nucleotide biosynthesis. Fbp2, which was also increased in expression in Ta, converts F1,6-BP back into fructose 6-phosphate. These combined findings may indicate that Ta initiates gluconeogenesis at the PFK2/FBPase2 step of the glucose metabolism process to reduce energy loss.

4.4. Contraction

Several components that contribute to sarcomere formation were identified as DE between So and Ta, including MyHCs, myosin light chains, troponin subunits, and NMJ-associated proteins (Figure 5a–e). The presence of type I (Myh7) and type II(B) (Myh4)-associated MyHC isoforms among DE genes provided proof of the concept for this analysis. The fold change in differential expression of Myh7 and Myh4, 34.7- and 57.7-fold, respectively (Table 3), reflected the appropriate percentage of types I and II(B) myofibers in each muscle. So and Ta both have type II(A) and II(D) myofibers, so differential expression of their specific MyHCs was predictably absent from our analysis. Our data also support a previous analysis correlating Myh4 and Myh7 with other genes associated with myofiber types I and II(B) [54], respectively (Table 3 and Tables S1 and S2). The myosin ELCs and RLCs mainly correlated to the appropriate muscle, except for Myl4 and Myl6 (Figure 5c, Table 3). The troponin complex subunits encoded by Tnnt1, Tnni1, and Tnnc1 were unique and highly expressed in So (Figure 5b, Table 3), suggesting that the troponin complex in So might be modified relative to Ta to keep the actin filament in an open position longer, facilitating longer durations of contraction.
Factors affecting NMJ specification and calcium handling have been theorized to be among the non-myofibrillar, muscle-specific systems that allow for increased maximal isometric stress after P28 [60]. Synaptic modeling affects the motility or twitching of the skeletal muscle. Muscle identity is determined before innervation, yet innervation maintains muscle identity. During regeneration, only innervated So muscles can upregulate slow isoform mRNAs [61]. Yet, muscles do not completely change into another phenotype when neural cues are altered [62]. Intrinsic signaling from nerves and other sources, and environmental cues, dictate adult muscle phenotype [61]. Genes are associated with neurotransmission and synaptic vesicles [43] at the presynaptic side of NMJs, in addition to clathrin-mediated endocytosis [47,49], and cell adhesion [48] of the postsynaptic side was highly expressed in So. The release and retrieval of the synaptic vesicles and their contents may play an important role in the slow-twitch mechanism. On the contrary, Ta had increased expression of genes associated with the presynaptic side of NMJs only, which modulate neuronal excitability [44], broaden the action potential [45], and slow down the degeneration of NMJs [46]. These genes may explain how Ta and its associated neurons handle the stimulation necessary for rapid contraction. The DE genes related to the ion transport pathways may be associated with the postsynaptic side of NMJs in Ta myofibers, such as the voltage-dependent calcium channel subunits (Table S1). However, all of this information and conjuncture is based on a small number of differentially expressed genes.

5. Conclusions

Two hindlimb muscles, So and Ta, which perform distinct locomotion functions, are characterized by unique metabolic and contractile properties. This focused investigation identified several genes uniquely tied to either So or Ta when these two muscle groups are compared. Approximately 10% of the mouse transcriptome was differentially expressed, including genes involved in lipid, glucose, and glycogen metabolism. Very highly differentially expressed genes such as Cyp4f39 highlight previously undescribed potential differences in fatty acids between the SM that may underlie susceptibility to pathologies. The high DE of Pfkm and Gys2 point to known differences in glucose and glycogen metabolism between the SM but add molecular details to the picture. Among contraction-related genes, RLCs, ELCs, and troponin complex subunits were higher expressed in So. DE of genes involved in NMJs suggests that electrochemical signaling is different between the two muscle groups, supporting the contractile differences.
Overall, this detailed study into the gene expression differences between these two muscle groups adds to the abundance of data on molecular differences among SM. Some similar results were found in single nuclei studies, which also identified unique expressions among individual skeletal muscle groups [54,63]. This study also highlights the continuing need to maintain and update gene and pathway ontology databases, especially adding more information on gene expression in different cell types.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/genes12111718/s1, Figure S1: Evaluation of Reactome overrepresented pathway analysis, Table S1. Differentially Expressed Genes in So and Ta involved in Ion Transport and Glycosaminoglycan biosynthesis pathways, Table S2. Differentially Expressed Genes in So and Ta organized by Immune Response, Signaling, and Cellular function pathways.

Author Contributions

C.K. conceived the study and designed experiments; C.-N.C. performed all experiments; K.B. provided biocomputational mentorship and developed the computational pipeline for RNA-Seq analysis; S.I.A. performed RNA-Seq analysis; T.M.F. provided project mentorship; S.I.A., T.M.F. and C.K. wrote the manuscript with contributions from all authors. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to thank NIH for the grant NIH-NIAMS 1R01AR054406 (CK) and the College of Pharmacy at Oregon State University for funding this investigation.

Institutional Review Board Statement

All experimental protocols were approved by the Environmental Health and Safety Committee at Oregon State University. All animal experiments were performed in accordance with institutional and National Health and Medical Research Council guidelines. The experimental protocol was approved by the Institutional Animal Care and Use Committee at Oregon State University.

Informed Consent Statement

Not applicable.

Data Availability Statement

FASTQ sequences were deposited to the NCBI gene expression omnibus (GEO) sequence read archive (SRA) using the accessions SRP127367 and SRP145066.

Acknowledgments

The authors want to thank Arun Singh and Thai Tran for their assistance during the data analysis stages.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Schiaffino, S.; Rossi, A.; Smerdu, V.; Leinwand, L.; Reggiani, C. Developmental Myosins: Expression Patterns and Functional Significance. Skelet. Muscle 2015, 5, 22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Schiaffino, S.; Reggiani, C. Fiber Types in Mammalian Skeletal Muscles. Physiol. Rev. 2011, 91, 1447–1531. [Google Scholar] [CrossRef] [Green Version]
  3. Fujita, N.; Fujimoto, T.; Tasaki, H.; Arakawa, T.; Matsubara, T.; Miki, A. Influence of Muscle Length on Muscle Atrophy in the Mouse Tibialis Anterior and Soleus Muscles. Biomed. Res. 2009, 30, 39–45. [Google Scholar] [CrossRef] [Green Version]
  4. Kammoun, M.; Cassar-Malek, I.; Meunier, B.; Picard, B. A Simplified Immunohistochemical Classification of Skeletal Muscle Fibres in Mouse. Eur. J. Histochem. 2014, 58, 2254. [Google Scholar] [CrossRef] [Green Version]
  5. Singh, A.J.; Chang, C.-N.; Ma, H.-Y.; Ramsey, S.A.; Filtz, T.M.; Kioussi, C. FACS-Seq Analysis of Pax3-Derived Cells Identifies Non-Myogenic Lineages in the Embryonic Forelimb. Sci. Rep. 2018, 8, 7670. [Google Scholar] [CrossRef] [PubMed]
  6. Gao, L.; Fang, Z.; Zhang, K.; Zhi, D.; Cui, X. Length Bias Correction for RNA-Seq Data in Gene Set Analyses. Bioinformatics 2011, 27, 662–669. [Google Scholar] [CrossRef] [Green Version]
  7. Robinson, M.D.; McCarthy, D.J.; Smyth, G.K. EdgeR: A Bioconductor Package for Differential Expression Analysis of Digital Gene Expression Data. Bioinformatics 2009, 26, 139–140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Croft, D.; O’Kelly, G.; Wu, G.; Haw, R.; Gillespie, M.; Matthews, L.; Caudy, M.; Garapati, P.; Gopinath, G.; Jassal, B.; et al. Reactome: A Database of Reactions, Pathways and Biological Processes. Nucleic Acids Res. 2011, 39, D691–D697. [Google Scholar] [CrossRef]
  9. Yu, G.; He, Q.-Y. ReactomePA: An R/Bioconductor Package for Reactome Pathway Analysis and Visualization. Mol. BioSyst. 2016, 12, 477–479. [Google Scholar] [CrossRef]
  10. Shih, D.M.; Gu, L.; Xia, Y.R.; Navab, M.; Li, W.F.; Hama, S.; Castellani, L.W.; Furlong, C.E.; Costa, L.G.; Fogelman, A.M.; et al. Mice Lacking Serum Paraoxonase Are Susceptible to Organophosphate Toxicity and Atherosclerosis. Nature 1998, 394, 284–287. [Google Scholar] [CrossRef]
  11. Fan, J.; Rone, M.; Papadopoulos, V. Translocator Protein 2 Is Involved in Cholesterol Redistribution during Erythropoiesis. J. Biol. Chem. 2009, 284, 30484–30497. [Google Scholar] [CrossRef] [Green Version]
  12. Yang, M.; Liu, Y.; Dai, J.; Li, L.; Ding, X.; Xu, Z.; Mori, M.; Miyahara, H.; Sawashita, J.; Higuchi, K. Apolipoprotein A-II Induces Acute-Phase Response Associated AA Amyloidosis in Mice through Conformational Changes of Plasma Lipoprotein Structure. Sci. Rep. 2018, 8, 5620. [Google Scholar] [CrossRef]
  13. Puig-Oliveras, A.; Revilla, M.; Castelló, A.; Fernández, A.I.; Folch, J.M.; Ballester, M. Expression-Based GWAS Identifies Variants, Gene Interactions and Key Regulators Affecting Intramuscular Fatty Acid Content and Composition in Porcine Meat. Sci. Rep. 2016, 6, 31803. [Google Scholar] [CrossRef]
  14. Bianco, A.C.; Salvatore, D.; Gereben, B.; Berry, M.J.; Larsen, P.R. Biochemistry, Cellular and Molecular Biology, and Physiological Roles of the Iodothyronine Selenodeiodinases. Endocr. Rev. 2002, 23, 38–89. [Google Scholar] [CrossRef]
  15. Panicker, V.; Cluett, C.; Shields, B.; Murray, A.; Parnell, K.; Perry, J.; Weedon, M.; Singleton, A.; Hernandez, D.; Evans, J.; et al. A Common Variation in Deiodinase 1 Gene DIO1 Is Associated with the Relative Levels of Free Thyroxine and Triiodothyronine. J. Clin. Endocrinol. Metab. 2008, 93, 3075–3081. [Google Scholar] [CrossRef]
  16. Salvatore, D.; Simonides, W.S.; Dentice, M.; Zavacki, A.M.; Larsen, P.R. Thyroid Hormones and Skeletal Muscle-New Insights and Potential Implications. Nat. Rev. Endocrinol. 2014, 10, 206–214. [Google Scholar] [CrossRef] [Green Version]
  17. Devalla, H.D.; Gélinas, R.; Aburawi, E.H.; Beqqali, A.; Goyette, P.; Freund, C.; Chaix, M.; Tadros, R.; Jiang, H.; le Béchec, A.; et al. TECRL, a New Life-threatening Inherited Arrhythmia Gene Associated with Overlapping Clinical Features of Both LQTS and CPVT. EMBO Mol. Med. 2016, 8, 1390–1408. [Google Scholar] [CrossRef] [PubMed]
  18. Meyer, S.U.; Sass, S.; Mueller, N.S.; Krebs, S.; Bauersachs, S.; Kaiser, S.; Blum, H.; Thirion, C.; Krause, S.; Theis, F.J.; et al. Integrative Analysis of MicroRNA and MRNA Data Reveals an Orchestrated Function of MicroRNAs in Skeletal Myocyte Differentiation in Response to TNF-α or IGF1. PLoS ONE 2015, 10, e0135284. [Google Scholar] [CrossRef] [PubMed]
  19. Miyamoto, M.; Itoh, N.; Sawai, M.; Sassa, T.; Kihara, A. Severe Skin Permeability Barrier Dysfunction in Knockout Mice Deficient in a Fatty Acid ω-Hydroxylase Crucial to Acylceramide Production. J. Investig. Dermatol. 2020, 140, 319–326. [Google Scholar] [CrossRef] [PubMed]
  20. Sokolowska, E.; Blachnio-Zabielska, A. The Role of Ceramides in Insulin Resistance. Front. Endocrinol. 2019, 10, 577. [Google Scholar] [CrossRef] [Green Version]
  21. Piekorz, R.P.; Gingras, S.; Hoffmeyer, A.; Ihle, J.N.; Weinstein, Y. Regulation of Progesterone Levels during Pregnancy and Parturition by Signal Transducer and Activator of Transcription 5 and 20α-Hydroxysteroid Dehydrogenase. Mol. Endocrinol. 2005, 19, 431–440. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  22. Parsons, W.H.; Kolar, M.J.; Kamat, S.S.; Cognetta, A.B., III; Hulce, J.J.; Saez, E.; Kahn, B.B.; Saghatelian, A.; Cravatt, B.F. AIG1 and ADTRP Are Atypical Integral Membrane Hydrolases That Degrade Bioactive FAHFAs. Nat. Chem. Biol. 2016, 12, 367–372. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Terao, M.; Romão, M.J.; Leimkühler, S.; Bolis, M.; Fratelli, M.; Coelho, C.; Santos-Silva, T.; Garattini, E. Structure and Function of Mammalian Aldehyde Oxidases. Arch. Toxicol. 2016, 90, 753–780. [Google Scholar] [CrossRef] [PubMed]
  24. Alnouti, Y.; Klaassen, C.D. Tissue Distribution, Ontogeny, and Regulation of Aldehyde Dehydrogenase (Aldh) Enzymes MRNA by Prototypical Microsomal Enzyme Inducers in Mice. Toxicol. Sci. 2008, 101, 51–64. [Google Scholar] [CrossRef] [Green Version]
  25. Zhang, Y.; Gu, J.; Wang, L.; Zhao, Z.; Pan, Y.; Chen, Y. Ablation of PPP1R3G Reduces Glycogen Deposition and Mitigates High-Fat Diet Induced Obesity. Mol. Cell. Endocrinol. 2017, 439, 133–140. [Google Scholar] [CrossRef]
  26. Irimia, J.M.; Meyer, C.M.; Peper, C.L.; Zhai, L.; Bock, C.B.; Previs, S.F.; McGuinness, O.P.; DePaoli-Roach, A.; Roach, P.J. Impaired Glucose Tolerance and Predisposition to the Fasted State in Liver Glycogen Synthase Knock-out Mice. J. Biol. Chem. 2010, 285, 12851–12861. [Google Scholar] [CrossRef] [Green Version]
  27. Ganesh, S.; Lal Agarwala, K.; Ueda, K.; Akagi, T.; Shoda, K.; Usui, T.; Hashikawa, T.; Osada, H.; Delgado-Escueta, A.V.; Yamakawa, K. Defective in the Progressive Myoclonus Epilepsy of Lafora Type, Is a Dual-Specificity Phosphatase Associated with Polyribosomes. Hum. Mol. Genet. 2000, 9, 2251–2261. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Garyali, P.; Siwach, P.; Singh, P.K.; Puri, R.; Mittal, S.; Sengupta, S.; Parihar, R.; Ganesh, S. The Malin-Laforin Complex Suppresses the Cellular Toxicity of Misfolded Proteins by Promoting Their Degradation through the Ubiquitin-Proteasome System. Hum. Mol. Genet. 2009, 18, 688–700. [Google Scholar] [CrossRef]
  29. Roach, P.J. Glycogen and Its Metabolism. Curr. Mol. Med. 2002, 2, 101–120. [Google Scholar] [CrossRef]
  30. Chesney, J.; Telang, S.; Yalcin, A.; Clem, A.; Wallis, N.; Bucala, R. Targeted Disruption of Inducible 6-Phosphofructo-2-Kinase Results in Embryonic Lethality. Biochem. Biophys. Res. Commun. 2005, 331, 139–146. [Google Scholar] [CrossRef]
  31. García, M.; Pujol, A.; Ruzo, A.; Riu, E.; Ruberte, J.; Arbós, A.; Serafín, A.; Albella, B.; Felíu, J.E.; Bosch, F. Phosphofructo-1-Kinase Deficiency Leads to a Severe Cardiac and Hematological Disorder in Addition to Skeletal Muscle Glycogenosis. PLoS Genet. 2009, 5, 615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Hook, J.; Lemckert, F.; Qin, H.; Schevzov, G.; Gunning, P. Gamma Tropomyosin Gene Products Are Required for Embryonic Development. Mol. Cell. Biol. 2004, 24, 2318–2323. [Google Scholar] [CrossRef] [Green Version]
  33. Rethinasamy, P.; Muthuchamy, M.; Hewett, T.; Boivin, G.; Wolska, B.M.; Evans, C.; Solaro, R.J.; Wieczorek, D.F. Molecular and Physiological Effects of Alpha-Tropomyosin Ablation in the Mouse. Circ. Res. 1998, 82, 116–123. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Wei, B.; Lu, Y.; Jin, J.-P. Deficiency of Slow Skeletal Muscle Troponin T Causes Atrophy of Type I Slow Fibres and Decreases Tolerance to Fatigue. J. Physiol. 2014, 592, 1367–1380. [Google Scholar] [CrossRef] [Green Version]
  35. Wei, B.; Jin, J.-P. TNNT1, TNNT2, and TNNT3: Isoform Genes, Regulation, and Structure–Function Relationships. Gene 2016, 582, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Ju, Y.; Li, J.; Xie, C.; Ritchlin, C.T.; Xing, L.; Hilton, M.J.; Schwarz, E.M. Troponin T3 Expression in Skeletal and Smooth Muscle Is Required for Growth and Postnatal Survival: Characterization of Tnnt3tm2a(KOMP)Wtsi Mice. Genesis 2013, 51, 667–675. [Google Scholar] [CrossRef] [Green Version]
  37. Acakpo-Satchivi, L.J.R.; Edelmann, W.; Sartorius, C.; Lu, B.D.; Wahr, P.A.; Watkins, S.C.; Metzger, J.M.; Leinwand, L.; Kucherlapati, R. Growth and Muscle Defects in Mice Lacking Adult Myosin Heavy Chain Genes. J. Cell Biol. 1997, 139, 1219–1229. [Google Scholar] [CrossRef] [PubMed]
  38. Woolner, S.; O’Brien, L.L.; Wiese, C.; Bement, W.M. Myosin-10 and Actin Filaments Are Essential for Mitotic Spindle Function. J. Cell Biol. 2008, 182, 77–88. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  39. Wang, Y.; Szczesna-Cordary, D.; Craig, R.; Diaz-Perez, Z.; Guzman, G.; Miller, T.; Potter, J.D. Fast Skeletal Muscle Regulatory Light Chain Is Required for Fast and Slow Skeletal Muscle Development. FASEB J. 2007, 21, 2205–2214. [Google Scholar] [CrossRef] [PubMed]
  40. McNamara, J.W.; Sadayappan, S. Skeletal Myosin Binding Protein-C: An Increasingly Important Regulator of Striated Muscle Physiology. Arch. Biochem. Biophys. 2018, 660, 121–128. [Google Scholar] [CrossRef]
  41. Markert, C.D.; Meaney, M.P.; Voelker, K.A.; Grange, R.W.; Dalley, H.W.; Cann, J.K.; Ahmed, M.; Bishwokarma, B.; Walker, S.J.; Yu, S.X.; et al. Functional Muscle Analysis of the Tcap Knockout Mouse. Hum. Mol. Genet. 2010, 19, 2268–2283. [Google Scholar] [CrossRef] [PubMed]
  42. MacArthur, D.G.; Seto, J.T.; Raftery, J.M.; Quinlan, K.G.; Huttley, G.A.; Hook, J.W.; Lemckert, F.A.; Kee, A.J.; Edwards, M.R.; Berman, Y.; et al. Loss of ACTN3 Gene Function Alters Mouse Muscle Metabolism and Shows Evidence of Positive Selection in Humans. Nat. Genet. 2007, 39, 1261–1265. [Google Scholar] [CrossRef] [PubMed]
  43. Missler, M.; Zhang, W.; Rohlmann, A.; Kattenstroth, G.; Hammer, R.E.; Gottmann, K.; Südhof, T.C. Alpha-Neurexins Couple Ca2+ Channels to Synaptic Vesicle Exocytosis. Nature 2003, 423, 939–948. [Google Scholar] [CrossRef]
  44. Delfini, M.-C.; Mantilleri, A.; Gaillard, S.; Hao, J.; Reynders, A.; Malapert, P.; Alonso, S.; François, A.; Barrere, C.; Seal, R.; et al. TAFA4, a Chemokine-like Protein, Modulates Injury-Induced Mechanical and Chemical Pain Hypersensitivity in Mice. Cell Rep. 2013, 5, 378–388. [Google Scholar] [CrossRef] [Green Version]
  45. Rowan, M.J.M.; Christie, J.M. Rapid State-Dependent Alteration in Kv3 Channel Availability Drives Flexible Synaptic Signaling Dependent on Somatic Subthreshold Depolarization. Cell Rep. 2017, 18, 2018–2029. [Google Scholar] [CrossRef]
  46. Taetzsch, T.; Tenga, M.J.; Valdez, G. Muscle Fibers Secrete FGFBP1 to Slow Degeneration of Neuromuscular Synapses during Aging and Progression of ALS. J. Neurosci. 2017, 37, 70–82. [Google Scholar] [CrossRef] [Green Version]
  47. Seet, L.-F.; Hong, W. Endofin Recruits Clathrin to Early Endosomes via TOM1. J. Cell. Sci. 2005, 118, 575–587. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Karayannis, T.; Au, E.; Patel, J.C.; Kruglikov, I.; Markx, S.; Delorme, R.; Heron, D.; Salomon, D.; Glessner, J.; Restituito, S.; et al. Cntnap4 Differentially Contributes to GABAergic and Dopaminergic Synaptic Transmission. Nature 2014, 511, 236–240. [Google Scholar] [CrossRef] [Green Version]
  49. Wang, S.; Hu, Y.; Yang, J.; Smith, C.E.; Nunez, S.M.; Richardson, A.S.; Pal, S.; Samann, A.C.; Hu, J.; Simmer, J.P. Critical Roles for WDR72 in Calcium Transport and Matrix Protein Removal during Enamel Maturation. Mol. Genet. Genom. Med. 2015, 3, 302–319. [Google Scholar] [CrossRef]
  50. Nishi, M.; Takeshima, H.; Hashimoto, K.; Kano, M.; Hashimoto, K.; Kano, M.; Kuriyama, K.; Komazaki, S.; Shibata, S. Motor Discoordination in Mutant Mice Lacking Junctophilin Type 3. Biochem. Biophys. Res. Commun. 2002, 292, 318–324. [Google Scholar] [CrossRef] [PubMed]
  51. Zhang, B.; Kaufman, R.J.; Ginsburg, D. LMAN1 and MCFD2 Form a Cargo Receptor Complex and Interact with Coagulation Factor VIII in the Early Secretory Pathway. J. Biol. Chem. 2005, 280, 25881–25886. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Vagnozzi, R.; Gatto, G.; Kallander, L.; Hoffman, N.; Mallilankaraman, K.; Ballard, V.; Lawhorn, B.; Stoy, P.; Philp, J.; Graves, A.; et al. Inhibition of the Cardiomyocyte-Specific Kinase TNNI3K Limits Oxidative Stress, Injury, and Adverse Remodeling in the Ischemic Heart. Sci. Transl. Med. 2013, 5, 207ra141. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Gallicano, G.I.; Kouklis, P.; Bauer, C.; Yin, M.; Vasioukhin, V.; Degenstein, L.; Fuchs, E. Desmoplakin Is Required Early in Development for Assembly of Desmosomes and Cytoskeletal Linkage. J. Cell Biol. 1998, 143, 2009–2022. [Google Scholar] [CrossRef]
  54. Terry, E.E.; Zhang, X.; Hoffmann, C.; Hughes, L.D.; Lewis, S.A.; Li, J.; Wallace, M.J.; Riley, L.; Douglas, C.M.; Gutierrez-Montreal, M.A.; et al. Transcriptional Profiling Reveals Extraordinary Diversity among Skeletal Muscle Tissues. eLife 2018, 7, e34613. [Google Scholar] [CrossRef]
  55. Błachnio-Zabielska, A.; Zabielski, P.; Baranowski, M.; Gorski, J. Effects of Streptozotocin-Induced Diabetes and Elevation of Plasma FFA on Ceramide Metabolism in Rat Skeletal Muscle. Horm. Metab. Res. 2010, 42, 1–7. [Google Scholar] [CrossRef]
  56. Chintala, S.; Tan, J.; Gautam, R.; Rusiniak, M.E.; Guo, X.; Li, W.; Gahl, W.A.; Huizing, M.; Spritz, R.A.; Hutton, S.; et al. The Slc35d3 Gene, Encoding an Orphan Nucleotide Sugar Transporter, Regulates Platelet-Dense Granules. Blood 2007, 109, 1533–1540. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Kinnear, N.P.; Boittin, F.-X.; Thomas, J.M.; Galione, A.; Evans, A.M. Lysosome-Sarcoplasmic Reticulum Junctions A Trigger Zone for Calcium Signaling by Nicotinic Acid Adenine Dinucleotide Phosphate and Endothelin-1. J. Biol. Chem. 2004, 279, 54319–54326. [Google Scholar] [CrossRef] [Green Version]
  58. Tomatsu, S.; Orii, K.O.; Vogler, C.; Nakayama, J.; Levy, B.; Grubb, J.H.; Gutierrez, M.A.; Shim, S.; Yamaguchi, S.; Nishioka, T.; et al. Mouse Model on N-Acetylgalactosamine-6-Sulfate Sulfatase Deficiency (Galns-/-) Produced by Targeted Disruption of the Gene Defective in Morquio a Disease. Hum. Mol. Genet. 2003, 12, 3349–3358. [Google Scholar] [CrossRef] [Green Version]
  59. Yamamoto, S.; Oka, S.; Inoue, M.; Shimuta, M.; Manabe, T.; Takahashi, H.; Miyamoto, M.; Asano, M.; Sakagami, J.; Sudo, K.; et al. Mice Deficient in Nervous System-Specific Carbohydrate Epitope HNK-1 Exhibit Impaired Synaptic Plasticity and Spatial Learning. J. Biol. Chem. 2002, 277, 27227–27231. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Gokhin, D.S.; Ward, S.R.; Bremner, S.N.; Lieber, R.L. Quantitative Analysis of Neonatal Skeletal Muscle Functional Improvement in the Mouse. J. Exp. Biol. 2008, 211, 837–843. [Google Scholar] [CrossRef] [Green Version]
  61. Esser, K.; Gunning, P.; Hardeman, E. Nerve-Dependent and -Independent Patterns of MRNA Expression in Regenerating Skeletal Muscle. Dev. Biol. 1993, 159, 173–183. [Google Scholar] [CrossRef] [PubMed]
  62. Pette, D.; Düsterhöft, S. Altered Gene Expression in Fast-Twitch Muscle Induced by Chronic Low-Frequency Stimulation. Am. J. Physiol. 1992, 262, R333–R338. [Google Scholar] [CrossRef] [PubMed]
  63. Dos Santos, M.; Backer, S.; Saintpierre, B.; Izac, B.; Andrieu, M.; Letourneur, F.; Relaix, F.; Sotiropoulos, A.; Maire, P. Single-Nucleus RNA-Seq and FISH Identify Coordinated Transcriptional Activity in Mammalian Myofibers. Nat. Commun. 2020, 11, 5102. [Google Scholar] [CrossRef] [PubMed]
Genes 12 01718 g001 550
Figure 1. RNA-Seq workflow and data quality assessments. (a) Biplot of sample-standardized RNA counts (all genes) from the So (red) and Ta (blue) samples. Principal Component 1 (PC1) constitutes 70% of the variability among samples, and Principal Component 2 (PC2) constitutes 11% of the variability. (b) Volcano plot of all 24,421 annotated genes. Genes with an FDR-corrected p-value (p_c) less than 0.05 c are gray dots, and those with a p-value higher than 0.05 c are black. A select set of highly significant DE genes with large fold changes for So (log FC > 3, −log(p_c) > 50) are shown in red and for Ta (log FC < −3, −log(p_c) > 75) in blue. (c) Heatmap of z-scored log-transformed read counts. Positive z-scores are in orange, and negative z-scores are in blue. Each sample is a column, and each row represents a statistically significant DE gene from the 2481 set. (d) Workflow used to reduce the number of genes for further analysis. A subset of the 2481 DE genes with extremely high differential expression was selected; 119 DE genes had a log-fold expression greater than twenty (purple). Separately, 97 of the 2481 genes mapped to sixteen statistically overrepresented molecular pathways using the Reactome database R package (p < 0.05 c) and had a log-fold expression change >3 (orange).
Figure 1. RNA-Seq workflow and data quality assessments. (a) Biplot of sample-standardized RNA counts (all genes) from the So (red) and Ta (blue) samples. Principal Component 1 (PC1) constitutes 70% of the variability among samples, and Principal Component 2 (PC2) constitutes 11% of the variability. (b) Volcano plot of all 24,421 annotated genes. Genes with an FDR-corrected p-value (p_c) less than 0.05 c are gray dots, and those with a p-value higher than 0.05 c are black. A select set of highly significant DE genes with large fold changes for So (log FC > 3, −log(p_c) > 50) are shown in red and for Ta (log FC < −3, −log(p_c) > 75) in blue. (c) Heatmap of z-scored log-transformed read counts. Positive z-scores are in orange, and negative z-scores are in blue. Each sample is a column, and each row represents a statistically significant DE gene from the 2481 set. (d) Workflow used to reduce the number of genes for further analysis. A subset of the 2481 DE genes with extremely high differential expression was selected; 119 DE genes had a log-fold expression greater than twenty (purple). Separately, 97 of the 2481 genes mapped to sixteen statistically overrepresented molecular pathways using the Reactome database R package (p < 0.05 c) and had a log-fold expression change >3 (orange).
Genes 12 01718 g001
Genes 12 01718 g002 550
Figure 2. Lipid metabolism associated genes in So and Ta myofibers. Substrates and products of the enzymatic pathways for cholesterol transport, β-oxidation, FA elongation, and ceramide de-novo synthesis are shown based on their cellular location. The gene transcripts involved an enzymatic process and had increased expression in So or Ta, labeled in red or blue font, respectively. Font size correlates to relative FC.
Figure 2. Lipid metabolism associated genes in So and Ta myofibers. Substrates and products of the enzymatic pathways for cholesterol transport, β-oxidation, FA elongation, and ceramide de-novo synthesis are shown based on their cellular location. The gene transcripts involved an enzymatic process and had increased expression in So or Ta, labeled in red or blue font, respectively. Font size correlates to relative FC.
Genes 12 01718 g002
Genes 12 01718 g003 550
Figure 3. Glycogen metabolism associated genes in So and Ta myofibers. (a) Glucose is converted into glycogen during glycogenesis. Glycogenolysis converts glycogen back into glucose molecules. Represented are magnesium ions (+Mg2+), adenosine diphosphate (ADP), uridine-5′-triphosphate (UTP), pyrophosphate (PPi), and inorganic phosphate (Pi). (b) Two subunits of the protein phosphorylase 1 were increased in So. Represented are an encircled ‘P’ indicating phosphorylation and uridine-diphosphate (UDP). (c) Three isoform subunits of phosphorylase kinase were increased in Ta. Represented are the addition (+Ca2+) and the removal of calcium ions (−Ca2+). Substrates and products of the enzymatic steps in glycogen metabolism include the interconversion of inactive to active states of regulatory enzymes. Enzymes involved at specific stages are encircled and adjacent to the corresponding reaction arrow. Gene symbols located underneath enzyme names in red or blue font were increased in So or Ta, respectively. Font size correlates to relative FC.
Figure 3. Glycogen metabolism associated genes in So and Ta myofibers. (a) Glucose is converted into glycogen during glycogenesis. Glycogenolysis converts glycogen back into glucose molecules. Represented are magnesium ions (+Mg2+), adenosine diphosphate (ADP), uridine-5′-triphosphate (UTP), pyrophosphate (PPi), and inorganic phosphate (Pi). (b) Two subunits of the protein phosphorylase 1 were increased in So. Represented are an encircled ‘P’ indicating phosphorylation and uridine-diphosphate (UDP). (c) Three isoform subunits of phosphorylase kinase were increased in Ta. Represented are the addition (+Ca2+) and the removal of calcium ions (−Ca2+). Substrates and products of the enzymatic steps in glycogen metabolism include the interconversion of inactive to active states of regulatory enzymes. Enzymes involved at specific stages are encircled and adjacent to the corresponding reaction arrow. Gene symbols located underneath enzyme names in red or blue font were increased in So or Ta, respectively. Font size correlates to relative FC.
Genes 12 01718 g003
Genes 12 01718 g004 550
Figure 4. Glucose metabolism associated genes in So and Ta myofibers. Enzymes involved at specific stages are encircled and adjacent to the corresponding reaction arrow. (a) Aldoa, Aldoc, Tpi1, Pgk1, Eno3, Pfkm, Pkm, and Fbp2 were increased in Ta. Isoforms of enolase (Eno2) and PEP carboxykinase (Pck1) were increased in So. Represented are guanosine triphosphate (GTP), guanosine diphosphate (GDP), and inorganic phosphate (P). (b) Subunits of the PFK2/FBPase2 complex, Pfkfb3, Pfkfb1, and Pfkfb4, were increased in Ta. The red or blue font of gene symbols underneath enzyme names represent increased expression in So or Ta, respectively. Font size correlates to relative FC.
Figure 4. Glucose metabolism associated genes in So and Ta myofibers. Enzymes involved at specific stages are encircled and adjacent to the corresponding reaction arrow. (a) Aldoa, Aldoc, Tpi1, Pgk1, Eno3, Pfkm, Pkm, and Fbp2 were increased in Ta. Isoforms of enolase (Eno2) and PEP carboxykinase (Pck1) were increased in So. Represented are guanosine triphosphate (GTP), guanosine diphosphate (GDP), and inorganic phosphate (P). (b) Subunits of the PFK2/FBPase2 complex, Pfkfb3, Pfkfb1, and Pfkfb4, were increased in Ta. The red or blue font of gene symbols underneath enzyme names represent increased expression in So or Ta, respectively. Font size correlates to relative FC.
Genes 12 01718 g004
Genes 12 01718 g005 550
Figure 5. Contraction-related genes in So and Ta. (a) Graphic representation of the sarcomere. The M-line, A-zone, I-zone, and Z-line are depicted. The red band represents the thick myosin filament. The blue spheres represent the thin filament. Within the I-zone, resides an elastic protein called titin. (b) Graphic representation of thin filament structure. Tpm3, Tnnc1, Tnni1, Tnnt1, Tnni3, Tnnt2, and Tnnt3 were increased in So. Actc1, Tmod1, Tpm1, and Tnnt3 were increased in Ta. (c) Thick filament close-up. Myh7, Myl6b, Myl3, Myl4, Myl2, and Myl10 were increased in So. Myh4, Mylpf, Mybpc2, and Mybph were increased in Ta. (d) Z-line close-up. The Z-line region of the sarcomere is where actin and titin are connected to an antiparallel complex, including titin cap (Tcap) and actinin (Actn3). (e) Myofiber structure. Jph3, Lman1l, Phactr3, Tnni3k, Dsp, Grp, Tom1, Cntnap4, Slc30a3, Wdr72, Lypd1, Asphd1, Nrxn3, and Nrsn2 were increased in So. Pak1, Wwp1, Fam19a4, Kcnc4, and Fgfbp1 were increased in Ta. Red or blue gene symbols underneath structure names represent increased expression in So or Ta, respectively. Font size correlates to relative FC.
Figure 5. Contraction-related genes in So and Ta. (a) Graphic representation of the sarcomere. The M-line, A-zone, I-zone, and Z-line are depicted. The red band represents the thick myosin filament. The blue spheres represent the thin filament. Within the I-zone, resides an elastic protein called titin. (b) Graphic representation of thin filament structure. Tpm3, Tnnc1, Tnni1, Tnnt1, Tnni3, Tnnt2, and Tnnt3 were increased in So. Actc1, Tmod1, Tpm1, and Tnnt3 were increased in Ta. (c) Thick filament close-up. Myh7, Myl6b, Myl3, Myl4, Myl2, and Myl10 were increased in So. Myh4, Mylpf, Mybpc2, and Mybph were increased in Ta. (d) Z-line close-up. The Z-line region of the sarcomere is where actin and titin are connected to an antiparallel complex, including titin cap (Tcap) and actinin (Actn3). (e) Myofiber structure. Jph3, Lman1l, Phactr3, Tnni3k, Dsp, Grp, Tom1, Cntnap4, Slc30a3, Wdr72, Lypd1, Asphd1, Nrxn3, and Nrsn2 were increased in So. Pak1, Wwp1, Fam19a4, Kcnc4, and Fgfbp1 were increased in Ta. Red or blue gene symbols underneath structure names represent increased expression in So or Ta, respectively. Font size correlates to relative FC.
Genes 12 01718 g005
Table
Table 1. DE genes associated with lipid, amino acid, and one-carbon metabolic pathways.
Table 1. DE genes associated with lipid, amino acid, and one-carbon metabolic pathways.
Gene SymbolGene NameFold Change
(So/Ta)		Localization and Biochemical Properties
Metabolism
Dio1deiodinase, iodothyronine, type I69.1de-iodinate thyroid hormone interaction; FA oxidation; oxidative phosphorylation uncoupling; causes mitochondrial heat production [15]
Pon1paraoxonase 168.83protection against oxidation for HDLs and LDLs [10]
Tecrltrans-2,3-enoyl-CoA reductase-like55FA elongation in polyunsaturated FA biosynthesis [17]
Tspo2translocator protein 251.12free cholesterol trafficking in erythroid cells [11]
Acsm5acyl-CoA synthetase medium-chain family member 530.67FA β-oxidation [13]
Apoa2apolipoprotein A-II27.64modulator of reverse cholesterol transport [12]
Akr1c18aldo-keto reductase family 1, member C1824.68catalyzes progesterone into 20-α-dihydroprogesterone (20-α-OHP) [21]
Adtrpandrogen dependent TFPI regulating protein21.66hydrolyzes FA esters of hydroxy FAs [22]
Cyp4f39cytochrome P450,
family 4, subfamily f, polypeptide 39		−269.231ω and 2ω diesters types and cholesteryl (O-Acyl)-w-hydroxy FAs (OAHFAs) production; FA ω-hydroxylase; acyl-ceramide production [19]
Amino Acid Metabolism
Aox3aldehyde oxidase 330.48potentially linked to amino acid and retinol metabolisms; oxidizes aliphatic or aromatic aldehydes into carboxylic acid [23]
One-Carbon Metabolism
Gnmtglycine N-methyltransferase26.87Methylation of DNA, RNA, proteins, and lipids via the methionine cycle
Miscellaneous
Aldh1a7aldehyde dehydrogenase family 1, subfamily A7−20.59Oxidoreductase, NAD/NADP acceptor [24]
Positive or negative fold change represents increased gene expression in So or Ta, respectively.
Table
Table 2. Differentially expressed genes associated with glucose and glycogen metabolic pathways.
Table 2. Differentially expressed genes associated with glucose and glycogen metabolic pathways.
Gene SymbolGene NameFold Change
(So/Ta)		Localization and Biochemical Properties
Glycogenolysis and Glycogenesis
Ppp1r1cprotein phosphatase 1, regulatory inhibitor subunit 1C75.78indirect regulation of glycogenesis and glycogenolysis
Ppp1r3gprotein phosphatase 1, regulatory subunit 3G36.99PP1 control; indirect regulation of glycogenesis and glycogenolysis [25]
Gys2glycogen synthase 2−90.2converts glucose 1-phosphate into glycogen [26]
Phkg1phosphorylase kinase γ 1−8.1activates glycogen phosphorylase [29]
Phka1phosphorylase kinase α 1−8activates glycogen phosphorylase [29]
Phkbphosphorylase kinase β−7.5activates glycogen phosphorylase [29]
Pgm2phosphoglucomutase 2−6.7converts of glucose 6-phosphate and glucose 1-phosphate
Pygmmuscle glycogen phosphorylase−5.8converts glycogen into glucose 1-phosphate
Aglamylo-1,6-glucosidase, 4-α-glucanotransferase−4.4glucosidase activity that converts α-1,6-linked branches into α-1,4-linked branches
Epm2aepilepsy, progressive myoclonic epilepsy, type 2 gene α−3.5facilitates the de-phosphorylation of glycogen to promote branching [27]
Pgm2l1phosphoglucomutase 2-like 1−3conversion of glucose 6-phosphate and glucose 1-phosphate
Glycolysis and Gluconeogenesis
Pck1phosphoenolpyruvate carboxykinase 1, cytosolic9.3converts oxaloacetate to phosphoenolpyruvate
Eno2enolase 2, γ neuronal3converts 3-PG into pyruvate
Pfkfb36-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 3−9.3converts fructose-6-phosphate to F2,6-BP; indirectly regulates glycolytic flux [30]
Eno3enolase 3, β muscle−5.4converts 3-PG into pyruvate in glycolysis
Pgk1phosphoglycerate kinase 1−5.1converts 1,2-bisphosphoglycerate and 3-phosphoglycerate
Tpi1triose-phosphate isomerase 1−4.5converts dihydroxyacetone phosphate to D-glyceraldehyde 3-phosphate
Fbp2fructose bisphosphatase 2−4.2converts F1,6-BP to fructose 6-phosphate
Pfkfb46-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 4−4.1converts F2,6-BP to fructose-6-phosphate; indirectly regulates glycolytic flux
Pkmpyruvate kinase, muscle−3.9converts PEP to pyruvate
Aldocaldolase C, fructose-bisphosphate−3.9converts F1,6-BP to dihydroxyacetone phosphate
Aldoaaldolase A, fructose-bisphosphate−3.8converts F1,6-BP to dihydroxyacetone phosphate; present in developing embryo and adult muscle
Pfkfb16-phosphofructo-2-kinase/fructose-2,6-bisphosphatase 1−3.1regulates glycolytic flux
Pfkmphosphofructokinase, muscle−3regulates glycolytic flux; converts fructose 6-phosphate to F1,6-BP [31]
Positive or negative FC represents increased gene expression in So or Ta, respectively.
Table
Table 3. Differentially expressed genes associated with contractile related structures.
Table 3. Differentially expressed genes associated with contractile related structures.
Gene SymbolGene NameFold Change
(So/Ta)		Localization and Biochemical Properties
Neuromuscular Junction
Grpgastrin releasing peptide149.61modulates autonomic system; regulates male sexual function; conveys itch sensation
Tom1target of myb1 trafficking protein84.21recruits clathrin to endosomal structures [47]
Jph3junctophilin 366.32formation of junctional membrane complexes; motor coordination [50]
Cntnap4contactin associated protein-like 458.74cell adhesion molecule and receptor [48]
Slc30a3
(Znt3)		solute carrier family 30 (zinc transporter), member 334.8transport of zinc in synaptic vesicles
Wdr72WD repeat domain 7233.63endocytosis, protein reabsorption, and calcium excretion [49]
Lypd1Ly6/Plaur domain containing 126.17regulates neuronal nicotinic receptors
Asphd1aspartate β-hydroxylase domain-containing protein 125.7possibly neurotransmission and synaptic vesicle location/function
Nrxn3neurexin III24.74synapse organization; regulating neurotransmitter release [43]
Nrsn2neurensin 224.72transporting small vesicles to perinuclear region to exit towards organelles
Fam19a4 (Tafa4)TAFA chemokine like family member 4−199.82modulates neuronal excitability and controls somatic sensation threshold [44]
Kcnc4 (Kv3.4)potassium voltage gated channel, Shaw-related subfamily, member 4−33.68broadens action potential to prevent somatic depolarization [45]
Fgfbp1fibroblast growth factor binding protein 1−28.02slows neuromuscular junctions (NMJs) degeneration [46]
Sarcolemma, Sarcoplasm, and Sarcoplasmic Reticulum
Lman1llectin, mannose-binding 1 like61.63binds to glycoproteins transported from endoplasmic reticulum (ER) to ER-Golgi intermediate [51]
Tnni3kTNNI3 interacting kinase27.58promotes oxidative stress and myocyte death; part of costamere attached to sarcolemma [52]
Phactr3phosphatase and actin regulator 326.09binds to cytoplasmic actin and regulates PP1
Dspdesmoplakin20.92assembles functional desmosomes; maintains cytoskeletal architecture [53]
Pak1p21 (RAC1) activated kinase 1−7.3regulates muscle-specific kinase to maintain NMJs, actin remodeling, and glucose uptakes in skeletal muscle
Wwp1WW domain containing E3 ubiquitin protein ligase 1−3.5cell proliferation and apoptosis
Thin Filament
Actin
Actc1actin, α, cardiac muscle 1−20.7assembles into filaments that are involved in muscle contraction, cell motility, cell signaling, and vesicle movement; associated with fetal skeletal muscle
Tropomyosin and Tropomodulin
Tpm3tropomyosin 3, γ29.17controls the actin filament with tropomodulin [32]
Tmod1tropomodulin 1−3.6controls the Ca2+-regulated thin filament end with tropomyosin
Tpm1tropomyosin 1, α−3.4controls the Ca2+-regulated thin filament end with tropomodulin [33]
Troponin Complex
Tnnc1troponin C, cardiac/slow skeletal47.7binds to calcium and exposes the myosin head binding sites; slow skeletal muscle
Tnni1troponin I, skeletal, slow 146.2binds to actin and inhibits ATPase activity; slow skeletal muscle
Tnnt1troponin T1, skeletal, slow36.6anchors to tropomodulin; slow skeletal muscle [34]
Tnni3troponin I, cardiac 314.3binds to actin and inhibits ATPase activity; cardiac muscle
Tnnt2troponin T2, cardiac3.3anchors to tropomodulin; cardiac, embryonic, and neonatal skeletal muscles [35]
Tnnt3troponin T3, skeletal, fast−3anchors to the tropomodulin; fast skeletal muscle [36]
Thick Filament
Myosin Heavy Chains
Myh7myosin, heavy polypeptide 7, cardiac muscle, β34.71isoform of myosin present in slow (type I) skeletal muscle fibers [1]
Myh4myosin, heavy polypeptide 4, skeletal muscle−57.73isoform of myosin present in adult type IIB skeletal muscle fibers [37]
Essential Myosin Light Chains
Myl6bmyosin, light polypeptide 6B54.1controls cell adhesion, cell migration, tissue architecture, cargo transport, and endocytosis; promotes p53 protein ubiquitination and degradation
Myl3myosin, light polypeptide 327involved in force development and fine-scale coordinated muscle contraction
Myl4myosin, light polypeptide 43.6binds to Ca2+; embryonic skeletal muscle and atrial myocardium [2]
Regulatory Myosin Light Chains
Myl2myosin, light polypeptide 2, regulatory, cardiac, slow41.7stiffens myosin neck and aids in myosin head movement
Myl10myosin, light chain 10, regulatory4stiffens myosin neck and aids in myosin head movement [38]
Mylpfmyosin light chain, phosphorylatable, fast skeletal muscle−4regulates myofilament activation by phosphorylation [39]
Myosin Binding Proteins
Mybpc2myosin binding protein C, fast type−40.7regulates of myofilament contraction [40]
Mybphmyosin binding protein H−42.77sarcomere contraction; maturation process of auto-phagosomal membranes; inhibition of non-muscle RLC MYL2A and MYH2A
Z-line
Tcaptitin-cap3.5binds to titin in an anti-parallel complex and stabilizes Z-line [41]
Actn3actinin α 3−33.3anchors actin filaments and supports sarcomere [42]
Positive or negative FC represents increased gene expression in So or Ta, respectively.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Alto, S.I.; Chang, C.-N.; Brown, K.; Kioussi, C.; Filtz, T.M. Gene Expression Profiling of Skeletal Muscles. Genes 2021, 12, 1718. https://doi.org/10.3390/genes12111718

AMA Style

Alto SI, Chang C-N, Brown K, Kioussi C, Filtz TM. Gene Expression Profiling of Skeletal Muscles. Genes. 2021; 12(11):1718. https://doi.org/10.3390/genes12111718

Chicago/Turabian Style

Alto, Sarah I., Chih-Ning Chang, Kevin Brown, Chrissa Kioussi, and Theresa M. Filtz. 2021. "Gene Expression Profiling of Skeletal Muscles" Genes 12, no. 11: 1718. https://doi.org/10.3390/genes12111718

APA Style

Alto, S. I., Chang, C. -N., Brown, K., Kioussi, C., & Filtz, T. M. (2021). Gene Expression Profiling of Skeletal Muscles. Genes, 12(11), 1718. https://doi.org/10.3390/genes12111718

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