Next Article in Journal
The Secretome Analysis of Activated Human Renal Fibroblasts Revealed Beneficial Effect of the Modulation of the Secreted Peptidyl-Prolyl Cis-Trans Isomerase A in Kidney Fibrosis
Next Article in Special Issue
The Diversity of Muscles and Their Regenerative Potential across Animals
Previous Article in Journal
Engineering CD4+ T Cells to Enhance Cancer Immunity
Previous Article in Special Issue
“The Social Network” and Muscular Dystrophies: The Lesson Learnt about the Niche Environment as a Target for Therapeutic Strategies
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Cells
Volume 9
Issue 7
10.3390/cells9071723
cells-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

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

High-Dimensional Single-Cell Quantitative Profiling of Skeletal Muscle Cell Population Dynamics during Regeneration

by
Lucia Lisa Petrilli
1,2,†,
Filomena Spada
1,†,
Alessandro Palma
1,
Alessio Reggio
1,
Marco Rosina
1,
Cesare Gargioli
1,
Luisa Castagnoli
1,
Claudia Fuoco
1,* and
Gianni Cesareni
1,3
1
Department of Biology, University of Rome “Tor Vergata”, 00133 Rome, Italy
2
Department of Onco-hematology, Gene and Cell Therapy—Bambino Gesù Children’s Hospital—IRCCS, 00146 Rome, Italy
3
Fondazione Santa Lucia Istituto di Ricovero e Cura a Carattere Scientifico (IRCCS), 00143 Rome, Italy
*
Author to whom correspondence should be addressed.
†
These authors contributed equally to the study.
Cells 2020, 9(7), 1723; https://doi.org/10.3390/cells9071723
Submission received: 17 April 2020 / Revised: 3 June 2020 / Accepted: 14 July 2020 / Published: 18 July 2020
(This article belongs to the Special Issue Muscle Homeostasis and Regeneration: From Molecular Mechanisms to Therapeutic Opportunities)
Browse Figures
Versions Notes

Abstract

:
The interstitial space surrounding the skeletal muscle fibers is populated by a variety of mononuclear cell types. Upon acute or chronic insult, these cell populations become activated and initiate finely-orchestrated crosstalk that promotes myofiber repair and regeneration. Mass cytometry is a powerful and highly multiplexed technique for profiling single-cells. Herein, it was used to dissect the dynamics of cell populations in the skeletal muscle in physiological and pathological conditions. Here, we characterized an antibody panel that could be used to identify most of the cell populations in the muscle interstitial space. By exploiting the mass cytometry resolution, we provided a comprehensive picture of the dynamics of the major cell populations that sensed and responded to acute damage in wild type mice and in a mouse model of Duchenne muscular dystrophy. In addition, we revealed the intrinsic heterogeneity of many of these cell populations.

1. Introduction

In physiological conditions, the adult skeletal muscle has a relatively low cell turnover [1]. However, physical activity, trauma, or muscle pathologies, undermining tissue integrity, trigger a tightly controlled regeneration process. Although satellite cells (SCs) are the main actors of myofiber regeneration after damage [2,3,4,5,6], successful muscle healing requires the participation of additional cell types that directly or indirectly contribute to this process. In this context, immune cells and fibro/adipogenic progenitors (FAPs) play a prominent role in supporting the clearance of the damaged tissue, while assisting SCs in their regenerative role [7,8,9,10,11]. However, the orchestrated crosstalk of the regeneration machinery gradually fails in patients affected by muscle-related disorders, such as dystrophies [12]. Here, the accumulation of intrinsic cell defects and the changes in the stem cell niche lead to infiltrations of fat and fibrotic deposition, compromising muscle functions [9,13,14,15,16].
Over the past decades, the complex cell crosstalk occurring during muscle regeneration has been studied in detail [17,18,19]. However, to date, most studies have mainly relied on the analysis of bulk cell populations identified by the expression of a few specific markers and sorted for ex vivo analysis. As a consequence, due to the lack of technologies suitable to address this issue, little is known about muscle cell population heterogeneity. Only recently, the development of technologies to determine the transcriptome of single-cells or their exposed antigen repertoires has permitted to reveal the extent of this heterogeneity and its possible implication in muscle physiology and pathology [20,21,22,23,24]. Moreover, the dynamic changes and the relative abundance of muscle cell populations, upon acute or chronic damage, still remain largely uncharacterized. To fill this gap, we explored, via single-cell mass cytometry, the changes in the multidimensional antigen repertoires of the main players, colonizing the muscle stem cell niche after injury. Here, we described, at single-cell resolution, the time-dependent changes of muscle population dynamics upon myotoxin-induced damage in wild type (wt) and in a mouse model of Duchenne muscular dystrophy (the mdx model).

2. Materials and Methods

2.1. Mouse Strains and Animal Procedures

C57BL/6J (RRID:IMSR_JAX:000664) and C57BL/10ScSn-Dmdmdx/J mice (RRID:IMSR_JAX:001801), hereafter referred to as wt and mdx mice, respectively, were purchased from the Jackson Laboratory.
Mice were bred respecting the standard animal facility procedures, and all the procedures were conducted in accordance with rules of good animal experimentation I.A.C.U.C. n°432 of 12 March 2006 and under ethical approval released on 23/October/2017 from the Italian Ministry of Health, protocol #820/2017-PR.
For muscle injury, 45-day-old wt and mdx mice were anesthetized with an intramuscular injection of saline solution containing ketamine (5 mg/mL) and xylazine (1 mg/mL) prior to the intramuscular administration of 20 μL of 10 μM cardiotoxin solution, isolated from Naja Pallida (Latoxan L8102, Portes les valence, France), into tibialis anterior, quadriceps, and gastrocnemius muscles.

2.2. Histological Analysis

Tibialis anterior (TA) muscles were collected, embedded in optimal cutting temperature compound (Killik—O.C.T., Bio Optica, Milan, Italy), and snap-frozen in liquid nitrogen for 10 s. Embedded muscles were stored at −80 °C for transverse cryo-sectioning with a Leica cryostat. Cryosections (10 μm thickness) were collected on Superfrost glass slides (Thermo Fisher Scientific, Monza, Italy), and tissue slides were stained with hematoxylin and eosin (H&E).
For the H&E, cryosections were fixed with 4% paraformaldehyde (PFA, Santa Cruz Biotechnology, D.B.A. Italia S.r.l., Segrate Milan, Italy) for 15 min at room temperature (RT). After washing in 1X PBS, tissue slides were incubated in the hematoxylin solution for 15 min and rinsed for 5 min in tap water. Cryosections were then counterstained with an alcoholic solution of eosin for 30 min. Following the eosin staining, cryosections were dehydrated in increasing concentrations of alcohol, clarified with the histo-clear solution (Agar Scientific Ltd, Stansted, UK), and finally mounted on coverslips, using the resinous Eukitt mounting medium (Electron Microscopy Sciences, Hatfield Township, PA, USA).
H&E images were captured using the Zeiss Lab A1 AX10 microscope at the 20× magnification in the bright field.

2.3. Skeletal Muscle Mononuclear Cell Purification

Isolation of mononuclear cell populations was performed as in Spada et al. [25]. Mice were sacrificed by cervical dislocation, and the hind limbs were washed with 70% ethanol. Mice hind limbs were then dissected and finely minced in Hank’s balanced salt solution (HBSS) with calcium and magnesium (Gibco- Thermo Fisher Scientific, Monza, Italy) supplemented with 0.2% bovine serum albumin (BSA) (AppliChem, Cinisello Balsamo, Milan, Italy) and 1% penicillin-streptomycin (P/S) (Life Technologies, Monza, Italy, 10,000 U/mL) (HBSS+) under a sterile hood. The homogenized tissue preparation was centrifuged at 70× g for 10 min at 4 °C to separate fat and subjected to enzymatic digestion for 1 h at 37 °C, with gentle mixing in a solution containing 2 μg/μL collagenase A (Roche- Merck KGaA, Darmstadt, Germany), 2.4 U/mL dispase II (Roche- Merck KGaA, Darmstadt, Germany), and 10 μg/mL DNase I (Roche- Merck KGaA, Darmstadt, Germany) diluted in Dulbecco’s phosphate-buffered saline (D-PBS) with calcium and magnesium (Gibco-Thermo Fisher Scientific, Monza, Italy). The reaction was inactivated with HBSS+, and the cell suspension was subjected to three sequential filtrations through 100 μm, 70 μm, and 40 μm cell strainers (BD Falcon, BD Italia, Milan, Italy) and centrifugations at 700× g for 5 min. The lysis of red blood cells was performed by incubating with RBC Lysis Buffer (Santa Cruz Biotechnology, D.B.A. Italia S.r.l., Segrate, Milan, Italy) for 150 s on ice, prior to the 40 µm filtration step.

2.4. Single-Cell Mass Cytometry

