Next Article in Journal
Mining Proteome Research Reports: A Bird’s Eye View
Previous Article in Journal / Special Issue
Phosphorylation-Dependent Interactome of Ryanodine Receptor Type 2 in the Heart
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Advanced Fiber Type-Specific Protein Profiles Derived from Adult Murine Skeletal Muscle

1
Medizinisches Proteom-Center, Medical Faculty, Ruhr-University Bochum, 44801 Bochum, Germany
2
Medical Proteome Analysis, Center for Protein Diagnostics (PRODI), Ruhr-University Bochum, 44801 Bochum, Germany
3
Institute of Neuropathology, University Hospital Erlangen, Friedrich-Alexander University Erlangen-Nürnberg, 91054 Erlangen, Germany
4
German Aerospace Center, Institute of Aerospace Medicine, 51147 Cologne, Germany
5
Center for Physiology and Pathophysiology, Institute of Vegetative Physiology, Medical Faculty, University of Cologne, 50931 Cologne, Germany
*
Authors to whom correspondence should be addressed.
Submission received: 14 April 2021 / Revised: 1 June 2021 / Accepted: 2 June 2021 / Published: 8 June 2021
(This article belongs to the Special Issue Striated Muscle Proteomics II)

Abstract

:
Skeletal muscle is a heterogeneous tissue consisting of blood vessels, connective tissue, and muscle fibers. The last are highly adaptive and can change their molecular composition depending on external and internal factors, such as exercise, age, and disease. Thus, examination of the skeletal muscles at the fiber type level is essential to detect potential alterations. Therefore, we established a protocol in which myosin heavy chain isoform immunolabeled muscle fibers were laser microdissected and separately investigated by mass spectrometry to develop advanced proteomic profiles of all murine skeletal muscle fiber types. All data are available via ProteomeXchange with the identifier PXD025359. Our in-depth mass spectrometric analysis revealed unique fiber type protein profiles, confirming fiber type-specific metabolic properties and revealing a more versatile function of type IIx fibers. Furthermore, we found that multiple myopathy-associated proteins were enriched in type I and IIa fibers. To further optimize the assignment of fiber types based on the protein profile, we developed a hypothesis-free machine-learning approach, identified a discriminative peptide panel, and confirmed our panel using a public data set.

1. Introduction

Skeletal muscle is a heterogeneous tissue composed of muscle fibers, connective tissue, nerve fibers, and blood vessels. The individual protein expression pattern of a single muscle fiber is of fundamental interest to improve our current understanding of specific molecular mechanisms, including the selective loss of different fiber subtypes in response to a variety of factors. Fiber types are classically subdivided into slow type I and fast type IIa, IIb, and IIx fibers [1,2]. Type IIb fibers are exclusively expressed in rodents [3]. All fiber types have defined tasks, e.g., a holding function, continuous movements, or abrupt movements, and therefore, they have unique proteomic compositions. Addressing the latter, it is thus crucial to study individual muscle fiber types in order to determine their specific proteomic patterns in response to, e.g., training, aging, denervation, or intrinsic myopathic processes.
In general, three classification approaches are widely used for the distinction of different skeletal muscle fibers: (1) traditional classification based on an ATPase activity staining technique; (2) differentiation via energetic properties. Slow fibers have a high level of oxidative phosphorylation activity, while fast fibers mainly rely on glycolytic events to sustain their energy supply [4,5]. The third and most widely used approach is differentiation via the specific myosin heavy chain (MYH) isoform expression pattern [6,7,8]. Type I fibers (slow) mainly contain myosin-7 (MYH7), type IIa fibers (fast) contain myosin-2 (MYH2), and type IIx fibers (fast) contain myosin-1 (MYH1), whereas myosin-4 (MYH4) is mainly found in type IIb fibers (very fast), as summarized in Table 1 [3,5]. Notably, skeletal muscles are highly dynamic with a remarkable level of plasticity that allows individual fibers to alter their MYH isoform pattern. In addition to pure slow and fast fiber types, there are intermediate forms that contain multiple MYH isoforms [2,3,9,10]. This plasticity allows fast adaptation in response to various conditions, including exercise and training, physical inactivity, aging, and myogenic, neurogenic, inflammatory, and metabolic disease processes [5,11,12,13,14,15,16,17]. For example, aging studies have been associated with selective wasting and a loss of type II fibers [13,18,19,20,21], whereas type I fiber predominance or the selective affection of fast fiber types has been described in a wide variety of primary muscle diseases, e.g., in the group containing myofibrillar myopathies (MFM) and Duchenne muscular dystrophy (DMD) [22,23,24,25,26]. Until now, the associated molecular processes, inducing fiber type shift, have been poorly understood, but knowledge of these processes could lead to novel therapeutic options. Hence, fiber type analysis at the omics level is emerging as a major field of study to determine the selective involvement of fiber types in processes related to age, exercise, and disease.
Thus, further progress in our current understanding of the protein diversity of individual muscle fiber types was provided by transcriptomic and proteomic analyses. For instance, the transcriptome profiles of total muscle tissue in comparison with derived type I and type IIa fibers provided insight into the levels of mitochondrial, cytosolic, and muscle structure coding genes as well as the distinct expression of different protein isoforms, e.g., caveolin, synaptotagmin, and tropomodulin, within the specific fiber subtypes [27]. Later proteomic studies of muscle fiber subtypes succeeded in revealing further specific features such as fiber type-specific protein changes during muscle atrophy [28], different mitochondrial adaptations of the fiber types as well as age-dependent differences in energetic profiles between young and old muscle fibers derived from healthy human subjects [12,13]. In the latter studies, single fiber proteome profiles were post-analytically classified into their respective types (type I, type IIa, type IIx, type IIb, and hybrid fibers) depending on their myosin isoform abundance.
In the context of focally restricted skeletal muscle alterations, the laser microdissection (LMD) technique provides the possibility to isolate pathogenically relevant tissue or even fiber areas for further proteomic analyses. Using this method, we previously characterized protein aggregates derived from genetically different myofibrillar myopathies [11,29]. Furthermore, this approach led to the identification of proteins that had not previously been described as components of disease-specific protein aggregates. Moreover, we used this approach for the isolation of type I fibers in R349P desmin knock-in mice and delineated a distinct mitochondrial pathology with this desminopathy model [30]. The combination of the immunohistochemical identification of fiber subtypes using MYH isoform-specific antibodies prior to LMD and mass spectrometry (MS) processing allowed us to obtain fiber type-specific proteomes covering up to 1620 proteins as well as to quantify even less abundant proteins. With our method, we were able to analyze the actual proteomic composition of all subfiber types without retrospective annotation.
Our study yielded fiber type-specific proteomic profiles, supporting well-known fiber type characteristics that are commonly reported in the literature. We also provided novel information on fiber type composition and functions. With our protein profiling, we demonstrated that type IIx fibers are not primarily glycolytic, but they do inherit a large number of mitochondrial proteins that are essential for oxidative energetic processes. Our analysis provided us with further evidence that type I and type IIa fibers contain enhanced amounts of proteins known to cause neuromuscular disorders when altered. Our results also suggest that proteins that are essential for accurate cellular maintenance, such as heat shock proteins and the UPS system, have a higher abundance in type I fibers, indicating that there are potential superior rescue mechanisms in this fiber type, as previously reported in [13,31]. With these results, we were able to form a cautious hypothesis on the potential cause of the frequently observed type I fiber predominance. Finally, we developed a hypothesis-free machine learning-based approach, the Random Forest-based Muscle Fiber Peptide Selection Pipeline (FiPSPi), to optimize fiber type assignment by first identifying discriminative peptide panels and then confirming our discriminative peptide panel using a public data set. We are now able to offer our confirmatory R-Script for other researchers, and we plan to standardize our discriminative R-Script to enable the identification of other peptide panels.

2. Materials and Methods

2.1. Animals

The WT mice used in this study were housed under specific and opportunistic pathogen-free (SOPF) conditions with free access to water and food in a standard environment. Health monitoring was carried out as recommended by the Federation of European Laboratory Animal Science Associations (FELASA). Mice were handled in accordance with the German Animal Welfare Act and the German Regulation for the protection of animals used for experimental purposes or other scientific purposes. All investigations were approved by the responsible governmental animal care and use the office (North Rhine-Westphalia State Agency for Nature, Environment and Consumer Protection (LANUV), Recklinghausen, Germany; reference number 84-02.05.40.14.057). Detailed information on sex, genotype, and age can be found in Table 2.

2.2. Preparation of Skeletal Muscle Cryosections, Immunohistochemistry, and Immunofluorescence Staining

Mice were sacrificed, and skeletal muscle specimens were collected and immediately frozen in isopentane precooled in liquid nitrogen. Soleus and tibialis anterior muscles were cryosected into 10 µm thick slices (Cyrostat Microm HM550, Thermo Fisher Scientific, Schwerte, Germany) and placed on a PET membrane frame slide for laser microdissection (LMD). Staining was carried out as described in [32], except that primary antibody incubation was carried out overnight. Antibodies directed against the specific MYH isoforms were purchased from the Developmental Studies Hybridoma Bank (DSHB Iowa City, IA, USA) and used as described in Table 3. PET membrane slides were completely dried prior to LMD.

2.3. Laser Microdissection

Laser microdissection (LMD; 6500, Leica Microsystems, Wetzlar, Germany) was carried out as described in [32] except that a total fiber area of 1,000,000 µm2 (corresponding to approximately 500 excised fibers) was excised for further analysis with mass spectrometry. Fiber types were excised at 20× magnification. LMD parameters, such as laser intensity, cutting speed, and fluorescence intensity, were adjusted for each PET membrane slide separately to ensure an optimal sample outcome. After the LMD process, samples were covered with 40 µL of formic acid and incubated for 20 min at RT. Sonication of samples for 5 min led to tissue and PET membrane disruption. Samples were stored at −80 °C until further processing.

2.4. Lysis of Whole Soleus and Tibialis Anterior Muscles for Spectral Library Generation

Complete murine soleus and tibialis anterior muscles were pulverized in liquid nitrogen, homogenized on ice, and then resuspended in urea buffer (7 M urea, 2 M thiourea, 20 mM Tris base, pH 8.5) as described in [33]. Resuspended samples were sonicated six times for ten seconds with ten seconds rest on ice to support the lysis. The protein concentration was determined by Bradford assay.

2.5. In-Solution Tryptic Digestion

In-solution tryptic digestion was performed as described in [30]. In brief, formic acid was completely removed by vacuum vaporization (Vacuum concentrator RCV2-25 CD Plus, Martin Christ Gefriertrocknungsanlagen, Osterode, Germany), and the total digestion volume was set to 50 µL using 50 mM ammonium bicarbonate and 5 mM 1,4-dithiothreitol (DTT) (end concentration). Reduction with DTT was performed for 20 min at 56 °C followed by alkylation for 15 min in the dark at room temperature with 15 mM IAA (end concentration). Tryptic digestion was carried out overnight using 0.1 µg of trypsin per sample, and the reaction was stopped after 16 h by adding 0.5% trifluoroacetic acid (TFA) (end concentration). The sample volume was completely vaporized in a vacuum concentrator, and peptides were solved in 50 µL of 0.1% TFA.

2.6. High pH Reversed-Phase Fractionation

