Abstract Pathogenic factors associated with aging, such as oxidative stress and hormone depletion converge on mitochondria and impair their function via opening of the mitochondrial permeability transition pore (MPTP). The MPTP is a large non-selective pore regulated by cyclophilin D (CypD) that disrupts mitochondrial membrane integrity. MPTP involvement has been firmly established in degenerative processes in heart, brain, and muscle. Bone has high energy demands and is therefore expected to be highly sensitive to mitochondrial dysfunction. Despite this, the role of mitochondria and the MPTP in bone maintenance and bone pathology has not been elucidated. Our goal was to determine whether mitochondria are impaired in aging bone and to see if protecting mitochondria from MPTP opening via CypD deletion protects against bone loss. We found that bone mass, strength, and formation progressively decline over the course of 18 months in C57BL/6J mice. Using metabolomics and electron microscopy, we determined that oxidative metabolism is impaired in aging bone leading to a glycolytic shift, imbalance in nucleotides, and decreased NAD+/NADH ratio. Mitochondria in osteocytes appear swollen which is a major marker of MPTP opening. CypD deletion by CypD knockout mouse model (CypD KO) protects against bone loss in 13- and 18-month-old mice and prevents decline in bone formation and mitochondrial changes observed in wild type C57BL/6J mice. Together, these data demonstrate that mitochondria are impaired in aging bone and that CypD deletion protects against this impairment to prevent bone loss. This implicates CypD-regulated MPTP and mitochondrial dysfunction in the impairment of bone cells and in aging-related bone loss. Our findings suggest mitochondrial metabolism as a new target for bone therapeutics and inhibition of CypD as a novel strategy against bone loss.

Citation: Shum LC, White NS, Nadtochiy SM, Bentley KLdM, Brookes PS, Jonason JH, et al. (2016) Cyclophilin D Knock-Out Mice Show Enhanced Resistance to Osteoporosis and to Metabolic Changes Observed in Aging Bone. PLoS ONE 11(5): e0155709. https://doi.org/10.1371/journal.pone.0155709 Editor: Yufeng Dong, Louisiana State University, UNITED STATES Received: March 17, 2016; Accepted: May 3, 2016; Published: May 16, 2016 Copyright: © 2016 Shum et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting Information files. Funding: LS was supported by National Center for Advancing Translational Sciences TL1-TR000096 grant. Work in the lab of PSB was supported by National Institutes of Health R01-HL071158 grant. RE was supported by National Institute of Arthritis and Musculoskeletal and Skin Diseases K01 AR064610, R03 AR061515, and P30 AR061307, and National Center for Advancing Translational Sciences UL1-TR000042 grants. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction Aging is a complex phenomenon marked by deterioration in form and function of many tissues. It is a multifactorial process with well-established hallmarks, such as loss of hormonal stimulation and increased oxidative stress. These and other age-associated factors converge on mitochondria, and disrupt their function [1]. As mitochondria deteriorate, they lose capacity for oxidative phosphorylation (OxPhos), biosynthesis, and ion transport. The relationship between mitochondrial dysfunction and aging has been studied in many tissues, including the brain, muscle, and cardiovascular system [2–7], but not bone. Since bone tissue has high energy demands due to continuous remodeling, intact and functioning mitochondria are important to bone health. The most common cause of mitochondrial dysfunction is increased activity of the mitochondrial permeability transition pore (MPTP), a large, non-selective mitochondrial pore that opens in response to stresses such as excessive calcium, reactive oxygen species, and hormone depletion [8–10]. MPTP opening dissipates the mitochondrial membrane potential, inducing mitochondrial swelling [11]. MPTP-mediated mitochondrial dysfunction has been associated with cell death during aging and other pathologies in many tissues [11–13]. The MPTP is opened by its protein regulator, cyclophilin D (CypD), a chaperone that also assists in protein folding. Deletion of CypD protects mitochondrial function and prevents necrotic cell death, by preventing opening of the MPTP [14]. In fact, CypD deletion via global knockout (KO) in mice showed beneficial results in models of Alzheimer's disease and cardiac reperfusion injuries [15,16]. Despite this, CypD deletion has not been investigated in bone pathologies. Reduced bone mass and mineral density are common characteristics in aging bone, leading to a significant increase of fracture risk in the elderly. In addition to hormone depletion, oxidative stress has come into focus as a second major pathogenic factor in age-associated bone loss [17,18]. Mitochondria are both a major source and major sensor of oxidative stress and are sensitive to changes in hormonal stimuli; however their role in aging bone has not been thoroughly studied. Our recent work and studies from other groups highlight the importance of mitochondrial oxidative phosphorylation in osteogenic cells [19–21]. Therefore, osteogenically differentiating bone marrow stem/stromal cells (BMSCs) and their progeny, osteoblasts (OB) and osteocytes (OT), are expected to be especially sensitive to mitochondrial dysfunction, leading to decreased bone formation and increased bone loss. Considering the multifactorial nature of aging and the fact that aging-related stresses converge on mitochondria via the CypD-mediated MPTP, we tested the hypothesis that mitochondrial dysfunction and the MPTP are involved in bone loss. Using metabolomics, electron microscopy and a global CypD knockout mouse model (CypD KO), we found that CypD deletion protects mitochondrial function and is beneficial for aging bone.

