Impaired mitochondrial oxidative metabolism in skeletal progenitor cells leads to musculoskeletal disintegration

Lin, Chujiao; Yang, Qiyuan; Guo, Dongsheng; Xie, Jun; Yang, Yeon-Suk; Chaugule, Sachin; DeSouza, Ngoc; Oh, Won-Taek; Li, Rui; Chen, Zhihao; John, Aijaz A.; Qiu, Qiang; Zhu, Lihua Julie; Greenblatt, Matthew B.; Ghosh, Sankar; Li, Shaoguang; Gao, Guangping; Haynes, Cole; Emerson, Charles P.; Shim, Jae-Hyuck

doi:10.1038/s41467-022-34694-8

Download PDF

Article
Open access
Published: 11 November 2022

Impaired mitochondrial oxidative metabolism in skeletal progenitor cells leads to musculoskeletal disintegration

Nature Communications volume 13, Article number: 6869 (2022) Cite this article

4464 Accesses
10 Citations
5 Altmetric
Metrics details

Subjects

Abstract

Although skeletal progenitors provide a reservoir for bone-forming osteoblasts, the major energy source for their osteogenesis remains unclear. Here, we demonstrate a requirement for mitochondrial oxidative phosphorylation in the osteogenic commitment and differentiation of skeletal progenitors. Deletion of Evolutionarily Conserved Signaling Intermediate in Toll pathways (ECSIT) in skeletal progenitors hinders bone formation and regeneration, resulting in skeletal deformity, defects in the bone marrow niche and spontaneous fractures followed by persistent nonunion. Upon skeletal fracture, Ecsit-deficient skeletal progenitors migrate to adjacent skeletal muscle causing muscle atrophy. These phenotypes are intrinsic to ECSIT function in skeletal progenitors, as little skeletal abnormalities were observed in mice lacking Ecsit in committed osteoprogenitors or mature osteoblasts. Mechanistically, Ecsit deletion in skeletal progenitors impairs mitochondrial complex assembly and mitochondrial oxidative phosphorylation and elevates glycolysis. ECSIT-associated skeletal phenotypes were reversed by in vivo reconstitution with wild-type ECSIT expression, but not a mutant displaying defective mitochondrial localization. Collectively, these findings identify mitochondrial oxidative phosphorylation as the prominent energy-driving force for osteogenesis of skeletal progenitors, governing musculoskeletal integrity.

Multimodal cell atlas of the ageing human skeletal muscle

Article Open access 22 April 2024

Human skeletal muscle aging atlas

Article Open access 15 April 2024

Fibroblast and myofibroblast activation in normal tissue repair and fibrosis

Article 08 April 2024

Introduction

Mitochondrial diseases are associated with decreased energy production in organs with high energy requirements and are mainly due to defects in the oxidative phosphorylation (OXPHOS) machinery^1,2,3,4. Releasing energy through OXPHOS requires the orchestrated action of five multi-heteromeric enzyme complexes (complexes I–V), located in the inner mitochondrial membrane⁵. Mitochondrial complex 1 (CI) is the largest complex comprised of at least 45 different proteins and is the primary entry point that generates the proton motive force for electrons into the electron transport chain driving ATP production^6,7,8. CI-deficiency is the most frequent cause of defects of the OXPHOS system that cause multisystem dysfunction^{9,10,11,12,13}, including skeletal disorders^14,15.

Emerging epidemiological evidence^16,17,18 and data from animal models^3,19 suggest an association of mitochondrial dysfunction and high oxidative stress with skeletal disorders, including a clinically significant decrease in bone strength, reduced bone mineral density, premature aging of bone or a high risk or fragility fractures. For example, Kearns Sayre Syndrome (KSS), a rare mitochondrial DNA deletion syndrome, shows clinical manifestations of muscle weakness and wasting, accompanied with multiple bone deformities and severe fractures, causing loss of the ability to walk and early lethality²⁰. In addition, mitochondrial respiratory chain deficiency is related to spontaneous bone fractures and myopathy²¹. However, it is unclear how mitochondrial defects contribute to skeletal dysfunction associated with myopathy.

Glycolysis is the primary metabolic process that generates ATP for energy and intermediate metabolites, including active lipids, nucleotides, and amino acids, for biosynthesis^22,23. It serves as the major energy and carbon source for committed osteoblast-lineage cells during osteogenic differentiation^22,24,25,26. In particular, metabolic intermediates fuel the synthesis of amino acids to produce extracellular matrix proteins and hydroxyapatite in bone²². In addition to glycolysis, glutamine metabolism also provides osteoblast progenitors with an additional energy source as well as a carbon and nitrogen source²⁶. Glutamine is converted to citrate through the tricarboxylic acid (TCA) cycle in the mitochondria, producing ATP through the OXPHOS pathway, while providing metabolic intermediates for the synthesis of amino acids, nucleotides, glutathione, and hexoamine²⁷. Recent studies demonstrated an in vivo role for glutamine and glucose metabolism in skeletal progenitors using conditional deletion of the enzyme glutaminase (Gls) and the glucose transporter Glut1, respectively. Gls-deficiency in Prx1⁺ skeletal progenitors impaired,cell proliferation and osteogenesis during bone remodeling while Glut1-deficiency in Prx1⁺ skeletal progenitors abolished chondrogenesis, not osteogenesis, during skeletal development²⁸. Notably, glutamine is likely to serve as major substrates for biosynthetic and antioxidant demands, not an energetic substrate, in skeletal progenitors since GLS inhibition did not induce any energetic alteration²⁹. Thus, other energy metabolic pathways controlling osteogenesis of skeletal progenitors remain to be elucidated.

To investigate in vivo roles for mitochondrial oxidative metabolism in musculoskeletal system, we developed a mouse model for mitochondrial CI deficiency by conditionally deleting Ecsit in multiple skeletal and muscle cells. Herein, we highlight the importance of mitochondrial OXPHOS in determining the osteogenic commitment and differentiation of skeletal progenitors as a critical rheostat of skeletal and muscular development. Notably, restoration of impaired OXPHOS pathway using recombinant adeno-associated virus (rAAV) may represent an untapped therapeutic avenue for mitochondrial disorders in the musculoskeletal system.

Results

Ecsit-deficiency in skeletal progenitors impairs early skeletal development

The CI intermediate assembly (MCIA) complex, containing the core subunits (ECSIT, NDUFAF1, ACAD9, and TMEM126B) is required for the formation of mitochondrial complex I (CI)^{6,30,31,32,33}. Among these subunits, ECSIT is highly expressed in skeletal cells, including endosteal osteoblasts, osteocytes, periosteal cells, and osteoclasts (Fig. 1a), and its deficiency reduces the expression of other MCIA and CI components^31,32,33. Mitochondrial localization of ECSIT in human bone marrow-derived stromal cells (BMSCs) and mouse osteoblasts residing on the surface of trabecular and cortical bones was also confirmed by immunofluorescence and subcellular fractionation analyses (Fig. 1b, c). The role of mitochondrial activity in bone formation was examined in vivo using mice with targeted deletion of Ecsit in skeletal progenitors (Ecsit^Prx1). Deletion efficiency of Ecsit in these cells was confirmed by fluorescence microscopy and immunoblotting analyses (Supplemental Fig. 1a, Fig. 2c). Ecsit^Prx1 mice developed limb deformities along with spontaneous fractures in long bones right after birth (Fig. 1d, right). At postnatal day 21 (P21), they displayed severe osteoporosis and delayed formation of primary and secondary ossification centers in long bones (Fig. 1d, e, Supplemental Fig. 1b, c). Severe hypomineralization and delayed formation of primary ossification centers was also evident in Ecsit^Prx1 embryos and neonates (Fig. 1f, Supplemental Fig. 1e, f, Supplemental Table 1). Consistent with immunofluorescence analysis showing nearly complete absence of osteocalcin (OCN)-expressing mature osteoblasts in the femur of P21 Ecsit^Prx1 mice (Fig. 1f, right), cKO Prx1⁺ skeletal progenitors, isolated from P0 Ecsit^Prx1 hindlimbs as Prx1⁺CD45^-TIE2^-TER119^- cells, failed to differentiate into mature osteoblasts (Fig. 1g). Of note, the bone marrow niche function of Ecsit^Prx1 skeletal cells was also disrupted, as shown by a significant reduction in total bone marrow cell numbers and a low frequency of hematopoietic stem cells (HSCs), B cells, macrophages, granulocytes, and T cells (Fig. 1h, Supplemental Fig. 2). These results suggest that ECSIT is essential for osteogenesis of Prx1⁺ skeletal progenitors and for the formation of bone marrow HSC niche. Notably, GLUT1-mediated glucose metabolism mediates chondrogenesis, not osteogenesis of Prx1⁺ skeletal progenitors²⁸ while GLS-mediated glutamine metabolism in Prx1⁺ skeletal progenitors is dispensable for skeletal development²⁹. Thus, ECSIT-mediated mitochondrial regulation, not GLUT1-mediated glucose metabolism or GLS-mediated glutamine metabolism, plays a critical role in the osteogenic activity of Prx1-lineage skeletal progenitors.

**Fig. 1: *Ecsit*^Prx1 mice show defects in skeletogenesis, fracture healing, and the bone marrow niche.**

Prx1-lineage skeletal cells are major cell populations (Fig. 1i) that initiate callus formation and vascularization for successful bone fracture repair^34,35,36. Ecsit deletion in Prx1⁺ skeletal progenitors resulted in nearly complete arrest of cortical bone and periosteal formation despite expansion of Ecist-deficient periosteal cells (Fig. 1j, Supplemental Fig. 1f). Accordingly, Ecsit^Prx1 femurs show spontaneous fractures occurring after birth. The resulting fracture sites displayed persistent nonunion with a lack of callus mineralization and mature osteoblasts and the persistence of an unmineralized fibrous callus with cartilaginous elements (Fig. 1j, k, Supplemental Fig. 1d). To define where ECSIT is required in the osteoblast maturation process, skeletal cell populations were isolated from P0 Ecsit^Prx1 limbs and subjected to flow cytometry³⁷, demonstrating an increased frequency of skeletal stem cells (SSCs) and pre-bone, cartilage, stromal progenitors (pre-BCSPs) and a decreased frequency of BCSPs (Fig. 1l, m, Supplemental Fig. 3). Thus, Ecsit deletion in Prx1⁺ skeletal progenitors impairs the differentiation of SSCs and pre-BCSPs to BCSPs, resulting in defective commitment and differentiation of skeletal progenitors to osteoblasts.