Fractionation of the highly complex muscle sample enabled higher proteome coverage and, through that, high protein identification rates in the MS analysis. High pH reversed-phase fractionation of digested muscle samples was performed using the Pierce™ High pH Reversed-Phase Peptide Fractionation Kit (Thermo Fisher Scientific, Bleiswijk, The Netherlands), as recommended in the manual. Briefly, 50 µg of digested muscle lysate was loaded onto the equilibrated spin column, and an increasing organic solvent gradient led to the fractionation of peptides. In total, 8 fractions were collected and used for spectral library generation.

2.7. MS Analysis

The mass spectrometry proteomics data have been deposited into the ProteomeXchange Consortium via the PRIDE [34] partner repository with the data set identifier PXD025359.

2.7.1. Data Dependent Acquisition (DDA) for Spectral Library Generation

Nano HPLC analysis was performed on an UltiMate 3000 RSLC nano LC system (Thermo Fisher Scientific, Bremen, Germany), as described in [11]. The HPLC system was online-coupled to the nano ESI source of a QExactive HF mass spectrometer (Thermo Fisher Scientific, Bremen, Germany). In the ESI-MS/MS analysis, full MS spectra were acquired in the range from 350 to 1400 m/z with a resolution of 60,000 at 200 m/z for the detection of precursor ions (AGC target 3 × 106, 80 ms maximum injection time). The m/z values initiating MS/MS were set on a dynamic exclusion list for 30 s, and the ten most intensive ions (charge state +2, +3, +4) were selected for fragmentation. MS/MS fragments were generated by high-energy collision-induced dissociation (HCD). The normalized collision energy (NCE) was either set to a fixed value of 27, or a stepped NCE was applied (25.5, 27, 30). The fragments were analyzed in an Orbitrap analyzer with a 30,000 resolution at 200 m/z (AGC 1 × 106, maximum injection time 120 ms).

2.7.2. Spectral Library Generation

Data-independent acquisition (DIA) regularly uses a spectral library to ensure the correct assignment of fragment ion spectra to their respective precursor ions. To generate an optimal DIA method, a sample-specific library was created using laser microdissected muscle tissue samples as well as whole soleus and tibialis anterior muscles (Spectral library Supplementary Tables S1 and S2). All samples were lysed and partially fractionated by high pH fractionation prior to DDA analysis. Prior to analysis, all samples were labeled by using internal retention time (iRT) peptides (Biognosys, Schlieren, Switzerland), as recommended by the vendor’s protocol. Different gradients, as well as NCE settings, were included in the library to extend the protein identification rates, resulting in a total of 56 DDA runs (for annotated raw files, see Supplementary Table S3). Spectral library generation was carried out using SpectronautTM Pulsar and the Pulsar search engine (v. 12.0.20491.7.17149, Biognosys, Schlieren, Switzerland). Data were searched against the Uniprot KB mus musculus reference proteome set, including iRT peptides (version 2018_06, 53,560 entries) and a contaminant database, resulting in a library size of 28,004 peptides and 5206 proteins. Biognosys factory settings were applied, and trypsin was chosen as the digestion enzyme. Specific modifications were set according to the sample treatment: carbamidomethylation (C) (fixed modification) and oxidation of (M), deamidation (NQ), and carbamidomethylation (N-term) (variable modifications) were included.

2.7.3. Data-Independent Acquisition (DIA) for Differential Proteome Analysis

Prior to DIA, labeling of samples by internal retention time (iRT) calibration was performed for the evaluation of data with SpectronautTM Pulsar. iRT peptides were prepared as recommended by Biognoysis, and a volume of 1 µL per sample was spiked in the prior analysis. All steps of DIA were again carried out on the QExactive HF using an 8 µL sample volume. In the ESI-MS/MS analysis, full MS spectra were scanned in a range from 350 to 1400 m/z with a resolution of 120,000 at 200 m/z for the detection of precursor ions (AGC target 3 × 106, 20 ms maximum injection time). MS/MS fragments were generated by high-energy collision-induced dissociation (HCD) in which ion dissociation was performed at an NCE of 27%, a fixed first mass of 130.0 m/z, and 24 isolation windows of 45 m/z. The fragments were analyzed in an Orbitrap analyzer with 30,000 resolution at 200 m/z (AGC 1 × 106, maximum injection time 120 ms).

2.8. Absolute Protein Concentration Estimation Using aLFQ

After exporting the non-normalized protein and peptide intensities out of Spectronaut Pulsar (Supplementary Tables S4 and S5), locally weighted scatter plot smoothing normalization (LOESS) was applied using the limma R package, version 3.36.5 [35,36]. Calculation of protein quantities was performed using the aLFQ package (version 1.3.5 [37]) within R (version 3.6.1). First, for each protein, an iBAQ value [38] was calculated from the LOESS-normalized peptide intensities (Supplementary Table S6). Absolute protein quantities were obtained using the known total protein concentration of 200 ng per sample. In aLFQ, this is achieved by dividing all iBAQ values from all samples by the total sum of iBAQ values and multiplying this with the specified total protein concentration. We modified the aLFQ workflow and conducted the protein quantification step separately for each sample to ensure there was a total normalized protein intensity of 200 ng in each sample. The modified R-script (aLFQ) is provided in the Supplementary Material. Percentage values were calculated determined by calculating the sum of all protein aLFQ values for each sample separately. The resulting sum value was equal to 100%. Subsequently, single protein percentage values were calculated for each protein and each sample.

2.9. Hierarchical Clustering and Subsequent go Term Enrichment and Pathway Enrichment

Non-normalized data obtained from Spectronaut were normalized as described above and used for unsupervised hierarchical clustering with Perseus (v. 1.6.1.3) [39]. For this, the normalized values were log transformed. Prior to clustering, the data set was filtered by setting a minimum percentage of valid values of 50%. Existing missing values were replaced using a normal distribution by imputing with a width of 0.3 and a downshift of 1.8. Groups were averaged using the medians, and data were Z-scored prior to clustering. For unsupervised hierarchical clustering, the Euclidian distance was chosen with an average linkage and no constraint. K-means was enabled, and a value of 300 clusters was chosen with a maximal number of 10 iterations. The resulting clusters were exported and used for subsequent enrichment analysis.
Pathway and Gene annotation enrichment analyses were carried out using DAVID Bioinformatics Resources 6.8 [40,41,42]. For this, the functional annotation tool was used. Uniprot accession of proteins was uploaded. Here, only the first accession for protein group results was used to enable accurate analysis. The Uniprot ID identifier was set, and the list type used was the gene list. As a background, the whole mus musculus genome was chosen. The KEGG pathway was used for the determination of enriched pathways in the data set. The term “cellular compartment” was chosen for GO term enrichment studies. Fold enrichment scores, as well as FDR, were used as additional options. Threshold factory settings were maintained. The resulting tables were exported into excel and ranked using p-value corrected values or fold enrichment scores.

2.10. Random Forest-Based Muscle Fiber Peptide Selection Pipeline (FiPSPi)