For single-cell mass cytometry experiments, 3 × 106 cells were used for each condition. Each time point was analyzed in triplicate, starting from mononuclear cells purified from three different mice. Cells were centrifuged at 600× g for 5 min and washed in D-PBS w/o calcium and magnesium (BioWest- VWR INTERNATIONAL PBI S.r.l., Milan, Italy). To minimize the inter-sample antibody staining variation, we applied a mass-tag barcoding protocol on fixed cells. Cells were fixed with 1 mL of Fix I Buffer (Fluidigm, South San Francisco, CA, USA) and then incubated for 10 min at RT. The fixation was quenched with Barcode Perm Buffer (Fluidigm, South San Francisco, CA, USA). The different samples were barcoded by individually incubating the cell suspensions with the appropriate combination of palladium isotopes from the Cell-IDTM 20-Plex Pd Barcoding Kit (Fluidigm, South San Francisco, CA, USA) in Barcode Perm Buffer for 30 min at RT. The staining was quenched with MaxPar Cell Staining Buffer (Fluidigm, South San Francisco, CA, USA).
The antibody staining with metal-tagged antibodies that target surface and intracellular antigens was performed on the samples pooled after mass-tag barcoding. Samples were collected in a single tube, and the surface antibody staining protocol was performed according to manufacturers’ instructions for 30 min at RT. Surface-stained cells were then washed twice with MaxPar Cell Staining Buffer (Fluidigm South San Francisco, CA, USA) and permeabilized with ice-cold methanol for 10 min on ice. Membrane-permeabilized cells were washed twice with MaxPar Cell Staining Buffer (Fluidigm, South San Francisco, CA, USA) and incubated with antibodies against intracellular antigens for 30 min at RT according to manufacturers’ instructions. The full list of antibodies is detailed in Table 1. All the antibodies listed were purchased from Fluidigm (South San Francisco, CA, USA). After intracellular antibody staining, cells were washed twice with MaxPar Cell Staining Buffer and stained for 1 h at RT with the intercalation solution, composed of Cell-ID Intercalator-Ir (191Ir and 193Ir, Fluidigm South San Francisco, CA, USA) in MaxPar Fix and Perm Buffer (Fluidigm, South San Francisco, CA, USA) at a final concentration of 125 nM. Cells were washed twice with MaxPar Cell Staining Buffer and MaxPar Water.
For mass cytometry analysis, cells were resuspended at the final concentration of 2.5 × 105 cells/mL in MaxPar Water containing 10% of EQTM Four Element Calibration Beads (Fluidigm, South San Francisco, CA, USA) and filtered through a 30-μm filter-cap FACS tube. Samples were kept on ice prior to the acquisition by using the mass cytometry platform CyTOF2 System (Fluidigm, South San Francisco, CA, USA).

2.5. CyTOF Data Analysis

Following data acquisition, channel intensity was normalized using calibration beads [26], and the normalized fcs file was de-barcoded by using the Debarcoder software (Fluidigm, South San Francisco, CA, USA). Data have been pre-processed using the Cytobank software platform [27]. Cells were manually gated from debris on the basis of DNA content monitored by the incorporation of the iridium (Ir) intercalator. Doublets were then excluded according to the event length parameter, and single live cells were finally manually gated by using the cisplatin (Pt) intercalator signal. Manually gated singlet (191Ir+ 193Ir+), viable (195Pt−) events were imported into Cytofkit for further analysis [28]. Cytofkit [28] parameters were set as follows: 14 biomarkers were included for clustering all the detected live cells per sample (“all” merge method); transformation method: cytofAsinh; FlowSOM was used as clustering algorithm with k = 15, tSNE perplexity set to 30; 2000 iterations and seed: 42.
Data were analyzed with Cytofkit shiny app [28] https://github.com/JinmiaoChenLab/cytofkit) and R scripts. Gating for marker-positive cells was performed by setting the mean +/- standard deviation as a threshold, depending on the expression value distribution of the specific marker.

2.6. Statistical Analysis

The experiments were performed at least in biological triplicates, that is, from at least 3 independent mononuclear cell preparations for each experiment. Only for the time point at day 5 in the wild type time series, we had just two biological repeats. Results were presented as mean ± SEM unless otherwise mentioned. Statistical evaluation was done by using One-way or Two-way ANOVA. Comparisons were considered statistically significant at * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001. All statistical analysis was performed using Prism 6 (GraphPad, San Diego, CA, USA).

3. Results

3.1. Histological Profiling of Skeletal Muscle Tissue from wt and mdx Muscles Following Acute Damage

To gain insights into the skeletal muscle repair process, we aimed at describing the dynamics of muscle cell populations following acute damage. To induce muscle injury, we used a well-established protocol based on the injection of the snake (Naja pallida) myotoxin (cardiotoxin, CTX) into the hind limb muscles of wt and mdx dystrophic mice (Figure 1A) [29,30]. Cardiotoxin, by inhibiting protein kinase C (PKC), induced the increase of cytosolic calcium, causing myofiber myolysis that, in turn, triggers regeneration [31,32,33]. First, we monitored, by hematoxylin and eosin staining, the changes in the skeletal muscle architecture at five different time points that were chosen to monitor the key events of the muscle healing process after damage: necrosis, inflammation, regeneration, and remodeling (Figure 1B,C).
Sections of uninjured wt muscles were characterized by polygonal fibers of uniform size containing peripheral nuclei (Figure 1B) [34]. Upon CTX injection, the skeletal muscle underwent degeneration, setting in motion the regeneration process. The degeneration of the muscle architecture was clearly observable at day 1 after injury, while interstitial cells, either resident cell populations [8] or infiltrating immune cells [35], became conspicuous at day 3 after damage. The regeneration process was completed after 20 days, as highlighted by the presence of multinucleated regenerated myofibers.
The mdx skeletal muscle, on the other hand, even in the absence of acute insult, was characterized by infiltrating inflammatory interstitial cells and centrally nucleated myofibers of different sizes, hallmarks of the dystrophic pathology (Figure 1C) [12,36,37]. Following CTX injection, the injured tissue underwent a regeneration process, without being apparently impacted further by extensive necrosis. If anything, the mdx muscle seemed to be more resilient to the myotoxin-induced damage and did not undergo the massive structural damage that was observed early after cardiotoxin injection in the wt muscle. Altogether, the histological analysis showed that the muscles from the two genetic backgrounds responded differently to CTX-induced injury.

3.2. Single-Cell Quantitative Profiling of Skeletal Muscle Populations Following Myotoxin-Induced Injury