To test whether Ecsit deletion in skeletal progenitors is a primary cause of impaired fracturing healing in Ecsit^Prx1 mice, Ecsit^fl/fl mice were crossed with Prx1-Cre^ERT2-EGFP mice (Ecsit^Prx1-ERT/GFP) to delete ECSIT expression in skeletal progenitors in an inducible manner. Fracture surgery was performed on the femurs of 6-week-old Ecsit^Prx1-ERT/GFP and Prx1-ERT/GFP control mice three days after tamoxifen injection (Supplemental Fig. 4a). One month later, tamoxifen-induced expression of Cre recombinase and deletion of Ecsit were validated by fluorescence microscopy and RT-PCR, respectively (Supplemental Fig. 4b, c). MicroCT analysis revealed that basal bone mass in control and Ecsit^Prx1-ERT/GFP femurs was comparable, as shown by equivalent trabecular bone mass and cortical bone thickness of non-fractured femurs (Supplemental Fig. 4d, e). Unlike Ecsit^Prx1 mice, tamoxifen-treated Ecsit^Prx1-ERT/GFP mice displayed normal periosteal and callus formation and fracture unionization in the fractured sites (Supplemental Fig. 4b, f, g). These results demonstrate that inducible deletion of Ecsit in Prx1⁺ skeletal progenitors in adult mice does not affect bone fracture healing and homeostasis. Thus, impaired fracture healing seen in Ecsit^Prx1 mice may result from pre-existing skeletal phenotypes and/or alteration of skeletal progenitor population during early skeletal development. However, it cannot be fully excluded that inadequately complete deletion of Ecsit in this system contributes to the lack of a phenotype observed. Taken together, ECSIT-mediated mitochondrial functions play critical roles in osteogenic commitment and differentiation of skeletal progenitors during early skeletal development.

ECSIT controls mitochondrial OXPHOS in skeletal progenitors

Despite the original identification of ECSIT as a key regulator of NF-κB signaling downstream of pro-inflammatory receptors, Toll-like receptor (TLR) and interleukin (IL)-1β^38,39, transcriptome analysis showed an enrichment of genes associated with metabolic pathways and mitochondrial function, not the NF-κB pathway in cKO Prx1⁺ skeletal cells (Fig. 2a, b). Specifically, expression of the genes involved in mitochondrial CI biogenesis and OXPHOS was markedly downregulated in the absence of Ecsit while there was little to no altered expression in glucose and glutamine metabolism (Fig. 2b). This is consistent with previous reports demonstrating ECSIT as a key component of the MCIA complex that controls the stability of CI subunits^31,32,33. To directly test the role of ECSIT in mitochondrial OXPHOS, the expression of mitochondrial CI, II, III, IV, and V subunit proteins in the OXPHOS pathway was examined in WT and cKO Prx1⁺ skeletal progenitors, demonstrating that expression of core CI subunit proteins, including ND6, NDUFB8 and NDUFS3, was ablated in the absence of Ecsit while Ecsit deletion minimally affected the expression of proteins in other mitochondrial complexes (Fig. 2c, left). Similarly, Ecsit-deficiency significantly reduced expression of NDUFAF1, a key component of the MCIA complex^32,33 (Fig. 2c, right), suggesting the importance of ECSIT in stabilizing the components of CI and MICA complexes in Prx1⁺ skeletal progenitors. Accompanying our transcriptome analysis that show little to no enrichment of genes associated with mitochondrial biogenesis (Fig. 2b), mRNA and protein levels of Atf5^40,41 and Hsp90^42,43 were largely intact in the absence of Ecsit (Fig. 2c, Supplemental Fig. 5a). Of note, expression of Lc3b and Pink1⁴⁴, key regulators of mitophagy, was markedly increased in cKO Prx1⁺ skeletal progenitors (Fig. 2c), suggesting enhanced mitophagy activity. These results suggest that Ecsit deletion in Prx1⁺ skeletal progenitors destabilizes the components of the MCIA and CI complexes while enhancing mitophagy activity. Accordingly, cKO Prx1⁺ skeletal progenitors displayed a significant decrease in mitochondrial numbers and transmembrane potential (Fig. 2d), in the ratio of mitochondrial to nuclear DNA levels (Fig. 2e), and in ATP production (Fig. 2f). Along with a decrease in basal and maximal mitochondrial oxygen consumption (Fig. 2g), cKO Prx1⁺ skeletal progenitors failed to respond to mitochondrial respiratory modulators, demonstrating the nearly complete absence of oxygen consumption rates (OCR, Fig. 2h). Of note, the glycolytic proton efflux rate was markedly upregulated in these cells, as shown by greater extracellular acidification rate (ECAR, Fig. 2i). However, enhanced glycolysis was incapable of compensating for the reduced ATP production by the damaged mitochondrial OXPHOS. In contrast to the involvement of glucose metabolism in RUNX2-mediated osteogenesis⁴⁵, expression and transcription activity of RUNX2 were largely intact in cKO Prx1⁺ skeletal progenitors (Supplemental Fig. 5b, c), suggesting that ECSIT-mediated mitochondrial oxidative metabolism is dispensable for RUNX2-mediated osteogenesis.

**Fig. 2: *Ecsit*-deficiency in skeletal progenitors impairs mitochondrial OXPHOS.**

Since ECSIT is localized on the membrane of mitochondria via an N-terminal mitochondria localization sequence (MLS, 1–48 aa)^32,46, to investigate mechanistic actions of ECSIT in regulating mitochondrial OXPHOS, cKO Prx1⁺ skeletal progenitors were reconstituted with FLAG-tagged ECSIT full length (ECSIT-FL) or a MLS-deletion mutant (ECSIT-ΔMLS) that failed to localize to mitochondria via rAAV-mediated delivery (Fig. 2j). Enforced expression of ECSIT-FL, not the ECSIT-ΔMLS mutant, reversed alterations in expression of the CI subunit NDUFS3, mitochondrial numbers and transmembrane activity, and osteogenic potential in cKO Prx1⁺ skeletal progenitors (Fig. 2k–m). Thus, mitochondrial ECSIT regulates CI assembly and mitochondrial OXPHOS in Prx1⁺ skeletal progenitors, which in turn determines the osteogenic commitment and differentiation of skeletal progenitors.

To confirm these findings in vivo, a bone-targeting rAAV was used to express FLAG-tagged ECSIT-FL and ECSIT-ΔMLS proteins in Ecsit^Prx1 mice (Fig. 3a). Given our previous report showing high transduction efficiency of the AAV9 serotype in osteoblast-lineage cells in vivo⁴⁷, systemic delivery of AAV9 vector expressing mCherry to P1 Prx1;Rosa^mT/mG neonates was effective in transducing GFP-expressing Prx1⁺ osteoblast-lineage cells, including skeletal progenitors, osteoprogenitors, and mature osteoblasts, in both trabecular and cortical compartments of long bones in addition to lung, heart, liver, kidney, and skeletal muscle (Fig. 3b, c). AAV-mediated expression of FLAG-tagged ECSIT-FL and ECSIT-ΔMLS proteins in the femur and liver was confirmed (Fig. 3f, Supplemental Fig. 6). Remarkably, enforced expression of ECSIT-FL, but not a vector control or ECSIT-ΔMLS, markedly improved the survival rates of Ecsit^Prx1 mice (Fig. 3d). The osteogenic activity of Ecsit^Prx1 skeletal progenitors was recovered by ECSIT-FL expression (Fig. 3f–h), leading to recovery of bone formation and regeneration activity at multiple anatomic compartments. This led to an almost complete rescue of skeletal deformities in Ecsit^Prx1 forelimbs and hindlimbs (Fig. 3e) and abnormal bone marrow compartment (Fig. 3f, left) and impaired fracture healing and periosteal bone formation of Ecsit^Prx1 long bones (Fig. 3g, h). However, none of these phenotypes were improved by AAV-mediated expression of ECSIT-ΔMLS (Fig. 3d–h). As expression of ECSIT-FL, not ECSIT-ΔMLS, restored OCN⁺ mature osteoblasts and osteoblast-supporting CD31⁺EDMC⁺ type H endothelium⁴⁸ in the fracture callus of Ecsit^Prx1 long bones, fracture healing progressed normally in these mice (Fig. 3g, h). Restoring ECSIT levels in Prx1⁺ osteoblast-lineage cells, including skeletal progenitors, osteoprogenitors, and mature osteoblasts, reverses skeletal defects and retrieves bone regeneration during fracture in Ecsit^Prx1 mice. Thus, ECSIT-mediated regulation of mitochondrial OXPHOS in skeletal progenitors is crucial for bone accrual, fracture repair, and bone marrow niche maintenance during early skeletal development.

**Fig. 3: AAV-mediated expression of ECSIT reverses *Ecsit*^Prx1 skeletal phenotypes.**

ECSIT is dispensable of mitochondrial OXPHOS in committed skeletal cells