Materials and Methods Materials Chemicals were from Sigma unless otherwise noted. Cell culture media was from Invitrogen or Stem Cell Technologies. Animals C57BL/6J wild type mice were from the National Institute of Aging and were acclimated at our vivarium facility. CypD KO mice (C57BL/6J background) were a kind gift from Dr. George Porter (University of Rochester) [22]. CypD KO mice were backcrossed to C57BL/6J for five generations. CypD KO mice were not significantly different in size compared to the C57BL/6J mice. Male mice were used for experiments. Animal husbandry and experiments were performed in accordance with the Division of Laboratory Animal Medicine, University of Rochester, state and federal law, and National Institutes of Health policy. University of Rochester Institutional Animal Care and Use Committee (IACUC) specifically approved this study. Conditions of animal husbandry, i.e. food type, light schedule, cage density, etc., at NIA were similar to the conditions at our animal facility. Bone microCT and biomechanical testing After sacrifice, femurs, tibiae, and spine were isolated and cleaned of soft tissue. Right femurs were wrapped in phosphate buffered saline soaked gauze and stored at -80°C to prevent drying until biomechanical testing was performed. For micro computed tomography (microCT), bones were fixed in 10% neutral buffered formalin (NBF) for 72 hours and imaged using a VivaCT 40 tomograph (Scanco Medical). A calibrating phantom was used to standardize radiodensities among scans. Volume quantification was performed using Scanco analysis software. Bone vs total volume (BV/TV), trabecular number (Tb. Number), trabecular thickness (Tb. Thickness), trabecular separation (Tb. Separation) and cortical thickness (Cort. Thickness) were determined for tibiae, femurs, and the third lumbar vertebrae of the spine (L3). For biomechanical testing, right femurs were subjected to torsion testing. Briefly, samples were held in bone cement (DePuy Orthopaedics) in aluminum holders and tested using an EnduraTec TestBench system (Bose) at a displacement rate of 1°/s until failure. The torque data were plotted against rotational deformation to determine the maximum torque, torsional rigidity, and energy to maximum. Histology NBF-fixed samples were processed for histology via decalcification in EDTA for two weeks followed by paraffin embedding. Sections were cut to 5 μm and stained with Alcian Blue/Hematoxylin and Orange G (ABH/OG). Histomorphometry ABH/OG stained slides were visualized using an Axioscope 40 (Zeiss) microscope equipped with an Olympus DP72 camera (Olympus) and evaluated with Osteomeasure software (OsteoMetrics). Slides were analyzed to measure trabecular bone area to total area in proximal tibia, distal femur, and L3. Three different slides were counted per mouse and averaged. There were at least five mice per age group. Bone formation rate Mice were labeled via intraperitoneal injection with Alizarin Red (Day 0), followed by calcein (Day 7). Seven days after calcein injection (Day 14), mice were sacrificed. Long bones were collected, stripped of soft tissue, and either fixed in 70% ethanol and shipped to the Orthopaedic Histology and Histomorphometry Laboratory at the Yale School of Medicine for further processing, embedding into methylmethacrylate, and sectioning, or processed for frozen sectioning as this technique became available in the Center for Musculoskeletal Research at the University of Rochester. Sections were visualized using a fluorescence Axioscope 40 microscope (Zeiss) equipped with an Olympus DP72 camera (Olympus) and analyzed using ImageJ software to calculate bone formation rate (BFR) as follows: BFR = MAR x (MS/BS) where MAR = Ir.L.Th/Ir.L.t. and MS = (dLS+sLS/2)/BS (MAR, mineral apposition rate; BS, bone surface; dLS, double labelled surface; sLS, single labelled surface; Ir.L.Th, distance between labels; and Ir.L.t, time between labels). Five mice per group were analyzed. Osteoclast activity Bone sections were stained for tartrate-resistant acid phosphatase (TRAP), counter-stained with FastGreen and scanned in an Olympus VS120 whole slide imager. TRAP+ OCs were assessed using Visiopharm software that calculated osteoclast surface per bone surface (Oc.S/B.S.). In addition, blood was collected and serum prepared and stored at -80°C until assayed. The bone resorption marker, CTX-I, was assessed in serum samples using the Ratlaps CTX-I EIA kit (Immune Diagnostic Systems) according to the manufacturer’s instructions. Absorbance was measured at 450 nm and 650 nm. Data were analyzed using GraphPad Prism 5 (GraphPad Software, Inc) to generate a standard curve and determine concentrations in pg/mL. Five mice per group were analyzed. Metabolomics Samples for metabolomic experiments were prepared as follows: bone tissue was cleaned of soft tissue, bone marrow, cartilage, and periosteum, then flash frozen in liquid nitrogen. Equal quantities of tissue were pulverized in liquid nitrogen and metabolites were extracted in 5 ml of 80% methanol. Extracts were dried under nitrogen stream and reconstituted in 200 μl of 50% methanol, and then analyzed using reverse phase liquid chromatography (LC) with an ion-pairing reagent in a Shimadzu HPLC coupled to a Thermo Quantum triple-quad mass spectrometer (MS). Data were analyzed using Mzrock machine learning tool kit (http://code.google.com/p/mzrock/), which automates analysis of targeted metabolomics data based on chromatographic retention time, whole molecule mass, collision energy, and resulting fragment mass. Data are presented as fold change over 3-mo-old mice. Three mice per group were analyzed. Electron microscopy Mice were perfused with 2.5% glutaraldehyde in 0.1 M sodium cacodylate buffer and bones were dissected, fixed for 24 hours, and decalcified as described above. Samples were post-fixed in 1.0% OsO4, dehydrated in a graded series of ethanol to 100%, and transitioned into propylene oxide for infiltration and embedding into EPON/Araldite epoxy resin. Thin sections (70nm) cut using a diamond knife were placed onto carbon/formvar slot grids, stained with aqueous uranyl acetate and lead citrate, examined and photographed using a Hitachi H-7650 transmission electron microscope with an attached Gatan Erlangshen 11 megapixel digital camera. The number of mitochondria in normal orthodox conformation and in swollen conformation was blindly counted by three investigators. An orthodox conformation is characterized by an electron-dense and crista-rich matrix, and a swollen conformation is characterized by an increased surface area and a less electron-dense matrix with fewer cristae. A total of 15 cells/sample were analyzed. Three mice per group were analyzed. Statistical analysis Data were analyzed using Prism 5.01 (GraphPad Software). Mean values and standard deviation were calculated, and the statistical significance (p < 0.05) established using ANOVA for >2 groups of variables or Student t-test for 2 groups of variables.

Discussion Here we show that bone loss in aging C57BL/6J mice coincides with metabolic changes in bone tissue indicative of mitochondrial dysfunction. We also observed increased mitochondrial swelling in osteocytes, indicating increased MPTP activity. Importantly, in mice with MPTP loss-of-function (CypD KO), we observed increased resistance against age-associated bone loss. CypD deletion also prevented the decline in bone formation, glycolytic shift and other metabolic changes in bone tissue indicative of impaired oxidative metabolism as well as mitochondrial swelling. There was also an effect of CypD deletion on bone resorbing OCs and a more pronounced decline in TRAP positive OCs in 13-mo old CypD KO mice than in 13-mo old wild type C57BL/6J mice. This effect requires further detailed investigation, however previous studies [29] already showed inverse correlation between mitochondrial activity and OC function. Thus, improved mitochondrial function due to CypD deletion could be a reason for accelerated OC decline in aging in CypD KO mice. While cells have multiple ways to produce energy, the most efficient is via mitochondrial oxidative phosphorylation. Hence, BMSCs upregulate oxidative phosphorylation during differentiation into the osteogenic lineage [19–21]. Glycolytic metabolism was found to be equally important for the osteogenic lineage [30,31], however, these findings by no means imply that mitochondria are not critical to bone physiology. Glycolysis dominance over OxPhos may be required at the initial proliferative stages of osteogenic differentiation, while activation of mitochondria may be needed for massive protein biosynthesis at later stages. These different findings indicate the metabolic plasticity of bone cells and suggest that further detailed studies of metabolic mechanisms in bone cells are required. There is still much to learn about the underlying mechanisms of bone loss during aging. Age-associated bone loss was previously thought to be primarily due to hormone depletion, such as the decline in estrogen production after menopause in females. More current opinions hold that this may not be the only cause, and that oxidative stress may also have strong impact [17,18]. Both of these elements are detrimental to bone, but their respective contributions have not been clarified. Importantly, both estrogen depletion and oxidative stress have harmful effects on mitochondria by inducing the MPTP [12,32–35]. Since this study was performed on male mice of moderate age, it is likely that oxidative stress was more relevant than hormone depletion. Pathological levels of oxidative stress (ROS) can activate CypD and, thus, promote the MPTP opening [16]. Oxidative stress found in aging bone tissue may damage mitochondria via the MPTP in bone cells, leading to decreased cell function and disruption of tissue homeostasis. The cause and effect relationship between ROS and MPTP in bone require detailed investigation and is a subject of our ongoing studies. We assessed changes in aging bone using multiple techniques, such as microCT, histomorphometry, and biomechanical testing. These assays clearly show the age-associated decrease in bone and impaired bone formation in wild type C57BL/6J mice, and lack thereof in CypD KO mice. Currently, the most informative technique that allows a snapshot of energy metabolism in whole tissue is metabolomics. We therefore used this technique to assess the metabolic profile of bone and detected a glycolytic shift indicative of mitochondrial dysfunction in aging wild type C57BL/6J, but not in CypD KO mice. To our knowledge, this is first time this technique has been applied to bone tissue. OTs comprise more than 70% of the cellular component of the bone tissue that was used in our metabolomic assay, and therefore, it is reasonable to assume that the metabolic changes observed in 13-mo-old wild type C57BL/6J mice were due to changes in OTs. This is further indicated by detection of mitochondrial swelling by EM in wild type C57BL/6J OTs and not in OBs. As evident from previously published works from our group and others [19–21], mitochondria are activated late in the course of osteogenic differentiation. Therefore while more work is needed to elucidate the contribution of different cell types, it is likely that mitochondria in OTs are most active and, thus, most vulnerable to the MPTP opening. We have not yet examined mice of advanced age, i.e. older than 24 mo, and thus cannot claim if MPTP inhibition via CypD KO can protect bone during advanced aging. Our data in 18-mo-old mice show that the protective effect of CypD deletion is somewhat diminished compared to 13-mo-old mice, and so it is possible that CypD deletion cannot completely protect bone in aging, and instead delay the onset of osteoporosis. Despite the fact that our CypD KO mice are of C57BL/6J background, our mice were not littermates and were bred in different facilities. This was a major limitation of our study that did not allow us to directly compare CypD KO mice to their wild type counterparts. Additionally, this work was done in global KO mice, so the observed effects could be caused by systemic changes induced by CypD deletion. However, a global KO is clinically relevant as a pharmacological CypD inhibitor would produce a similar systemic effect. To our knowledge, this is the first study using metabolomics in bone tissue and showing that protecting mitochondrial function is linked to improved bone health in aging. This opens a broad and overlooked target for new bone therapeutics: mitochondrial metabolism.

Acknowledgments We would like to thank Drs. Regis O’Keefe, Edward Schwarz, Matthew Hilton, Hani Awad, Brendan Boyce, Edward Puzas, and Cheryl Ackert-Bicknell for fruitful discussions and sharing resources, Dr. Eric Beier for help with aging mice data interpretation, Dr. Longze Zhang for help with Visiopharm, Dr. Tzong-jen Sheu for help with mouse perfusions, and Dr. George Porter for supplying CypD KO mice. We also thank Center for Musculoskeletal Research Core facilities: Histology, Biochemistry and Molecular Imaging (HBMI) core for technical assistance with histology, and the Biomechanics and Multimodal Tissue Imaging (BMTI) Core for technical assistance with microCT and the biomechanical testing. We thank Gayle Schneider for specimen sectioning in the URMC EM Shared Resource Lab, and the laboratory of Joshua Munger (Biochemistry, Rochester) for help with metabolomics analysis by LC-MS.

Author Contributions Conceived and designed the experiments: LS RE. Performed the experiments: LS NW SN KB JJ. Analyzed the data: LS NW SN KB JJ RE. Contributed reagents/materials/analysis tools: SN KB PB JJ RE. Wrote the paper: LS NW SN RE.