Next, we sought to characterize the different mononuclear cell populations in the two mouse models and monitor their abundance changes during the regeneration process. To this end, we resorted to using mass cytometry [8].
In wt muscles, the number of isolated mononuclear cells increased after damage and peaked at day 3 to return to almost baseline levels at day 10 (Figure 1D). In contrast, consistent with the resilience of the mdx muscle observed in the histological analysis, this response was not detected in the injured muscles of mdx mice, where the mononuclear cell number remained constant over the whole regeneration process (Figure 1E).
The mononuclear cell samples from both animal models were separately barcoded and labeled with a panel of 23 metal-tagged antibodies (Table 1) targeting antigens expressed by muscle resident cells and/or by cell populations from the hematopoietic compartment. Mononuclear cells were purified from the uninjured and injured muscles at five time-points and, after barcoding and labeling, analyzed with a CyTOF2 mass cytometer in a single-run experiment [38]. Signals were debarcoded, and live cells were identified using the cisplatin (Pt) intercalator signal. Live/dead cell analysis highlighted that dead wt cells significantly increased at day 3 (Figure S1A), while for mdx muscles, the live/dead cell ratio remained constant (Figure S1B).
Single-cell data were analyzed by applying, as a dimensionality reduction method, the t-distributed stochastic neighbor embedding (t-SNE) algorithm implemented in Cytofkit [28,39].
As the readouts of phospho-antibodies were barely above the background signal, albeit cell-specific (Figure S2), they were not considered in this analysis. The readouts of the 14 antigens in Figure 2A were used as input for the t-SNE algorithm. This approach yielded a two-dimensional map of the antigenic expression profiles of mononuclear cell populations in the wt skeletal muscle (Figure 2A) and led to the identification of 15 different cell clusters (Figure 2B).
The 15 identified clusters (Figure 2B) were further grouped into eight cell types by matching their expression profile to that of cell types already described in the literature [20,40]. More specifically, we were able to identify populations expressing antigens typical of immune cells (CD45+) (clusters 7, 8, 12, 13, 14, and 15), macrophages (CD45+ and F4/80+) (clusters 7, 13, and 14), myogenic progenitors (MPs) (α7-integrin+) (cluster 1), fibro/adipogenic progenitors (FAPs) (SCA1+, CD34+, CD140α+, CD90.2+, and vimentin+) (cluster 10), endothelial progenitor cells (CD31+) (clusters 3 and 4), pericyte-like cells (CD146+ and CD140β+) (cluster 11), and mesenchymal-like cells (CD90.2+) (cluster 5). The expression profiles of three remaining non-abundant clusters (clusters 2, 6, and 9) could not be matched to any of the already described muscle cell types and were collectively dubbed “others” (Figure 2B,C).
As we analyzed samples at different time points, we could also monitor the dynamic of cell populations after CTX damage. The picture in Figure 2C, in fact, was not static as the relative abundance of the different cell populations changed along the regeneration process. These could be best appreciated by comparing the bidimensional t-SNE maps at different time points (Figure 2D). Here the colors, relating to cell density, from blue (low density) to red (high density), allowed to appreciate the significant changes in the cell population proportions during regeneration.
To quantitatively describe the dynamics of these changes, we first classified the observed cell populations into two main clusters according to CD45 expression: (i) the immune (CD45+) and the (ii) non-immune (CD45−) subpopulations of mononuclear cells. The data were shown as population relative abundance in the mononuclear cell preparations analyzed in the CyTOF. However, this data representation could be easily transformed into changes in the absolute numbers of each cell population in the mouse muscle as the total number of mononuclear cells in each condition was known (Figures S3 and S4). The curve trends in the two representations were similar.
We first focused on CD45+ hematopoietic cells as they play a critical role in the regeneration process by sending regulatory signals by removing the damaged-fiber debris and by stimulating proliferation and differentiation of myogenic progenitors [19,41,42,43,44,45]. In a homeostatic muscle, approximately 60% of the mononuclear cells exposed the CD45 antigen. This cell compartment increased in number, after damage, to peak at days 1 and 3 (Figure 3A).
In physiological conditions, a fraction of the hematopoietic cells, less than 2% of the total recorded events, also expressed the F4/80+ antigen, a pan-macrophage marker (Figure 3B). However, already 3 days after injury, the macrophage population significantly increased, reaching a maximum of approximately 40%. At later times, macrophage abundance gradually decreased, returning close to baseline levels at day 20. The macrophage population is not homogeneous as it contains pro-inflammatory (M1) and anti-inflammatory (M2) macrophages that differ in the expression of the CD206 antigen and play a different role in muscle regeneration [46,47].
Noteworthy, CD206+ M2 macrophages, which are responsible for the resolution of the inflammatory response [48,49,50], increased significantly on day 5 at the expense of inflammatory macrophages M1 that peaked at day 3 (Figure 3B,C).
We next focused on the cell populations that did not express the CD45 antigen. This group included the two main players of the regeneration process: myogenic progenitors (MPs) and fibro/adipogenic progenitors (FAPs) [5,6,8,51,52]. Our single-cell analysis revealed that both cell types more than doubled in number from day 3 to 10 (Figure 3D,E).
We could not obtain mass cytometry grade antibodies, specifically recognizing the paired box protein PAX7 antigen, which labels satellite cells. Thus, we resorted to using the α7-integrin antigen as a marker of the myogenic progenitor (MP) cluster (Figure 3D), including both satellite cells and myoblasts. Within the population expressing α7-integrin, we could identify two smaller clusters expressing different levels of vimentin (Figure S1C). During regeneration, the MP cell population, after an initial decrease, became more populated (Figure 3D). The relative abundance of the two subpopulations changed in time with the vimentin expressing MP, significantly increasing in number at day 3 (Figure S1C). At day 3 post-injury, the MPs accounted for about 6% of the total recorded cell events.
The FAP compartment was defined in our t-SNE map by a cluster of cells expressing SCA1, CD34, CD140α, CD90.2, and vimentin. Similar to the MP population, albeit with a different trend, the FAP population became more numerous as the regeneration process progressed, reaching the maximum expansion between 5–10 days after CTX injury (about 11.3% of total events), to return to almost control levels at day 20 (Figure 3E).
Our analysis also allowed us to characterize the kinetic of vessel-associated populations during regeneration [53,54]. These were identified, among CD45 negative cells, as they expressed the CD31 antigen. We were able to distinguish two subpopulations of cells expressing additional markers of endothelial populations at different levels (SCA1, CD146, and CD34) (Figure S1D). CD31 cluster 4 was considered a myoendothelial cell subpopulation, containing cells that also express high levels of the α7-integrin antigen [55]. Overall, we observed that the whole endothelial progenitor pool followed a kinetic that was different from the populations examined so far. They decreased significantly in number in the first few days after damage and then increased over the homeostatic level toward the end of the regeneration process (Figure 3F).
We also considered a population of cells whose expression profile was reminiscent of that of pericytes [56,57]. They were characterized by the expression of CD146, CD140β, and α7-integrin (Figure 3G). These pericyte-like cells followed a trend that was very similar to that of endothelial cells, with a sharp decrease early after damage and a rapid increase in the late regeneration phase, when vascularization took place (Figure 3G). We also looked at a cluster of cells that we were not able to match to any of the cell types described so far. We dubbed this population as mesenchymal-like cells, as they were highly positive for CD90.2 and CD140β markers. The abundance of this population followed a kinetic that seemed to be governed by the regeneration process and was also very similar to that of the classical endothelial cells (Figure 3H).

3.3. Single-Cell Profiling of the mdx Muscle

We further aimed at characterizing the response to acute damage and the ensuing regeneration process in the mdx dystrophic muscle. Mononuclear cells were isolated from the uninjured and CTX-injured skeletal muscles of mdx mice, following the same procedure described for the wt in the previous section (Figure 4A–D). Mononuclear cells were separated from muscle fibers, barcoded, labeled with the same antibody panel, and prepared for mass cytometry. The resulting single-cell antigen expression profile was processed, as described for the wt cells, in order to obtain a two-dimensional representation of antigen expression in the different cell types. As the wt and mdx mass cytometry analyses were performed at different times, the resulting t-SNE maps could not be directly compared (Figure 2A and Figure 4A). However, by comparing the antigen expression profiles in the different clusters, we were able to associate each cluster to one of the main muscle mononuclear cell types and match it to those observed in the wt (Figure 4B,C). As observed in injured-recovering wt muscles, the relative numerosity of the different cell populations also changed after damage and during regeneration in the mdx model (Figure 4D).
The number of cells in the CD45+ compartment (clusters 9, 10, 11, 12, 13, 14, and 15), as a whole, did not change significantly along the regeneration process (Figure 5A). However, the distribution of cells in the different populations in the compartment was found to be highly dynamic. In particular, macrophages (clusters 11, 14, and 15), which were already abundant in the mdx muscle before acute damage, increased significantly in the early days after CTX treatment, as observed in the wt model, to return to mdx baseline levels, approximately 30% of total mononuclear cells, at day 20 (Figure 5B). M2 macrophages, on the other hand, remained rather constant early after an injury to increase only late in the regeneration process (Figure 5C).
When compared to wt, myogenic progenitors (MPs) (cluster 4) and fibro/adipogenic progenitors (FAPs) (cluster 2) followed a different trend. As the mdx muscle was under chronic stress, we observed that both progenitor cell populations were more populated in the absence of acute damage. However, while MPs decreased in number as early as one day after acute damage to return to unperturbed levels at later times, FAPs remained at a constant high level along the whole regeneration process (Figure 5D,E).
As already observed for the wt, also for the mdx mouse model, the CD31 expressing cells (clusters 1 and 6), defining endothelial progenitors, comprised two sub-populations differing for the expression of α7-integrin and other markers (CD90.2, CD140β, and CXCR4) (Figure S1E). The population expressing higher levels of α7-integrin was twice as abundant as the other. Nevertheless, both sub-populations reacted similarly to acute damage and first dropped in abundance by a factor of approximately three, to recover at the late stages of the regeneration process (Figure 5F). However, consistently, the cluster enriched in α7-integrin expressing cells (cluster 6) had not fully recovered at day 20, suggesting that the vascularization was still ongoing at times when histology seemed to indicate completion of the regeneration process. This differed from what was observed in the wt. Consistent with this consideration, and differently from what observed in the wt, also the pericyte-like population (cluster 5) behaved similarly, dropping in abundance immediately after damage and then slowly recovering without, however, reaching full recovery at day 20 (Figure 5G).

4. Discussion

