An ITGA11 expressing subpopulation as predictor for the donor-specific osteogenic capacity of stromal cells

Abstract

Stromal progenitor cells of bone marrow origin are non-hematopoietic cells that give rise to osteoblasts and adipocytes in the postnatal organism. Marrow stromal cells (also known as mesenchymal stem cells - MSCs) are currently being employed in a large number of clinical trials for regenerative purposes post in vitro expansion. However, the clinical outcome has been variable, which might in part be due to the heterogeneity of the cells and the lack of a defined cell product with a molecular signature that favors tissue regeneration. In this study, we determined the cellular heterogeneity of primary stromal cultures and examined how inter-donor variation in subpopulation composition contributes to the differentiation potential of primary cultures. We profiled 136 014 stromal progenitors from 26 donors and identified 5 subpopulations that were linked to distinct bone-related pathways and genetic traits of bone mineral density and morphology. Abundance of one cluster characterized by high expression of ITGA11 (integrin alpha-11) and genes related to matrix function, collagen organization, and elevated expression upon lineage commitment was positively correlated with osteoblastic differentiation capacity in vitro. In addition, ITGA11 protein expression in progenitor cells was a predictive marker for matrix mineralization in vitro and ectopic bone formation in vivo. Sorting stromal progenitors into ITGA11^high and ITGA11^low cells established cultures with high and low osteoblastic differentiation potential and revealed transcriptional differences reflective of the subpopulation-specific signature, which was not affected by siRNA-mediated knockdown of ITGA11 expression. Our findings corroborate the presence of an extensive donor-dependent cellular heterogeneity that persists in cultured stromal cells, and that ITGA11 can be employed as a marker for isolating cells with high bone-forming potential, a feature likely to benefit clinical trials of bone regeneration.

References

[1] Kassem M, Bianco P. Skeletal stem cells in space and time. Cell, 2015, 160: 17-19.

[2] Shi Y, et al.. How mesenchymal stem cells interact with tissue immune responses. Trends Immunol., 2012, 33: 136-143.

[3] Garcia-Bernal D, et al.. The current status of mesenchymal stromal cells: controversies, unresolved issues and some promising solutions to improve their therapeutic efficacy. Front. Cell Dev. Biol., 2021, 9: 650664.

[4] Miao C, et al.. A brief review: the therapeutic potential of bone marrow mesenchymal stem cells in myocardial infarction. Stem Cell Res. Ther., 2017, 8: 242.

[5] Zaher W, et al.. An update of human mesenchymal stem cell biology and their clinical uses. Arch. Toxicol., 2014, 88: 1069-1082.

[6] Trounson A, McDonald C. Stem cell therapies in clinical trials: progress and challenges. Cell Stem Cell, 2015, 17: 11-22.

[7] Mizoguchi T, et al.. Osterix marks distinct waves of primitive and definitive stromal progenitors during bone marrow development. Dev. Cell, 2014, 29: 340-349.

[8] Mendez-Ferrer S, et al.. Mesenchymal and haematopoietic stem cells form a unique bone marrow niche. Nature, 2010, 466: 829-834.

[9] Mizuhashi K, et al.. Resting zone of the growth plate houses a unique class of skeletal stem cells. Nature, 2018, 563: 254-258.

[10] Worthley DL, et al.. Gremlin 1 identifies a skeletal stem cell with bone, cartilage, and reticular stromal potential. Cell, 2015, 160: 269-284.

[11] Zhou BO, et al.. Leptin-receptor-expressing mesenchymal stromal cells represent the main source of bone formed by adult bone marrow. Cell Stem Cell, 2014, 15: 154-168.

[12] Rennert RC, et al.. High-resolution microfluidic single-cell transcriptional profiling reveals clinically relevant subtypes among human stem cell populations commonly utilized in cell-based therapies. Front. Neurol., 2016, 7: 41.

[13] Chan CKF, et al.. Identification of the human skeletal stem cell. Cell, 2018, 175: 43-56. e21.

[14] Khong SML, et al.. Single-cell transcriptomics of human mesenchymal stem cells reveal age-related cellular subpopulation depletion and impaired regenerative function. Stem Cells, 2019, 37: 240-246.

[15] Chan CK, et al.. Identification and specification of the mouse skeletal stem cell. Cell, 2015, 160: 285-298.

[16] Ambrosi TH, et al.. Human skeletal development and regeneration are shaped by functional diversity of stem cells across skeletal sites. Cell Stem Cell, 2025, 32: 811-823.

[17] To K, et al.. A multi-omic atlas of human embryonic skeletal development. Nature, 2024, 635: 657-667.

