CoChat
Sign Up
Cited by 323Single-cell and spatial transcriptomicsMuscle Physiology and DisordersAdipose Tissue and Metabolism

The DOI matches its Crossref record; 72 of the 72 references it cites were found in a registry or an index.

Single-cell transcriptional profiles in human skeletal muscle

Open access
Rubenstein, Aliza B., Smith, Gregory R., Raue, Ulrika, Begue, Gwénaëlle, Minchev, Kiril, Ruf-Zamojski, Frederique, Nair, Venugopalan D., Wang, Xingyu, Zhou, Lan, Zaslavsky, Elena, Trappe, Todd A., Trappe, Scott, Sealfon, Stuart C.·2020·Scientific Reports
DOI 10.1038/s41598-019-57110-6PMC PMC6959232
Open this paper in CoChat — free accountOpen PDF

Abstract

Skeletal muscle is a heterogeneous tissue comprised of muscle fiber and mononuclear cell types that, in addition to movement, influences immunity, metabolism and cognition. We investigated the gene expression patterns of skeletal muscle cells using RNA-seq of subtype-pooled single human muscle fibers and single cell RNA-seq of mononuclear cells from human vastus lateralis, mouse quadriceps, and mouse diaphragm. We identified 11 human skeletal muscle mononuclear cell types, including two fibro-adipogenic progenitor (FAP) cell subtypes. The human FBN1+ FAP cell subtype is novel and a corresponding FBN1+ FAP cell type was also found in single cell RNA-seq analysis in mouse. Transcriptome exercise studies using bulk tissue analysis do not resolve changes in individual cell-type proportion or gene expression. The cell-type gene signatures provide the means to use computational methods to identify cell-type level changes in bulk studies. As an example, we analyzed public transcriptome data from an exercise training study and revealed significant changes in specific mononuclear cell-type proportions related to age, sex, acute exercise and training. Our single-cell expression map of skeletal muscle cell types will further the understanding of the diverse effects of exercise and the pathophysiology of muscle disease.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.1038/s41598-019-57110-6

Cites