The skeletal muscle has a remarkable capacity to self-repair if damaged [7,8,9,10,49,58]. However, this healing process may fail, owing to excessive damage, aging, or genetic disorders [12,59,60,61,62]. As a consequence of this failure in the repair process, as in muscle dystrophies, the tissue undergoes degeneration, leading to progressive muscle wasting and weakness characterized by chronic inflammation and, at later stages, fat and fibrotic tissue infiltrations [9,13,14].
Here, we exploited the resolution power of mass cytometry to characterize the modulation of the profile of the mononuclear cell population following chronic or acute damage [63,64,65]. To this end, we assembled a panel of 23 metal-tagged antibodies and characterized, at the single-cell level, the dynamics of muscle mononuclear cell populations after acute damage in wild type (wt) and in a mouse model of Duchenne muscular dystrophy (mdx). The regeneration process was monitored by examining samples of muscle mononuclear cells at 1, 3, 5, 10, and 20 days after cardiotoxin injection in wt or mdx mice [29,30].
This approach yielded a reach multiparametric dataset, disclosing the details of how the composition and heterogeneity of mononuclear cell populations changed in time as the muscle healing process proceeded. By applying a dimensionality reduction technique, such as the t-distributed stochastic neighbor embedding (t-SNE) algorithm, we generated bidimensional maps, providing a visual description of the regeneration process. Furthermore, this representation contributed to revealing subtle differences in the expression of specific markers in subpopulations within the major clusters that identified “classical” muscle populations. This population heterogeneity could not be only explained by experimental variability, even if any functional implication in muscle physiology or pathology remains to be established. In this report, we only dwelt upon a few of these t-SNE map features, while the dataset remains as a resource for additional analysis for the community.
A time-dependent variation in the abundance of muscle mononuclear cells was observed in wt mice, starting at day 1, indicating that the system sensed the damage and promptly responded to restore muscle tissue homeostasis. Inflammatory cells were the predominant population 3 days after the injury, accounting for over 40% of the total mononuclear cells. The activation of the inflammatory compartment was limited to the first few days after the injury as macrophages returned to baseline levels after 10 days.
Fibro/adipogenic progenitors (FAPs), here identified as SCA1+, CD34+, CD140α+, CD90.2+, and vimentin+ [66], stimulate satellite cell activation and differentiation, thus playing a positive role in muscle regeneration [8]. FAPs are quiescent in intact muscles, while they proliferate in response to injury [8,67,68]. Consistently, we observed that FAPs, after CTX injection, rapidly expanded, in wt muscles, reaching a peak between days 5 and 10, to return eventually to baseline level, having accomplished their function of providing a transient production of pro-differentiation signals and of depositing extracellular matrix for muscle remodeling.
In response to the secreted inflammatory and fibro/adipogenic stimulating signals, the myogenic compartment also promptly activated at day 3, as confirmed by vimentin expression, which marked activated MPs or myoblasts [69]. After an initial drop in concentration, α7-integrin+ MPs, including satellite cells and myoblasts, expanded on day 3, still remaining high on day 10, when they started to decrease.
Different kinetics was observed for endothelial progenitor cells, here, characterized as CD31, SCA1, CD34, and CD146 expressing cells, and for pericyte-like cells, mainly identified as CD146 and CD140β-positive cells. The endothelial and pericyte-like cell contribution to skeletal muscle recovery only took place at the end of the regeneration process; once the inflammatory cells were removed, FAPs decreased in number, and MPs differentiated.
Due to the incomplete nature of our panel, some cell clusters in our t-SNE maps remained loosely defined. One cluster included cells with an antigen repertoire reminiscent of mesenchymal cells, as they were positive for the mesenchymal markers CD90.2 and CD140β [70]. This population had a clear kinetic, suggesting a late intervention of the mesenchymal population in the muscle regeneration process.
We also investigated the response to acute injury of a muscle environment that was already chronically perturbed as in the mdx mice. Mononuclear cells extracted from the muscle of mdx mice did not experience any significant increase in number following cardiotoxin injury. However, by looking into the details of the population profiles at each time point, we observed a significant modulation of the population distributions as the increase in the number of the cells in one population was counterbalanced by the decrease in another one. For instance, while we observed that the cells in the inflammatory compartment, essentially macrophages, expanded early after damage, the endothelial and pericyte-like clusters dropped in numerosity to recover the initial values only late in the process. On the other hand, differently from wt, the myogenic and fibro/adipogenic compartments, the two main players in muscle regeneration, showed little variation during the regeneration process, probably because they were already chronically activated.
Overall, our multiparametric analysis offered a comprehensive description of both muscle tissue homeostasis and the rearrangements induced in the mononuclear cell population profile by a perturbation of the muscle system, be it a chronic condition, as in the case of mdx mice, or acute stress, as that triggered by cardiotoxin injection.

Supplementary Materials

The following are available online at https://www.mdpi.com/2073-4409/9/7/1723/s1. Figures S1–S4 and “FCS folder” with two subfolders (“MDX_FCS” and “WT_FCS”) containing the mass cytometry raw data (.fcs files).

Author Contributions

Conceptualization, L.L.P., F.S., C.G., L.C., C.F., and G.C.; methodology, L.L.P., F.S., A.P., A.R., M.R., C.G., and C.F.; software, L.L.P., F.S., and A.P.; validation, L.L.P., F.S., A.P., A.R., M.R., C.G., and C.F.; formal analysis, L.L.P., F.S., A.P., A.R., and M.R.; investigation, L.L.P., F.S., A.P., A.R., and M.R; resources, C.G., L.C., C.F., and G.C.; data curation, L.L.P., F.S., and A.P.; writing—original draft preparation, L.L.P. and F.S.; writing—review and editing, L.L.P., F.S., C.G., L.C., C.F., and G.C.; visualization, L.L.P., F.S., C.G., L.C., C.F., and G.C.; supervision, C.G, L.C., C.F., and G.C.; project administration, L.C., C.F., and G.C.; funding acquisition, G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by a grant of the European Research Council Grant N. 322749 (G.C.).

Acknowledgments