In order to identify a discriminative peptide panel for optimized fiber type assignment by a hypothesis-free machine learning-based method, the Random Forest-based Muscle Fiber Peptide Selection Pipeline (FiPSPi) was implemented as part of this work. FiPSPi was implemented in the Medizinisches Proteom-Center and BioInfra.Prot [43] in R programming language [44] using the R packages randomForest [45] and rpart [46]. Currently, FiPSPi consists of two R scripts, one to identify panels of discriminative peptides and one to confirm them with independent data. The original FiPSPi code is available as supplementary material for this paper, and further developed versions can be downloaded from GitHub (https://github.com/mpc-bioinformatics/FiPSPi, accessed on 11 November 2019).

2.11. Validation by Normalized Quantitative Western Blotting Based on Standardized Fluorescent Labeling

Due to the very limited amount of LMD tissue, whole murine soleus and tibialis anterior muscles were used for the Western Blot analysis to validate the findings of the global protein analysis [47,48]. Soleus and tibialis anterior muscles were lysed as described in [33], and 50 µg of protein lysate was separated by 1D-PAGE (polyacrylamide gel electrophoresis) (50 V for 15 min and 180 V for 1 h) using MES running buffer (50 mM MES, 50 mM Tris, 2 mM EDTA). The Quantitative Smart Protein Layers (SPL) Western Blotting system (NH DyeAGNOSTICS GmbH, Halle, Germany) was used for the detection of quantitative differences between selected candidate proteins, and it was used as recommended in the vendors’ protocol. In short, the total protein was pre-labeled with a red fluorescent fluorophore (700 nm detection wavelength), and a 12.5 kDa green fluorescence-labeled standard protein was spiked in (800 nm detection wavelength), enabling error correction of differing sample loadings and data normalization between experiments. The resulting gel was scanned at wavelengths of 700 and 800 nm prior to blotting to determine the total protein and spike in protein quantities and identify inconsistent loss during the blotting procedure. Semi-dry Western blotting was performed at 128 mA for one hour using two different buffer systems (Anode buffer: 300 mM Tris, 100 mM Tricine pH: 8.7, Cathode buffer: 300 mM6-Aminocaproic acid, 30 mM Tris pH: 8.7). Antibodies for desmin (DAKO, Agilent, Frankfurt, Germany 1:100) and myotilin (Santa Cruz, Heidelberg, Germany, 1:200) were used for validation. Incubation was carried out overnight at 8 °C. Blots were washed three times with 1× TBS and incubated with secondary antibody (IRDye 800CW, LI-COR Biosciences, Bad Homburg, Germany, Goat anti-mouse, 1: 10,000) for one hour at room temperature. After washing, Western blots were scanned at wavelengths of 700 and 800 nm using the Odyssey® imaging system (LI-COR Biosciences, Bad Homburg, Germany). Quantification of results was carried out using the vendor’s software (SPL LabImage Software), as described by the vendor’s protocol, resulting in quantitative protein volumes (SPL normalized volume).

3. Results

3.1. Study Overview and Proteomics Workflow

Recent fiber type classification methods rely on immunolabeling of myosin heavy chain (MYH) isoforms in skeletal muscle tissue specimens. We used this approach for the precise excision of skeletal muscle fiber types out of complex muscle tissue using laser microdissection (LMD) followed by data-independent acquisition (DIA) mass spectrometry (MS) to uncover differences in four murine skeletal muscle fiber type proteomes. For this, 15 mice were euthanized, and their soleus and tibialis anterior muscles were excised, cryosected, and stained for the collection of type I, type IIa, type IIb, and type IIx fibers. MYH isoform-specific antibodies enabled the distinction between fiber types and their selective isolation (Supplementary Figure S1). Using an MS-based approach, whole muscle tissue and LMD muscle tissue were analyzed, resulting in the identification of 5206 muscle-specific proteins. This in-depth proteome served as a reference for the subsequent fiber type characterization. The analysis of excised fiber types (total number of fiber type samples n = 56) resulted in the identification of over 1600 (1037 quantified) different fiber-specific proteins using a minimal amount of sample (approximately 150 ng of peptide extract). The fiber-specific proteome analysis revealed 889 proteins in type I, 875 proteins in type IIa, 738 proteins in type IIb, and 800 proteins in type IIx fibers.

3.2. Fiber Types Are Distinguishable at the Proteomic Level

In the first step, we assessed whether the different fiber subtypes could be generally distinguished based on their protein profiles. For this, sample correlation was carried out by principal component analysis (PCA) for each fiber type sample (Figure 1a). A significant separation of type IIb and type IIx fibers was observed, while type I and type IIa fibers formed a combined, but significantly separated, cluster from type IIx and type IIb fibers with moderate segregation into type I and type IIa fibers. From these results, we concluded that type I and type IIa fibers have similar protein profiles. Type IIb and type IIx fibers, on the other hand, harbor unique proteomic profiles that are clearly distinguishable from those of all other fiber types. In order to identify the strongest separating factors, PCA loadings were examined (Figure 1 and Supplementary Table S7). We confirmed that the four fiber type characteristic MYH isoforms were the major separators (Figure 1b), supporting their fiber type specificity. However, well-known slow and fast fiber type marker proteins, such as specific isoforms of troponin (troponin I, troponin T) [49], myomesin (myomesin-1, myomesin-2), and myozenin (myozenin-1, myozenin-2) [27,50] were found to be essential for fiber type separation. Moreover, a fiber type-specific shift in the distribution of heat shock and proteasomal proteins was observed (segregation toward type I and type IIa fibers), suggesting there is fiber type-specific protein expression for a variety of functional groups (for further details, see Supplementary Table S8).
The current gold standard for fiber type discrimination is based on the detection of MYH isoforms with the following assignments: MYH7 for type I, MYH2 for type IIa, MYH4 for type IIb, and MYH1 for type IIx fibers. Thus, we decided to examine the quantity distribution for each MYH isoform in all fiber types using representative peptides with high-intensity and suitable detectability (Table 4, Supplementary Table S9).
Overall, for each fiber type, the characteristic MYH isoform was detected with the highest abundance but was not the only isoform present. MYH4 showed the most distinct pattern in type IIb fibers with neglectable levels of other isoforms. Type I fibers instead displayed high intensities not only for MYH7 but also for MYH2. At this stage, our results show that the expected MYH isoforms are all present in the respective fibers, but the presence of other forms is also indicated. We further focused our attention on other known marker proteins, especially sarcomeric proteins for slow and fast fibers, which have been verified in several studies [27,49,51,52,53] (Table 5, Supplementary Table S8). Type I and type IIb fibers could unequivocally be identified by having a high abundance of either slow (type I) or fast (type IIb) sarcomeric protein markers, such as different isoforms of myomesin and myozenin (Table 4). Notably, the traditionally classified fast type IIa fibers showed an enhanced concentration of slow fiber type protein markers, such as alpha-actinin-2, myozenin-2, and troponin C. Type IIx fibers could not be unambiguously classified as either slow or fast, as they displayed a high expression of marker proteins that are indicative of fast but also slow fibers.
In the next step, we performed a cluster analysis, resulting in 4 distinct protein expression clusters, (see Figure 2) and examined them by GO (gene ontology) term and pathway enrichment analyses (Supplementary Table S10, Supplementary Figures S2–S9). The formation of clusters is based on combining proteins with similar expression profiles into a common cluster. Each protein is assigned only once, resulting in unique clusters. This can provide information about the individual protein profiles of individual fiber types but also about similarities between them. As described in the literature, fast fiber types are highly glycolytic, as they need a fast but nonpersistent supply of energy. Accordingly, members of cluster 1, which were mainly annotated as enzymes involved in glycolysis, sarcomere, and cytoskeleton-associated proteins and proteins involved in calcium release and storage, were found to be associated with type IIb and type IIx fibers. Notably, cluster 1 was found to be the only cluster displaying GO terms associated with the sarcoplasmic reticulum, the central storage area for calcium within the muscle (Supplementary Table S10 and Supplementary Figures S2 and S6), underlying the unique twitching ability of fast type IIx and type IIb fibers. Members of cluster 2 (137 proteins) were found to be mainly associated with type IIa and type IIx fibers. This cluster primarily comprised proteins involved in oxidative processes. The GO term enrichment analysis revealed a high level of enrichment of mitochondria-associated terms, such as the “respiratory chain” and the “mitochondrial proton-transporting ATP synthase complex”, indicating that the oxidative energy supply is of high importance in both fiber types (Supplementary Table S10, Supplementary Figures S3 and S7). In cluster 3 (77 proteins), proteins involved in protein homeostasis (e.g., ribosomes and proteasomes) and proteins associated with the Z-disk were found (Supplementary Table S10, Supplementary Figures S4 and S8). Here, type I, type IIa, and type IIb fibers showed comparably high abundances, while type IIx fibers displayed lower values. Type I fibers known to contain a large number of mitochondria displayed a large number of proteins involved in mitochondrial processes assigned to cluster 4, such as essential energetic processes like the respiratory chain and fatty acid beta oxidation (cluster 4, 271 protein members (Supplementary Table S10, Supplementary Figures S5 and S9)). Type IIa fibers showed comparable expression patterns in this cluster, which is in accordance with their high oxidative energetic potential. Notably, cluster 4 contained a high number of GO terms and pathway annotations associated with structural elements of the sarcomeres, cytoskeleton, and extracellular matrix. Among others, the GO terms “dystroglycan complex” and “sarcoglycan complex” showed fold enrichment scores of more than 45-fold. Additionally, the terms “M-band”, “I-band”, “A-band”, and “Z-disk” were found to be common with fold enrichments of more than 10-fold.
Although the GO term and pathway annotations of all clusters showed similarities either with respect to oxidative-energy-providing processes and/or sarcomere structure, proteins within each cluster differed, as all proteins were assigned to a cluster only once. Nevertheless, the GO term and pathway analysis using DAVID resulted in multiple annotations for several proteins. Therefore, we decided to manually annotate each protein into a single group based on pathway or localization, and we calculated their percentage values with respect to cluster size (Table 6 and Supplementary Table S10a). Our manual annotation of cluster members revealed that cluster 2 exhibited the highest number of proteins associated with respiratory chain complexes (53 proteins in cluster 2 (37.23%), 4 proteins in cluster 3 (5.19%), and 32 proteins in cluster 4 (7.38%)). Cluster 4, instead, displayed a remarkably high number of proteins associated with fatty acid beta oxidation (14 proteins, (5.17%)). Although the content of proteins associated with structural components of the sarcomere and the cytoskeleton were similar in cluster 1 and cluster 4, our manual annotation revealed that cluster 1 mainly exhibited proteins associated with fast-twitch muscles, such as fast troponin isoforms and myosin-4. Cluster 4, instead, displayed enhanced intensities for slow fiber type markers and components of the extracellular matrix (Supplementary Table S10a).
In addition to this, previously unknown characteristics were revealed, such as increased mitochondrial protein content in type IIx fibers. In comparison with type IIb fibers, all other fiber types showed a significantly greater abundance of proteins that are essential for the oxidative supply of ATP, a feature that the literature mainly attributes to slow type I fibers and, to a lesser extent, to fast type IIa fibers. According to our data, type IIx fibers behave more as intermediate fibers: they contain the highest numbers of both oxidative and glycolytic enzymes (Supplementary Table S8). Type IIb fibers, on the other hand, showed a high percentage of glycolytic enzymes and proteins involved in calcium influx and efflux, which is typical for fast fibers. Compared with all other groups, type IIb fibers exhibited the clearest fast fiber phenotype. This proteomic classification indicates a high level of similarity between slow-twitch type I and fast-twitch type IIa fibers and provides novel information on the intermediate nature of type IIx fibers (see Table 7).

3.3. Definition of a Discriminative Peptide Panel for Optimized Fiber Type Assignment

To further optimize the clear allocation of the four different fiber types, we used a hypothesis-free machine learning-based feature selection and classification approach with random forests [54] (i.e., random forest-based recursive feature elimination [43,55]) to clearly define the discriminatory protein/peptide panel. To do so, we established the Random Forest-based Muscle Fiber Peptide Selection Pipeline (FiPSPi) script to identify peptides that could differentiate between the four fiber types based on our data set. The best discrimination panel consisted of one alpha-actinin-3 (TINEVENQVLTR) and one myosin-1 (IAEQELLDASER) peptide, two well-known fiber type-specific markers for fast fiber types. With this panel, an average classification accuracy value of 99.94% was obtained as well as improved fiber type clustering, giving clear separation of type I and type IIa fibers (see Figure 3b). Moreover, the alpha-actinin-3 peptide was found to differentiate well between type IIa, type IIb, and type IIx fibers, and myosin-1 peptide differentiates well between type I and all other fiber types, especially type IIa (Figure 3c,d).
To further validate the discriminatory potential of these two peptides, we developed a second confirmatory script for peptide validation. To validate our discriminative peptides, an independent public data set, in which single muscle fibers were mechanically dissected and retrospectively annotated, was used [12]. Detailed information about the criteria applied to the public data set and the subsequent samples taken for validation can be found in Appendix A Table A1 and Supplementary Table S11. As the fiber type isolation technique and mass spectrometric methods used differed from our approach, we expected significantly decreased classification results for this data set, where a classification accuracy of about 75–80% would be sufficient to validate the panel. However, an average classification accuracy of 87.2% was achieved using the machine learning-based classification approach described in this study. Here, the alpha-actinin-3 peptide was able to distinguish between type I and the other fiber types, and the myosin-1 peptide could differentiate between type IIx and all other fiber types (Supplementary Figure S11a,b). To make the panel more intuitively applicable for the reader, we further constructed a decision tree for both data sets (Supplementary Figures S10 and S11c), which uses optimized threshold values between the fiber type distributions of the two peptides, which are indicated by the well-separable boxplots presented in Figures S10c,d and S11a,b as classification rules. As can be seen in Supplementary Figure S10, the decision tree-based approach works well for our data set since no ambiguities can be observed in the terminal nodes. Our classification rules could be generalized to other data sets, and this may provide a valuable tool for further studies and clinical applications. To allow this, we commented on the associated R code and made it available in the Supplementary Materials (FiPSPi.zip), as well as in a GitHub repository (https://github.com/mpc-bioinformatics/FiPSPi, accessed on 11 November 2019), in which we aim to generalize the code for other fiber type-based proteomic data sets and thereby enable the further identification of discriminative peptides.

3.4. Higher Expression Levels of Muscular Disease Associated Proteins in Type I and Type IIa Fibers

Skeletal muscle fibers react differently to changes in life, for instance, aging and exercise, but also in the context of neuromuscular disorders. In order to assess the abundance of various myopathy-related proteins in the four fiber types, we investigated the expression patterns of dystrophin and several myofibrillar myopathy-associated proteins in more detail. The latter group comprised desmin [56], αB-Crystallin [57], myotilin [57,58], Z-band alternatively spliced PDZ motif-containing protein (ZASP) [59,60], filamin-C [61], BAG family molecular chaperone regulator 3 [62], the four-and-a-half-LIM protein 1 [63], plectin, alpha actin 1, heat shock protein B8, and DnaJ heat shock protein family member B6 [64]. Our detailed analyses of the profiles of these proteins revealed that dystrophin displayed enhanced expression levels in type I and type IIa fibers. Additionally, proteins involved in forming the dystrophin-associated glycoprotein complex, such as dystrophin, beta-sarcoglycan, delta-sarcoglycan, and gamma-sarcoglycan, were all found to be over-represented in slow type I fibers when compared with fast type IIb fibers (Supplementary Table S8). With respect to myofibrillar myopathy-associated proteins, alpha-crystallin B, desmin, myotilin, BAG family molecular chaperone regulator 3, and heat shock protein beta-8 showed enhanced expression in type I and type IIa fibers. Filamin C and plectin showed comparable expression levels for type I, type IIa, and type IIx fibers (Table 8). Alpha actin instead displayed the highest level of expression in type I fibers, and the four-and-a-half-LIM protein 1, as well as the Z-band, alternatively spliced PDZ motif-containing protein, could not be identified in the data set. We further verified the MS results for desmin and myotilin by a quantitative Western Blot approach (Figure 4) that compared the soleus (solely composed of type I and type IIa fibers) with the tibialis anterior muscle (mainly composed of type IIb and type IIx fibers). In accordance with our proteomic results, myotilin was found to have 3.7 fold higher expression in the soleus muscle, and desmin even showed enrichment of eight-fold in the soleus muscles of healthy WT mice.
In addition, sarcomere-damage-associated components showed increased expression in type I and type IIa fibers. Among these is Xin-actin-binding-repeat-containing protein 2 (Supplementary Table S8). Xin-actin-binding-repeat-containing protein 2 is not only known to selectively stabilize filamin-C [65] but has also recently been identified as a marker of muscular damage [66]. Additionally, a large number of subunits of the 26S proteasome were found to be over-represented in type I and type IIa fibers (Supplementary Table S8). These subunits are essential for the degradation of sarcomeric proteins, as the ubiquitin-proteasome-dependent proteolytic pathway is mainly responsible for the breakdown of myofibrillar proteins [67].
In summary, our analysis broadens our insight into the specific protein composition of muscle fiber types. Our analysis highlights that type I and type IIa fibers, apart from type IIx and type IIb fibers, share similar protein profiles. Further, type IIx fibers, which are traditionally classified as highly glycolytic, showed enhanced expression of proteins associated with an oxidative energy supply. Moreover, our study revealed that type I and type IIa contain higher concentrations of myopathy-related proteins. Additionally, a variety of proteins involved in protein degradation and protein folding were found to be more highly expressed in these two fiber types.

4. Discussion

In this study, we set out to characterize the proteomic compositions of different subfibers in the skeletal muscle tissue of healthy mice. Our applied analysis workflow has the following advantages: (i) the protocol requires very little sample material, which makes it highly applicable for small samples, such as those obtained with skeletal muscle biopsies; (ii) our method is not dependent on postanalytical assignment of fiber types and enables ad priori distinction of the four fiber types by immunohistochemical staining and selective enrichment by laser microdissection. This strategy serves as a powerful tool that allows in-depth proteome analysis of fiber types, giving a more distinct view of metabolic differences and susceptibility to pathophysiological processes in different types of skeletal muscle fiber. Type IIa and type IIx fibers showed less distinct MYH isoform profiles, leading us to the hypothesis that these two types of fibers may take on variable functions and tasks according to their composition and function. Although both fiber types are referred to as fast, the literature states that type IIa fibers mainly sustain their energetic needs using oxidative pathways and, therefore, by containing high numbers of mitochondria, whereas type IIx fibers rely on glycolysis to produce ATP rapidly [5,10,68]. Our further analysis identified a more versatile function of type IIx fibers in the body. With our technique, we proved that type IIx fibers inherit both oxidative and glycolytic properties while exhibiting solely fast fiber type marker proteins. Although the adaptive behavior of type IIx fibers and their high content of proteins essential for the oxidative energy supply is rarely documented in the literature, a similar phenomenon has already been observed in other mammalian species, such as springbok and deer. In these species, the oxidative capacity of type IIx fibers could prevent early fatigue of muscles after sprinting [69]. Therefore, the observed dual nature of type IIx fibers might have originated from adaptive muscle requirements, stressing the versatile function of type IIx fibers. Although type IIx fibers were shown to have a dual nature, differences between type IIx, type I, and type IIa fibers in terms of their oxidative energy supply became clear during the course of our analysis. In contrast to type IIx fibers, type I and type IIa fibers displayed increased intensities in proteins involved in fatty acid beta oxidation, confirming that type I and type IIa fibers have a greater ability to provide a long-lasting supply of ATP than other fiber types by breaking down fatty acids. Thus, our results allow the conclusion that the traditional and/or immunohistochemical classification based on ATPase staining and/or MYH isoform expression reflects only some of the molecular characteristics of the different fiber types. Based on our proteomic data, the current classification of fiber types can further be refined by our new insights into the abundances of proteins involved in the metabolic processes of muscle fibers.
Furthermore, we aimed to develop an unbiased approach to reclassify our fiber types. Our machine-learning strategy FiPSPi that was developed in-house allowed us to optimize fiber type assignment, resulting in the selection of a discriminative two-peptide panel. Using this tool, we are able to determine every fiber type at the peptide level. FiPSPi further showed the ability to identify type I, type IIx, and TB fibers correctly in an independent publicly available proteomic single fiber type data set that focused on a different approach of single fiber dissection and annotation. Unfortunately, the group of type IIa fibers could not be included in the validation study because only three samples remained after the application of our stringent criteria (for further information, see Appendix A), which was not sufficient for subsequent statistical testing. However, for type I, type IIx, and type IIb fibers, high classification accuracy was reached, validating our identified discriminative peptide panel. We have deposited the complete FiPSPi code consisting of two R Scripts with a detailed description in the Supplementary Material and in a GitHub repository, which can be reached via https://github.com/mpc-bioinformatics/FiPSPi, accessed on 11 November 2019. The first script is used to identify suitable panels of discriminative peptides between fiber types and is currently built on our data set. The second script is used to validate the discrimination of fiber types by previously identified peptides based on independent proteomic data sets and can already be used by other researchers. Furthermore, we plan to generalize the first script to enable other researchers to analyze their proteomic fiber type-based data sets and thereby potentially identify further discriminative peptides.
We further aimed to investigate our proteomic fiber type data with respect to currently unknown unique fiber type features, as the literature regularly indicates the unique properties of skeletal muscle fiber types and their high adaptability to lifestyle changes [19,28,70,71,72,73,74] and their differential levels of susceptibility in muscular disorders [24,26,75,76].
Strikingly, proteins are known to be involved or mutated in neuromuscular disorders such as myofibrillar myopathies (MFM) type I, type IIa [23,24,60,62,64,77,78], or Duchenne muscular dystrophy (DMD) [79] show a high level of enrichment in type I and type IIa fibers. In our study, we identify five MFM-causing proteins as being enriched in type I fibers (desmin, myotilin, filamin-C, BAG family molecular chaperone regulator 3, and heat shock protein beta-8). This observation seems somehow counterintuitive, as a higher abundance of disease-relevant proteins in type I fibers leads to the conclusion that these fiber types should be potentially affected more severely. Additionally, type I and type IIa fibers likewise have higher mitochondrial contents. Mitochondrial pathology has already been identified for a large number of neuromuscular diseases, which are characterized by extensive mitochondrial abnormalities underlying their potential higher vulnerability during disease progression [30,80]. Indeed, the reason for the frequently observed phenomenon of type I fiber predominance still remains largely unexplained, and studies focusing on this phenomenon have largely concluded that there is selective wasting of fast T2 fibers [26,81].
Nevertheless, our further in-depth proteomic profiling of fibers also revealed an enrichment of proteins essential for the proper folding of sarcomeric components as well as their degradation in type I and type IIa fibers. This leads us to the vague hypothesis that type I and type IIa fibers might inherit superior compensation or advanced rescue mechanisms, preventing the accumulation of misfolded proteins present in several neuromuscular diseases. Indeed, at this level, we cannot say with certainty why type I fiber predominance regularly occurs. Further research is needed, e.g., combining omics techniques and physiological/functional studies, to understand why the phenomenon of type I fiber predominance is observed in numerous diseases.
Finally, it should be noted that our approach also works optimally for clinical samples, as it requires only minimal amounts of tissue and protein and may be further optimized by establishing label-free-based detection of fiber types by coherent anti-Stokes Raman scattering as described in [82]. Needle biopsies of muscles are one of the most common procedures used to study the histological effects of lifestyle changes on skeletal muscles. Deciphering the proteomic composition has emerged as a valid tool to identify the molecular mechanisms involved in skeletal muscles during life progression, and the information gained could be used in clinical applications [11,32]. Our approach may enhance human fiber type analysis by revealing distinct fiber type behaviors.

Supplementary Materials

The following are available online at https://0-www-mdpi-com.brum.beds.ac.uk/article/10.3390/proteomes9020028/s1: Figure S1: MYH Isoform Staining. Figure S2: Bubbleplot_GO Term CC Cluster 1. Figure S3: Bubbleplot_GO Term CC Cluster 2. Figure S4: Bubbleplot_GO Term CC Cluster 3. Figure S5: Bubbleplot_GO Term CC Cluster 4. Figure S6: Bubbleplot_Pathway Cluster 1. Figure S7: Bubbleplot_Pathway Cluster 2. Figure S8: Bubbleplot_Pathway Cluster 3. Figure S9: Bubbleplot_Pathway Cluster 4. Figure S10: Decision Tree. Figure S11: Peptide_Panel_Decision_Tree_Public_Data_Set. Table S1: Spectral library information. Table S2: Spectral library. Table S3: Study annotation. Table S4: Protein_Report_Fiber_types_Spectronaut_Pulsar_no_norm. Table S5: Peptide_Report_Fiber_types_Spectronaut_Pulsar_no_norm. Table S6: Peptide_loess_originalscale_for_aLFQ. Table S7: PCA_loading. Table S8: aLFQ. Table S9: Peptides_MYH_isoform_distribution. Table S10: GO_Term_Pathway_analysis_DAVID. Table S10a: Functional_annotation_Clusters. Table S11: Public Data set_Samples. FiPSPi.zip. Rcode aLFQ.

Author Contributions

Conceptualization, B.E., K.B., R.S., C.S.C. and K.M.; Data curation, B.E. and M.T.; Formal analysis, B.E. and K.S.; Funding acquisition, M.E., R.S., C.S.C. and K.M.; Investigation, B.E., K.S. and M.T.; Methodology, B.E., K.S., M.T., R.S. and C.S.C.; Project administration, R.S., C.S.C. and K.M.; Resources, M.E., R.S. and C.S.C.; Software, B.E., K.S., M.T. and M.E.; Supervision, R.S., C.S.C. and K.M.; Validation, B.E.; Visualization, B.E., K.S. and M.T.; Writing—original draft, B.E., K.S. and M.T.; Writing—review and editing, B.E., M.T., K.B., R.S., C.S.C. and K.M. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by de.NBI, a project of the German Federal Ministry of Education and Research (BMBF) (grant number FKZ 031 A 534A) and P.U.R.E. (Protein Research Unit Ruhr within Europe) and Center for Protein Diagnostics (ProDi) grants, both from the Ministry of Innovation, Science and Research of North Rhine-Westphalia, Germany.

Institutional Review Board Statement

Mice were handled in accordance with the German Animal Welfare Act and the German Regulation for the protection of animals used for experimental purposes or other scientific purposes. All the investigations were approved by the responsible governmental animal care and use office (North Rhine-Westphalia State Agency for Nature, Environment and Consumer Protection (LANUV), Recklinghausen, Germany; reference number 84-02.05.40.14.057).

Data Availability Statement

The mass spectrometry proteomics data have been deposited into the ProteomeXchange Consortium via the PRIDE [34] partner repository with the data set identifier PXD025359.

Acknowledgments

The monoclonal antibodies BA-F8, BF-F3 and SC-71, deposited by S. Schiaffino, and the antibody 6H1, deposited by C. Lukas, were obtained from the Developmental Studies Hybridoma Bank, created by the NICHD of the NIH, and maintained at The University of Iowa, Department of Biology, Iowa City, IA 52242, USA.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

FiPSPi Validation of Public Data Set Exclusion Criteria

As an independent validation data set for the identified two-peptide panel, search files from Murgia et al. [12] were downloaded from the PRIDE repository (PXD001641). The unnormalized peptide intensities were normalized using the same method as that described in Section 2.8 of the Materials and Methods Section in the main manuscript (i.e., LOESS normalization). Due to the metadata available as supplementary information for the original paper by Murgia et al., not all samples from the original data set were appropriate for panel confirmation within this work. Consequently, exclusion criteria were formulated to exclude particular samples from the validation data set that were not appropriate and/or useful for the confirmation of the identified two-peptide panel described in the main manuscript.
The following exclusion criteria were applied:
  • Samples without a predominant MYH isoform and/or without an unambiguous fiber type in the column “Fiber type, renamed” in the Supplementary Material for Murgia et al. (Supplementary Table S1, file embr201439757-sup-0002-supp_tables1.xlsx) were excluded;
  • Samples with at least one missing value in our reanalysis results for peptide validation were excluded;
  • After application of the first two exclusion criteria, group “type IIa” had only three samples. Since this group was too small for the training of a machine learning-based model, the whole group was excluded from validation;
  • Finally, samples with a weakly predominant MYH isoform and that were outliers in the principal component analysis were excluded.
As a result, a total of 26 samples were excluded due to the above criteria. The remaining validation data set includes 29 samples and three groups of muscle fibers: 9 in group “type I”, 12 in group “T2b”, and 8 in group “type IIx”. The following table shows all 29 samples from Murgia et al. that were included in the validation data set.
Table A1. Samples in the validation data set. In this table, all 29 samples that were included in the validation data set for the identified two-peptide panel for the discrimination of muscle fiber types are listed. In the column “Original sample ID”, the original sample IDs used by Murgia et al. and for the raw files in PRIDE are shown. In the column “Sample ID within this work”, the standardized sample IDs used within this work are shown.
Table A1. Samples in the validation data set. In this table, all 29 samples that were included in the validation data set for the identified two-peptide panel for the discrimination of muscle fiber types are listed. In the column “Original sample ID”, the original sample IDs used by Murgia et al. and for the raw files in PRIDE are shown. In the column “Sample ID within this work”, the standardized sample IDs used within this work are shown.
SampleOriginal Sample IDSample ID Within This Work
11107fib_E34Type IIb_2
21107fib_S1Type I_1
3Fib_1011E26Type IIx_1
4Fib_1011E8aType IIx_2
5Fib_1011E8bType IIx_3
6Fib_1011Sol2Type I_2
7Fib_1011Sol3Type IIx_4
8SA_0430_e10Type IIb_4
9SA_0430_e5Type IIb_5
10SA_0430_S11Type I_3
11SA_1211E111Type IIb_7
12SA_1211E112Type IIb_8
13SA_1211E92Type IIx_6
14SA_1211S15Type I_4
15SA_24Ott_25EType IIb_9
16SA_24Ott_e14Type IIb_10
17SA_fib_0108_EDL1_2Type IIb_11
18SA_fib_0108_EDL2_1Type IIx_8
19SA_fib_0108_EDL2_2Type IIx_9
20SA_Fib_0716_E13_130718153850Type IIb_12
21SA_Fib_0716_E9Type IIb_13
22SA_fib_0716_S61Type I_5
23SA_fib_0716_S62Type I_6
24SA_fib_0719_14EType IIb_14
25SA_fib_0719_18EType IIb_15
26SA_fib_0719_S1Type I_7
27SA_Fib_1128_E10Type IIx_10
28SA_Fib_1128_S21Type I_8
29SA_Fib0114_S8Type I_10

References

  1. Burke, R.E.; Levine, D.N.; Zajac, F.E. Mammalian motor units: Physiological-histochemical correlation in three types in cat gastrocnemius. Science 1971, 174, 709–712. [Google Scholar] [CrossRef] [PubMed]
  2. Schiaffino, S.; Reggiani, C. Myosin isoforms in mammalian skeletal muscle. J. Appl. Physiol. 1994, 77, 493–501. [Google Scholar] [CrossRef]
  3. Schiaffino, S.; Gorza, L.; Sartore, S.; Saggin, L.; Ausoni, S.; Vianello, M.; Gundersen, K.; Lømo, T. Three myosin heavy chain isoforms in type 2 skeletal muscle fibres. J. Muscle Res. Cell Motil. 1989, 10, 197–205. [Google Scholar] [CrossRef] [PubMed]
  4. Mishra, P.; Varuzhanyan, G.; Pham, A.H.; Chan, D.C. Mitochondrial Dynamics Is a Distinguishing Feature of Skeletal Muscle Fiber Types and Regulates Organellar Compartmentalization. Cell Metab. 2015, 22, 1033–1044. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Schiaffino, S.; Reggiani, C. Fiber types in mammalian skeletal muscles. Physiol. Rev. 2011, 91, 1447–1531. [Google Scholar] [CrossRef] [Green Version]
  6. DeNardi, C.; Ausoni, S.; Moretti, P.; Gorza, L.; Velleca, M.; Buckingham, M.; Schiaffino, S. Type 2X-myosin heavy chain is coded by a muscle fiber type-specific and developmentally regulated gene. J. Cell Biol. 1993, 123, 823–835. [Google Scholar] [CrossRef] [PubMed]
  7. Pette, D.; Staron, R.S. Myosin isoforms, muscle fiber types, and transitions. Microsc. Res. Tech. 2000, 50, 500–509. [Google Scholar] [CrossRef]
  8. Stuart, C.A.; Stone, W.L.; Howell, M.E.A.; Brannon, M.F.; Hall, H.K.; Gibson, A.L.; Stone, M.H. Myosin content of individual human muscle fibers isolated by laser capture microdissection. Am. J. Physiol. Cell Physiol. 2016, 310, C381-9. [Google Scholar] [CrossRef] [Green Version]
  9. Staron, R.S.; Pette, D. The continuum of pure and hybrid myosin heavy chain-based fibre types in rat skeletal muscle. Histochemistry 1993, 100, 149–153. [Google Scholar] [CrossRef]
  10. Pette, D.; Staron, R.S. Mammalian Skeletal Muscle Fiber Type Transitions. Int. Rev. Cytol. 1997, 170, 143–223. [Google Scholar] [CrossRef]
  11. Maerkens, A.; Olivé, M.; Schreiner, A.; Feldkirchner, S.; Schessl, J.; Uszkoreit, J.; Barkovits, K.; Güttsches, A.K.; Theis, V.; Eisenacher, M.; et al. New insights into the protein aggregation pathology in myotilinopathy by combined proteomic and immunolocalization analyses. Acta Neuropathol. Commun. 2016, 4, 8. [Google Scholar] [CrossRef] [Green Version]
  12. Murgia, M.; Nagaraj, N.; Deshmukh, A.S.; Zeiler, M.; Cancellara, P.; Moretti, I.; Reggiani, C.; Schiaffino, S.; Mann, M. Single muscle fiber proteomics reveals unexpected mitochondrial specialization. EMBO Rep. 2015, 16, 387–395. [Google Scholar] [CrossRef] [PubMed]
  13. Murgia, M.; Toniolo, L.; Nagaraj, N.; Ciciliot, S.; Vindigni, V.; Schiaffino, S.; Reggiani, C.; Mann, M. Single Muscle Fiber Proteomics Reveals Fiber-Type-Specific Features of Human Muscle Aging. Cell Rep. 2017, 19, 2396–2409. [Google Scholar] [CrossRef] [Green Version]
  14. Albers, P.H.; Pedersen, A.J.T.; Birk, J.B.; Kristensen, D.E.; Vind, B.F.; Baba, O.; Nøhr, J.; Højlund, K.; Wojtaszewski, J.F.P. Human muscle fiber type-specific insulin signaling: Impact of obesity and type 2 diabetes. Diabetes 2015, 64, 485–497. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Bagley, J.R.; McLeland, K.A.; Arevalo, J.A.; Brown, L.E.; Coburn, J.W.; Galpin, A.J. Skeletal Muscle Fatigability and Myosin Heavy Chain Fiber Type in Resistance Trained Men. J. Strength Cond. Res. 2017, 31, 602–607. [Google Scholar] [CrossRef] [PubMed]
  16. Bassel-Duby, R.; Olson, E.N. Signaling pathways in skeletal muscle remodeling. Annu. Rev. Biochem. 2006, 75, 19–37. [Google Scholar] [CrossRef] [Green Version]
  17. Carroll, K.M.; Bazyler, C.D.; Bernards, J.R.; Taber, C.B.; Stuart, C.A.; DeWeese, B.H.; Sato, K.; Stone, M.H. Skeletal Muscle Fiber Adaptations Following Resistance Training Using Repetition Maximums or Relative Intensity. Sports 2019, 7, 169. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  18. Frontera, W.R.; Zayas, A.R.; Rodriguez, N. Aging of human muscle: Understanding sarcopenia at the single muscle cell level. Phys. Med. Rehabil. Clin. N. Am. 2012, 23, 201–207. [Google Scholar] [CrossRef] [Green Version]
  19. Lexell, J.; Henriksson-Larsén, K.; Winblad, B.; Sjöström, 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]
  20. Piec, I.; Listrat, A.; Alliot, J.; Chambon, C.; Taylor, R.G.; Bechet, D. Differential proteome analysis of aging in rat skeletal muscle. FASEB J. 2005, 19, 1143–1145. [Google Scholar] [CrossRef]
  21. Purves-Smith, F.M.; Sgarioto, N.; Hepple, R.T. Fiber typing in aging muscle. Exerc. Sport Sci. Rev. 2014, 42, 45–52. [Google Scholar] [CrossRef] [PubMed]
  22. Kyriakides, T.; Silberstein, J.M.; Jongpiputvanich, S.; Silberstein, E.P.; Walsh, P.J.; Gubbay, S.S.; Kakulas, B.A. The clinical significance of type 1 fiber predominance. Muscle Nerve 1993, 16, 418–423. [Google Scholar] [CrossRef] [PubMed]
  23. Claeys, K.G.; van der Ven, P.F.M.; Behin, A.; Stojkovic, T.; Eymard, B.; Dubourg, O.; Laforêt, P.; Faulkner, G.; Richard, P.; Vicart, P.; et al. Differential involvement of sarcomeric proteins in myofibrillar myopathies: A morphological and immunohistochemical study. Acta Neuropathol. 2009, 117, 293–307. [Google Scholar] [CrossRef] [PubMed]
  24. Kley, R.A.; Hellenbroich, Y.; van der Ven, P.F.M.; Fürst, D.O.; Huebner, A.; Bruchertseifer, V.; Peters, S.A.; Heyer, C.M.; Kirschner, J.; Schröder, R.; et al. Clinical and morphological phenotype of the filamin myopathy: A study of 31 German patients. Brain 2007, 130, 3250–3264. [Google Scholar] [CrossRef] [PubMed]
  25. Kley, R.A.; Serdaroglu-Oflazer, P.; Leber, Y.; Odgerel, Z.; van der Ven, P.F.M.; Olivé, M.; Ferrer, I.; Onipe, A.; Mihaylov, M.; Bilbao, J.M.; et al. Pathophysiology of protein aggregation and extended phenotyping in filaminopathy. Brain 2012, 135, 2642–2660. [Google Scholar] [CrossRef]
  26. Webster, C.; Silberstein, L.; Hays, A.P.; Blau, H.M. Fast muscle fibers are preferentially affected in Duchenne muscular dystrophy. Cell 1988, 52, 503–513. [Google Scholar] [CrossRef]
  27. Chemello, F.; Bean, C.; Cancellara, P.; Laveder, P.; Reggiani, C.; Lanfranchi, G. Microgenomic analysis in skeletal muscle: Expression signatures of individual fast and slow myofibers. PLoS ONE 2011, 6, e16807. [Google Scholar] [CrossRef] [Green Version]
  28. 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]
  29. Kley, R.A.; Maerkens, A.; Leber, Y.; Theis, V.; Schreiner, A.; van der Ven, P.F.M.; Uszkoreit, J.; Stephan, C.; Eulitz, S.; Euler, N.; et al. A combined laser microdissection and mass spectrometry approach reveals new disease relevant proteins accumulating in aggregates of filaminopathy patients. Mol. Cell. Proteom. 2013, 12, 215–227. [Google Scholar] [CrossRef] [Green Version]
  30. Winter, L.; Wittig, I.; Peeva, V.; Eggers, B.; Heidler, J.; Chevessier, F.; Kley, R.A.; Barkovits, K.; Strecker, V.; Berwanger, C.; et al. Mutant desmin substantially perturbs mitochondrial morphology, function and maintenance in skeletal muscle tissue. Acta Neuropathol. 2016, 132, 453–473. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Deshmukh, A.S.; Steenberg, D.E.; Hostrup, M.; Birk, J.B.; Larsen, J.K.; Santos, A.; Kjøbsted, R.; Hingst, J.R.; Schéele, C.C.; Murgia, M.; et al. Deep muscle-proteomic analysis of freeze-dried human muscle biopsies reveals fiber type-specific adaptations to exercise training. Nat. Commun. 2021, 12, 304. [Google Scholar] [CrossRef]
  32. Maerkens, A.; Kley, R.A.; Olivé, M.; Theis, V.; van der Ven, P.F.M.; Reimann, J.; Milting, H.; Schreiner, A.; Uszkoreit, J.; Eisenacher, M.; et al. Differential proteomic analysis of abnormal intramyoplasmic aggregates in desminopathy. J. Proteom. 2013, 90, 14–27. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Oertzen-Hagemann, V.; Kirmse, M.; Eggers, B.; Pfeiffer, K.; Marcus, K.; de Marées, M.; Platen, P. Effects of 12 Weeks of Hypertrophy Resistance Exercise Training Combined with Collagen Peptide Supplementation on the Skeletal Muscle Proteome in Recreationally Active Men. Nutrients 2019, 11, 1072. [Google Scholar] [CrossRef] [Green Version]
  34. Perez-Riverol, Y.; Csordas, A.; Bai, J.; Bernal-Llinares, M.; Hewapathirana, S.; Kundu, D.J.; Inuganti, A.; Griss, J.; Mayer, G.; Eisenacher, M.; et al. The PRIDE database and related tools and resources in 2019: Improving support for quantification data. Nucleic Acids Res. 2019, 47, D442–D450. [Google Scholar] [CrossRef] [PubMed]
  35. Ritchie, M.E.; Phipson, B.; Wu, D.; Hu, Y.; Law, C.W.; Shi, W.; Smyth, G.K. Limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015, 43, e47. [Google Scholar] [CrossRef]
  36. Bolstad, B.M.; Irizarry, R.A.; Astrand, M.; Speed, T.P. A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics 2003, 19, 185–193. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Rosenberger, G.; Ludwig, C.; Röst, H.L.; Aebersold, R.; Malmström, L. aLFQ: An R-package for estimating absolute protein quantities from label-free LC-MS/MS proteomics data. Bioinformatics 2014, 30, 2511–2513. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  38. Schwanhäusser, B.; Busse, D.; Li, N.; Dittmar, G.; Schuchhardt, J.; Wolf, J.; Chen, W.; Selbach, M. Global quantification of mammalian gene expression control. Nature 2011, 473, 337–342. [Google Scholar] [CrossRef] [Green Version]
  39. Tyanova, S.; Temu, T.; Sinitcyn, P.; Carlson, A.; Hein, M.Y.; Geiger, T.; Mann, M.; Cox, J. The Perseus computational platform for comprehensive analysis of (prote)omics data. Nat. Methods 2016, 13, 731–740. [Google Scholar] [CrossRef]
  40. Huang, D.W.; Sherman, B.T.; Tan, Q.; Collins, J.R.; Alvord, W.G.; Roayaei, J.; Stephens, R.; Baseler, M.W.; Lane, H.C.; Lempicki, R.A. The DAVID Gene Functional Classification Tool: A novel biological module-centric algorithm to functionally analyze large gene lists. Genome Biol. 2007, 8, R183. [Google Scholar] [CrossRef] [Green Version]
  41. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Bioinformatics enrichment tools: Paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009, 37, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Huang, D.W.; Sherman, B.T.; Lempicki, R.A. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat. Protoc. 2009, 4, 44–57. [Google Scholar] [CrossRef] [PubMed]
  43. Turewicz, M.; Kohl, M.; Ahrens, M.; Mayer, G.; Uszkoreit, J.; Naboulsi, W.; Bracht, T.; Megger, D.A.; Sitek, B.; Marcus, K.; et al. BioInfra.Prot: A comprehensive proteomics workflow including data standardization, protein inference, expression analysis and data publication. J. Biotechnol. 2017, 261, 116–125. [Google Scholar] [CrossRef] [PubMed]
  44. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2018. [Google Scholar]
  45. Liaw, A.; Wiener, M. Classification and Regression by randomForest. R News 2002, 2, 18–22. [Google Scholar]
  46. Therneau, T.; Atkinson, B. rpart: Recursive Partitioning and Regression Trees. R Package Version 4.1-15. Available online: https://CRAN.R-project.org/package=rpart (accessed on 11 November 2019).
  47. Taylor, S.C.; Posch, A. The design of a quantitative western blot experiment. Biomed Res. Int. 2014, 2014, 361590. [Google Scholar] [CrossRef]
  48. Faden, F.; Eschen-Lippold, L.; Dissmeyer, N. Normalized Quantitative Western Blotting Based on Standardized Fluorescent Labeling. Methods Mol. Biol. 2016, 1450, 247–258. [Google Scholar] [CrossRef]
  49. Dhoot, G.K.; Perry, S.V. Distribution of polymorphic forms of troponin components and tropomyosin in skeletal muscle. Nature 1979, 278, 714–718. [Google Scholar] [CrossRef]
  50. Lynch, C.J.; Xu, Y.; Hajnal, A.; Salzberg, A.C.; Kawasawa, Y.I. RNA sequencing reveals a slow to fast muscle fiber type transition after olanzapine infusion in rats. PLoS ONE 2015, 10, e0123966. [Google Scholar] [CrossRef]
  51. Gahlmann, R.; Wade, R.; Gunning, P.; Kedes, L. Differential expression of slow and fast skeletal muscle troponin C. Slow skeletal muscle troponin C is expressed in human fibroblasts. J. Mol. Biol. 1988, 201, 379–391. [Google Scholar] [CrossRef]
  52. Hallauer, P.L.; Bradshaw, H.L.; Hastings, K.E. Complex fiber-type-specific expression of fast skeletal muscle troponin I gene constructs in transgenic mice. Development 1993, 119, 691–701. [Google Scholar] [CrossRef]
  53. Deshmukh, A.S.; Murgia, M.; Nagaraj, N.; Treebak, J.T.; Cox, J.; Mann, M. Deep proteomics of mouse skeletal muscle enables quantitation of protein isoforms, metabolic pathways, and transcription factors. Mol. Cell. Proteom. 2015, 14, 841–853. [Google Scholar] [CrossRef] [Green Version]
  54. Breiman, L. Random Forests. Mach. Learn. 2001, 45, 5–32. [Google Scholar] [CrossRef] [Green Version]
  55. Jiang, H.; Deng, Y.; Chen, H.-S.; Tao, L.; Sha, Q.; Chen, J.; Tsai, C.-J.; Zhang, S. Joint analysis of two microarray gene-expression data sets to select lung adenocarcinoma marker genes. BMC Bioinform. 2004, 5, 81. [Google Scholar] [CrossRef] [Green Version]
  56. Goldfarb, L.G.; Vicart, P.; Goebel, H.H.; Dalakas, M.C. Desmin myopathy. Brain 2004, 127, 723–734. [Google Scholar] [CrossRef] [PubMed]
  57. Vicart, P.; Caron, A.; Guicheney, P.; Li, Z.; Prévost, M.C.; Faure, A.; Chateau, D.; Chapon, F.; Tomé, F.; Dupret, J.M.; et al. A missense mutation in the alphaB-crystallin chaperone gene causes a desmin-related myopathy. Nat. Genet. 1998, 20, 92–95. [Google Scholar] [CrossRef]
  58. Salmikangas, P.; Mykkänen, O.M.; Grönholm, M.; Heiska, L.; Kere, J.; Carpén, O. Myotilin, a novel sarcomeric protein with two Ig-like domains, is encoded by a candidate gene for limb-girdle muscular dystrophy. Hum. Mol. Genet. 1999, 8, 1329–1336. [Google Scholar] [CrossRef]
  59. Faulkner, G.; Pallavicini, A.; Formentin, E.; Comelli, A.; Ievolella, C.; Trevisan, S.; Bortoletto, G.; Scannapieco, P.; Salamon, M.; Mouly, V.; et al. ZASP: A new Z-band alternatively spliced PDZ-motif protein. J. Cell Biol. 1999, 146, 465–475. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Selcen, D.; Ohno, K.; Engel, A.G. Myofibrillar myopathy: Clinical, morphological and genetic studies in 63 patients. Brain 2004, 127, 439–451. [Google Scholar] [CrossRef] [PubMed]
  61. Vorgerd, M.; van der Ven, P.F.M.; Bruchertseifer, V.; Löwe, T.; Kley, R.A.; Schröder, R.; Lochmüller, H.; Himmel, M.; Koehler, K.; Fürst, D.O.; et al. A mutation in the dimerization domain of filamin c causes a novel type of autosomal dominant myofibrillar myopathy. Am. J. Hum. Genet. 2005, 77, 297–304. [Google Scholar] [CrossRef] [Green Version]
  62. Selcen, D. Myofibrillar myopathies. Neuromuscul. Disord. 2011, 21, 161–171. [Google Scholar] [CrossRef] [Green Version]
  63. Quinzii, C.M.; Vu, T.H.; Min, K.C.; Tanji, K.; Barral, S.; Grewal, R.P.; Kattah, A.; Camaño, P.; Otaegui, D.; Kunimatsu, T.; et al. X-linked dominant scapuloperoneal myopathy is due to a mutation in the gene encoding four-and-a-half-LIM protein 1. Am. J. Hum. Genet. 2008, 82, 208–213. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  64. Kley, R.A.; Olivé, M.; Schröder, R. New aspects of myofibrillar myopathies. Curr. Opin. Neurol. 2016, 29, 628–634. [Google Scholar] [CrossRef] [PubMed]
  65. Leber, Y.; Ruparelia, A.A.; Kirfel, G.; van der Ven, P.F.M.; Hoffmann, B.; Merkel, R.; Bryson-Richardson, R.J.; Fürst, D.O. Filamin C is a highly dynamic protein associated with fast repair of myofibrillar microdamage. Hum. Mol. Genet. 2016, 25, 2776–2788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Nilsson, M.I.; Nissar, A.A.; Al-Sajee, D.; Tarnopolsky, M.A.; Parise, G.; Lach, B.; Fürst, D.O.; van der Ven, P.F.M.; Kley, R.A.; Hawke, T.J. Xin is a marker of skeletal muscle damage severity in myopathies. Am. J. Pathol. 2013, 183, 1703–1709. [Google Scholar] [CrossRef]
  67. Ulbricht, A.; Gehlert, S.; Leciejewski, B.; Schiffer, T.; Bloch, W.; Höhfeld, J. Induction and adaptation of chaperone-assisted selective autophagy CASA in response to resistance exercise in human skeletal muscle. Autophagy 2015, 11, 538–546. [Google Scholar] [CrossRef] [Green Version]
  68. Pette, D.; Staron, R.S. Cellular and molecular diversities of mammalian skeletal muscle fibers. Rev. Physiol. Biochem. Pharmacol. 1990, 116, 1–76. [Google Scholar]
  69. Curry, J.W.; Hohl, R.; Noakes, T.D.; Kohn, T.A. High oxidative capacity and type IIx fibre content in springbok and fallow deer skeletal muscle suggest fast sprinters with a resistance to fatigue. J. Exp. Biol. 2012, 215, 3997–4005. [Google Scholar] [CrossRef] [Green Version]
  70. Komi, P.V.; Karlsson, J. Skeletal muscle fibre types, enzyme activities and physical performance in young males and females. Acta Physiol. Scand. 1978, 103, 210–218. [Google Scholar] [CrossRef]
  71. Larsson, L.; Biral, D.; Campione, M.; Schiaffino, S. An age-related type IIB to IIX myosin heavy chain switching in rat skeletal muscle. Acta Physiol. Scand. 1993, 147, 227–234. [Google Scholar] [CrossRef]
  72. Larsson, L.; Sjödin, B.; Karlsson, J. Histochemical and biochemical changes in human skeletal muscle with age in sedentary males, age 22–65 years. Acta Physiol. Scand. 1978, 103, 31–39. [Google Scholar] [CrossRef]
  73. Lexell, J.; Taylor, C.C.; Sjöström, M. What is the cause of the ageing atrophy? Total number, size and proportion of different fiber types studied in whole vastus lateralis muscle from 15- to 83-year-old men. J. Neurol. Sci. 1988, 84, 275–294. [Google Scholar] [CrossRef]
  74. Lievens, E.; Klass, M.; Bex, T.; Derave, W. Muscle fiber typology substantially influences time to recover from high-intensity exercise. J. Appl. Physiol. 2020, 128, 648–659. [Google Scholar] [CrossRef] [PubMed]
  75. Ljubicic, V.; Burt, M.; Jasmin, B.J. The therapeutic potential of skeletal muscle plasticity in Duchenne muscular dystrophy: Phenotypic modifiers as pharmacologic targets. FASEB J. 2014, 28, 548–568. [Google Scholar] [CrossRef]
  76. Fröhlich, T.; Kemter, E.; Flenkenthaler, F.; Klymiuk, N.; Otte, K.A.; Blutke, A.; Krause, S.; Walter, M.C.; Wanke, R.; Wolf, E.; et al. Progressive muscle proteome changes in a clinically relevant pig model of Duchenne muscular dystrophy. Sci. Rep. 2016, 6, 33362. [Google Scholar] [CrossRef]
  77. Schröder, R.; Schoser, B. Myofibrillar myopathies: A clinical and myopathological guide. Brain Pathol. 2009, 19, 483–492. [Google Scholar] [CrossRef] [PubMed]
  78. Bär, H.; Fischer, D.; Goudeau, B.; Kley, R.A.; Clemen, C.S.; Vicart, P.; Herrmann, H.; Vorgerd, M.; Schröder, R. Pathogenic effects of a novel heterozygous R350P desmin mutation on the assembly of desmin intermediate filaments in vivo and in vitro. Hum. Mol. Genet. 2005, 14, 1251–1260. [Google Scholar] [CrossRef] [PubMed]
  79. Hoffman, E.P.; Brown, R.H.; Kunkel, L.M. Dystrophin: The protein product of the duchenne muscular dystrophy locus. Cell 1987, 51, 919–928. [Google Scholar] [CrossRef]
  80. Lindgren, U.; Roos, S.; Hedberg Oldfors, C.; Moslemi, A.-R.; Lindberg, C.; Oldfors, A. Mitochondrial pathology in inclusion body myositis. Neuromuscul. Disord. 2015, 25, 281–288. [Google Scholar] [CrossRef] [PubMed]
  81. Miike, T.; Ohtani, Y.; Tamari, H.; Ishitsu, T.; Une, Y. Muscle fiber type transformation in nemaline myopathy and congenital fiber type disproportion. Brain Dev. 1986, 8, 526–532. [Google Scholar] [CrossRef]
  82. Niedieker, D.; Großerüschkamp, F.; Schreiner, A.; Barkovits, K.; Kötting, C.; Marcus, K.; Gerwert, K.; Vorgerd, M. Label-free identification of myopathological features with coherent anti-Stokes Raman scattering. Muscle Nerve 2018, 58, 456–459. [Google Scholar] [CrossRef]