Since ECSIT is highly expressed in osterix (OSX)⁺ committed osteoprogenitors and dentin matrix acidic phosphoprotein 1 (DMP1)⁺ mature osteoblasts and osteocytes (Figs. 1a and 4a), the role of ECSIT in these cells were examined using conditional deletion with OSX-cre (Ecsit^Osx) or DMP1-cre (Ecsit^Dmp1) mice. mRNA and protein levels of Ecsit in these mice were examined by RT-PCR, immunoblot, and immunohistochemistry analyses, demonstrating lack of Ecsit expression (Fig. 4b, d, g). Remarkably, unlike Ecsit^Prx1 mice, little to no gross skeletal phenotypes were seen in two-month-old Ecsit^Osx and Ecsit^Dmp1 mice. Bone morphology, trabecular bone mass, and cortical bone thickness were all normal in the long bones of these mice (Fig. 4c, e, Supplemental Fig. 7). In addition, osteogenic differentiation and bone fracture healing were largely intact in the absence of Ecsit (Fig. 4f, Supplemental Fig. 8). These results suggest that ECSIT function is dispensable for the development of committed osteoblast-lineage cells, including Osx⁺ osteoprogenitors and Dmp1⁺ mature osteoblasts and osteocytes.

**Fig. 4: *Ecsit*^Osx and *Ecsit*^Dmp1 mice show no skeletal phenotype and intact mitochondrial function.**

To gain insights on this mechanism, mitochondrial OXPHOS was examined in Ecsit^Osx osteoprogenitors and Ecsit^Dmp1 mature osteoblasts. Unlike Ecsit^prx1 skeletal progenitors, expression of the CI protein NDUFS3, ND6, and NDUFB8, and the MCIA subunit NDUFAF1 (Fig. 4g), the ratio of mitochondrial DNA to nuclear DNA (Fig. 4h), mitochondrial numbers and transmembrane potential (Fig. 4i) were all intact in these cells, suggesting that ECSIT function is dispensable for mitochondrial OXPHOS in Osx⁺ committed osteoprogenitors and Dmp1⁺ mature osteoblasts and osteocytes. Notably, while a similar decrease in ECSIT expression was observed among Ecsit^Prx1 skeletal progenitors, Ecsit^Osx osteoprogenitors, and Ecsit^Dmp1 osteoblasts, only Ecsit^Prx1 cells showed a significant decrease in NDUFS3 expression (Supplemental Fig. 9). However, it cannot be fully excluded that inadequately deletion of Ecsit in osteoblasts may contribute to the lack of an observed phenotype. Thus, the role of ECSIT in controlling mitochondrial oxidative metabolism is context- and tissue-dependent.

Ecsit deficiency in skeletal progenitors induces muscle atrophy

Accompanying the onset spontaneous fractures in Ecsit^Prx1 mice is an onset of myopathy and muscle atrophy. Specifically, the size and weight of skeletal muscles, including tibialis anterior (TA) and gastrocnemius (GA) muscles, were dramatically decreased (Fig. 5a, b, Supplemental Fig. 10a). While muscle structure of Ecsit^Prx1 neonates is grossly normal after birth (Fig. 5c, Supplemental Fig. 10b), Ecsit^Prx1 mice gradually developed severe degeneration of skeletal muscle, as evident by significant decrease in mature myofibers, extracellular matrix in the interstitial space, and the expression of mature myoblast genes (Fig. 5d, e, Supplemental Fig. 10c). Despite a substantial decrease in total myocyte numbers in Ecsit^Prx1 mice, the frequency of skeletal muscle-resident satellite cells⁴⁹ (β-integrin⁺CXCR4⁺Sca1^-CD45^-CD11b^-TER119^-) and mesenchymal progenitors⁵⁰ (Sca1⁺integrin α7^-CD45^-CD31^-) was dramatically increased in Ecsit^Prx1 muscle (Fig. 5f, Supplemental Fig. 11), indicating a defect in myoblast differentiation, but not the maintenance of muscle-resident stem cells. Of note, this phenotype does not result from the intrinsic effects of Ecsit deficiency in muscle cells, as shown by intact expression of Ecsit mRNA and protein in Ecsit^Prx1 muscle (Fig. 5d, e). This is accompanied by fluorescence microscopy showing little to no Prx1-cre driven reporter activity in skeletal muscle cells of Prx1-cre;Rosa26^mT/mG mice (Fig. 5g). To directly test this, mice with targeted deletion of Ecsit in myoblast precursors (Ecsit^Myf5) were generated, demonstrating little to no gross phenotypes in the skeletal muscle, bone, and bone marrow (Fig. 5h). Likewise, the expression of mature myoblast genes in the skeletal muscle of these mice was also normal (Fig. 5i). Additionally, primary human myoblast progenitors (17Ubic)⁵¹ with Ecsit knockdown were cultured under myogenic conditions, demonstrating no effect of Ecsit-deficiency on myoblast differentiation (Fig. 5j). As seen in Osx⁺ committed osteoprogenitors and Dmp1⁺ mature osteoblasts (Fig. 4g), deletion of Ecsit in myoblast progenitors did not affect expression of CI protein NDUFS3, the ratio of mitochondrial DNA to nuclear DNA, mitochondrial numbers and transmembrane potential while these cells showed a mild reduction in the MCIA subunit Ndufaf1 (Fig. 5k, l). These results suggest that ECSIT function is dispensable for mitochondrial OXPHOS for ATP production during myoblast differentiation. Thus, Ecsit deletion in Prx1⁺ skeletal progenitors, not muscle cells, is responsible for muscle atrophy in Ecsit^Prx1 mice.

**Fig. 5: *Ecsit*^Prx1 mice display skeletal muscle atrophy.**

To determine how skeletal fracture could lead to muscle atrophy, the migration of GFP-expressing Ecsit^Prx1 skeletal progenitors into skeletal muscle was monitored at fracture sites. Wild-type mice undergoing fracture did not experience muscle degeneration, showing normal callus formation at 7- and 14-days post-fracture, and restriction of wild-type Prx1-lineage skeletal cells to the callus and not adjacent skeletal muscle (Fig. 6a). On the other hand, Ecsit^Prx1 mice showed abnormal callus development and fibrotic tissue at fracture sites and failed to form periosteum, where Ecsit^Prx1 skeletal progenitors migrated into the skeletal muscle and induced degeneration of myofibers (Fig. 6b). In non-fractured long bone showing normal muscle structure, Ecsit^Prx1 skeletal progenitors were only detected inside the bone (Fig. 6c). Thus, impaired fracture healing and periosteal formation allow migration of Ecsit^Prx1 skeletal progenitors to adjacent skeletal muscle, resulting in muscle atrophy. To directly test this hypothesis, conditioned media (CM) harvested from WT or cKO Prx1⁺ skeletal progenitors were added to cultures of human myoblast progenitors (17Ubic). These studies demonstrated a significant reduction in myoblast differentiation in the presence of cKO-CM compared to WT-CM (Fig. 6d). Transcriptome and RT-PCR analyses of FACS-sorted WT and cKO Prx1⁺ skeletal progenitors and SSCs revealed elevated levels of the soluble myogenic suppressor TGF-β1^52,53 in cKO cells (Fig. 6e, f). Since TGF-β1 expression was unchanged in satellite cells (MuSCs) isolated from the skeletal muscle of Ecsit^Prx1 mice and Ecsit-deficient myoblast progenitors (Fig. 6f, g), cKO Prx1⁺ skeletal progenitors and SSCs, not MuSCs and myoblast progenitors, are likely to be major contributors of TGF-β1 expression in skeletal muscle. Accordingly, the TGF-β1 inhibitor SB-431542 almost completely reversed the ability of cKO-CM to suppress myoblast differentiation (Fig. 6d). Since the TGF-β pathway plays a role in the pathogenesis of skeletal muscle myopathies and its inhibition is considered for the treatment of skeletal muscle fibrosis⁵⁴, TGF-β1 is likely to be an important factor mediating the myoblast suppressive activity of Ecsit^Prx1 skeletal progenitors, though other pathways may contribute as well.

**Fig. 6: *Ecsit*^Prx1 skeletal progenitors induce skeletal muscle degeneration.**

Intravenous administration of P1 Ecsit^Prx1 neonates with bone-targeting rAAV9 expressing FLAG-tagged ECSIT-FL, not ECSIT-ΔMLS mutant, almost completely recovered degeneration of skeletal muscle (Fig. 6h). This was accompanied with little to no migration of Ecsit^Prx1 skeletal progenitors expressing ECSIT-FL, not ECSIT-ΔMLS mutant, to the skeletal muscle (Fig. 6i). Finally, enforced expression of ECSIT-FL markedly decreased TGF-β transcription in FACS-sorted cKO Prx1⁺ skeletal progenitors while TGF-β expression was partially reversed by ECSIT-ΔMLS mutant (Fig. 6j), suggesting the presence of additional functional domain(s) within Ecsit regulating TGF-β expression. Thus, impaired OXPHOS in skeletal progenitors is likely to be responsible for muscle atrophy seen in Ecsit^Prx1 mice. Taken together, ECSIT-mediated regulation of mitochondrial oxidative metabolism is a key determinant of the osteogenic commitment and differentiation of Prx1⁺ skeletal progenitors, governing musculoskeletal development and bone marrow niche functions.

Discussion

Our study provides in vivo evidence showing the importance of mitochondrial oxidative metabolism in determining osteogenic commitment and differentiation of skeletal progenitors during early skeletal development. Despite emerging evidence suggesting an association of mitochondrial dysfunction and high oxidative stress in skeletal cells with skeletal disorders¹⁷, only few genetic models with mitochondrial OXPHOS defects have been successfully generated due to either embryonic lethality or the absence of any relevant phenotypes. Herein, we generated a mouse model with CI deficiency by conditionally deleting the key CI assembly factor Ecsit in various musculoskeletal cells. Deletion of Ecsit in Prx1⁺ skeletal progenitors disrupts mitochondrial OXPHOS for ATP production while enhancing glucose metabolism as a compensatory mechanism. Lack of OXPHOS activity ablates the osteogenic potential of skeletal progenitors, resulting in defects in skeletal development, bone marrow niche maintenance, and bone fracture repair. Surprisingly, breakdown of this regulatory circuit leads to the infiltration of Prx1⁺ skeletal progenitors into adjacent muscle, which induces myopathy mediated by TGF-β1 (Supplemental Fig. 12). This unexpected finding indicates that OXPHOS defects or TGF-β1 may be implicated in paired skeletal fragility/myopathy disorders, and there may be potential therapeutic targets in this context. Intriguingly, despite the requirement of mitochondrial localization of Ecsit in regulating mitochondrial oxidative metabolism in Prx1⁺ skeletal progenitors, our findings that cytosolic Ecsit (Ecsit-ΔMLS) can partially improve survival of Ecsit^Prx1 mice and TGF-β1 expression in Ecsit^Prx1 skeletal progenitors suggest the presence of additional functional domain(s) within Ecsit regulating survival and TGF-β1 expression. Further studies will be necessary to identify these functional domains in Prx1⁺ skeletal progenitors.