[18] Boulais PE, et al.. The majority of CD45(-) Ter119(-) CD31(-) bone marrow cell fraction is of hematopoietic origin and contains erythroid and lymphoid progenitors. Immunity, 2018, 49: 627-639.e626.

[19] Bandyopadhyay S, et al.. Mapping the cellular biogeography of human bone marrow niches using single-cell transcriptomics and proteomic imaging. Cell, 2024, 187: 3120-3140.e3129.

[20] Li H, et al.. Identification of phenotypically, functionally, and anatomically distinct stromal niche populations in human bone marrow based on single-cell RNA sequencing. eLife, 2023, 12: e8165

[21] Kowal JM, et al.. Single-cell high-content imaging parameters predict functional phenotype of cultured human bone marrow stromal stem cells. Stem Cells Transl. Med., 2020, 9: 189-202.

[22] Kowal JM, et al.. Identification of a clinical signature predictive of differentiation fate of human bone marrow stromal cells. Stem Cell Res. Ther., 2021, 12: 265.

[23] Dominici M, et al.. Minimal criteria for defining multipotent mesenchymal stromal cells. The International Society for Cellular Therapy position statement. Cytotherapy, 2006, 8: 315-317.

[24] Omatsu Y, et al.. The essential functions of adipo-osteogenic progenitors as the hematopoietic stem and progenitor cell niche. Immunity, 2010, 33: 387-399.

[25] Groza T, et al.. The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease. Nucleic Acids Res., 2023, 51: D1038-D1045.

[26] Morris JA, et al.. An atlas of genetic influences on osteoporosis in humans and mice. Nat. Genet., 2019, 51: 258-266.

[27] Rauch A, et al.. Osteogenesis depends on commissioning of a network of stem cell transcription factors that act as repressors of adipogenesis. Nat. Genet., 2019, 51: 716-727.

[28] Simonsen JL, et al.. Telomerase expression extends the proliferative life-span and maintains the osteogenic potential of human bone marrow stromal cells. Nat. Biotechnol., 2002, 20: 592-596.

[29] Bausch-Fluck D, et al.. A mass spectrometric-derived cell surface protein atlas. PLoS ONE, 2015, 10: e0121314.

[30] Shen B, et al.. Integrin alpha11 is an Osteolectin receptor and is required for the maintenance of adult skeletal bone mass. eLife, 2019, 8: e42274.

[31] Galgaro BC, et al.. The adenosinergic pathway in mesenchymal stem cell fate and functions. Med. Res. Rev., 2021, 41: 2316-2349.

[32] Hemler ME. Tetraspanin functions and associated microdomains. Nat. Rev. Mol. Cell Biol., 2005, 6: 801-811.

[33] Abdallah BM, Ditzel N, Kassem M. Assessment of bone formation capacity using in vivo transplantation assays: procedure and tissue analysis. Methods Mol. Biol., 2008, 455: 89-100.

[34] Boraldi F, et al.. Mineralization by mesenchymal stromal cells is variously modulated depending on commercial platelet lysate preparations. Cytotherapy, 2018, 20: 335-342.

[35] Shanbhag S, et al.. Influence of platelet storage time on human platelet lysates and platelet lysate-expanded mesenchymal stromal cells for bone tissue engineering. Stem Cell Res Ther., 2020, 11: 351.

[36] Chen L, et al.. KIAA1199 deficiency enhances skeletal stem cell differentiation to osteoblasts and promotes bone regeneration. Nat. Commun., 2023, 14. ArticleID: 2016

[37] Gao Y, et al.. Multi-omics analysis of human mesenchymal stem cells shows cell aging that alters immunomodulatory activity through the downregulation of PD-L1. Nat. Commun., 2023, 14. ArticleID: 4373

[38] Wang Z, et al.. Single-cell transcriptome atlas of human mesenchymal stem cells exploring cellular heterogeneity. Clin. Transl. Med., 2021, 11. ArticleID: e650

[39] Xie Z, et al.. Single-cell RNA sequencing analysis of human bone-marrow-derived mesenchymal stem cells and functional subpopulation identification. Exp. Mol. Med., 2022, 54: 483-492.

[40] Zhang C, et al.. Single-cell transcriptomic analysis reveals the cellular heterogeneity of mesenchymal stem cells. Genomics Proteom. Bioinform., 2022, 20: 70-86.

[41] Post S, et al.. Demonstration of the presence of independent pre-osteoblastic and pre-adipocytic cell populations in bone marrow-derived mesenchymal stem cells. Bone, 2008, 43: 32-39.