Figure 1. Differentiation among the four fiber types. (a) Principal component analysis (PCA) of fiber types (n = 56). Significant separation of type IIb (purple) and type IIx fibers (green) can be observed in principle component (PC2), while type I (blue) and type IIa fibers (red) formed a combined, but significantly separated, cluster from type IIx and type IIb fibers in PC2 with moderate segregation into type I and type IIa fibers. (b) PCA loadings, visualizing the proteins separating the four groups. The four fiber type characteristic myosin heavy chain (MYH) isoforms are indicated in red and were identified as major separators.
Figure 1. Differentiation among the four fiber types. (a) Principal component analysis (PCA) of fiber types (n = 56). Significant separation of type IIb (purple) and type IIx fibers (green) can be observed in principle component (PC2), while type I (blue) and type IIa fibers (red) formed a combined, but significantly separated, cluster from type IIx and type IIb fibers in PC2 with moderate segregation into type I and type IIa fibers. (b) PCA loadings, visualizing the proteins separating the four groups. The four fiber type characteristic myosin heavy chain (MYH) isoforms are indicated in red and were identified as major separators.
Proteomes 09 00028 g001
Figure 2. Hierarchical clustering: high protein intensities are shown in red, and low protein intensities are shown in green. The profile plots indicate trends of protein expression intensity in the different fiber types. Cluster 1 displays high intensities in type IIx and type IIb fibers. Proteins of cluster 1 were annotated to glycolysis, calcium signaling, and the sarcoplasmic reticulum. High intensities in cluster 2 were annotated to type IIa and type IIx fibers and could be associated with mitochondrial processes such as the tricarboxylic acid (TCA) cycle the oxidative phosphorylation (OXPHOS) pathway and the mitochondrial respiratory chain complexes. Cluster 3 shows high-intensity values in type I, type IIa, and type IIb and is composed of ribosomal and proteasomal proteins. Cluster 4 can be annotated to structural components of the sarcomere as well as mitochondrial processes, such as the OXPHOS pathway, the TCA cycle, and the respiratory chain but also has a unique annotation to fatty acid beta oxidation and high intensities in GO terms and pathways associated with structural components of the sarcomere and the cytoskeleton.
Figure 2. Hierarchical clustering: high protein intensities are shown in red, and low protein intensities are shown in green. The profile plots indicate trends of protein expression intensity in the different fiber types. Cluster 1 displays high intensities in type IIx and type IIb fibers. Proteins of cluster 1 were annotated to glycolysis, calcium signaling, and the sarcoplasmic reticulum. High intensities in cluster 2 were annotated to type IIa and type IIx fibers and could be associated with mitochondrial processes such as the tricarboxylic acid (TCA) cycle the oxidative phosphorylation (OXPHOS) pathway and the mitochondrial respiratory chain complexes. Cluster 3 shows high-intensity values in type I, type IIa, and type IIb and is composed of ribosomal and proteasomal proteins. Cluster 4 can be annotated to structural components of the sarcomere as well as mitochondrial processes, such as the OXPHOS pathway, the TCA cycle, and the respiratory chain but also has a unique annotation to fatty acid beta oxidation and high intensities in GO terms and pathways associated with structural components of the sarcomere and the cytoskeleton.
Proteomes 09 00028 g002
Figure 3. (a) Principal component analysis (PCA) for all fiber types (n = 56) based on our complete proteomic data. Clear discrimination of type IIb and T2X fibers is possible, whereas type I and type IIa fibers cluster together. (b) PCA based on our alpha-actinin-3 and myosin-1 peptide panel. Discrimination of type I and type IIa fibers is especially enhanced. (c,d) Boxplots of the intensity value distributions on a log2-scale of the four fiber types, type I, type IIa, type IIb, and type IIx, for alpha-actinin-3 peptide (c) and the myosin-1 peptide (d), which were used in our panel. The alpha-actinin-3 peptide TINEVENQVLTR differentiates well between the type IIa, type IIb, and type IIx fibers, and the myosin-1-peptide IAEQELLDASER differentiates well between type I and all other fiber types, especially type IIa. Thus, only these two peptides are needed to differentiate all four fiber types.
Figure 3. (a) Principal component analysis (PCA) for all fiber types (n = 56) based on our complete proteomic data. Clear discrimination of type IIb and T2X fibers is possible, whereas type I and type IIa fibers cluster together. (b) PCA based on our alpha-actinin-3 and myosin-1 peptide panel. Discrimination of type I and type IIa fibers is especially enhanced. (c,d) Boxplots of the intensity value distributions on a log2-scale of the four fiber types, type I, type IIa, type IIb, and type IIx, for alpha-actinin-3 peptide (c) and the myosin-1 peptide (d), which were used in our panel. The alpha-actinin-3 peptide TINEVENQVLTR differentiates well between the type IIa, type IIb, and type IIx fibers, and the myosin-1-peptide IAEQELLDASER differentiates well between type I and all other fiber types, especially type IIa. Thus, only these two peptides are needed to differentiate all four fiber types.
Proteomes 09 00028 g003
Figure 4. (a,b) Quantitative Smart Protein Layers (SPL) Western blot analysis of whole murine soleus and tibialis anterior muscle. SPL Western blotting system (NH DyeAGNOSTICS GmbH) was used for the detection of quantitative differences between selected candidate proteins. Total protein was pre-labeled with a red fluorescent fluorophore (700 nm), and a 12.5 kDa green fluorescence-labeled standard protein (SMA S, 800 nm) was spiked in, enabling error correction of differing sample loading and data normalization between experiments. Primary antibodies against myotilin (a), myotilin (E-10), sc-393957), and desmin (b), DAKO, M0760 Clone D33) were visualized by infrared secondary antibodies (IRDye goat anti-mouse 800 nm). SPL system analysis resulted in quantitative protein volumes (SPL normalized volume), determining a higher concentration in the soleus muscle for both proteins.
Figure 4. (a,b) Quantitative Smart Protein Layers (SPL) Western blot analysis of whole murine soleus and tibialis anterior muscle. SPL Western blotting system (NH DyeAGNOSTICS GmbH) was used for the detection of quantitative differences between selected candidate proteins. Total protein was pre-labeled with a red fluorescent fluorophore (700 nm), and a 12.5 kDa green fluorescence-labeled standard protein (SMA S, 800 nm) was spiked in, enabling error correction of differing sample loading and data normalization between experiments. Primary antibodies against myotilin (a), myotilin (E-10), sc-393957), and desmin (b), DAKO, M0760 Clone D33) were visualized by infrared secondary antibodies (IRDye goat anti-mouse 800 nm). SPL system analysis resulted in quantitative protein volumes (SPL normalized volume), determining a higher concentration in the soleus muscle for both proteins.
Proteomes 09 00028 g004
Table 1. Detailed information on myosin isoforms (MYH): accession (Uniprot), protein name, protein name abbreviation, gene name, and corresponding fiber type.
Table 1. Detailed information on myosin isoforms (MYH): accession (Uniprot), protein name, protein name abbreviation, gene name, and corresponding fiber type.
Accession (Uniprot)Protein Name Protein Name AbbreviationGene NameCorresponding
Fiber Type
Q91Z83Myosin-7MYH7Myh7Type I
G3UW82Myosin, heavy polypeptide 2, skeletal muscle, adult *MYH2Myh2Type IIa
Q5SX39Myosin-4MYH4Myh4Type IIb
Q5SX40Myosin-1MYH1Myh1Type IIx
* further referred to as myosin-2.
Table 2. Sample information. Mouse identifier (ID), sex, genotype, and age in months and days for each mouse used in the present study are displayed.
Table 2. Sample information. Mouse identifier (ID), sex, genotype, and age in months and days for each mouse used in the present study are displayed.
Mouse IDSexGenotypeAge in MonthsAge in Days
1fwt21657
2fwt21657
3mwt15466
4mwt14448
5mwt14448
6mwt14448
7mwt13401
8mwt11361
9mwt11361
10mwt11361
11mwt395
12mwt395
13mwt396
14fwt396
15fwt396
Table 3. Primary antibodies (product numbers Developmental Studies Hybridoma Bank (DSHB), Iowa) and dilutions used for the detection of fiber types (type I, type IIa, type IIb, and type IIx) using the fiber type-specific myosin (MYH) isoform. Secondary antibodies and dilutions necessary for the detection of primary antibody signals.
Table 3. Primary antibodies (product numbers Developmental Studies Hybridoma Bank (DSHB), Iowa) and dilutions used for the detection of fiber types (type I, type IIa, type IIb, and type IIx) using the fiber type-specific myosin (MYH) isoform. Secondary antibodies and dilutions necessary for the detection of primary antibody signals.
Primary Antibody/DilutionSecondary Antibody/DilutionFiber TypeMYH Isoform
BA-F8/1:20TexasRed/1:60Type IMYH7
SC-71/1:20Cy2/1:200Type IIaMYH2
BF-F3/1:20AlexaFluor594/1:500Type IIbMYH4
6H1/1:15AlexaFluor594/1:500Type IIxMYH1
Table 4. Summed intensities (Sum) of unique peptides (proteotypicity = TRUE) detected for each myosin heavy chain (MYH) isoform that is indicative of fiber types (type I, type IIa, type IIb, and type IIx).
Table 4. Summed intensities (Sum) of unique peptides (proteotypicity = TRUE) detected for each myosin heavy chain (MYH) isoform that is indicative of fiber types (type I, type IIa, type IIb, and type IIx).
ProteinPeptide SequenceProteotypicitySum Type I (10 × 106)Sum Type IIa (10 × 106)Sum Type IIb (10 × 106)Sum Type IIx (10 × 106)
MYH7DTQIQLDDAVRTRUE226.5745.220.000.00
MYH7IEDEQALGSQLQKTRUE162.6534.330.050.19
MYH7ANDDLKENIAIVERTRUE137.4127.700.010.01
MYH7GQNVQQVSYAIGALAKTRUE80.4615.950.040.03
MYH7SLQSLLKTRUE74.6715.450.020.03
ProteinPeptide SequenceProteotypicitySum Type I (10 × 106)Sum Type IIa (10 × 106)Sum Type IIb (10 × 106)Sum Type IIx (10 × 106)
MYH2LINDLTTQRTRUE95.23216.188.3214.98
MYH2AAYLQGLNSADLLKTRUE77.60179.828.1711.76
MYH2VLNASAIPEGQYIDSKTRUE74.62173.666.2311.20
MYH2GQTVEQVTNAVGALAKTRUE39.80106.673.767.45
MYH2IEDEQAIGIQLQKTRUE25.5871.062.535.65
ProteinPeptide SequenceProteotypicitySum Type I (10 × 106)Sum Type IIa (10 × 106)Sum Type IIb (10 × 106)Sum Type IIx (10 × 106)
MYH4VAEQELLDASERTRUE3.783.82221.7149.27
MYH4LINELSTQKTRUE4.174.55175.1137.76
MYH4AAYLTSLNSADLLKTRUE1.872.98129.6824.67
MYH4LQDAEEHVEAVNSKTRUE1.421.3475.1017.27
MYH4RVAEQELLDASERTRUE1.522.4874.1819.12
ProteinPeptide SequenceProteotypicitySum Type I (10 × 106)Sum Type IIa (10 × 106)Sum Type IIb (10 × 106)Sum Type IIx (10 × 106)
MYH1AAYLQNLNSADLLKTRUE14.0220.5425.48122.96
MYH1SALAHALQSSRTRUE14.1920.2817.7697.12
MYH1DSLVSQLSRTRUE8.3711.4213.7949.63
MYH1NAYEESLDHLETLKRTRUE4.606.718.6841.85
MYH1QLDEKDSLVSQLSRTRUE4.596.828.7741.53
Table 5. Percentage values of sarcomeric protein markers and their indication for either slow or fast fibers in all four skeletal muscle fiber types (type I, type IIa, type IIb, and type IIx).
Table 5. Percentage values of sarcomeric protein markers and their indication for either slow or fast fibers in all four skeletal muscle fiber types (type I, type IIa, type IIb, and type IIx).
Protein NamesIndicative For% Type I% Type IIa% Type IIb% Type IIx
Alpha-actinin-2Slow2.601.890.391.15
Alpha-actinin-3Fast0.100.100.930.63
Myomesin 2Fast0.050.110.360.34
Myomesin-1Slow0.040.020.000.00
Myosin-binding protein C. slow typeSlow0.090.090.040.06
Myosin-binding protein C. fast typeFast0.080.110.640.33
Myozenin-1Fast0.100.190.300.31
Myozenin-2Slow0.200.140.010.02
Troponin C. slow skeletal and cardiac muscles (TN-C)Slow0.450.170.000.00
Troponin I. fast skeletal muscle (troponin I. fast-twitch isoform)Fast0.060.120.230.26
Table 6. Percentages of manually annotated proteins into 13 groups (Sarcomere/Cytoskeleton, Respiratory chain, Tricarboxylic acid cycle, Glycolysis, Glycogenolysis, Mitochondria, Proteasome, Heat shock, Nucleus, Fatty acid beta oxidation, Calcium regulation/transport, Ribosome, Others) for each cluster.
Table 6. Percentages of manually annotated proteins into 13 groups (Sarcomere/Cytoskeleton, Respiratory chain, Tricarboxylic acid cycle, Glycolysis, Glycogenolysis, Mitochondria, Proteasome, Heat shock, Nucleus, Fatty acid beta oxidation, Calcium regulation/transport, Ribosome, Others) for each cluster.
AnnotationCluster 1 (%)Cluster 2 (%)Cluster 3 (%)Cluster 4 (%)
Sarcomere/Cytoskeleton18.245.1114.2919.56
Respiratory chain3.3837.235.197.38
Tricarboxylic acid cycle0.687.301.301.48
Glycolysis4.730.000.000.74
Glycogenolysis5.410.730.000.00
Mitochondria0.009.491.305.90
Proteasome0.000.009.094.43
Heat shock0.000.002.603.69
Nucleus0.680.002.603.32
Fatty acid beta oxidation2.030.731.305.17
Calcium regulation/transport8.780.730.001.85
Ribosome8.784.3828.574.43
Other47.3034.3133.7742.07
Table 7. Classification of fiber types according to the traditional classification based on ATPase staining, the oxidative phosphorylation capacity, the immunohistochemical classification based on myosin heavy chain (MYH) isoforms, and the classification based on the proteomic data.
Table 7. Classification of fiber types according to the traditional classification based on ATPase staining, the oxidative phosphorylation capacity, the immunohistochemical classification based on myosin heavy chain (MYH) isoforms, and the classification based on the proteomic data.
Fiber TypeTraditional Classification (Based on ATPase Staining)Oxidative Phosphorylation CapacityImmunohistochemical
Classification via MYH Isoform
Proteomic Classification
Type I SlowVery highMYH7Highly oxidative
Type IIaFastHighMYH2Highly oxidative
Type IIbFastVery lowMYH4Highly glycolytic
Type IIxFastLow MYH1Oxidative and glycolytic
Table 8. Mass spectrometry-based calculated percentages of disease-related proteins and their distribution (in%) in skeletal muscle fiber types (type I, type IIa, type IIb, and type IIx). All proteins do show a higher abundance in type I and type IIa fibers.
Table 8. Mass spectrometry-based calculated percentages of disease-related proteins and their distribution (in%) in skeletal muscle fiber types (type I, type IIa, type IIb, and type IIx). All proteins do show a higher abundance in type I and type IIa fibers.
Protein Names% Type I% Type IIa% Type IIb% Type IIx
Actin. alpha skeletal muscle30.4924.2320.1421.98
Alpha-crystallin B chain0.060.060.010.02
BAG family molecular chaperone regulator 30.010.010.000.01
Desmin0.370.400.140.20
Dystrophin0.080.070.020.02
Filamin-C0.120.150.080.13
Heat shock protein beta-80.010.010.000.00
Myotilin0.380.410.130.25
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Eggers, B.; Schork, K.; Turewicz, M.; Barkovits, K.; Eisenacher, M.; Schröder, R.; Clemen, C.S.; Marcus, K. Advanced Fiber Type-Specific Protein Profiles Derived from Adult Murine Skeletal Muscle. Proteomes 2021, 9, 28. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020028

AMA Style

Eggers B, Schork K, Turewicz M, Barkovits K, Eisenacher M, Schröder R, Clemen CS, Marcus K. Advanced Fiber Type-Specific Protein Profiles Derived from Adult Murine Skeletal Muscle. Proteomes. 2021; 9(2):28. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020028

Chicago/Turabian Style

Eggers, Britta, Karin Schork, Michael Turewicz, Katalin Barkovits, Martin Eisenacher, Rolf Schröder, Christoph S. Clemen, and Katrin Marcus. 2021. "Advanced Fiber Type-Specific Protein Profiles Derived from Adult Murine Skeletal Muscle" Proteomes 9, no. 2: 28. https://0-doi-org.brum.beds.ac.uk/10.3390/proteomes9020028

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