Glucose and glutamine metabolism and mitochondrial OXPHOS have been known as major pathways of energy metabolism regulating the production of ATP in osteoblast-lineage cells²⁴. In particular, glucose and glutamine metabolism are prominent metabolic features of committed osteoblast-lineage cells^55,56,57,58 and integral to collagen production and mineralizing activity by generating ATP and intermediate metabolites^59,60. In contrast to our findings that ECSIT-mediated mitochondrial OXPHOS provides an energy source crucial to the osteogenic commitment and differentiation of Prx1⁺ skeletal progenitors during skeletal development, glutamine metabolism mainly contributed to biosynthesis and redox homeostasis, but not energy production, during bone remodeling. Of note, since amino-acid-transaminase-derived αKG is the critical downstream glutamine metabolite regulating proliferation of skeletal progenitors²⁹, cell proliferation rates were largely not altered in Ecsi^Prx1 skeletal progenitors with normal glutamine metabolism. Intriguingly, despite extensive studies showing that glucose is a major nutrient for osteoblast development, mice lacking the glucose transporter Glut1 in Prx1⁺ skeletal progenitors displayed impaired cartilage development without any alteration in skeletal mineralization²⁸. This suggests that glycolysis is critical for chondrogenesis, not osteogenesis, of Prx1⁺ skeletal progenitors during skeletal development. Notably, Ecsit-deficiency in Prx1⁺ skeletal progenitors enhanced glucose metabolism as a compensatory mechanism while cartilage development occurred normally in Ecsi^Prx1 mice. Thus, our data implies that mitochondrial OXPHOS, not glucose or glutamine metabolism, is the major energy-driving force for osteogenesis in skeletal progenitors and that it is, therefore, critical for skeletal development and the maintenance of musculoskeletal integrity. We note that reduced mitochondrial mass in Ecsi^Prx1 skeletal progenitors also potentially results in a broader set of cellular abnormalities than just decreased ATP synthesis, as many TCA cycle metabolites, such as acetyl-CoA and α-ketoglutaric acid, have direct roles in lipid synthesis, protein acetylation, histone demethylation, and other processes^29,61,62.

Since ECSIT function is dispensable for Osx⁺ committed osteoprogenitors and Dmp1⁺ mature osteoblasts and terminally differentiated osteocytes, glycolysis may offset this as glycolysis is the primary metabolic pathway for ATP production during the differentiation of committed osteoblast-lineage cells⁵⁹. Further studies will be necessary to determine how ECSIT controls mitochondrial oxidative metabolism and osteogenic development in a context and tissue-dependent manner. Additional further studies are warranted to define the precise mechanisms by which mitochondrial OXPHOS impacts skeletal progenitors through other mitochondrial regulatory mechanisms, such as mitochondrial biogenesis, the degradation of damaged mitochondria via mitophagy, reactive oxygen species (ROS) generation, and mtDNA replication and/or transcription.

The unexpected finding that mitochondrial defects lead Ecsit^Prx1 skeletal progenitors to infiltrate into muscle and mediate pathogenic TGF-β1 secretion, provides a new potential explanation for a number of clinical disorders that pair prominent or primary skeletal phenotypes with myopathy. For instance, Kearns Sayre Syndrome displays linked skeletal and myopathy phenotypes due to mitochondrial dysfunction, raising the possibility that the myopathy may be in part secondary to pathogenic infiltration of skeletal progenitors into adjacent muscle²⁰. rAAVs have been considered the most promising viral vector for gene therapy⁶³, as they have a clinical track record of being well tolerated and have the potential to efficiently mediate genetic modifications that can persist for years after a single treatment. Remarkably, AAV-mediated expression of full length Ecsit almost completely restored impaired mitochondrial OXPHOS and osteogenic potentials of Ecsit^Prx1 progenitors, preventing musculoskeletal deformities of Ecsit^Prx1 mice. Thus, our findings provide proof-of-concept demonstration for gene therapy of developmental skeletal disorders. Collectively, this study reveals a previously unappreciated role of mitochondrial oxidative metabolism in determining osteogenesis of skeletal progenitors and demonstrates a new gene therapy approach for mitochondrial diseases in the musculoskeletal system.

Methods

Mice

Mice were housed in a constant environment under a standard mouse chow diet (up to 5 mice per cage); ambient temperature of 21 ± 2 °C, circulating air, and constant humidity of 50 ± 10%, in a 12 h light, 12 h dark cycle. Mice were monitored every three days for the amount of their food and water intake and signs for distress. For signs of severe distress, including general malaise, severe cachexia, or more than 20% loss of body weight, humane euthanasia was performed in consultation with veterinary staff. Mice were euthanized in a carbon dioxide chamber, followed by cervical dislocation. Neonates were euthanized by decapitation. Ecsit^fl/fl mice were generated as previously described and maintained on a C57BL/6J background⁴⁶. Cre deleter mice (C57BL/6J) that express Cre recombinase under the control of the prx1 promoter (Prx1-cre), the osterix promoter (Osx-cre), and the dmp1 promoter (Dmp1-cre) were purchased from The Jackson Laboratory. To label Cre-expressing cells, Ecsit;^fl/flPrx1-cre mice were further crossed with Rosa^mT/mG cre reporter mice (C57BL/6J)⁶⁴. Mouse genotypes were determined by PCR on tail genomic DNA and primer sequences are available upon request. Control littermates were used and analyzed in all experiments. A phenotypic summary of Ecsit^Prx1 mice is described in Supplemental Table 1. All animals were used in accordance with the NIH Guide for the Care and Use of Laboratory Animals and were handled according to protocols approved by the University of Massachusetts Medical School Institutional Animal Care and Use Committee (IACUC).

Generation of rAAV vectors

DNA sequences for ECSIT-FL and ECSIT-ΔMSL were synthesized as gBlocks, cloned into the intronic region of the pAAVsc-CB6-Egfp plasmid at the restriction enzyme sites (PstI and BglII)⁶⁵, and packaged into an AAV9 capsid. rAAV production was performed by transient transfection of HEK293 cells, purified by CsCl sedimentation, titered by droplet digital PCR (ddPCR) on a QX200 ddPCR system (Bio-Rad) using the Egfp prime/probe set as previously described^47,65. The sequences of gBlocks and oligonucleotides for ddPCR and are listed in Supplemental Table 2.

MicroCT, radiography, and skeletal preparation

MicroCT was used for qualitative and quantitative assessment of trabecular and cortical bone microarchitecture and performed by an investigator blinded to the genotypes of the animals under analysis. Femurs excised from the indicated mice were scanned using a microCT 35 (Scanco Medical) with a spatial resolution of 7 μm. For trabecular bone analysis of the distal femur, an upper 2.1 mm region beginning 280 μm proximal to the growth plate was contoured. For cortical bone analysis of the femur and tibia, a midshaft region of 0.6 mm in length was used. MicroCT scans of forelimbs and hindlimbs were performed using isotropic voxel sizes of 12 μm. 3D reconstruction images were obtained from contoured 2D images by methods based on distance transformation of the binarized images. All images presented are representative of the respective genotypes (n = 5–8).

X-radiographic images of whole body and hindlimbs were taken by the Trident Specimen Radiography System (Hologic, Inc.).

Skeletons were prepared for gross morphology analysis using the McLeod method, as previously described⁶⁶. Briefly, mice were sacrificed, skinned, eviscerated, and fixed in 95% ethanol for one day. Then, skeletons were stained by alizarin red s and alcian blue solutions (Sigma, A3157) and sequentially cleared in 1% potassium hydroxide. All images presented are representative of the respective genotypes (n = 3–5).

Histology, immunofluorescence, and immunohistochemistry

Histological analysis was performed on liver and hindlimbs, and skeletons and skeletal muscles were dissected from mice fixed in 4% PFA for one day, and decalcified by daily changes of 15% tetrasodium EDTA for one to two weeks. Tissues were dehydrated by passage through an ethanol series, cleared twice in xylene, embedded in paraffin, and sectioned at 5 μm thickness along the coronal plate from anterior to posterior. Decalcified femoral sections were stained with hematoxylin and eosin (H&E) or safranin O. Alternatively, intact hindlimbs containing bones and skeletal muscles were immediately fixed in ice-cold 4% paraformaldehyde solution for one day for histological analysis of skeletal muscles. Semi-decalcification was carried out for five days in 0.5 M EDTA pH 7.4 at 4 °C with constant shaking (age ≥1 week), and infiltration was followed with a mixture of 20% sucrose phosphate buffer for one day. All samples were embedded in a 50/50 mixture of 25% sucrose solution and OCT compound (Sakura) and cut into 10 μm thick cross sections using a cryostat (Leica Microsystems, USA).