[42] Larsen KH, et al.. Identifying a molecular phenotype for bone marrow stromal cells with in vivo bone-forming capacity. J. Bone Miner. Res., 2010, 25: 796-808.

[43] Sundelacruz S, Levin M, Kaplan DL. Comparison of the depolarization response of human mesenchymal stem cells from different donors. Sci. Rep., 2015, 5. ArticleID: 18279

[44] O’Keefe RJ, et al.. American society for bone and mineral research-orthopaedic research society joint task force report on cell-based therapies. J. Bone Min. Res., 2020, 35: 3-17.

[45] Oryan A, et al.. Role of mesenchymal stem cells in bone regenerative medicine: what is the evidence?. Cells Tissues Organs, 2017, 204: 59-83.

[46] Gomez-Barrena E, et al.. Bone regeneration: stem cell therapies and clinical studies in orthopaedics and traumatology. J. Cell Mol. Med., 2011, 15: 1266-1286.

[47] Halfon S, et al.. Markers distinguishing mesenchymal stem cells from fibroblasts are downregulated with passaging. Stem Cells Dev., 2011, 20: 53-66.

[48] Perez-Silos V, Camacho-Morales A, Fuentes-Mera L. Mesenchymal stem cells subpopulations: application for orthopedic regenerative medicine. Stem Cells Int., 2016, 2016: 3187491.

[49] Chung MT, et al.. CD90 (Thy-1)-positive selection enhances osteogenic capacity of human adipose-derived stromal cells. Tissue Eng. Part A, 2013, 19: 989-997.

[50] Abdallah BM, et al.. CD34 defines an osteoprogenitor cell population in mouse bone marrow stromal cells. Stem Cell Res., 2015, 15: 449-458.

[51] Harkness L, et al.. CD146/MCAM defines functionality of human bone marrow stromal stem cell populations. Stem Cell Res. Ther., 2016, 7: 4.

[52] Wagner W, et al.. Replicative senescence of mesenchymal stem cells: a continuous and organized process. PLoS One, 2008, 3: e2213.

[53] Pittenger MF, et al.. Multilineage potential of adult human mesenchymal stem cells. Science, 1999, 284: 143-147.

[54] Simmons PJ, Torok-Storb B. Identification of stromal cell precursors in human bone marrow by a novel monoclonal antibody, STRO-1. Blood, 1991, 78: 55-62.

[55] Gronthos S, et al.. Molecular and cellular characterisation of highly purified stromal stem cells derived from human bone marrow. J. Cell Sci., 2003, 116: 1827-1835.

[56] Gothard D, et al.. Prospective isolation of human bone marrow stromal cell subsets: A comparative study between Stro-1-, CD146- and CD105-enriched populations. J. Tissue Eng., 2014, 5: 2041731414551763.

[57] Lee RH, et al.. The CD34-like protein PODXL and alpha6-integrin (CD49f) identify early progenitor MSCs with increased clonogenicity and migration to infarcted heart in mice. Blood, 2009, 113: 816-826.

[58] Costa MA, et al.. On the role of subtype selective adenosine receptor agonists during proliferation and osteogenic differentiation of human primary bone marrow stromal cells. J. Cell Physiol., 2011, 226: 1353-1366.

[59] Takedachi M, et al.. CD73-generated adenosine promotes osteoblast differentiation. J. Cell Physiol., 2012, 227: 2622-2631.

[60] Shih YV, et al.. Dysregulation of ectonucleotidase-mediated extracellular adenosine during postmenopausal bone loss. Sci. Adv., 2019, 5: eaax1387.

[61] Canepa DD, et al.. Identification of ALP⁺/CD73⁺ defining markers for enhanced osteogenic potential in human adipose-derived mesenchymal stromal cells by mass cytometry. Stem Cell Res. Ther., 2021, 12: 7.

[62] Yue, R., Shen, B, Morrison, S. J. Clec11a/osteolectin is an osteogenic growth factor that promotes the maintenance of the adult skeleton. eLife5, e18782 (2016).

[63] Dadras M, et al.. Comparative proteomic analysis of osteogenic differentiated human adipose tissue and bone marrow-derived stromal cells. J. Cell Mol. Med., 2020, 24: 11814-11827.

[64] Fregni G, et al.. Reciprocal modulation of mesenchymal stem cells and tumor cells promotes lung cancer metastasis. eBioMedicine, 2018, 29: 128-145.

[65] Kaltz N, et al.. Novel markers of mesenchymal stem cells defined by genome-wide gene expression analysis of stromal cells from different sources. Exp. Cell Res., 2010, 316: 2609-2617.