We acknowledge the Umberto Veronesi Foundation for awarding C. Fuoco with post-doctoral fellowship 2019.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Decary, S.; Mouly, V.; Hamida, C.B.; Sautet, A.; Barbet, J.P.; Butler-Browne, G.S. Replicative potential and telomere length in human skeletal muscle: Implications for satellite cell-mediated gene therapy. Hum. Gene Ther. 1997, 8, 1429–1438. [Google Scholar] [CrossRef] [PubMed]
  2. Mauro, A. Satellite cell of skeletal muscle fibers. J. Biophys. Biochem. Cytol. 1961, 9, 493–495. [Google Scholar] [CrossRef] [PubMed]
  3. Yin, H.; Price, F.; Rudnicki, M.A. Satellite cells and the muscle stem cell niche. Physiol. Rev. 2013, 93, 23–67. [Google Scholar] [CrossRef] [Green Version]
  4. Lepper, C.; Partridge, T.A.; Fan, C.-M. An absolute requirement for Pax7-positive satellite cells in acute injury-induced skeletal muscle regeneration. Development 2011, 138, 3639–3646. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Sambasivan, R.; Yao, R.; Kissenpfennig, A.; Van Wittenberghe, L.; Paldi, A.; Gayraud-Morel, B.; Guenou, H.; Malissen, B.; Tajbakhsh, S.; Galy, A. Pax7-expressing satellite cells are indispensable for adult skeletal muscle regeneration. Development 2011, 138, 3647–3656. [Google Scholar] [CrossRef] [Green Version]
  6. Seale, P.; Sabourin, L.A.; Girgis-Gabardo, A.; Mansouri, A.; Gruss, P.; Rudnicki, M.A. Pax7 is required for the specification of myogenic satellite cells. Cell 2000, 102, 777–786. [Google Scholar] [CrossRef] [Green Version]
  7. Heredia, J.E.; Mukundan, L.; Chen, F.M.; Mueller, A.A.; Deo, R.C.; Locksley, R.M.; Rando, T.A.; Chawla, A. Type 2 innate signals stimulate fibro/adipogenic progenitors to facilitate muscle regeneration. Cell 2013, 153, 376–388. [Google Scholar] [CrossRef] [Green Version]
  8. Joe, A.W.B.; Yi, L.; Natarajan, A.; Le Grand, F.; So, L.; Wang, J.; Rudnicki, M.A.; Rossi, F.M.V. Muscle injury activates resident fibro/adipogenic progenitors that facilitate myogenesis. Nat. Cell Biol. 2010, 12, 153–163. [Google Scholar] [CrossRef] [Green Version]
  9. Uezumi, A.; Ito, T.; Morikawa, D.; Shimizu, N.; Yoneda, T.; Segawa, M.; Yamaguchi, M.; Ogawa, R.; Matev, M.M.; Miyagoe-Suzuki, Y.; et al. Fibrosis and adipogenesis originate from a common mesenchymal progenitor in skeletal muscle. J. Cell Sci. 2011, 124, 3654–3664. [Google Scholar] [CrossRef] [Green Version]
  10. Murphy, M.M.; Lawson, J.A.; Mathew, S.J.; Hutcheson, D.A.; Kardon, G. Satellite cells, connective tissue fibroblasts and their interactions are crucial for muscle regeneration. Development 2011, 138, 3625–3637. [Google Scholar] [CrossRef] [Green Version]
  11. Malecova, B.; Gatto, S.; Etxaniz, U.; Passafaro, M.; Cortez, A.; Nicoletti, C.; Giordani, L.; Torcinaro, A.; Bardi, M.; Bicciato, S.; et al. Dynamics of cellular states of fibro-adipogenic progenitors during myogenesis and muscular dystrophy. Nat. Commun. 2018, 9, 1–12. [Google Scholar] [CrossRef] [PubMed]
  12. Emery, A.E.H. The muscular dystrophies. Lancet 2002, 359, 687–695. [Google Scholar] [CrossRef]
  13. 1Marden, F.A.; Connolly, A.M.; Siegel, M.J.; Rubin, D.A. Compositional analysis of muscle in boys with Duchenne muscular dystrophy using MR imaging. Skelet. Radiol. 2004, 34, 140–148. [Google Scholar] [CrossRef]
  14. Contreras, O.; Rebolledo, D.L.; Oyarzún, J.E.; Olguin, H.C.; Brandan, E. Connective tissue cells expressing fibro/adipogenic progenitor markers increase under chronic damage: Relevance in fibroblast-myofibroblast differentiation and skeletal muscle fibrosis. Cell Tissue Res. 2016, 364, 647–660. [Google Scholar] [CrossRef] [PubMed]
  15. Hogarth, M.W.; Defour, A.; Lazarski, C.; Gallardo, E.; Diaz-Manera, J.; Partridge, T.A.; Nagaraju, K.; Jaiswal, J.K. Fibroadipogenic progenitors are responsible for muscle loss in limb girdle muscular dystrophy 2B. Nat. Commun. 2019, 10, 1–13. [Google Scholar] [CrossRef] [Green Version]
  16. Madaro, L.; Passafaro, M.; Sala, D.; Etxaniz, U.; Lugarini, F.; Proietti, D.; Alfonsi, M.V.; Nicoletti, C.; Gatto, S.; De Bardi, M.; et al. Denervation-activated STAT3-IL-6 signalling in fibro-adipogenic progenitors promotes myofibres atrophy and fibrosis. Nat. Cell Biol. 2018, 20, 917–927. [Google Scholar] [CrossRef] [PubMed]
  17. Sirabella, D.; De Angelis, L.; Berghella, L. Sources for skeletal muscle repair: From satellite cells to reprogramming. J. Cachexia Sarcopenia Muscle 2013, 4, 125–136. [Google Scholar] [CrossRef]
  18. Dey, D.; Goldhamer, D.J.; Yu, P.B. Contributions of muscle-resident progenitor cells to homeostasis and disease. Curr. Mol. Biol. Rep. 2015, 1, 175–188. [Google Scholar] [CrossRef] [Green Version]
  19. Pillon, N.J.; Bilan, P.J.; Fink, L.N.; Klip, A. Cross-talk between skeletal muscle and immune cells: Muscle-derived mediators and metabolic implications. Am. J. Physiol. Endocrinol. Metab. 2013, 304, E453–E465. [Google Scholar] [CrossRef]
  20. Giordani, L.; He, G.J.; Negroni, E.; Sakai, H.; Law, J.Y.C.; Siu, M.M.; Wan, R.; Corneau, A.; Tajbakhsh, S.; Cheung, T.H.; et al. High-dimensional single-cell cartography reveals novel skeletal muscle-resident cell populations. Mol. Cell 2019, 74, 609–621.e6. [Google Scholar] [CrossRef]
  21. Gatto, S.; Puri, P.L.; Malecova, B. Single cell gene expression profiling of skeletal muscle-derived cells. Methods Mol. Biol. 2017, 1556, 191–219. [Google Scholar] [PubMed]
  22. Rubenstein, A.B.; Smith, G.R.; Raue, U.; Begue, G.; Minchev, K.; Ruf-Zamojski, F.; Nair, V.D.; Wang, X.; Zhou, L.; Zaslavsky, E.; et al. Single-cell transcriptional profiles in human skeletal muscle. Nat. Publ. Group 2020, 10, 229. [Google Scholar] [CrossRef] [Green Version]
  23. Dell’Orso, S.; Juan, A.H.; Ko, K.-D.; Naz, F.; Perovanovic, J.; Gutierrez-Cruz, G.; Feng, X.; Sartorelli, V. Single cell analysis of adult mouse skeletal muscle stem cells in homeostatic and regenerative conditions. Development 2019, 146, dev174177. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Marinkovic, M.; Fuoco, C.; Sacco, F.; Cerquone Perpetuini, A.; Giuliani, G.; Micarelli, E.; Pavlidou, T.; Petrilli, L.L.; Reggio, A.; Riccio, F.; et al. Fibro-adipogenic progenitors of dystrophic mice are insensitive to NOTCH regulation of adipogenesis. Life Sci. Alliance 2019, 2, e201900437. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Spada, F.; Fuoco, C.; Pirrò, S.; Paoluzi, S.; Castagnoli, L.; Gargioli, C.; Cesareni, G. Characterization by mass cytometry of different methods for the preparation of muscle mononuclear cells. New Biotechnol. 2016, 33, 514–523. [Google Scholar] [CrossRef] [Green Version]
  26. Finck, R.; Simonds, E.F.; Jager, A.; Krishnaswamy, S.; Sachs, K.; Fantl, W.; Pe’er, D.; Nolan, G.P.; Bendall, S.C. Normalization of mass cytometry data with bead standards. Cytom. Part 2013, 83, 483–494. [Google Scholar] [CrossRef]
  27. Kotecha, N.; Krutzik, P.O.; Irish, J.M. Web-based analysis and publication of flow cytometry experiments. Curr. Protoc. Cytom. 2010. [Google Scholar] [CrossRef]
  28. Chen, H.; Lau, M.C.; Wong, M.T.; Newell, E.W.; Chen, J. Cytofkit: A bioconductor package for an integrated mass cytometry data analysis pipeline. PLoS Comput. Biol. 2016, 12, e1005112. [Google Scholar] [CrossRef]
  29. Ramadasan-Nair, R.; Gayathri, N.; Mishra, S.; Sunitha, B.; Mythri, R.B.; Nalini, A.; Subbannayya, Y.; Harsha, H.C.; Kolthur-Seetharam, U.; Srinivas Bharath, M.M. Mitochondrial alterations and oxidative stress in an acute transient mouse model of muscle degeneration: Implications for muscular dystrophy and related muscle pathologies. J. Biol. Chem. 2014, 289, 485–509. [Google Scholar] [CrossRef] [Green Version]
  30. Duchen, L.W.; Excell, B.J.; Patel, R.; Smith, B. Changes in motor end-plates resulting from muscle fibre necrosis and regeneration. A light and electron microscopic study of the effects of the depolarizing fraction (cardiotoxin) of Dendroaspis jamesoni venom. J. Neurol. Sci. 1974, 21, 391–417. [Google Scholar] [CrossRef]
  31. Wang, H.-X.; Lau, S.-Y.; Huang, S.-J.; Kwan, C.-Y.; Wong, T.-M. Cobra venom cardiotoxin induces perturbations of cytosolic calcium homeostasis and hypercontracture in adult rat ventricular myocytes. J. Mol. Cell. Cardiol. 1997, 29, 2759–2770. [Google Scholar] [CrossRef] [PubMed]
  32. Raynor, R.L.; Zheng, B.; Kuo, J.F. Membrane interactions of amphiphilic polypeptides mastoparan, melittin, polymyxin B, and cardiotoxin. Differential inhibition of protein kinase C, Ca2+/calmodulin-dependent protein kinase II and synaptosomal membrane Na,K-ATPase, and Na+ pump and differentiation of HL60 cells. J. Biol. Chem. 1991, 266, 2753–2758. [Google Scholar] [PubMed]
  33. Suh, B.C.; Song, S.K.; Kim, Y.K.; Kim, K.T. Induction of cytosolic Ca2+ elevation mediated by Mas-7 occurs through membrane pore formation. J. Biol. Chem. 1996, 271, 32753–32759. [Google Scholar] [CrossRef] [Green Version]
  34. Bentzinger, C.F.; Wang, Y.X.; Dumont, N.A.; Rudnicki, M.A. Cellular dynamics in the muscle satellite cell niche. EMBO Rep. 2013, 14, 1062–1072. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. St Pierre, B.A.; Tidball, J.G. Differential response of macrophage subpopulations to soleus muscle reloading after rat hindlimb suspension. J. Appl. Physiol. 1994, 77, 290–297. [Google Scholar] [CrossRef]
  36. Sacco, A.; Mourkioti, F.; Tran, R.; Choi, J.; Llewellyn, M.; Kraft, P.; Shkreli, M.; Delp, S.; Pomerantz, J.H.; Artandi, S.E.; et al. Short telomeres and stem cell exhaustion model Duchenne muscular dystrophy in mdx/mTR mice. Cell 2010, 143, 1059–1071. [Google Scholar] [CrossRef] [Green Version]
  37. Anderson, J.E.; Ovalle, W.K.; Bressler, B.H. Electron microscopic and autoradiographic characterization of hindlimb muscle regeneration in the mdx mouse. Anat. Rec. 1987, 219, 243–257. [Google Scholar] [CrossRef]
  38. Zunder, E.R.; Finck, R.; Behbehani, G.K.; Amir, E.-A.D.; Krishnaswamy, S.; Gonzalez, V.D.; Lorang, C.G.; Bjornson, Z.; Spitzer, M.H.; Bodenmiller, B.; et al. Palladium-based mass tag cell barcoding with a doublet-filtering scheme and single-cell deconvolution algorithm. Nat. Protoc. 2015, 10, 316–333. [Google Scholar] [CrossRef]
  39. Amir, E.-A.D.; Davis, K.L.; Tadmor, M.D.; Simonds, E.F.; Levine, J.H.; Bendall, S.C.; Shenfeld, D.K.; Krishnaswamy, S.; Nolan, G.P.; Pe’er, D. VISNE enables visualization of high dimensional single-cell data and reveals phenotypic heterogeneity of leukemia. Nat. Biotechnol. 2013, 31, 545–552. [Google Scholar] [CrossRef] [Green Version]
  40. Palma, A.; Cerquone Perpetuini, A.; Ferrentino, F.; Fuoco, C.; Gargioli, C.; Giuliani, G.; Iannuccelli, M.; Licata, L.; Micarelli, E.; Paoluzi, S.; et al. Myo-REG: A portal for signaling interactions in muscle regeneration. Front. Physiol. 2019, 10, 1216. [Google Scholar] [CrossRef] [Green Version]
  41. Arnold, L.; Henry, A.; Poron, F.; Baba-Amer, Y.; Van Rooijen, N.; Plonquet, A.; Gherardi, R.K.; Chazaud, B. Inflammatory monocytes recruited after skeletal muscle injury switch into antiinflammatory macrophages to support myogenesis. J. Exp. Med. 2007, 204, 1057–1069. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. McLennan, I.S. Degenerating and regenerating skeletal muscles contain several subpopulations of macrophages with distinct spatial and temporal distributions. J. Anat. 1996, 188, 17–28. [Google Scholar]
  43. Lolmede, K.; Campana, L.; Vezzoli, M.; Bosurgi, L.; Tonlorenzi, R.; Clementi, E.; Bianchi, M.E.; Cossu, G.; Manfredi, A.A.; Brunelli, S.; et al. Inflammatory and alternatively activated human macrophages attract vessel-associated stem cells, relying on separate HMGB1- and MMP-9-dependent pathways. J. Leukoc. Biol. 2009, 85, 779–787. [Google Scholar] [CrossRef]
  44. Chazaud, B.; Sonnet, C.; Lafuste, P.; Bassez, G.; Rimaniol, A.-C.; Poron, F.; Authier, F.-J.; Dreyfus, P.A.; Gherardi, R.K. Satellite cells attract monocytes and use macrophages as a support to escape apoptosis and enhance muscle growth. J. Cell Biol. 2003, 163, 1133–1143. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Dumont, N.; Frenette, J. Macrophages protect against muscle atrophy and promote muscle recovery in vivo and in vitro: A mechanism partly dependent on the insulin-like growth factor-1 signaling molecule. Am. J. Pathol. 2010, 176, 2228–2235. [Google Scholar] [CrossRef] [PubMed]
  46. Martinez, F.O.; Gordon, S. The M1 and M2 paradigm of macrophage activation: Time for reassessment. F1000Prime Rep. 2014, 6, 13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Kharraz, Y.; Guerra, J.; Mann, C.J.; Serrano, A.L.; Muñoz-Cánoves, P. Macrophage plasticity and the role of inflammation in skeletal muscle repair. Mediat. Inflamm. 2013. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Mantovani, A.; Biswas, S.K.; Galdiero, M.R.; Sica, A.; Locati, M. Macrophage plasticity and polarization in tissue repair and remodelling. J. Pathol. 2013, 229, 176–185. [Google Scholar] [CrossRef] [PubMed]
  49. Saclier, M.; Yacoub-Youssef, H.; Mackey, A.L.; Arnold, L.; Ardjoune, H.; Magnan, M.; Sailhan, F.; Chelly, J.; Pavlath, G.K.; Mounier, R.; et al. Differentially activated macrophages orchestrate myogenic precursor cell fate during human skeletal muscle regeneration. Stem Cells 2013, 31, 384–396. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  50. Stout, R.D.; Jiang, C.; Matta, B.; Tietzel, I.; Watkins, S.K.; Suttles, J. Macrophages sequentially change their functional phenotype in response to changes in microenvironmental influences. J. Immunol. 2005, 175, 342–349. [Google Scholar] [CrossRef]
  51. Fu, X.; Wang, H.; Hu, P. Stem cell activation in skeletal muscle regeneration. Cell. Mol. Life Sci. 2015, 72, 1663–1677. [Google Scholar] [CrossRef] [Green Version]
  52. Nguyen, H.X.; Tidball, J.G. Interactions between neutrophils and macrophages promote macrophage killing of rat muscle cells in vitro. J. Physiol. (Lond.) 2003, 547, 125–132. [Google Scholar] [CrossRef] [PubMed]
  53. Musarò, A. The basis of muscle regeneration. Adv. Biol. 2014, 2014, 1–16. [Google Scholar] [CrossRef] [Green Version]
  54. Hansen-Smith, F.M.; Hudlicka, O.; Egginton, S. In vivo angiogenesis in adult rat skeletal muscle: Early changes in capillary network architecture and ultrastructure. Cell Tissue Res. 1996, 286, 123–136. [Google Scholar] [CrossRef] [PubMed]
  55. Zheng, B.; Cao, B.; Crisan, M.; Sun, B.; Li, G.; Logar, A.; Yap, S.; Pollett, J.B.; Drowley, L.; Cassino, T.; et al. Prospective identification of myogenic endothelial cells in human skeletal muscle. Nat. Biotechnol. 2007, 25, 1025–1034. [Google Scholar] [CrossRef] [PubMed]
  56. Dellavalle, A.; Sampaolesi, M.; Tonlorenzi, R.; Tagliafico, E.; Sacchetti, B.; Perani, L.; Innocenzi, A.; Galvez, B.G.; Messina, G.; Morosetti, R.; et al. Pericytes of human skeletal muscle are myogenic precursors distinct from satellite cells. Nat. Cell Biol. 2007, 9, 255–267. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Birbrair, A.; Zhang, T.; Wang, Z.-M.; Messi, M.L.; Enikolopov, G.N.; Mintz, A.; Delbono, O. Role of pericytes in skeletal muscle regeneration and fat accumulation. Stem Cells Dev. 2013, 22, 2298–2314. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Christov, C.; Chrétien, F.; Abou-Khalil, R.; Bassez, G.; Vallet, G.; Authier, F.-J.; Bassaglia, Y.; Shinin, V.; Tajbakhsh, S.; Chazaud, B.; et al. Muscle satellite cells and endothelial cells: Close neighbors and privileged partners. Mol. Biol. Cell 2007, 18, 1397–1409. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Almada, A.E.; Wagers, A.J. Molecular circuitry of stem cell fate in skeletal muscle regeneration, ageing and disease. Nat. Rev. Mol. Cell Biol. 2016, 17, 267–279. [Google Scholar] [CrossRef] [Green Version]
  60. Blau, H.M.; Cosgrove, B.D.; Ho, A.T.V. The central role of muscle stem cells in regenerative failure with aging. Nat. Med. 2015, 21, 854–862. [Google Scholar] [CrossRef] [Green Version]
  61. Kuswanto, W.; Burzyn, D.; Panduro, M.; Wang, K.K.; Jang, Y.C.; Wagers, A.J.; Benoist, C.; Mathis, D. Poor repair of skeletal muscle in aging mice reflects a defect in local, interleukin-33-dependent accumulation of regulatory T cells. Immunity 2016, 44, 355–367. [Google Scholar] [CrossRef] [Green Version]
  62. Pastoret, C.; Sebille, A. Age-related differences in regeneration of dystrophic (mdx) and normal muscle in the mouse. Muscle Nerve 1995, 18, 1147–1154. [Google Scholar] [CrossRef] [PubMed]
  63. Ornatsky, O.; Bandura, D.; Baranov, V.; Nitz, M.; Winnik, M.A.; Tanner, S. Highly multiparametric analysis by mass cytometry. J. Immunol. Methods 2010, 361, 1–20. [Google Scholar] [CrossRef] [PubMed]
  64. Bendall, S.C.; Nolan, G.P.; Roederer, M.; Chattopadhyay, P.K. A deep profiler’s guide to cytometry. Trends Immunol. 2012, 33, 323–332. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  65. Bandura, D.R.; Baranov, V.I.; Ornatsky, O.I.; Antonov, A.; Kinach, R.; Lou, X.; Pavlov, S.; Vorobiev, S.; Dick, J.E.; Tanner, S.D. Mass cytometry: Technique for real time single cell multitarget immunoassay based on inductively coupled plasma time-of-flight mass spectrometry. Anal. Chem. 2009, 81, 6813–6822. [Google Scholar] [CrossRef] [PubMed]
  66. Reggio, A.; Rosina, M.; Palma, A.; Cerquone Perpetuini, A.; Petrilli, L.L.; Gargioli, C.; Fuoco, C.; Micarelli, E.; Giuliani, G.; Cerretani, M.; et al. Adipogenesis of skeletal muscle fibro/adipogenic progenitors is affected by the WNT5a/GSK3/β-catenin axis. Cell Death Differ. 2020, 108, 1–21. [Google Scholar] [CrossRef] [PubMed]
  67. Uezumi, A.; Fukada, S.-I.; Yamamoto, N.; Takeda, S.; Tsuchida, K. Mesenchymal progenitors distinct from satellite cells contribute to ectopic fat cell formation in skeletal muscle. Nat. Cell Biol. 2010, 12, 143–152. [Google Scholar] [CrossRef]
  68. Lemos, D.R.; Babaeijandaghi, F.; Low, M.; Chang, C.-K.; Lee, S.T.; Fiore, D.; Zhang, R.-H.; Natarajan, A.; Nedospasov, S.A.; Rossi, F.M.V. Nilotinib reduces muscle fibrosis in chronic muscle injury by promoting TNF-mediated apoptosis of fibro/ adipogenic progenitors. Nat. Med. 2015, 21, 1–11. [Google Scholar] [CrossRef]
  69. Vater, R.; Cullen, M.J.; Harris, J.B. The expression of vimentin in satellite cells of regenerating skeletal muscle in vivo. Histochem. J. 1994, 26, 916–928. [Google Scholar] [CrossRef]
  70. Chirieleison, S.M.; Feduska, J.M.; Schugar, R.C.; Askew, Y.; Deasy, B.M. Human muscle-derived cell populations isolated by differential adhesion rates: Phenotype and contribution to skeletal muscle regeneration in Mdx/SCID mice. Tissue Eng. Part 2012, 18, 232–241. [Google Scholar] [CrossRef] [Green Version]