Immunofluorescence staining and analysis were performed, as described previously⁸. Briefly, after treatment with 0.2% Triton X-100 for 10 min, cryosections were blocked with 5% donkey serum at room temperature for 30 min and incubated overnight at 4 °C with primary antibody and then visualized with fluorescence-conjugated secondary antibody (1:400, Molecular Probes). Nuclei were counterstained with 4–6,diamidino-2-phenylindole (DAPI). An Olympus IX81 confocal microscope or Leica TCS SP5 II Zeiss LSM-880 confocal microscope or EVOS epifluorescence microscope was used to image samples. Osteocalcin (OCN, sc-365797, Santa Cruz Biotechnology, 1:100), ECSIT (HPA042979, Sigma,1:100), ECSIT (gifted from Dr. Sankar Gosh, Columbia University), NDUFS3 (43–9200, Invitrogen,1:100), CD31/PECAM-1 (AF3628, R&D systems, 1:50), Endomucin (sc-65495, Santa Cruz Biotechnology, 1:100), COL1A1 (A1352, ABclonal, 1:50), MF20 (AB_2147781, DSHB, 1:100) and anti-FLAG-tag (MA1-142-A488, Sigma-Aldrich, 1:100) were used for primary antibodies.

For immunohistochemistry, paraffin sections were dewaxed and stained according to the manufacturer’s directions, using the Discovery XT automated IHC stainer (Ventana Medical Systems, Inc., Tucson, AZ, USA). CC1 standard (pH 8.4 buffer contained Tris/Borate/EDTA) and inhibitor D (3% H₂O₂, Endogenous peroxidase) were used for antigen retrieval and blocking, respectively. Sections were incubated with antibodies specific to ECSIT (A7804, ABclonal,1:200) for 40 min at 37 °C, and a secondary antibody for 20 min at 37 °C. Subsequently, they were incubated in SA-HRP D for 16 min at 37 °C and then DAB + H₂O₂ substrate for 8 min, followed by hematoxylin and bluing reagent counterstain at 37 °C. Reaction buffer (pH 7.6 Tris buffer) was used as washing solution. Stained samples were visualized using an Aperio virtual microscope (Leica Microsystems, USA) and images of the sample were analyzed by the Aperio image scope program (ver. 12.3.2.8013, Leica Microsystems, USA).

Human myoblast and mouse osteoblast and osteoclast differentiation

Human primary myoblasts 17Ubic were kindly gifted from Dr. Charles P Emerson at UMass Chan Medical school. Informed consent was obtained from the patients who donated tissue for the production of cell lines. IRB protocols were approved by UMass Chan Medical School (H00006581-10 and H00006581-11) and by Kennedy Krieger Institute (B0410080117). For human myoblast differentiation, human primary myoblasts 17Ubic⁵¹ were cultured in growth medium (20% FBS,0.5% CEE and 1.2 mM CaCl₂ in Ham’s F10 medium) and plated into a collagen-coated 100 mm Petri-dish. To induce myogenic differentiation and fusion, cells were cultured in growth medium until >95% confluency and changed with differentiation medium (OPTI-MEM). Knockdown of ECSIT expression in 17Ubic cells was performed using lentivirus-mediated delivery of human ECSIT shRNA hairpins and knockdown efficiency was confirmed by immunoblotting. Alternatively, conditioned medium was harvested from cultured Ecsit;^+/+Prx1-Cre;Rosa^mT/mG or Ecsit;^fl/flPrx1-Cre;Rosa^mT/mG skeletal cells and added to 17Ubic cells. 6 days after myogenic culture (30% conditioned medium + 70% myogenic differentiation medium), myoblast differentiation was analyzed by immunofluorescence for MF20 and RT-PCR for myogenic gene expression. The ALK inhibitor SB-431542 (S4317) was purchased from Sigma-Aldrich and added to the culture at a final concentration of 1 µM.

For mouse osteoblast differentiation, GFP-expressing Prx1⁺ skeletal cells were FACS-sorted from P1 Ecsit;^+/+Prx1-Cre;Rosa^mT/mG and Ecsit;^fl/flPrx1-Cre;Rosa^mT/mG neonates using cell surface markers (CD45^-Tie2^-Ter119^-). 6 days after osteogenic culture, skeletal cells were fixed with10% neutral formalin buffer and stained with the solution containing Fast Blue (Sigma-Aldrich, FBS25) and Naphthol AS-MX (Sigma-Aldrich, 855) to assess alkaline phosphatase (ALP) activity. Subsequently, cells were washed and incubated with a solution containing 6.5 mM Na2CO3, 18.5 mM NaHCO₃, 2 mM MgCl₂, and phosphatase substrate (Sigma-Aldrich, S0942), and ALP activity was measured by spectrometer (BioRad). Finally, cell survival was analyzed by staining with Alamar Blue solution (DAL1100, Invitrogen).

For osteoclast differentiation, bone marrow cells were flushed from long bones of two-month-old Ecsit^fl/fl mice and cultured in the presence of 10 ng/ml M-CSF (R&D systems) to obtain bone marrow-derived macrophages (BMMs). Twenty-four hours later, the non-adherent cells were re-plated, cultured in the medium containing 20 ng/ml of M-CSF for 3 days, and then 10 ng/ml of RANKL were added to induce osteoclast differentiation. Twelve hours later, BMMs were treated with rAAV9 vectors expressing control vector or Cre recombinase, cultured in the presence of 20 ng/ml of M-CSF and 10 ng/ml of Rank ligand for 6 days, and osteoclast differentiation was analyzed by TRAP staining and RT-PCR analysis.

Quantitative RT-PCR and immunoblotting

Fresh tibialis anterior (TA) muscle and gastrocnemius (GA) muscle tissues were dissected and homogenized mechanically in QIAzol (QIAGEN) to isolate total RNA samples. cDNA was synthesized using the HighCapacity cDNA Reverse Transcription Kit from Applied Biosystems. Quantitative RT-PCR was performed using SYBR® Green PCR Master Mix (Applied Biosystems) with QuantStudio (TM) 6 Flex System (Applied Biosystems). The primers used for PCR are described in Supplemental Table 2.

For immunoblotting analysis, cells were lysed in TNT lysis buffer (50 mM Tris-HCl (pH 7.4), 150 mM NaCl, 1% Triton X-100, 1 mM EDTA, 1 mM EGTA, 50 Mm NaF, 1 mM Na₃VO₄, 1 mM PMSF, and protease inhibitor cocktail (Sigma-Aldrich), and protein amounts from cell lysates were measured using the DC protein assay (BioRad). Equivalent amounts of proteins were subjected to SDS-PAGE, transferred to Immunobilon-P membranes (Millipore), immunoblotted with antibodies specific to ECSIT (HPA042979, Sigma,1:1000), ECSIT (a gift from Dr. Sankar Gosh, Columbia University), NDUFS3 (43–9200, Invitrogen, 1:1000), PINK1 (BC100, Novus Biologicals, 1:1000), LC3B (2775S, Cell Signaling Technology, 1:1000), ND6 (NBP2–94464, Novus Biologicals, 1:1000), NDUFAF1 (PA5–59051, Thermofisher, 1:1000), HSP90 (4877s, Cell Signaling Technology, 1:1000), RUNX2 (sc-390351, Santa Cruz Biotechnology,1:1000) and developed with ECL (Thermo Fisher Scientific). Immunoblotting with anti-GAPDH antibody (CB1001, Emd millipore, 1:1000) was used as a loading control. Anti-OXPHOS Cocktail Kit was purchased from ThermoScientific (458099, 1:250).

RNA sequencing and whole transcriptome analysis

RNA-seq analysis was performed using OneStopRNAseq⁶⁷. Specifically, FastQC⁶⁸ and MultiQC⁶⁹ were used for raw read quality control and QoRTs⁷⁰ for post-alignment quality control. Paired-end reads were aligned to mouse genome mm10, with star_2.7.5a, annotated with gencode.vM25.primary_assembly annotation⁷¹. Aligned exon fragments with mapping quality higher than 20 were counted toward gene expression with featureCounts_2.0.0⁷² with default settings except for the following parameters ‘-Q 20–minOverlap 1–fracOverlap 0 -p -B -C’. Differential expression (DE) analysis was performed with DESeq2_1.28.1⁷³. Within DE analysis, ’ashr’ was used to create log2 Fold Change (LFC) shrinkage for each comparison. Significant DE genes (DEGs) were filtered with the criteria FDR <0.05 and absolute log2 fold change (|LFC|) >0.585. Gene set enrichment analysis was performed with GSEA⁷⁴.

Flow cytometry analysis

Bone marrow cells and splenocytes were isolated from p10 Ecsit^Prx1 and Ecsit^fl/fl mice. After incubation with red blood cells (RBC) lysis buffer containing NH₄Cl, KHCO₃, and EDTA), cells were passed through a 40 μm cell strainer, washed with FACS buffer (cold PBS (pH 7.2) containing 0.5% BSA (Fraction V) and 1 mM EDTA), and stained with mouse hematopoietic lineage antibody Cocktail (a mixture of antibodies against CD3, CD4, CD8, B220, Gr-1, CD11b and Ter119, 88–7772-72, Invitrogen, 1:100) and antibodies for cKit (12-1171-83, Invitrogen, 1:100) and Sca-1 (17-5981-82, Invitrogen, 1:100) for 30 min at 4 ˚C, washed with PBS, then 7-amino-actinomycin D (7AAD, 00-6993-50, Invitrogen, 1:1000) was added right before flow cytometry analysis to stain dead cells. Alternatively, cells were stained with Gr-1 (11-5931-85, Invitrogen, 1:100), B220 (48-0452-82, Invitrogen, 1:100), CD3ε (17-0031-82, Invitrogen, 1:100), CD11b (12-0112-82, Invitrogen, 1:100) to analyze immune cell populations. After washing three times, cells were re-suspended in cold PBS (pH 7.2) with 1 mM EDTA and analyzed with a LSRII (BD Biosciences) with the exclusion of 7AAD⁺ cells and doublets.