72 of 72 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1 Schematic showing workflow of gene signature determination and skeletal muscle tissue deconvolution. Top (1): a skeletal muscle biopsy from the vastus lateralis obtained from a young healthy male was prepared by dissociation and filtering of skeletal muscle fibers. The cells were subjected to scRNA-seq using 10X Chromium. Normalization and clustering was performed by Seurat and clusters were manually identified as cell types. Top right (2): Confirmation of selected human skeletal muscle scRNA-seq cell types by mouse skeletal muscle scRNA-seq assay. Middle (3): skeletal muscle biopsies from the vastus lateralis of nine old healthy males were each used to isolate 96 individual muscle fibers. A small portion of each muscle fiber was clipped and subjected to SDS-PAGE to determine the MHC isoform. For each subject, all Type I fibers were pooled and all Type IIa fibers were pooled, for a total of 18 fiber-type specific samples which were then subjected to RNA-seq. After quality control and normalization, gene expression was averaged over all Type I samples and all Type IIa samples separately to obtain fiber-type specific gene signatures. Bottom (4): a previously published dataset8, consisting of microarray profiling of skeletal muscle biopsies from the vastus lateralis of 28 healthy young and old males and females was deconvolved using the previously determined mononuclear and multinucleated cell-type gene signatures. Biopsies were obtained before and 4 h after an acute resistance exercise bout at the onset and end of 12 weeks of resistance training (3 d/wk).
Figure 1 Schematic showing workflow of gene signature determination and skeletal muscle tissue deconvolution. Top (1): a skeletal muscle biopsy from the vastus lateralis obtained from a young healthy male was prepared by dissociation and filtering of skeletal muscle fibers. The cells were subjected to scRNA-seq using 10X Chromium. Normalization and clustering was performed by Seurat and clusters were manually identified as cell types. Top right (2): Confirmation of selected human skeletal muscle scRNA-seq cell types by mouse skeletal muscle scRNA-seq assay. Middle (3): skeletal muscle biopsies from the vastus lateralis of nine old healthy males were each used to isolate 96 individual muscle fibers. A small portion of each muscle fiber was clipped and subjected to SDS-PAGE to determine the MHC isoform. For each subject, all Type I fibers were pooled and all Type IIa fibers were pooled, for a total of 18 fiber-type specific samples which were then subjected to RNA-seq. After quality control and normalization, gene expression was averaged over all Type I samples and all Type IIa samples separately to obtain fiber-type specific gene signatures. Bottom (4): a previously published dataset8, consisting of microarray profiling of skeletal muscle biopsies from the vastus lateralis of 28 healthy young and old males and females was deconvolved using the previously determined mononuclear and multinucleated cell-type gene signatures. Biopsies were obtained before and 4 h after an acute resistance exercise bout at the onset and end of 12 weeks of resistance training (3 d/wk).CC BY · PubMed Central
Figure 2 Single-cell RNA-seq cell-type analysis of human mononuclear muscle cells. (a) Cell-type labeled t-SNE plot of mononuclear cells from a combined set of four samples of a muscle biopsy. Cells are colored by their expression of top cell-type differentiating markers. Gray cells do not express any top cell-type differentiating markers, which may be due to transcript drop-out. (b) Bar graph of cell-type composition for each of the four muscle samples highlighting the heterogeneity of muscle cell-type composition. (c) Dot plot of gene expression of three top cell-type differentiating markers. Dot color reflects average gene expression and dot size represents percent of cells expressing the gene.
Figure 2 Single-cell RNA-seq cell-type analysis of human mononuclear muscle cells. (a) Cell-type labeled t-SNE plot of mononuclear cells from a combined set of four samples of a muscle biopsy. Cells are colored by their expression of top cell-type differentiating markers. Gray cells do not express any top cell-type differentiating markers, which may be due to transcript drop-out. (b) Bar graph of cell-type composition for each of the four muscle samples highlighting the heterogeneity of muscle cell-type composition. (c) Dot plot of gene expression of three top cell-type differentiating markers. Dot color reflects average gene expression and dot size represents percent of cells expressing the gene.CC BY · PubMed Central
Figure 3 Human mononuclear muscle cell-type clustering analysis. (a) Heatmap and dendrogram of mononuclear cell types in muscle, clustered by their expression of top cell-type differentiating marker genes. (b) Box plot of cell-cell pearson correlation coefficients for each cell type. Higher values for a given cell type suggest greater homogeneity of gene expression between the cells of that cell type.
Figure 3 Human mononuclear muscle cell-type clustering analysis. (a) Heatmap and dendrogram of mononuclear cell types in muscle, clustered by their expression of top cell-type differentiating marker genes. (b) Box plot of cell-cell pearson correlation coefficients for each cell type. Higher values for a given cell type suggest greater homogeneity of gene expression between the cells of that cell type.CC BY · PubMed Central
Figure 4 Single-cell RNA-seq cell-type analysis of mouse mononuclear muscle cells. (a) UMAP plot of FAP cells isolated from quadriceps and diaphragm, colored by their expression of FBN1 (green) and LUM (orange). (b) Correlation plot of the gene expression for 11 top markers for each human FAP cell subtype (FBN1+ and LUM+) over the mouse samples. Genes are arranged via hierarchical clustering with their human FAP subtype marker assignment labeled by the black bars on the left of the figure. The location of FBN1 and LUM are emphasized by boxes.
Figure 4 Single-cell RNA-seq cell-type analysis of mouse mononuclear muscle cells. (a) UMAP plot of FAP cells isolated from quadriceps and diaphragm, colored by their expression of FBN1 (green) and LUM (orange). (b) Correlation plot of the gene expression for 11 top markers for each human FAP cell subtype (FBN1+ and LUM+) over the mouse samples. Genes are arranged via hierarchical clustering with their human FAP subtype marker assignment labeled by the black bars on the left of the figure. The location of FBN1 and LUM are emphasized by boxes.CC BY · PubMed Central
Figure 5 Fiber-type gene signatures and fiber-type specific tissue deconvolution. (a) Heatmap of gene expression for twenty markers per fiber-type over eighteen fiber-type specific tissue samples. Heatmap values are regularized-log transformed gene expression values. (b) Correlation heatmap for twenty gene markers per fiber-type. Estimated cell-type proportions (SPVs) for each fiber-type delineated in black; SPVs correlate with gene markers for each fiber-type. (c) Box plots showing estimated proportions of Type I fibers (left plot) and Type IIa fibers (right plot) within Type I specific tissue samples (orange boxes) and Type IIa specific tissue samples (blue boxes).
Figure 5 Fiber-type gene signatures and fiber-type specific tissue deconvolution. (a) Heatmap of gene expression for twenty markers per fiber-type over eighteen fiber-type specific tissue samples. Heatmap values are regularized-log transformed gene expression values. (b) Correlation heatmap for twenty gene markers per fiber-type. Estimated cell-type proportions (SPVs) for each fiber-type delineated in black; SPVs correlate with gene markers for each fiber-type. (c) Box plots showing estimated proportions of Type I fibers (left plot) and Type IIa fibers (right plot) within Type I specific tissue samples (orange boxes) and Type IIa specific tissue samples (blue boxes).CC BY · PubMed Central
Figure 6 Fiber-type deconvolution of microarray dataset. (a) Correlation heatmap of twenty marker genes for each fiber-types. Estimated cell-type proportions (SPVs) for each fiber-type delineated in black; SPVs correlate with gene markers for each fiber-type. (b,c) Correlation of estimated Type I and Type IIa fiber-type proportions with biochemically measured slow-twitch and fast-twitch fiber-type proportions, respectively (slow-twitch: p = 0.02, fast-twitch: p = 0.03).
Figure 6 Fiber-type deconvolution of microarray dataset. (a) Correlation heatmap of twenty marker genes for each fiber-types. Estimated cell-type proportions (SPVs) for each fiber-type delineated in black; SPVs correlate with gene markers for each fiber-type. (b,c) Correlation of estimated Type I and Type IIa fiber-type proportions with biochemically measured slow-twitch and fast-twitch fiber-type proportions, respectively (slow-twitch: p = 0.02, fast-twitch: p = 0.03).CC BY · PubMed Central
Figure 7 PLIER deconvolution of all skeletal muscle cell types. (a) Heatmap of association of cell types with LVs. The heatmap scale is arbitrary and should only be used to compare one association with another. (b) Difference in LV59 (endothelial cells/pericytes) between trained and untrained subjects. Points represent means of all trained and all untrained samples and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old). (four-way ANOVA p = 4.5e–13) (c) Difference in LV45 (lymphocytes) between male and female subjects (four-way ANOVA p = 0.00056). Points represent means of all male and all female subjects and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old). (d) Difference in LV33 (myeloid cells [neutrophils]) between pre- and post-acute exercise bout (four-way ANOVA p = 2.3e−11). Points represent means of all pre- and all 4 hour post-acute exercise samples and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old).
Figure 7 PLIER deconvolution of all skeletal muscle cell types. (a) Heatmap of association of cell types with LVs. The heatmap scale is arbitrary and should only be used to compare one association with another. (b) Difference in LV59 (endothelial cells/pericytes) between trained and untrained subjects. Points represent means of all trained and all untrained samples and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old). (four-way ANOVA p = 4.5e–13) (c) Difference in LV45 (lymphocytes) between male and female subjects (four-way ANOVA p = 0.00056). Points represent means of all male and all female subjects and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old). (d) Difference in LV33 (myeloid cells [neutrophils]) between pre- and post-acute exercise bout (four-way ANOVA p = 2.3e−11). Points represent means of all pre- and all 4 hour post-acute exercise samples and error bars represent standard error. Both cohorts are depicted (orange for young and blue for old).CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Cited by