Cells 09 01723 g001 550
Figure 1. Cardiotoxin-induced injury on wild type (wt) and a mouse model of Duchenne muscular dystrophy (mdx) skeletal muscle tissue. (A) Experimental procedure. 45-day-old wt and mdx mice were injected intramuscularly with cardiotoxin (CTX) (10 μM), and the skeletal muscles were analyzed 1, 3, 5, 10, and 20 days (d) after injury. (B) Representative hematoxylin and eosin staining of uninjured wt tibialis anterior (TA) muscles and regenerating wt TA muscles at 1, 3, 5, 10, and 20 days after intramuscular CTX injection. Regenerating muscles were characterized by centrally located nuclei at day 5, but reconstituted multinucleated myofibers by day 10. (C) Representative hematoxylin and eosin staining on histological sections of uninjured and regenerating mdx TA muscles at 1, 3, 5, 10, and 20 days after intramuscular CTX injection. All along the considered time points, mdx muscles were characterized by infiltrating inflammatory interstitial cells and centrally nucleated myofibers of different sizes. (D,E) The number of cells (in millions) extracted from uninjured and CTX-injured wt (D) and mdx (E) mice (n=3; for 5d wt time point, n=2). All data were represented as mean ± SEM, and the statistical significance was estimated by two-way ANOVA (**** p < 0.0001). (B,C) 20× magnification; scale bar: 100 μm.
Figure 1. Cardiotoxin-induced injury on wild type (wt) and a mouse model of Duchenne muscular dystrophy (mdx) skeletal muscle tissue. (A) Experimental procedure. 45-day-old wt and mdx mice were injected intramuscularly with cardiotoxin (CTX) (10 μM), and the skeletal muscles were analyzed 1, 3, 5, 10, and 20 days (d) after injury. (B) Representative hematoxylin and eosin staining of uninjured wt tibialis anterior (TA) muscles and regenerating wt TA muscles at 1, 3, 5, 10, and 20 days after intramuscular CTX injection. Regenerating muscles were characterized by centrally located nuclei at day 5, but reconstituted multinucleated myofibers by day 10. (C) Representative hematoxylin and eosin staining on histological sections of uninjured and regenerating mdx TA muscles at 1, 3, 5, 10, and 20 days after intramuscular CTX injection. All along the considered time points, mdx muscles were characterized by infiltrating inflammatory interstitial cells and centrally nucleated myofibers of different sizes. (D,E) The number of cells (in millions) extracted from uninjured and CTX-injured wt (D) and mdx (E) mice (n=3; for 5d wt time point, n=2). All data were represented as mean ± SEM, and the statistical significance was estimated by two-way ANOVA (**** p < 0.0001). (B,C) 20× magnification; scale bar: 100 μm.
Cells 09 01723 g001
Cells 09 01723 g002 550
Figure 2. Dynamic changes of mononuclear cell subpopulations in wt muscles during regeneration. (A) The different mononuclear cell samples from uninjured and CTX-injured wt hind limb muscles were merged to create a single t-distributed stochastic neighbor embedding (t-SNE) map; colored according to expression levels of SCA1, CD90.2, α7-integrin, CD146, CD34, CD25, CD45, CD31, CD206, F4/80, CD140β, CXCR4, vimentin, and CD140α (blue: low expression; red: high expression of the selected marker). (B) FlowSOM heatmap of column normalized (Z-score) marker expression for each of the 15 identified clusters. Colors varied according to the expression level of the considered marker in a blue to red scale, indicating low to high expression, respectively. (C) Cell clusters defined by the FlowSOM analysis were assigned to arbitrary colors (yellow: immune cells; light blue: macrophages; pink: myogenic progenitors (MPs); dark green: fibro/adipogenic progenitors (FAPs); red: endothelial progenitors; blue: pericyte-like cells; light green: mesenchymal-like cells; purple: other). (D) Time course of the variation in cell subpopulation abundance upon CTX-induced injury. The dynamic changes were illustrated by density plots, colored according to cell density (blue: low density; red: high density).
Figure 2. Dynamic changes of mononuclear cell subpopulations in wt muscles during regeneration. (A) The different mononuclear cell samples from uninjured and CTX-injured wt hind limb muscles were merged to create a single t-distributed stochastic neighbor embedding (t-SNE) map; colored according to expression levels of SCA1, CD90.2, α7-integrin, CD146, CD34, CD25, CD45, CD31, CD206, F4/80, CD140β, CXCR4, vimentin, and CD140α (blue: low expression; red: high expression of the selected marker). (B) FlowSOM heatmap of column normalized (Z-score) marker expression for each of the 15 identified clusters. Colors varied according to the expression level of the considered marker in a blue to red scale, indicating low to high expression, respectively. (C) Cell clusters defined by the FlowSOM analysis were assigned to arbitrary colors (yellow: immune cells; light blue: macrophages; pink: myogenic progenitors (MPs); dark green: fibro/adipogenic progenitors (FAPs); red: endothelial progenitors; blue: pericyte-like cells; light green: mesenchymal-like cells; purple: other). (D) Time course of the variation in cell subpopulation abundance upon CTX-induced injury. The dynamic changes were illustrated by density plots, colored according to cell density (blue: low density; red: high density).
Cells 09 01723 g002
Cells 09 01723 g003 550
Figure 3. Characterization of the cell-population density rearrangements induced by CTX in the wt skeletal muscle after 1, 3, 5, 10, 20 days. (A) Bar plots, showing the CD45- (blue) / CD45+ (red) ratio in mononuclear cells from wt skeletal muscles at different time points after CTX injury. Data were represented as mean, while the statistical significance was estimated by two-way ANOVA. (B–H) Identification in the t-SNE maps of mononuclear cell populations. The different plots were color-coded according to the expression of surface antigens that characterize the relevant cell types. The bar plots quantitated the variation in population abundance in the wt limb muscles at different times during regeneration. Cell percentages were assessed on the total number of cells in each sample (n = 3; for 5d wt time point, n = 2). The statistical significance was estimated by one-way ANOVA. All data were represented as mean ± SEM, and the statistical significance was defined as * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Figure 3. Characterization of the cell-population density rearrangements induced by CTX in the wt skeletal muscle after 1, 3, 5, 10, 20 days. (A) Bar plots, showing the CD45- (blue) / CD45+ (red) ratio in mononuclear cells from wt skeletal muscles at different time points after CTX injury. Data were represented as mean, while the statistical significance was estimated by two-way ANOVA. (B–H) Identification in the t-SNE maps of mononuclear cell populations. The different plots were color-coded according to the expression of surface antigens that characterize the relevant cell types. The bar plots quantitated the variation in population abundance in the wt limb muscles at different times during regeneration. Cell percentages were assessed on the total number of cells in each sample (n = 3; for 5d wt time point, n = 2). The statistical significance was estimated by one-way ANOVA. All data were represented as mean ± SEM, and the statistical significance was defined as * p < 0.05; ** p < 0.01; *** p < 0.001; **** p < 0.0001.
Cells 09 01723 g003
Cells 09 01723 g004 550
Figure 4. Identification of the main mdx skeletal muscle cell populations upon CTX-induced injury. (A) A mononuclear cell suspension was purified from uninjured and CTX-injured mdx hind limb muscles. The different samples were combined, barcoded, and analyzed to create a single t-SNE map colored according to SCA1, CD90.2, α7-integrin, CD146, CD34, CD25, CD45, CD31, CD206, F4/80, CD140β, CXCR4, vimentin, and CD140α expression levels (blue: low expression; high expression of the selected marker). (B) FlowSOM heatmap of column normalized (Z-score) marker expression for each of the 15 clusters identified by Cytofkit analysis. Colors varied according to the expression level of each considered marker in a blue to red scale, indicating low and high expression, respectively. (C) Cell population clusters defined by overlapping the expression for each of the different analyzed antigens were projected onto t-SNE space and assigned to specific colors (yellow: immune cells; light blue: macrophages; pink: myogenic progenitors (MPs); dark green: fibro/adipogenic progenitors (FAPs); red: endothelial progenitors; blue: pericyte-like cells; purple: other). (D) Density plots colored by density (blue: low density; red: high density), showing cell abundance variations at different times during the regeneration process.
Figure 4. Identification of the main mdx skeletal muscle cell populations upon CTX-induced injury. (A) A mononuclear cell suspension was purified from uninjured and CTX-injured mdx hind limb muscles. The different samples were combined, barcoded, and analyzed to create a single t-SNE map colored according to SCA1, CD90.2, α7-integrin, CD146, CD34, CD25, CD45, CD31, CD206, F4/80, CD140β, CXCR4, vimentin, and CD140α expression levels (blue: low expression; high expression of the selected marker). (B) FlowSOM heatmap of column normalized (Z-score) marker expression for each of the 15 clusters identified by Cytofkit analysis. Colors varied according to the expression level of each considered marker in a blue to red scale, indicating low and high expression, respectively. (C) Cell population clusters defined by overlapping the expression for each of the different analyzed antigens were projected onto t-SNE space and assigned to specific colors (yellow: immune cells; light blue: macrophages; pink: myogenic progenitors (MPs); dark green: fibro/adipogenic progenitors (FAPs); red: endothelial progenitors; blue: pericyte-like cells; purple: other). (D) Density plots colored by density (blue: low density; red: high density), showing cell abundance variations at different times during the regeneration process.
Cells 09 01723 g004
Cells 09 01723 g005 550
Figure 5. Characterization of the population dynamics induced by CTX injury in the mdx skeletal muscle at day 1, 3, 5, 10, 20 after injury. (A) Bar plots, showing the CD45- (blue) / CD45+ (red) ratio in the mdx mononuclear cells at different time points after CTX injury. Data were represented as mean, and the statistical significance was estimated by two-way ANOVA. (B) t-SNE maps, showing the population gated as macrophage, colored for CD45 and F4/80 expression (red: high expression; blue: low expression) together with the macrophage trend observed in uninjured and CTX-injured mdx hind limb muscles. (C) M2 macrophages, expressing the CD206 antigen. (D) t-SNE maps of the myogenic progenitor (MP) population colored according to the levels of α7-integrin expression (red: high expression) and MP population dynamics during regeneration. (E) t-SNE maps, representing fibro/adipogenic progenitors (FAPs) identified among uninjured and CTX-injured mdx hind limb muscles, colored according to the expression of SCA1, CD34, CD140α, CD90.2, and vimentin (red: high expression; blue: low expression). The bar plot illustrated the population dynamics during regeneration. (F) Bar plots, showing the abundance of two different subclusters (cluster 1 and cluster 6) of endothelial progenitor cells at different time points. (G) t-SNE maps of pericyte-like cells identified in uninjured and CTX-injured mdx hind limb cell populations. Cell clusters were colored according to the expression of CD146, α7-integrin, and CD140β (red: high expression; blue: low expression). The bar plot illustrated the population dynamics during regeneration. Population abundance was assessed by calculating the percentage of cells in any given population over the total number of mononuclear cells in each sample (n = 3). The statistical significance was estimated by one-way ANOVA. All data were represented as mean ± SEM, and the statistical significance was defined as * p < 0.05; ** p < 0.01; **** p < 0.001.
Figure 5. Characterization of the population dynamics induced by CTX injury in the mdx skeletal muscle at day 1, 3, 5, 10, 20 after injury. (A) Bar plots, showing the CD45- (blue) / CD45+ (red) ratio in the mdx mononuclear cells at different time points after CTX injury. Data were represented as mean, and the statistical significance was estimated by two-way ANOVA. (B) t-SNE maps, showing the population gated as macrophage, colored for CD45 and F4/80 expression (red: high expression; blue: low expression) together with the macrophage trend observed in uninjured and CTX-injured mdx hind limb muscles. (C) M2 macrophages, expressing the CD206 antigen. (D) t-SNE maps of the myogenic progenitor (MP) population colored according to the levels of α7-integrin expression (red: high expression) and MP population dynamics during regeneration. (E) t-SNE maps, representing fibro/adipogenic progenitors (FAPs) identified among uninjured and CTX-injured mdx hind limb muscles, colored according to the expression of SCA1, CD34, CD140α, CD90.2, and vimentin (red: high expression; blue: low expression). The bar plot illustrated the population dynamics during regeneration. (F) Bar plots, showing the abundance of two different subclusters (cluster 1 and cluster 6) of endothelial progenitor cells at different time points. (G) t-SNE maps of pericyte-like cells identified in uninjured and CTX-injured mdx hind limb cell populations. Cell clusters were colored according to the expression of CD146, α7-integrin, and CD140β (red: high expression; blue: low expression). The bar plot illustrated the population dynamics during regeneration. Population abundance was assessed by calculating the percentage of cells in any given population over the total number of mononuclear cells in each sample (n = 3). The statistical significance was estimated by one-way ANOVA. All data were represented as mean ± SEM, and the statistical significance was defined as * p < 0.05; ** p < 0.01; **** p < 0.001.
Cells 09 01723 g005
Table
Table 1. List of the metal-tagged antibodies used in the mass cytometry experiments.
Table 1. List of the metal-tagged antibodies used in the mass cytometry experiments.
AntibodyMetal
Anti-mouse CD45147Sm
Anti-mouse Ly-6A/E (SCA1)164Dy
Anti-mouse CD90.2 (Thy-1.2)156Gd
Anti-mouse CD146141Pr
Anti-mouse F4/80146Nd
Anti-mouse CD140α148Nd
Anti-mouse CD140β151Eu
Anti-mouse α7-integrin161Dy
Anti-mouse CD206169Tm
Anti-mouse CD34144Nd
Anti-mouse CXCR4159Tb
Anti-mouse CD4172Yb
Anti-mouse CD25 (IL-2R)150Nd
Anti-vimentin154Sm
Anti-CD31 (PECAM-1)165Ho
Anti-pan-actin175Lu
Anti-mouse interleukin-6 (IL-6)167Er
Anti-phospho-Akt (S473)152Sm
Anti-phospho-Stat1 (Y701)153Eu
Anti-phospho-Erk1/2 (T202/Y204)171Yb
Anti-phospho-Stat3 (Y705)158Gd
Anti-cleaved caspase3142Nd
Anti-phospho-Creb176Yb
CD: Cluster Differentiation; SCA1: Stem Cell Antigen1; CXCR4: C-X-C Motif Chemokine Receptor 4; PECAM-1: Platelet Endothelial Cell adhesion-1; Akt: RAC alpha serine/threonine-protein kinase; Stat1: Signal Transducer and activator of transcription 1; Erk1/2: Extracellular signal-regulated kinases 1, 2; Stat3: Signal transducer and activator of transcription 3; cleaved caspase3: cysteine-aspartic proteases; Creb: cAMP response element-binding protein.