For flow cytometry analysis of muscle satellite cells, skeletal muscle tissues were dissected from P15 Ecsit^Prx1 and Ecsit^fl/fl hindlimbs and after 1 h enzymatic digestion at 37 °C, cells were stained with the following antibodies. CD45 (103103, BioLegend, 1:100), CD11b (101203, BioLegend, 1:100), Ter119 (116203, BioLegend, 1:100), Sca1 (108143, BioLegend, 1:100), integrin β-1 (102205, BioLegend, 1:100) and CXCR4 (146507, BioLegend, 1:100) for muscle-resident satellite cell staining while CD31 (102503, BioLegend, 1:100), CD45 (103103, BioLegend, 1:100), integrin-α7 (130-123-833, Miltenyi Biotec, 1:100), Sca1(108143, BioLegend, 1:100) for staining of skeletal muscle resident stem/progenitor cells. Ghost Dye (13–0865, Tonbo Biosciences, 1:100) was used to stain dead cells.

For flow cytometry analysis of skeletal stem cells, P1 Ecsit^Prx1 and Ecsit^fl/fl limbs were dissociated by mechanical and enzymatic digestion (1mg/ml of Collagenase D (11088866001, Sigma-Aldrich), 2 mg/ml of Dispase II (4942078001, Roche), 1 mg/ml of Hyaluronidase (H3506, Sigma-Aldrich) and 10000 unit/ml of DNase I (4716728001, Roche) for 1 hour at 37 °C under gentle agitation. After blocking with purified rat anti-mouse CD16/CD32 (101302, BioLegend, 1:100) for 1 hour on ice, cells were stained with Biotin-conjugated CD45 (103103, BioLegend, 1:100), Tie2 (124005, BioLegend, 1:100) and Ter119 (116203, BioLegend, 1:100) with BV421-conjuated streptavidin (405226, BioLegend, 1:100), PE-conjugated CD51(104105, BioLegend, 1:100), BV605-conjugated Thy1.1-2 (140317, BioLegend, 1:100), APC-conjugated CD200 (123809, BioLegend, 1:100), FITC-conjugated LY51 (108305, BioLegend, 1:100) and PE/Cy7-conjugated CD105 (120409, BioLegend, 1:100) for 45 min on ice. The grating strategy of flow cytometer to analyze skeletal progenitors is described in Supplemental Fig. 3. Alternatively, cells were isolated from P1 Ecsit;^+/+Prx1-Cre;Rosa^mT/mG, and Ecsit;^fl/flPrx1-Cre;Rosa^mT/mG limbs and stained with CD45, Tie2, and Ter119. After washing three times, cells were re-suspended in cold PBS (pH 7.2) with 1 mM EDTA and 1 µg/ml DAPI and sorted with a FACSAria II SORP cell sorter (Becton Dickinson) with exclusion of DAPI⁺ cells and doublets.

Bone fracture

Bone fracture surgery was performed on 10-week-old male Osx-cre, Ecsit^Osx, and Ecsit^fl/fl mice. A longitudinal skin incision was made along the lateral aspect of the thigh from the stifle joint to the hip. The lateral aspect of the femur was exposed by parting the vastus lateralis muscle and the rectus femoris muscle to expose the length of the femur while preserving the femoral nerve. The middle of the femoral shaft was excised with a surgical saw. Intramedullary fixation was performed with a 25G needle penetrating from the patella furrow of the distal femur to the greater trochanter tip of the femur. Both ends of the needle were bent and then cut with a wire cutter, leaving 1 mm. Fascia was sutured using a 4/0 Vicryl suture, and then skin was closed using a 4/0 Nylon suture. Radiographs of the injured legs were performed to monitor fracture healing 2 weeks post-surgery. Ten weeks later, microCT and histology were performed for skeletal analyses. Alternatively, 6 week-old Prx1;Rosa26^mTmG mice, while still under surgical plane of anesthesia, were placed into an Einhorn Device to create a closed femoral fracture⁷⁵, which was then confirmed by radiography. Nonunion fracture was defined as failure of bridging callus on anteroposterior and lateral fractured cortical bones by coronal and sagittal reconstruction view of microCT and radiography. Fracture unionization was defined as osseous consolidation evident on reconstruction view of microCT and radiography.

AAV-mediated expression of ECSIT proteins in mice

A single dose of rAAV9 vectors carrying control vector or FLAG-ECSIT constructs (2 x 10¹¹ GC, 50 μl) was randomly injected into P1 Ecsit^fl/fl and Ecsit^Prx1 neonates via the facial vein and 21 days later, skeletal and muscular phenotypes were assessed using X-ray, microCT, and histology. Alternatively, a single dose of rAAV9 vector carrying mCherry (2 x 10¹¹ GC, 50 μl) was injected into P1 Prx1;Rosa26^mTmG neonates via facial vein. Twenty-one days later, mCherry expression in individual tissues and the cryosectioned bone and its adjacent muscle were monitored by the IVIS-100 optical imaging and fluorescence microscopy, respectively.

Mitochondrial fractionation

Human bone marrow-derived stromal cells (BMSCs) were purchased from ScienCell Inc. (Cat #: 7500). Cells were resuspended in mitochondria buffer (70 mM sucrose, 1 mM EGTA, 210 mM sorbitol, 10 mM MOPS [pH 7.4])⁴¹ and then, the supernatant was collected after centrifugation for 10min (450g, 4 °C). For total fraction, 54 ul of the supernatant was collected after centrifugation for 10 min (550 g, 4 °C). The rest of supernatant was further centrifuged at 9500g for 10min at 4 °C and then, the supernatant was collected for the cytosol fraction. For mitochondria fraction, the pellet was resuspended in mitochondrial buffer, incubated on ice for 15 min, and centrifuged at 500 g for 10 min at 4 °C. The supernatant was further centrifuged at 9500g for 10min at 4 °C. Finally, the pellet was collected and resuspended in 54 ul mitochondria buffer as mitochondria fraction.

Statistical analysis

All data were presented as the mean ± SD. Sample sizes were calculated on the assumption that a 30% difference in the parameters measured would be considered biologically significant with an estimate of sigma of 10–20% of the expected mean. Alpha and Beta were set to the standard values of .05 and 0.8, respectively. No animals or samples were excluded from analysis, and animals were randomized to treatment versus control groups, where applicable. For relevant data analysis, where relevant, we first performed the Shapiro-Wilk normality test for checking normal distributions of the groups. If normality tests passed, two-tailed, unpaired Student’s t-test were used, but if normality tests failed, Mann-Whitney tests were used for the comparisons between two groups. For the comparisons of three or four groups, we used one-way ANOVA if normality tests passed, followed by Tukey`s multiple comparison test for all pairs of groups. If normality tests failed, Kruskal–Wallis test was performed and was followed by Dunn`s multiple comparison test. For the comparison of two proportions, we performed Fisher’s exact test. The GraphPad PRISM software (v.9.4.1, La Jolla, CA) and Microsoft Office Excel 2016 were used for statistical analysis. ns, not significant; P < 0.05 was considered statistically significant. *P < 0.05; **P < 0.01; ***P < 0.001; and ****P < 0.0001.

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.

Data availability

Data supporting the findings of this manuscript are available from the corresponding author upon request. The raw data are protected and are not available due to data privacy laws. However, the processed data generated in this study are provided in the Supplementary Information and Source Data file. In addition, the RNA-seq data generated in this study has been deposited in the NCBI database under accession code GSE214454. Source data are provided with this paper.

References

Suomalainen, A. & Battersby, B. J. Mitochondrial diseases: the contribution of organelle stress responses to pathology. Nat. Rev. Mol. cell Biol. 19, 77–92 (2018).
Article CAS PubMed Google Scholar
Trifunovic, A. et al. Premature ageing in mice expressing defective mitochondrial DNA polymerase. Nature 429, 417–423 (2004).
Article ADS CAS PubMed Google Scholar
Kujoth, G. C. et al. Mitochondrial DNA mutations, oxidative stress, and apoptosis in mammalian aging. Science 309, 481–484 (2005).
Article ADS CAS PubMed Google Scholar
Chinnery, P. & Turnbull, D. Mitochondrial DNA and disease. Lancet 354, S17–S21 (1999).
Article Google Scholar
Vercellino, I. & Sazanov, L.A. The assembly, regulation and function of the mitochondrial respiratory chain. Nat. Rev. Mol. Cell Biol. 1–21 (2021).
Guerrero-Castillo, S. et al. The assembly pathway of mitochondrial respiratory chain complex I. Cell Metab. 25, 128–139 (2017).
Article CAS PubMed Google Scholar
Hirst, J. Mitochondrial complex I. Annu. Rev. Biochem. 82, 551–575 (2013).
Article CAS PubMed Google Scholar
Jones, A. J., Blaza, J. N., Varghese, F. & Hirst, J. Respiratory complex I in Bos taurus and Paracoccus denitrificans pumps four protons across the membrane for every NADH oxidized. J. Biol. Chem. 292, 4987–4995 (2017).
Article CAS PubMed PubMed Central Google Scholar
Martín, M. A. et al. Leigh syndrome associated with mitochondrial complex I deficiency due to a novel mutation in the NDUFS1 gene. Arch. Neurol. 62, 659–661 (2005).
Article PubMed Google Scholar
Triepels, R., Van Den Heuvel, L., Trijbels, J. & Smeitink, J. Respiratory chain complex I deficiency. Am. J. Med. Genet. 106, 37–45 (2001).
Article CAS PubMed Google Scholar
Karamanlidis, G. et al. Mitochondrial complex I deficiency increases protein acetylation and accelerates heart failure. Cell Metab. 18, 239–250 (2013).
Article CAS PubMed PubMed Central Google Scholar
Kruse, S. E. et al. Mice with mitochondrial complex I deficiency develop a fatal encephalomyopathy. Cell Metab. 7, 312–320 (2008).
Article CAS PubMed PubMed Central Google Scholar
Koopman, W. J., Willems, P. H. & Smeitink, J. A. Monogenic mitochondrial disorders. N. Engl. J. Med. 366, 1132–1141 (2012).
Article CAS PubMed Google Scholar
Gorman, G. S. et al. Mitochondrial diseases. Nat. Rev. Dis. Prim. 2, 1–22 (2016).
MathSciNet Google Scholar
Manolagas, S. C. From estrogen-centric to aging and oxidative stress: a revised perspective of the pathogenesis of osteoporosis. Endocr. Rev. 31, 266–300 (2010).
Article CAS PubMed PubMed Central Google Scholar
Gandhi, S. S., Muraresku, C., McCormick, E. M., Falk, M. J. & McCormack, S. E. Risk factors for poor bone health in primary mitochondrial disease. J. Inherit. Metab. Dis. 40, 673–683 (2017).
Article CAS PubMed PubMed Central Google Scholar
Yang, S., Feskanich, D., Willett, W. C., Eliassen, A. H. & Wu, T. Association between global biomarkers of oxidative stress and hip fracture in postmenopausal women: a prospective study. J. Bone Miner. Res. 29, 2577–2583 (2014).
Article PubMed Google Scholar
Langdahl, J. H. et al. Mitochondrial point mutation m. 3243a> g associates with lower bone mineral density, thinner cortices, and reduced bone strength: a case‐control study. J. Bone Miner. Res. 32, 2041–2048 (2017).
Article CAS PubMed Google Scholar
Linnane, A., Ozawa, T., Marzuki, S. & Tanaka, M. Mitochondrial DNA mutations as an important contributor to ageing and degenerative diseases. Lancet 333, 642–645 (1989).
Article Google Scholar
Tzoufi, M. et al. A rare case report of simultaneous presentation of myopathy, Addison’s disease, primary hypoparathyroidism, and Fanconi syndrome in a child diagnosed with Kearns–Sayre syndrome. Eur. J. Pediatrics 172, 557–561 (2013).
Article Google Scholar
Cholley, F. et al. Mitochondrial respiratory chain deficiency revealed by hypothermia. Neuropediatrics 32, 104–106 (2001).
Article CAS PubMed Google Scholar
Flanagan, B. & Nichols, G. Jr Metabolic studies of bone in vitro: v. glucose metabolism and collagen biosynthesis. J. Biol. Chem. 239, 1261–1265 (1964).
Article CAS PubMed Google Scholar
Vander Heiden, M. G., Cantley, L. C. & Thompson, C. B. Understanding the Warburg effect: the metabolic requirements of cell proliferation. Science 324, 1029–1033 (2009).
Article ADS Google Scholar
Lee, W.-C., Guntur, A. R., Long, F. & Rosen, C. J. Energy metabolism of the osteoblast: implications for osteoporosis. Endocr. Rev. 38, 255–266 (2017).
Article PubMed PubMed Central Google Scholar
Guntur, A. R., Le, P. T., Farber, C. R. & Rosen, C. J. Bioenergetics during calvarial osteoblast differentiation reflect strain differences in bone mass. Endocrinology 155, 1589–1595 (2014).
Article PubMed PubMed Central Google Scholar
Esen, E. et al. WNT-LRP5 signaling induces Warburg effect through mTORC2 activation during osteoblast differentiation. Cell Metab. 17, 745–755 (2013).
Article CAS PubMed PubMed Central Google Scholar
Stegen, S. et al. HIF-1α promotes glutamine-mediated redox homeostasis and glycogen-dependent bioenergetics to support postimplantation bone cell survival. Cell Metab. 23, 265–279 (2016).
Article CAS PubMed PubMed Central Google Scholar
Lee, S.-Y., Abel, E. D. & Long, F. Glucose metabolism induced by Bmp signaling is essential for murine skeletal development. Nat. Commun. 9, 1–11 (2018).
Article ADS Google Scholar
Yu, Y. et al. Glutamine metabolism regulates proliferation and lineage allocation in skeletal stem cells. Cell Metab. 29, 966–978.e964 (2019).
Article CAS PubMed PubMed Central Google Scholar
Heide, H. et al. Complexome profiling identifies TMEM126B as a component of the mitochondrial complex I assembly complex. Cell Metab. 16, 538–549 (2012).
Article CAS PubMed Google Scholar
Formosa, L. E. et al. Dissecting the roles of mitochondrial complex I intermediate assembly complex factors in the biogenesis of complex I. Cell Rep. 31, 107541 (2020).
Article CAS PubMed Google Scholar
Vogel, R. O. et al. Cytosolic signaling protein Ecsit also localizes to mitochondria where it interacts with chaperone NDUFAF1 and functions in complex I assembly. Genes Dev. 21, 615–624 (2007).
Article CAS PubMed PubMed Central Google Scholar
Nouws, J. et al. Acyl-CoA dehydrogenase 9 is required for the biogenesis of oxidative phosphorylation complex I. Cell Metab. 12, 283–294 (2010).
Article CAS PubMed Google Scholar
de Lageneste, O. D. et al. Periosteum contains skeletal stem cells with high bone regenerative potential controlled by Periostin. Nat. Commun. 9, 1–15 (2018).
Google Scholar
Kawanami, A., Matsushita, T., Chan, Y. Y. & Murakami, S. Mice expressing GFP and CreER in osteochondro progenitor cells in the periosteum. Biochemical biophysical Res. Commun. 386, 477–482 (2009).
Article CAS Google Scholar
Moore, E. R., Yang, Y. & Jacobs, C. R. Primary cilia are necessary for Prx1-expressing cells to contribute to postnatal skeletogenesis. J. cell Sci. 131, jcs217828 (2018).
Article PubMed PubMed Central Google Scholar
Chan, C. F. et al. Identification and specification of the mouse skeletal stem cell. Cell 160, 285–298 (2015).
Kopp, E. et al. ECSIT is an evolutionarily conserved intermediate in the Toll/IL-1 signal transduction pathway. Genes Dev. 13, 2059–2071 (1999).
Article CAS PubMed PubMed Central Google Scholar
Wi, S. M. et al. TAK1-ECSIT-TRAF6 complex plays a key role in the TLR4 signal to activate NF-kappaB. J. Biol. Chem. 289, 35205–35214 (2014).
Article CAS PubMed PubMed Central Google Scholar
Fiorese, C. J. et al. The transcription factor ATF5 mediates a mammalian mitochondrial UPR. Curr. Biol. 26, 2037–2043 (2016).
Article CAS PubMed PubMed Central Google Scholar
Nargund, A. M., Fiorese, C. J., Pellegrino, M. W., Deng, P. & Haynes, C. M. Mitochondrial and nuclear accumulation of the transcription factor ATFS-1 promotes OXPHOS recovery during the UPRmt. Mol. Cell 58, 123–133 (2015).
Article CAS PubMed PubMed Central Google Scholar
Kang, B. H. et al. Regulation of tumor cell mitochondrial homeostasis by an organelle-specific Hsp90 chaperone network. Cell 131, 257–270 (2007).
Article CAS PubMed Google Scholar
Taipale, M., Jarosz, D. F. & Lindquist, S. HSP90 at the hub of protein homeostasis: emerging mechanistic insights. Nat. Rev. Mol. Cell Biol. 11, 515–528 (2010).
Article CAS PubMed Google Scholar
Ashrafi, G. & Schwarz, T. The pathways of mitophagy for quality control and clearance of mitochondria. Cell Death Differ. 20, 31–42 (2013).
Article CAS PubMed Google Scholar
Wei, J. et al. Glucose uptake and Runx2 synergize to orchestrate osteoblast differentiation and bone formation. Cell 161, 1576–1591 (2015).
Article CAS PubMed PubMed Central Google Scholar
Carneiro, F. R. G., Lepelley, A., Seeley, J. J., Hayden, M. S. & Ghosh, S. An essential role for ecsit in mitochondrial complex i assembly and mitophagy in macrophages. Cell Rep. 22, 2654–2666 (2018).
Article CAS PubMed PubMed Central Google Scholar
Yang, Y.-S. et al. Bone-targeting AAV-mediated silencing of Schnurri-3 prevents bone loss in osteoporosis. Nat. Commun. 10, 1–13 (2019).
ADS Google Scholar
Romeo, S. G. et al. Endothelial proteolytic activity and interaction with non-resorbing osteoclasts mediate bone elongation. Nat. Cell Biol. 21, 430–441 (2019).
Article CAS PubMed Google Scholar
Liu, L., Cheung, T. H., Charville, G. W. & Rando, T. A. Isolation of skeletal muscle stem cells by fluorescence-activated cell sorting. Nat. Protoc. 10, 1612–1624 (2015).
Article CAS PubMed PubMed Central Google Scholar
Joe, A. W. et al. Muscle injury activates resident fibro/adipogenic progenitors that facilitate myogenesis. Nat. Cell Biol. 12, 153–163 (2010).
Article CAS PubMed PubMed Central Google Scholar
Jones, T. I. et al. Facioscapulohumeral muscular dystrophy family studies of DUX4 expression: evidence for disease modifiers and a quantitative model of pathogenesis. Hum. Mol. Genet. 21, 4419–4430 (2012).
Article CAS PubMed PubMed Central Google Scholar
Girardi, F. et al. TGFβ signaling curbs cell fusion and muscle regeneration. Nat. Commun. 12, 1–16 (2021).
Article Google Scholar
Melendez, J. et al. TGFβ signalling acts as a molecular brake of myoblast fusion. Nat. Commun. 12, 1–11 (2021).
Article ADS Google Scholar
Ismaeel, A. et al. Role of transforming growth factor-beta in skeletal muscle fibrosis: a review. Int. J. Mol. Sci. 20, 2446 (2019).
Borle, A. B., Nichols, N. & Nichols, G. Jr Metabolic studies of bone in vitro: I. Normal bone. J. Biol. Chem. 235, 1206–1210 (1960).
Article CAS PubMed Google Scholar
Cohn, D. V. & Forscher, B. K. Aerobic metabolism of glucose by bone. J. Biol. Chem. 237, 615–618 (1962).
Article CAS PubMed Google Scholar
Peck, W. A., Birge, S. J. & Fedak, S. A. Bone cells: biochemical and biological studies after enzymatic isolation. Science 146, 1476–1477 (1964).
Article ADS CAS PubMed Google Scholar
Lee, W.-C., Ji, X., Nissim, I. & Long, F. Malic enzyme couples mitochondria with aerobic glycolysis in osteoblasts. Cell Rep. 32, 108108 (2020).
Article CAS PubMed PubMed Central Google Scholar
Peck, W. A., Birge, S. J. Jr & Brandt, J. Collagen synthesis by isolated bone cells: stimulation by ascorbic acid in vitro. Biochimica et. Biophysica Acta (BBA)-Nucleic Acids Protein Synth. 142, 512–525 (1967).
Article CAS Google Scholar
Gerstenfeld, L. C., Chipman, S. D., Glowacki, J. & Lian, J. B. Expression of differentiated function by mineralizing cultures of chicken osteoblasts. Developmental Biol. 122, 49–60 (1987).
Article CAS Google Scholar
Martínez-Reyes, I. & Chandel, N. S. Mitochondrial TCA cycle metabolites control physiology and disease. Nat. Commun. 11, 1–11 (2020).
Article Google Scholar
Martínez-Reyes, I. et al. TCA cycle and mitochondrial membrane potential are necessary for diverse biological functions. Mol. Cell 61, 199–209 (2016).
Article PubMed Google Scholar
Kuzmin, D. A. et al. The clinical landscape for AAV gene therapies. Nat. Rev. Drug Discov. 20, 173–174 (2021).
Muzumdar, M. D., Tasic, B., Miyamichi, K., Li, L. & Luo, L. A global double-fluorescent Cre reporter mouse. Genesis 45, 593–605 (2007).
Article CAS PubMed Google Scholar
Xie, J. et al. Short DNA hairpins compromise recombinant adeno-associated virus genome homogeneity. Mol. Ther. 25, 1363–1374 (2017).
Article CAS PubMed PubMed Central Google Scholar
McLeod, M. J. Differential staining of cartilage and bone in whole mouse fetuses by alcian blue and alizarin red S. Teratology 22, 299–301 (1980).
Article CAS PubMed Google Scholar
Li, R., Hu, K., Liu, H., Green, M. R. & Zhu, L. J. OneStopRNAseq: a web application for comprehensive and efficient analyses of RNA-Seq data. Genes 11, 1165 (2020).
Article CAS PubMed Central Google Scholar
Andrews, S. & Bioinformatics, B. FastQC: a quality control tool for high throughput sequence data, citeulike-article-id: 11583827. arXiv Prepr. arXiv 0906, 2747 (2010).
Google Scholar
Ewels, P., Magnusson, M., Lundin, S. & Käller, M. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics 32, 3047–3048 (2016).
Article CAS PubMed PubMed Central Google Scholar
Hartley, S. W. & Mullikin, J. C. QoRTs: a comprehensive toolset for quality control and data processing of RNA-Seq experiments. BMC Bioinforma. 16, 1–7 (2015).
Article Google Scholar
Harrow, J. et al. GENCODE: the reference human genome annotation for The ENCODE Project. Genome Res. 22, 1760–1774 (2012).
Article CAS PubMed PubMed Central Google Scholar
Liao, Y., Smyth, G. K. & Shi, W. featureCounts: an efficient general purpose program for assigning sequence reads to genomic features. Bioinformatics 30, 923–930 (2014).
Article CAS PubMed Google Scholar
Love, M. I., Huber, W. & Anders, S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol. 15, 1–21 (2014).
Article Google Scholar
Subramanian, A. et al. Gene set enrichment analysis: a knowledge-based approach for interpreting genome-wide expression profiles. Proc. Natl Acad. Sci. USA 102, 15545–15550 (2005).
Article ADS CAS PubMed PubMed Central Google Scholar
Bonnarens, F. & Einhorn, T. A. Production of a standard closed fracture in laboratory animal bone. J. Orthop. Res. 2, 97–101 (1984).
Article CAS PubMed Google Scholar

Download references

Acknowledgements

We would like to thank many individuals who provided reagents and the Flow Core for FACS sorting (S10 grant 1S10OD028576). This project was supported by a NIH NIAMS grant: R21AR077557 and AAVAA Therapeutics. M.B.G. holds a Career Award for Medical Scientists from the Burroughs Wellcome Fund, NIH grant: R01AR075585, and a Pershing Square Sohn Cancer Research Alliance award.

Author information

These authors contributed equally: Qiyuan Yang, Dongsheng Guo.

Authors and Affiliations

Department of Medicine/Division of Rheumatology, UMass Chan Medical School, Worcester, MA, USA
Chujiao Lin, Yeon-Suk Yang, Sachin Chaugule, Ngoc DeSouza, Won-Taek Oh, Zhihao Chen, Aijaz A. John, Qiang Qiu, Shaoguang Li & Jae-Hyuck Shim
Department of Molecular, Cell and Cancer Biology, UMass Chan Medical School, Worcester, MA, USA
Qiyuan Yang, Rui Li, Lihua Julie Zhu & Cole Haynes
Department of Neurology, UMass Chan Medical School, Worcester, MA, USA
Dongsheng Guo & Charles P. Emerson
Horae Gene Therapy Center, UMass Chan Medical School, Worcester, MA, USA
Jun Xie, Guangping Gao & Jae-Hyuck Shim
Department of Microbiology and Physiological Systems, UMass Chan Medical School, Worcester, MA, USA
Jun Xie & Guangping Gao
Viral Vector Core, UMass Chan Medical School, Worcester, MA, USA
Jun Xie & Guangping Gao
Department of Pathology and Laboratory Medicine, Weill Cornell Medical College, Cornell University, New York, NY, USA
Matthew B. Greenblatt
Research Divisions, Hospital for Special Surgery, New York, NY, USA
Matthew B. Greenblatt
Department of Microbiology and Immunology, Columbia University Vagelos College of Physicians and Surgeons, New York, NY, USA
Sankar Ghosh
Li Weibo Institute for Rare Diseases Research, UMass Chan Medical School, Worcester, MA, USA
Guangping Gao & Jae-Hyuck Shim
Wellstone Muscular Dystrophy Program, UMass Chan Medical School, Worcester, MA, USA
Charles P. Emerson

Authors

Chujiao Lin
View author publications
You can also search for this author in PubMed Google Scholar
Qiyuan Yang
View author publications
You can also search for this author in PubMed Google Scholar
Dongsheng Guo
View author publications
You can also search for this author in PubMed Google Scholar
Jun Xie
View author publications
You can also search for this author in PubMed Google Scholar
Yeon-Suk Yang
View author publications
You can also search for this author in PubMed Google Scholar
Sachin Chaugule
View author publications
You can also search for this author in PubMed Google Scholar
Ngoc DeSouza
View author publications
You can also search for this author in PubMed Google Scholar
Won-Taek Oh
View author publications
You can also search for this author in PubMed Google Scholar
Rui Li
View author publications
You can also search for this author in PubMed Google Scholar
Zhihao Chen
View author publications
You can also search for this author in PubMed Google Scholar
Aijaz A. John
View author publications
You can also search for this author in PubMed Google Scholar
Qiang Qiu
View author publications
You can also search for this author in PubMed Google Scholar
Lihua Julie Zhu
View author publications
You can also search for this author in PubMed Google Scholar
Matthew B. Greenblatt
View author publications
You can also search for this author in PubMed Google Scholar
Sankar Ghosh
View author publications
You can also search for this author in PubMed Google Scholar
Shaoguang Li
View author publications
You can also search for this author in PubMed Google Scholar
Guangping Gao
View author publications
You can also search for this author in PubMed Google Scholar
Cole Haynes
View author publications
You can also search for this author in PubMed Google Scholar
Charles P. Emerson
View author publications
You can also search for this author in PubMed Google Scholar
Jae-Hyuck Shim
View author publications
You can also search for this author in PubMed Google Scholar

Contributions

C.L. designed, executed, and interpreted the experiments. J.X. and G.G. generated the AAV vectors. Q.Y. and D.G. performed mitochondrial and skeletal muscle experiments, respectively. Y.Y., Z.C., and A.A.J. performed histology and microCT. S.C. N.D. and Q.Q. performed flow cytometer analysis of bone marrow. S.C. performed osteoclast experiments. W.T.O. performed bone fracture experiments. R.L. and L.J.Z. performed whole transcriptome analysis. S.G. generated Ecsit-floxed mice. M.G., S.L., C.H., and C.P.E. interpreted the experiments and helped draft the manuscript. J.H.S. supervised the research and prepared the manuscript.

Corresponding author

Correspondence to Jae-Hyuck Shim.

Ethics declarations

Competing interests

G.G. and J.H.S. are co-founders of AAVAA Therapeutics and hold equity of this company. G.G. is a scientific co-founder of Voyager Therapeutics and Aspa Therapeutics Inc., and an inventor on patents with potential royalties licensed to other biopharmaceutical companies, in which they hold equity. Other authors declare no competing interests.

Peer review

Peer review information

Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available.

Additional information

Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Peer Review File

Reporting Summary

Source data

Source Data

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

Reprints and permissions

About this article

Cite this article

Lin, C., Yang, Q., Guo, D. et al. Impaired mitochondrial oxidative metabolism in skeletal progenitor cells leads to musculoskeletal disintegration. Nat Commun 13, 6869 (2022). https://doi.org/10.1038/s41467-022-34694-8

Download citation

Received: 17 January 2022
Accepted: 27 October 2022
Published: 11 November 2022
DOI: https://doi.org/10.1038/s41467-022-34694-8

This article is cited by

Metabolic regulation of skeletal cell fate and function
- Steve Stegen
- Geert Carmeliet
Nature Reviews Endocrinology (2024)
Molecular mechanisms of cellular metabolic homeostasis in stem cells
- Xiaoyu Li
- Ou Jiang
- Songlin Wang
International Journal of Oral Science (2023)

Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.