Sorted by citations

See all 323

Related

See all

Metadata

Authors
Rubenstein, A. B., Smith, G. R., Raue, U., Begue, G., Minchev, K., Ruf-Zamojski, F., Nair, V. D., Wang, X., Zhou, L., Zaslavsky, E., Trappe, T. A., Trappe, S., & Sealfon, S. C.
Affiliations
  • Icahn School of Medicine at Mount Sinai (US)
  • Ball State University (US)
  • Boston University (US)
Corresponding author
Sealfon, Stuart C.
ORCID
Venue
Scientific Reports
Publisher
Springer Science and Business Media LLC
Published
2020-01-14
Volume
10
Issue
1
Pages
229-229
ISSN
2045-2322 (ISSN-L)
DOI
10.1038/s41598-019-57110-6
PMC
PMC6959232
OpenAlex
W2999612022
Cited by
323
References
72
Citations by year
2026: 41 · 2025: 53 · 2024: 57 · 2023: 59 · 2022: 52 · 2021: 37 · 2020: 24
FWCI
13.94
Type
journal article
Language
en
Open access
Licence
Locations
  • Scientific Reports (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion)
Funders
  • U.S. Department of Health & Human Services | NIH | National Institute of Arthritis and Musculoskeletal and Skin Diseases (1R01AR074428)
  • U.S. Department of Health & Human Services | National Institutes of Health (AG038576, AG184809, U24 DK112331)
  • Ball State University (Academic Excellence Award)
  • Empire State Development's Division of Science, Technology and Innovation
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Muscle Physiology and Disorders
  • Health Sciences › Medicine › Physiology (medical) › Adipose Tissue and Metabolism
Keywords
skeletal muscle, muscle fiber heterogeneity, fibro-adipogenic progenitor cells, exercise training, vastus lateralis, quadriceps, diaphragm, single-cell transcriptomics, cell type proportions
MeSH terms
Animals, Diaphragm, Female, Humans, Male, Biomarkers, Muscle, Skeletal, Adipogenesis, Mice, Quadriceps Muscle, Single-Cell Analysis, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review