[66] Popov C, et al.. Integrins alpha2beta1 and alpha11beta1 regulate the survival of mesenchymal stem cells on collagen I. Cell Death Dis., 2011, 2. ArticleID: e186

[67] Tsai CC, et al.. Hypoxia inhibits senescence and maintains mesenchymal stem cell properties through down-regulation of E2A-p21 by HIF-TWIST. Blood, 2011, 117: 459-469.

[68] Lennon DP, Edmison JM, Caplan AI. Cultivation of rat marrow-derived mesenchymal stem cells in reduced oxygen tension: effects on in vitro and in vivo osteochondrogenesis. J. Cell Physiol., 2001, 187: 345-355.

[69] Du G, et al.. TGF-beta induced key genes of osteogenic and adipogenic differentiation in human mesenchymal stem cells and MiRNA-mRNA regulatory networks. Front. Genet., 2021, 12: 759596.

[70] Saller MM, et al.. Increased stemness and migration of human mesenchymal stem cells in hypoxia is associated with altered integrin expression. Biochem. Biophys. Res Commun., 2012, 423: 379-385.

[71] Andreasen CM, et al.. Understanding age-induced cortical porosity in women: the accumulation and coalescence of eroded cavities upon existing intracortical canals is the main contributor. J. Bone Min. Res., 2018, 33: 606-620.

[72] Vedi S, Elkin SL, Compston JE. A histomorphometric study of cortical bone of the iliac crest in patients treated with glucocorticoids. Calcif. Tissue Int., 2005, 77: 79-83.

[73] Dalle Carbonare, et al.. Comparison of trabecular bone microarchitecture and remodeling in glucocorticoid-induced and postmenopausal osteoporosis. J. Bone Min. Res., 2001, 16: 97-103.

[74] Stenderup K, et al.. Aging is associated with decreased maximal life span and accelerated senescence of bone marrow stromal cells. Bone, 2003, 33: 919-926.

[75] Walzl A, et al.. The resazurin reduction assay can distinguish cytotoxic from cytostatic compounds in spheroid screening assays. J. Biomol. Screen., 2014, 19: 1047-1059.

[76] Parekh, S. et al. zUMIs - A fast and flexible pipeline to process RNA sequencing data with UMIs. Gigascience7, giy059 (2018).

[77] Lun ATL, et al.. EmptyDrops: distinguishing cells from empty droplets in droplet-based single-cell RNA sequencing data. Genome Biol., 2019, 20. ArticleID: 63

[78] McCarthy DJ, et al.. Scater: pre-processing, quality control, normalization and visualization of single-cell RNA-seq data in R. Bioinformatics, 2017, 33: 1179-1186.

[79] Lun AT, McCarthy DJ, Marioni JC. A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor. F1000Res, 2016, 5: 2122

[80] Butler A, et al.. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol., 2018, 36: 411-420.

[81] Korsunsky I, et al.. Fast, sensitive and accurate integration of single-cell data with Harmony. Nat. Methods, 2019, 16: 1289-1296.

[82] Dobin A, et al.. STAR: ultrafast universal RNA-seq aligner. Bioinformatics, 2013, 29: 15-21.

[83] Liao Y, Smyth GK, Shi W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics, 2014, 30: 923-930.

[84] Love MI, Huber W, Anders S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol., 2014, 15. ArticleID: 550

[85] Pico AR, et al.. WikiPathways: pathway editing for the people. PLoS Biol., 2008, 6: e184.

[86] Ashburner M, et al.. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat. Genet., 2000, 25: 25-29.

[87] Young MD, et al.. Gene ontology analysis for RNA-seq: accounting for selection bias. Genome Biol., 2010, 11. ArticleID: R14

[88] Lawrence M, et al.. Software for computing and annotating genomic ranges. PLoS Comput. Biol., 2013, 9: e1003118.

[89] Korotkevich G. et al. Fast gene set enrichment analysis. Preprint at bioRxivhttps://www.biorxiv.org/content/10.1101/060012v3 (2021).

Funding

Lundbeckfonden (Lundbeck Foundation)(R335-2019-2195)

Novo Nordisk Fonden (Novo Nordisk Foundation)(NNF22OC0078257)

Det Frie Forskningsråd (Danish Council for Independent Research)(0134-00292B)

Odense Universitetshospital (Odense University Hospital)(98-A4534)

Danish 3R-center - https://en.3rcenter.dk/

RIGHTS & PERMISSIONS

The Author(s)