Share and Cite

MDPI and ACS Style

Petrilli, L.L.; Spada, F.; Palma, A.; Reggio, A.; Rosina, M.; Gargioli, C.; Castagnoli, L.; Fuoco, C.; Cesareni, G. High-Dimensional Single-Cell Quantitative Profiling of Skeletal Muscle Cell Population Dynamics during Regeneration. Cells 2020, 9, 1723. https://doi.org/10.3390/cells9071723

AMA Style

Petrilli LL, Spada F, Palma A, Reggio A, Rosina M, Gargioli C, Castagnoli L, Fuoco C, Cesareni G. High-Dimensional Single-Cell Quantitative Profiling of Skeletal Muscle Cell Population Dynamics during Regeneration. Cells. 2020; 9(7):1723. https://doi.org/10.3390/cells9071723

Chicago/Turabian Style

Petrilli, Lucia Lisa, Filomena Spada, Alessandro Palma, Alessio Reggio, Marco Rosina, Cesare Gargioli, Luisa Castagnoli, Claudia Fuoco, and Gianni Cesareni. 2020. "High-Dimensional Single-Cell Quantitative Profiling of Skeletal Muscle Cell Population Dynamics during Regeneration" Cells 9, no. 7: 1723. https://doi.org/10.3390/cells9071723

APA Style

Petrilli, L. L., Spada, F., Palma, A., Reggio, A., Rosina, M., Gargioli, C., Castagnoli, L., Fuoco, C., & Cesareni, G. (2020). High-Dimensional Single-Cell Quantitative Profiling of Skeletal Muscle Cell Population Dynamics during Regeneration. Cells, 9(7), 1723. https://doi.org/10.3390/cells9071723

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