|Metabolomic Profiles Delineate Signature Metabolic Shifts during Estrogen Deficiency-Induced Bone Loss in Rat by GC-TOF/MS.|
|Jump to Full Text|
|PMID: 23408954 Owner: NLM Status: In-Data-Review|
|Postmenopausal osteoporosis is a complicated and multi-factorial disease. To study the metabolic profiles and pathways activated in osteoporosis, Eight rats were oophorectomized (OVX group) to represent postmenopausal osteoporosis and the other eight rats were sham operated (Sham group) to be the control. The biochemical changes were assessed with metabolomics using a gas chromatography/time-of-flight mass spectrometry. Metabolomic profile using serial blood samples obtained prior to and at different time intervals after OVX were analyzed by principal component analysis (PCA) and Partial least squares-discriminant analysis (PLS-DA). The conventional indicators (bone mineral density, serum Bone alkaline phosphatase (B-ALP) and N-telopeptide of type I collagen (NTx) of osteoporosis in rats were also determined simultaneously. In OVX group, the metabolomics method could describe the endogenous changes of the disease more sensitively and systematically than the conventional criteria during the progression of osteoporosis. Significant metabolomic difference was also observed between the OVX and Sham groups. The metabolomic analyses of rat plasma showed that levels of arachidonic acid, octadecadienoic acid, branched-chain amino acids (valine, leucine and isoleucine), homocysteine, hydroxyproline and ketone bodies (3-Hydroxybutyric Acid) significantly elevated, while levels of docosahexaenoic acid, dodecanoic acid and lysine significantly decreased in OVX group compared with those in the homeochronous Sham group. Considering such metabolites are closely related to the pathology of the postmenopausal osteoporosis, the results suggest that potential biomarkers for the early diagnosis or the pathogenesis of osteoporosis might be identified via metabolomic study.|
|Bo Ma; Jiannan Liu; Qi Zhang; Hanjie Ying; Jiye A; Jianguo Sun; Di Wu; Yonglu Wang; Jing Li; Yinhui Liu|
Related Documents :
|17713744 - Amino acids supply in culture media is not a limiting factor in the matrix synthesis of...
6208774 - Massive cutaneous hyalinosis. identification of the hyalin material as monoclonal kappa...
6431634 - Hormonal effects on glycosaminoglycans in thoracic aortas of rabbits.
8634804 - Comparison of the effect of various chemical peeling agents in a mini-pig model.
8097364 - Pharmacokinetics and pharmacodynamics of acid-suppressive agents in patients with gastr...
7748514 - Coagulation protein function. ii. influence of thiols upon acetaldehyde effects.
|Type: Journal Article Date: 2013-02-07|
|Title: PloS one Volume: 8 ISSN: 1932-6203 ISO Abbreviation: PLoS ONE Publication Date: 2013|
|Created Date: 2013-02-14 Completed Date: - Revised Date: -|
Medline Journal Info:
|Nlm Unique ID: 101285081 Medline TA: PLoS One Country: United States|
|Languages: eng Pagination: e54965 Citation Subset: IM|
|School of Pharmaceutical Sciences, Nanjing University of Technology, Nanjing, People's Republic of China.|
|APA/MLA Format Download EndNote Download BibTex|
Journal ID (nlm-ta): PLoS One
Journal ID (iso-abbrev): PLoS ONE
Journal ID (publisher-id): plos
Journal ID (pmc): plosone
Publisher: Public Library of Science, San Francisco, USA
Copyright: 2013 Ma et al
Received Day: 9 Month: 8 Year: 2012
Accepted Day: 17 Month: 12 Year: 2012
collection publication date: Year: 2013
Electronic publication date: Day: 7 Month: 2 Year: 2013
Volume: 8 Issue: 2
E-location ID: e54965
PubMed Id: 23408954
Publisher Id: PONE-D-12-23994
|Metabolomic Profiles Delineate Signature Metabolic Shifts during Estrogen Deficiency-Induced Bone Loss in Rat by GC-TOF/MS Alternate Title:Metabolomic Profiles in Bone Loss Rats|
|Julie A. Chowenedit1||
1School of Pharmaceutical Sciences, Nanjing University of Technology, Nanjing, People’s Republic of China
2Department of Respiration, Jiangsu Geriatric Hospital, Nanjing, People’s Republic of China
3School of Life Science and Pharmaceutical Engineering, Nanjing University of Technology, Nanjing, People’s Republic of China
4Key Laboratory of Drug Metabolism and Pharmacokinetics, China Pharmaceutical University, Nanjing, People’s Republic of China
5Laboratory for Systems Biology, Division of Clinical Pharmacology and Therapeutics, The Children’s Hospital of Philadelphia, Philadelphia, Pennsylvania, United States of America
6Kinetic Modeling and Simulation (KMAS) Core, Institute of Translational Medicine and Therapeutics, University of Pennsylvania, Philadelphia, Pennsylvania, United States of America
|Hosptial Infantil Universitario Niño Jesús, CIBEROBN, Spain
|Correspondence: * E-mail: firstname.lastname@example.org (QZ); email@example.com (HY)
[conflict] Competing Interests: The authors have declared that no competing interests exist.
Contributed by footnote: Conceived and designed the experiments: QZ HJY. Performed the experiments: BM JNL JGS YLW JL YHL. Analyzed the data: JYA. Contributed reagents/materials/analysis tools: QZ HJY. Wrote the paper: BM JYA DW.
Postmenopausal osteoporosis is a skeletal condition associated with reduced bone mineral and bone strength, involved in millions of people worldwide, especially those with pathological fracture. Osteoporosis is also called the “silent disease” in clinics because majority of people don’t know they have got osteoporosis until it has progressed and diagnosed at the point of fracture, most frequently occurred in the hip, wrist or spine, and the fracture often causes dangerous conditions and leads to deformity, and even death.
Bone mineral density (BMD) as a “gold standard” has been used in osteoporosis for a long time , . BMD test can indicate bone density at the normal, relatively low or osteoporotic levels, and predict the risk of fracture at the certain points. However, alterations in bone mineral density are slow in the dynamic disease progress of osteoporosis. Recently, Gourlay et al  attempted to standardize the bone-density testing interval (The BMD testing interval was defined as the estimated time for 10% of women to make the transition to osteoporosis before having a hip or clinical vertebral fracture, with adjustment for estrogen use and clinical risk factors.) and transition process to osteoporosis in elder women. Their data indicated that the bone-density testing interval for women with normal bone density or mild osteopenia as well as advanced osteopenia are 15, 5 and 1 year, respectively. As a clinical biomarker, bone mineral density has the disadvantages of slow change and low sensitivity, even frequent BMD testing is unlikely to improve the prediction of fracture and osteoporosis. For this reason, simple, sensitive and specific biomarkers are needed to be discovered, validated and applied for early diagnose of postmenopausal osteoporosis in clinic.
An association between an imbalance of bone formation and bone resorption was identified in pathological study on bone loss. Specific biochemical indicators for bone turnover, including bone formation markers (B-ALP; Osteocalcin et al), and bone resorption markers (NTx; Tartrate-resistant acid phosphatase-5b, TRCAP-5b; and Carboxy-terminal collagen crosslinks, CTX etc), might be used as index for disease progression of osteoporosis–. These sensitive and validated biochemical markers can offer an alternative to well-accepted BMD test to monitor disease progression of osteoporosis and therapeutic treatment , . The disadvantage of the biochemical markers is that they only reflect the alteration of bone formation or bone resorption, while the incidence of osteoporosis is attributed to the dual outcomes of bone formation and resorption.
Metabolomics as an important component of systems biology, including genomics, transcriptomics and proteomics, provide a wide spectrum of information on the biochemical finger print in cell, tissue or organism levels to elucidate novel mechanisms by detecting and comparing small-molecule metabolite profiles under difference conditions . Metabolomics is the endpoints of genotype functions and biochemical phenotype in body. Metabolic profiles detected by metabolomics in different conditions are linked closely to functions alteration in body . Biomarkers obtained by metabolomics are more sensitive to disease etiology and progression compared with those obtained by proteinomics and genomics , . Metabolomics has been used in the early detection and diagnosis of disease progression and provided prognostic biomarkers as novel therapeutic targets –. Postmenopausal osteoporosis is known as a complex disease, and many pathophysiologic factors involve in its occurrence and progression, including estrogen receptor , OPG/RANK/RANKL system , inflammatory factor  and oxidative stress . Considering there is no sensitive and specific biomarker indicating the pathogenesis of osteoporosis from a holistic point of view so far, metabolomics study might provide suitable approaches to investigate osteoporosis on disease etiology and progression, and identify a cluster of biomarkers of postmenopausal osteoporosis.
In our laboratory, metabolomic approach by using GC-MS or GC-TOF/MS has been successfully applied to the study of hypertension , hyperlipemia , chronic myeloid leukemia , colon carcinoma, and professional athletes  in serum or urine from rat or human. Metabolomic study on estrogen deficiency-induced obesity was performed recently and a series of biomarkers was identified to relate to estrogen deficiency . Based on our previous work, the metabolic alterations occurred in osteoporosis induced by estrogen deficiency are examined by using GC-TOF/MS. The goals of study are listed as following. 1) Using GC-TOF/MS to delineate the metabolic profile and monitor the endogenous metabolomics in rat plasma during the process of estrogen deficiency-induced bone loss. 2) Identifying and characterizing specific metabolic pathways served as targets for treatment intervention or as substrates for molecular imaging in postmenopausal osteoporosis. 3) Exploring mechanisms of bone loss induced by estrogen deficiency from the perspective of metabolomics.
Myristic-1,2-13C2 acid, 99 atom %13C, the stableisotopelabeled compound used as an internal standard (IS),Alkane series (C8–C40), MSTFA (N-methyl-N-trimethylsilyltrifluoroacetamid), and 1% TMCS (trimethylchlorosilane) were obtained from Sigma-Aldrich(St. Louis, MO, USA). HPLC-grade methanol and n-heptane were obtained from Merck (Damstadt, Germany). The reference substance octadecadienoic acid, arachidonic acid, docosahexaenoic acid, isoleucine, valine, leucine, homocysteine, hydroxyproline, tryptophan, cystine, succinic acid, citric acid, malic acid and lysine were purchased from Sigma-Aldrich (St. Louis, MO, USA). Purified water was produced by a Milli-Q system (Millipore, MA, USA).
Sixteen virgin Sprague–Dawley female rats with specific pathogen free (SPF) (6 months old, body weight 260±20.0 g) were obtained from the Animal Center of Nanjing Medical University in this study. Rats were housed under controlled conditions, including room temperature kept at 22±1°C, and standard solid food with water ad libitum provided during the experiment. This study was reviewed and approved by the Animal Ethical Committee of Nanjing University of Technology.
After acclimatization for 7 days, all animals were anesthetized with 30 mg/kg pentobarbital sodium (Fluka, Germany), intraperitoneal injection. Eight rats were oophorectomized (OVX group) and the remaining eight rats were sham operated (Sham group). Prior to the experiment and at 9∶00 am on the 3rd, 6th, 9th, 12th, 24th week after operation, 1 mL blood was collected and split equally into 2 tubes to separate serum and plasma. One aliquot whole blood was collected in the tube containing EDTA and placed for 15 min. Then plasma was immediately separated by centrifugation at 2000×g for 10 min at 4 °C. (no food but water ad libitum for 12 h). The serum was used for determinating bone alkaline phosphatase (B-ALP, an indicator of bone formation) and N-telopeptide of type I collagen (NTx) levels by ELISA method (Quantikine, R&D Systems, Minneapolis, MN, USA) following the manufacturer’s instruction with 6.7% intra- and 8.4% inter-assay variabilities for B-ALP as well as 7.2% intra- and 5.7% inter-assay variabilities for NTx, respectively. The plasma sample was used for metabonomic study. All samples were stored at −80°C until analysis. BMD of the lumbar vertebrae (L1–L5) and left femur from anesthetized rats was measured (isoflurane) by dual-energy X-ray absorptiometry (GE, Lunar ProdigyTM) at the 12th 24th week post-operation. BMD was calculated by dividing the measured area by the bone mineral content (BMC), as previously reported . Baseline of BMD was determined 1 week before surgery. All rats were weighed prior to surgical operation and monitored weekly post surgical operation. Uterine tissues were removed and weighed immediately to calculate uterus indexes by dividing the uterus weight by the body weight after rats were sacrificed. Meanwhile, the serum level of estradiol was analyzed immediately by ELISA method following the ELISA kit manufacturer’s instruction.
A 50 µL aliquot of plasma in an eppendorf tube was spiked with 200 uL methanol containing internal standard myristic-1, 2-13C2 acid (12.5 µg/ml) to deproteinize. After the mixture was vortexed for 3 min, the mixture was centrifuged for 10 min at 12000×g at 4°C. One hundred microliter of the supernatant was transferred into a GC vial, followed by being evaporated to dryness under vacuum in speedvac concentrator (Thermo Fisher Scientific, Asheville, USA). Subsequently, The residue was reconstituted in 30 µl methoxyamine in pyridine (15 mg/ml) and vortexed for 2 min. Methoximation reaction was performed at room temperature for 16 h, before trimethylsilylation for 1 h by adding 30 µl MSTFA with 1% TMCS as catalyst. Finally, the solution was vortexed for 30 s after methyl myristate in heptane (30 µg/ml), the external standard, was added to the GC vial. One microliter of solution was used for GC-TOF/MS analysis.
One microliter of the derivatized sample was injected into an Agilent 6980 GC equipped with an Agilent 7683 Series autosampler (Agilent, Atlanta, GA, USA).Chromatographic separation was conducted on a fused-silica capillary column (10 m×0.18 mm ID) chemically bonded with a 0.18 µm DB5-MS stationary phase (J&W Scientific, Folsom, CA, USA). Helium was served as carrier gas through the column. Mass spectra were detected and obtained by using Pegasus III TOF/MS (Leco).The gas flow rate through the column was 1 mL/min. The column initial temperature was kept at 70°C for 2 min. Temperature was increased from 70°C to 310°C at a rate of 20°C/min and held for 2 min. The transfer line temperature and ion source temperature were controlled at 250°C and 200°C, respectively. Ion source voltage and current: 70 eV at a current of 3.0 mA. The MS data were acquired in scan mode over the range between m/z 50 and 800 at a rate of 20 spectra/s, the acceleration voltage was turned on after a solvent delay of 170 s and the acceleration voltage was turned on at1,650 v , .
All GC/MS data was processed using ChromaTOF (version 3.25) software (Leco). Peak width in automatic peak detection and mass spectrum deconvolution were set to 2 s. Peaks with signal-to-noise (S/N) ratios lower than 30 were rejected. To obtain accurate peak areas for IS and specific peaks/compounds, one quant mass were specified for each peak and the data were reprocessed according to previous report, and each peak area was normalized using IS before multivariate data analysis. All compounds were identified and assigned by comparison of mass spectra and retention index of all detected compounds with the authentic reference standards and those in the NIST library 2.0 (2005), Wiley library, and in-house mass spectra library database established by the Key Laboratory of Drug Metabolism and Pharmacokinetics in China Pharmaceutical University , .
Before multivariate data analysis, the peak area ratio of each compound to a corresponding internal standard was calculated as the response. Then, the data were imported to SIMCA-P v11.5 software (Umetrics AB, Sweden) for multivariate data analysis. Partial least squares-discriminant analysis (PLS-DA) and Principal component analysis (PCA) were used for modeling metabonomic data from Sham group and OVX group.
Data obtained from BMD, body weight, commercial assay kits were presented as the mean ± SD. The data were analyzed using one-way analysis of variance followed by a least significant difference (LSD) test at P<0.01 and P<0.05.
The representative change tendency of body weight during the whole experiment period is listed in Table 1. The body weights of rats in OVX group were significantly higher than those in Sham group starting at 4th week post operation. Estrogen deficiency caused atrophy of uterine tissue, and serum E2 levels of rats in OVX group significantly deceased at the end of experiment compared with those in Sham group (shown in Table 1). Both data indicated the successful surgical procedure.
As shown in Table 1, at the12th week and 24th week after OVX, using DEXA analysis, significant decrease in BMD values for both the femur (P<0.05 for 12 weeks and P<0.05 for 24 weeks) and the lumbar vertebrae (P<0.05 for 12 weeks and P<0.01 for 24 weeks) were observed in rats in OVX group compared with those in Sham group.
Estrogen deficiency resulted in gradual increase in serum B-ALP and NTx in OVX rats after surgery. B-ALP and NTx of rats in OVX group significantly increased since the 12th week (P<0.05) and 9th week (P<0.05) post operation, respectively, compared with those in Sham group(shown in Figure 1A and 1B). The results indicated that OVX rat had the characteristic of bone high turnover rate according to the previous reports , .
Distinct differences of the major metabolites levels between the OVX and Sham groups were observed in GC-TOF/MS chromatogram, shown in Figure 2A and Figure 2B. PCA (data not shown) and PLS-DA were performed to further discern the presence of inherent similarities in chromatographic profiles. The PLS-DA score plots (Figure 3A) and 3D plots (Figure 3B) of the metabolites showed the metabolomic movement during the periods of osteoporotic progression including pre-operation, the 3rd, 6th, 9th and 12th week post operation (these five time points were marked as M0, M3, M6, M9, and M12 in the figures respectively). The distinguished separations among M0, M3, M6, M9, and M12 were observed clearly. Postoperative plasma metabolomic profiles were observed to be distinguished from preoperative samples starting at the 3rd week. However, classical pathological indexes including BMD, plasma B-ALP and NTx did not appear significant changes still the 9th weeks.
In PLS-DA, the R2X, R2Y and Q2(cum) parameters were used for the model evaluation, representing the explanation, fitness and prediction power respectively. (R2X is the percentage of all GC/TOF-MS response variables explained by the model. R2Y is the percentage of all observation or sample variables explained by the model. Q2 is the percentage of all observation or sample variables predicted by the model). In the PLS-DA score plots (Figure 3A), two principal components (PC1 and PC2) were calculated with the R2X, R2Y and Q2 parameters of 0.449, 0.848 and 0.56, respectively.
Kinetics of metabolic changes in rat plasma before and after bilateral oophorectomy was shown in Figure 4. Plasma levels of alanine, citric acid, succinic acid, and docosahexaenoic acid decreased significantly while homocysteine, branched-chain amino acids (leucine, isoleucine and valine), fatty acids (arachidonic acid and octadecadienoic acid) and cholesterol increased significantly after OVX operation compared with those before operation. It’s noted that the level of hydroxyproline was low at the 3rd week post operation, but unexpectedly and significantly increased at the 6th week. Figure 5C described the progression of the osteoporosis analyzed by PLS-DA, that plasma samples at the 12th week were well discriminated from those at the 24th week in OVX group.
The metabolic profiles of OVX and Sham groups at the same period were also analyzed and compared by PLS-DA. It’s shown in Figure 5A and Figure 5B that metabolic profiles of OVX and Sham groups separated distinctly at the 12th week and 24th week post operation. Two principal components were calculated with the set of R2X, R2Y and Q2 parameters as 0.335, 0.946 and 0.714, and 0.325, 0.996 and 0.582 respectively. Metabolic changes were listed in Figure 6 and Table 2. In detail, the plasma level of arachidonic acid in OVX group elevated at 12th week (2.5-fold, P = 0.0050) and 24th week (3.29-fold, P = 0.0023), compared with that in the homeochronous Sham group. Omega-6 fatty acids- octadecadienoic acid also elevated at 12th week (1.19-fold, P = 0.5302) and 24th week (1.85-fold, P = 0.0046). Homocysteine (1.53-fold at the 12th week, P = 0.042; 1.98-fold at the 24th week P = 0.0442) and Hydroxyproline (2.2-folds the 12th week, P = 0.0060; 2.43-folds at the 24th week P = 0.0053) concomitantly enhanced in plasma. However, the levels of omega-3 fatty acids (docosahexaenoic acid) in OVX group decreased (−1.22-fold P = 0.3508, −1.54-fold P = 0.0012) compared with those in the homeochronous Sham group. Elevated levels of branched-chain amino acids (leucine, isoleucine and valine), cholesterol, tryptophan, cystine and glycerol-3-phosphate were also found at the 12th week and 24thweek. Interestingly, a portion of metabolite changes were different in the12th week and 24th week. The levels of intermediate metabolites of TCA (succinic acid, citric acid and malic acid) significantly decreased at the 12th week, then gradually returned to the normal level at the 24th week.
A rat model of estrogen deficiency-induced bone loss was chosen to simulate postmenopausal osteoporosis. Comprehensive analysis of changes in the plasma metabolic profile was performed to determine step-wise disease progression from normal conditions to estrogen deficiency and eventually to bone loss.
As shown in Fig. 3A and Fig. 3B, distinct changes were observed from GC-TOF/MS chromatograms in OVX and Sham groups. Further investigations in biomarkers identification were performed and significant differences were shown in Table 2. The alteration of these molecules including amino acids, carbohydrates, fatty acids and cholesterol suggested that they might play certain roles in the etiology of postmenopausal osteoporosis. The early detection of metabolic changes linked to the onset of disease suggested that metabolic alterations might precede BMD and other related diagnostic indexes, and metabolomics might be a sensitive technology to diagnose osteoporosis superior to the classical pathology indexes.
Metabolite fingerprinting obtained from metabolomics could be applied not only to early diagnosis of diseases, but also to identify the causes of disease via biomarker index. Metabolic biomarkers identified by GC-TOF/MS would be beneficial to systemically understand the mechanisms of postmenopausal osteoporosis and to design the target points of anti-osteoporosis drugs. A great many endogenous molecules including fatty acids, cholesterol, amino acids, and carbohydrates exhibited significant difference in OVX group compared with those in Sham group. In the process of bone loss, the metabolism of fatty acids (octadecadienoic acid, arachidonic acid, oleic acid, docosahexaenoic acid, dodecanoic acid, aminomalonic acid, hexadecanoic acid), homocysteine, branched-chain amino acids (BCAA) and cholesterol were also significantly changed in our research.
Significant increase of the fatty acid levels in the process of bone loss are particularly noted in our metabolomic study. Arachidonic acid (AA, 20∶4n-6), the n-6 long chain polyunsaturated fatty acid (LCPUFA), is a precursor of prostaglandins, thromboxanes, and leukotrienes. AA play an important role in the regulation of the metabolism on lipid protein, the body’s inflammatory response, muscle growth, blood rheology, elasticity of blood vessels, white blood cell function and platelet activation–. Meanwhile, AA also affects the bone remodeling process. AA was found to promote osteoclastogenesis by stimulating RANK-L expression and inhibiting OPG secretion by osteoblasts. Its metabolites (prostaglandin and leukotrienes) also affect the bone formation and bone reabsorption. Coetzee et al, indicated that AA inhibited the secretion of OPG by osteoblasts, reduce the OPG/RANKL ratio and may increase osteoclastogenesis and prostaglandin E2 (PGE2), a pro-inflammatory lipid mediator derived from arachidonic acid by the activity of cyclooxygenase (COX), has similar effects on MG-63 and MC3T3-E1 osteoblast-like cells . Moreover, leukotrienes, the fatty signaling molecules derived from arachidonic acid, and particularly IL-1 and IL-6 (potent bone-resorbing cytokines), have been implicated in bone remodeling and disease–specifically in osteoporosis and rheumatoid arthritis–. Abnormal metabolism pathway of AA due to estrogen deficiency might act as an important factor in the pathogenesis of osteoporosis based on our study data.
In contrast to AA, the levels of oleic acid and docosahexaenoic acid (DHA, 22∶6n-3) decreased in plasma during bone loss induced by estrogen deficiency. An augmentation of the n-6/n-3 polyunsaturated fatty acid ratio was shown to lead to increased bone loss in both animal  and human . The elevated levels of the n-6 polyunsaturated fatty acids (AA, Octadecadienoic acid) and decreased levels of n-3 polyunsaturated fatty acids (Docosahexaenoic acid) also supported these results.
In our previous study, significant elevated levels of plasma cholesterol were found in estrogen deficiency-induced obesity in ovariectomized rats . Enhanced levels of cholesterol and increased body weight were also observed in bone loss rats induced by estrogen deficiency in the current study. Enhanced cholesterol uptake and biosynthesis or inhibition of cholesterol oxidation might explain the reason why estrogen deficiency increased total cholesterol accumulation in plasma. And, the negative relationship between estrogen and cholesterol might be an important index in the process of the menopause for the high incidence of cardiovascular disease, dementia, Alzheimer’s disease, type 2 diabetes and osteoporosis in menopausal women . Cholesterol and its metabolites would suppress the functional activity of osteoblasts and thereby induce reduced bone mineralization to decrease bone formation . By the same token, high cholesterol diet increases the risk of osteoporosis, possible via inhibiting the differentiation and proliferation of osteoblasts in rats, as You et al observed . At present, the drug treatment against osteoporosis can prevent bone loss, as well as lowering cholesterol levels. Arzoxifene, a new selective estrogen receptor modulator (SERM), maintained bone formation, prevented bone loss by decreasing osteoclast number following by reducing serum cholesterol in OVX rats in Ma and co-workers’ study . Statins, the enzyme HMG-CoA reductase inhibitors, initially used for treating hyperlipidemia, were recently observed not only to prevent fractures, but also trigger significant bone re-growth. From the mechanistic point of views, statins were found to inhibit mevalonate synthesis, prevent the synthesis of cholesterol and indirectly further alter osteoblasts and osteoclast activity, its direct biologic action also affect on osteoblast activity by stimulating gene expression of bone morphogenetic protein-2 , .
For the amino acids metabolism, above all, homocysteine levels significantly increased in OVX group compared to the Sham group. Homocysteine is a non-protein amino acid biosynthesized from methionine by the removal of its terminal Cε methyl group and recycled into methionine or converted into cysteine with the aid of B-vitamins. Homocysteine is still considered as part of a screen for people at the high risk of heart attack or stroke in clinics. A negative relationship between exogenous estrogen change and homocysteine levels was observed in post-menopause women, and estrogen blocking homocysteine-induced endothelial dysfunction in porcine coronary arteries was reported by Todd . In our study, a gradually increasing homocysteine level followed by decreasing estrogen concentrations was found in plasma after ovaries removed. A positive correlation of homocysteine and osteoporosis was also observed in metabolomic profiles. The research led by McLean and Jacques indicated homocysteine could be a predictive factor for osteoporotic fracture in older persons . For the correlation of homocysteine levels and osteoporosis, recent in vitro studies performed in mouse bone marrow cells have suggested that homocysteine directly stimulates osteoclast formation and activity (upregulated TRACP+ multinucleated cells and TRACP activity) by inducing p38 MAPK activity and intracellular reactive oxygen species (ROS) generation, which is a risk factor for osteoporosis and fracture . An increased homocysteine level in plasma was mainly attributed to the lack of estrogen caused in postmenopausal osteoporosis. In terms of specific biochemistry pathways, estrogen influenced metabolism transformation between methionine and homocysteine, especially interfering with the transsulfuration pathway. Estrogen could enhance cystathionine-β-synthase activity and directed homocysteine metabolism to form cysteine and glutathione to eliminate homocysteine and prevent its accumulation, that may act as a key aspect in elucidating this phenomenon .
It suggested that homocysteine could be a potential biomarker in diagnosis, and treatment against postmenopausal osteoporosis. Some anti-osteoporosis drugs such as strontium ranelate , raloxifene  etc, could inhibit increased homocysteine levels in plasma induced by estrogen deficiency and reduce the risk of osteoporosis. Therefore, homocysteine enhancing excretion of amino acids might be a strong indicator in both post-menopause and osteoporosis, and could be a potential target point or a key factor of upper stream controller for anti-osteoporosis treatment.
Hydroxyproline is a major component of the protein collagen and plays key role for collagen stability. Some research suggested that increased serum and urine levels of hydroxyproline had been demonstrated in Paget’s disease bone , . Meanwhile, mean plasma levels of branched-chain amino acids (BCAA, leucine, isoleucine, and valine) raised rapidly in the first three week after OVX and remained significantly elevated (p<0.05 versus sham-operated controls) for the entire 24 week. Overall, these observations indicated that increased homocysteine and hydroxyproline levels in OVX rats, besides increased accumulation of BCAA, may be potential biomarkers of postmenopausal osteoporosis.
Estrogen deficiency is the main reason of post-menopause osteoporosis. So, we used OVX rats, which are golden standard models of the post-menopause osteoporosis, to study the estrogen deficiency-induced bone loss and got the metabolomic profiles to delineate signature metabolic shifts during the progress of the estrogen deficiency. On the other hand, estrogen is involved in the regulation of many physiological functions of the body. Estrogen deficiency also induces other physical changes such as obesity. Therefore, on the basis of the current results, the further studies, that could link the metabolites found in this study to the osteoblasts or osteoclasts which are essential to the bone loss, should be carried out to confirm the specific biomarkers on estrogen deficiency-induced bone loss.
The global metabolomics approach applied in our research demonstrated that estrogen deficiency-induced bone loss resulted in a marked change in the plasma metabolic profile that appeared earlier than the known biomarkers of postmenopausal osteoporosis, such as BMD and serum B-ALP and NTx. In the preclinical osteoporosis model, significant changes of endogenous metabolites and specific metabolites were identified which not only regulated bone remodel effects, but also pointed to the mechanisms of disease etiology and progression. Further efforts should focuse on the validation of our findings and identification of genuine biomarkers from spurious biomarkers through preclinical and clinic study.
|1.||Villareal DT,, Binder EF,, Williams DB,, Schechtman KB,, Yarasheski KE,, et al. (Year: 2001) Bone mineral density response to estrogen replacement in frail elderly women: a randomized controlled trial. JAMA286: 815–82011497535|
|2.||Hillier TA,, Stone KL,, Bauer DC,, Rizzo JH,, Pedula KL,, et al. (Year: 2007) Evaluating the value of repeat bone mineral density measurement and prediction of fractures in older women: the study of osteoporotic fractures. Arch Intern Med167: 155–16017242316|
|3.||Gourlay ML,, Fine JP,, Preisser JS,, May RC,, Li C,, et al. (Year: 2012) Bone-density testing interval and transition to osteoporosis in older women. N Engl J Med366: 225–23322256806|
|4.||Kahl KG,, Rudolf S,, Stoeckelhuber BM,, Dibbelt L,, Gehl HB,, et al. (Year: 2005) Bone mineral density, markers of bone turnover, and cytokines in young women with borderline personality disorder with and without comorbid major depressive disorder. Am J Psychiatry162: 168–17415625216|
|5.||Bhattacharyya S,, Siegel ER,, Achenbach SJ,, Khosla S,, Suva LJ, (Year: 2008) Serum biomarker profile associated with high bone turnover and BMD in postmenopausal women. J Bone Miner Res23: 1106–111718302502|
|6.||Reginster JY,, Collette J,, Neuprez A,, Zegels B,, Deroisy R,, et al. (Year: 2008) Role of biochemical markers of bone turnover as prognostic indicator of successful osteoporosis therapy. Bone42: 832–83618316258|
|7.||Vasikaran SD, (Year: 2008) Utility of biochemical markers of bone turnover and bone mineral density in management of osteoporosis. Crit Rev Clin Lab Sci45: 221–25818415816|
|8.||Abe Y,, Ishikawa H,, Fukao A, (Year: 2008) Higher efficacy of urinary bone resorption marker measurements in assessing response to treatment for osteoporosis in postmenopausal women. Tohoku J Exp Med214: 51–5918212487|
|9.||Nicholson JK,, Lindon JC, (Year: 2008) Systems biology: Metabonomics. Nature455: 1054–105618948945|
|10.||Leslie RD,, Beyan H, (Year: 2011) Metabolomics makes a mark: early changes associated with autoimmune diabetes. Diabetes60: 2688–269022025776|
|11.||Patel VR,, Eckel-Mahan K,, Sassone-Corsi P,, Baldi P, (Year: 2012) CircadiOmics: integrating circadian genomics, transcriptomics, proteomics and metabolomics. Nat Methods9: 772–77322847108|
|12.||Patti GJ,, Yanes O,, Siuzdak G, (Year: 2012) Innovation: Metabolomics: the apogee of the omics trilogy. Nat Rev Mol Cell Biol13: 263–26922436749|
|13.||Sreekumar A,, Poisson LM,, Rajendiran TM,, Khan AP,, Cao Q,, et al. (Year: 2009) Metabolomic profiles delineate potential role for sarcosine in prostate cancer progression. Nature457: 910–91419212411|
|14.||Patti GJ,, Yanes O,, Shriver LP,, Courade JP,, Tautenhahn R,, et al. (Year: 2012) Metabolomics implicates altered sphingolipids in chronic pain of neuropathicorigin. Nat Chem Biol8: 232–23422267119|
|15.||Shi X,, Yao D,, Chen C, (Year: 2012) Identification of N-Acetyltaurine as a Novel Metabolite of Ethanol throughMetabolomics-guided Biochemical Analysis. J Biol Chem287: 6336–634922228769|
|16.||Chen C,, Shah YM,, Morimura K,, Krausz KW,, Miyazaki M,, et al. (Year: 2008) Metabolomics reveals that hepatic stearoyl-CoA desaturase 1 downregulation exacerbates inflammation and acute colitis. Cell Metab7: 135–14718249173|
|17.||Bonnelye E,, Aubin JE, (Year: 2005) Estrogen receptor-related receptor alpha: a mediator of estrogen response in bone. J Clin Endocrinol Metab90: 3115–312115713703|
|18.||Hofbauer LC,, Schoppet M, (Year: 2004) Clinical implications of the osteoprotegerin/RANKL/RANK system for bone and vascular diseases. JAMA292: 490–49515280347|
|19.||Pacifici R,, Vannice JL,, Rifas L,, Kimble RB, (Year: 1993) Monocytic secretion of interleukin-1 receptor antagonist in normal and osteoporotic women: effects of menopause and estrogen/progesterone therapy. J Clin Endocrinol Metab77: 1135–11418077304|
|20.||Wauquier F,, Leotoing L,, Coxam V,, Guicheux J,, Wittrant Y, (Year: 2009) Oxidative stress in bone remodelling and disease. Trends Mol Med15: 468–47719811952|
|21.||Lu Y,, A J,, Wang G,, Hao H,, Huang Q,, et al. (Year: 2008) Gas chromatography/time-of-flight mass spectrometry based metabonomic approach to differentiating hypertension- and age-related metabolic variation in spontaneously hypertensive rats. Rapid Commun Mass Spectrom22: 2882–288818720470|
|22.||Zhang Q,, Wang G,, A J,, Ma B,, Dua Y,, et al. (Year: 2010) Metabonomic profiling of diet-induced hyperlipidaemia in a rat model. Biomarkers15: 205–21620001218|
|23.||A J,, Qian S,, Wang G,, Yan B,, Zhang S,, et al. (Year: 2010) Chronic myeloid leukemia patients sensitive and resistant to imatinib treatment show different metabolic responses. PLOS ONE5: e1318620949032|
|24.||Yan B,, A J,, Wang G,, Lu H,, Huang X,, et al. (Year: 2009) Metabolomic investigation into variation of endogenous metabolites in professional athletes subject to strength-endurance training. J Appl Physiol106: 531–53819036890|
|25.||Ma B,, Zhang Q,, Wang GJ,, A JY,, Wu D,, et al. (Year: 2011) GC-TOF/MS-based metabolomic profiling of estrogen deficiency-induced obesity in ovariectomized rats. Acta Pharmacol Sin32: 270–27821293480|
|26.||Liu L,, Aa J,, Wang G,, Yan B,, Zhang Y,, et al. (Year: 2010) Differences in metabolite profile between blood plasma and serum. Anal Biochem406: 105–11220655863|
|27.||A J,, Trygg J,, Gullberg J,, Johansson AI,, Jonsson P,, et al. (Year: 2005) Extraction and GC/MS analysis of the human blood plasma metabolome. Anal Chem77: 8086–809416351159|
|28.||Jonsson P,, Johansson AI,, Gullberg J,, Trygg J,, A J,, et al. (Year: 2005) High-throughput data analysis for detecting and identifying differences between samples in GC/MS-based metabolomic analyses. Anal Chem77: 5635–564216131076|
|29.||Garnero P,, Sornay-Rendu E,, Chapuy MC,, Delmas PD, (Year: 1996) Increased bone turnover in late postmenopausal women is a major determinant of osteoporosis. J Bone Miner Res11: 337–3498852944|
|30.||Ma B,, Zhang Q,, Wu D,, Wang YL,, Hu YY,, et al. (Year: 2012) Strontium fructose 1,6-diphosphate prevents bone loss in a rat model of postmenopausal osteoporosis via the OPG/RANKL/RANK pathway. Acta Pharmacol Sin33: 479–48922426695|
|31.||Patel MI,, Kurek C,, Dong Q, (Year: 2008) The arachidonic acid pathway and its role in prostate cancer development and progression. J Urol179: 1668–167518343442|
|32.||Zordoky BN,, El-Kadi AO, (Year: 2010) Effect of cytochrome P450 polymorphism on arachidonic acid metabolism and their impact on cardiovascular diseases. Pharmacol Ther125: 446–46320093140|
|33.||Harizi H,, Corcuff JB,, Gualde N, (Year: 2008) Arachidonic-acid-derived eicosanoids: roles in biology and immunopathology. Trends Mol Med14: 461–46918774339|
|34.||Coetzee M,, Haag M,, Kruger MC, (Year: 2007) Effects of arachidonic acid, docosahexaenoic acid, prostaglandin E(2) and parathyroid hormone on osteoprotegerin and RANKL secretion by MC3T3-E1 osteoblast-like cells. J Nutr Biochem 200718: 54–63|
|35.||Tilg H,, Moschen AR,, Kaser A,, Pines A,, Dotan I, (Year: 2008) Gut, inflammation and osteoporosis: basic and clinical concepts. Gut57: 684–69418408105|
|36.||Pacifici R,, Rifas L,, Teitelbaum S,, Slatopolsky E,, McCracken R,, et al. (Year: 1987) Spontaneous release of interleukin 1 from human blood monocytes reflects bone formation in idiopathic osteoporosis. Proc Natl Acad Sci U S A84: 4616–46203496597|
|37.||Abrahamsen B,, Bonnevie-Nielsen V,, Ebbesen EN,, Gram J,, Beck-Nielsen H, (Year: 2000) Cytokines and bone loss in a 5-year longitudinal study–hormone replacement therapy suppresses serum soluble interleukin-6 receptor and increases interleukin-1-receptor antagonist: the Danish Osteoporosis Prevention Study. J Bone Miner Res15: 1545–155410934653|
|38.||Watkins BA,, Li Y,, Allen KG,, Hoffmann WE,, Seifert MF, (Year: 2000) Dietary ratio of (n-6)/(n-3) polyunsaturated fatty acids alters the fatty acid composition of bone compartments and biomarkers of bone formation in rats. J Nutr130: 2274–228410958824|
|39.||Weiss LA,, Barrett-Connor E,, von Mühlen D, (Year: 2005) Ratio of n-6 to n-3 fatty acids and bone mineral density in older adults: the Rancho Bernardo Study. Am J Clin Nutr81: 934–93815817874|
|40.||Omoigui S, (Year: 2005) Cholesterol synthesis is the trigger and isoprenoid dependent interleukin-6 mediated inflammation is the common causative factor and therapeutic target for atherosclerotic vascular disease and age-related disorders including osteoporosis and type 2 diabetes. Med Hypotheses65: 559–56915935563|
|41.||Parhami F,, Mody N,, Gharavi N,, Ballard AJ,, Tintut Y, (Year: 2002) Role of the cholesterol biosynthetic pathway in osteoblastic differentiation of marrow stromal cells. J Bone Miner Res17: 1997–200312412807|
|42.||You L,, Sheng ZY,, Tang CL,, Chen L,, Pan L,, et al. (Year: 2011) High cholesterol diet increases osteoporosis risk via inhibiting bone formation in rats. Acta Pharmacol Sin32: 1498–150422036861|
|43.||Ma YL,, Bryant HU,, Zeng Q,, Palkowitz A,, Jee WS,, et al. (Year: 2002) Long-term dosing of arzoxifene lowers cholesterol, reduces bone turnover, and preserves bone quality in ovariectomized rats. J Bone Miner Res17: 2256–226412469920|
|44.||Whitney C,, Warburton DE,, Frohlich J,, Chan SY,, McKay H,, et al. (Year: 2004) Are cardiovascular disease and osteoporosis directly linked. Sports Med34: 779–80715462612|
|45.||Jinno Y,, Ayukawa Y,, Ogino Y,, Atsuta I,, Tsukiyama Y,, et al. (Year: 2009) Vertical bone augmentation with fluvastatin in an injectable delivery system: a rat study. Clin Oral Implants Res20: 756–76019302395|
|46.||Spencer TA,, Chai H,, Fu W,, Ramaswami G,, Cox MW,, et al. (Year: 2004) Estrogen blocks homocysteine-induced endothelial dysfunction in porcine coronary arteries. J Surg Res118: 83–9015093721|
|47.||McLean RR,, Jacques PF,, Selhub J,, Tucker KL,, Samelson EJ,, et al. (Year: 2004) Homocysteine as a predictive factor for hip fracture in older persons. N Engl J Med350: 2042–204915141042|
|48.||Koh JM,, Lee YS,, Kim YS,, Kim DJ,, Kim HH,, et al. (Year: 2006) Homocysteine enhances bone resorption by stimulation of osteoclast formation and activity through increased intracellular ROS generation. J Bone Miner Res21: 1003–101116813521|
|49.||Dimitrova KR,, DeGroot K,, Myers AK,, Kim YD, (Year: 2002) Estrogen and homocysteine. Cardiovasc Res53: 577–58811861028|
|50.||Bayhan I,, Uygur D,, Ugurlu N,, Ozaksit G, (Year: 2009) Strontium ranelate decreases plasma homocysteine levels in postmenopausal osteoporotic women. Rheumatol Int29: 263–26618818926|
|51.||De Leo V,, la Marca A,, Morgante G,, Lanzetta D,, Setacci C,, et al. (Year: 2001) Randomized control study of the effects of raloxifene on serum lipids and homocysteine in older women. Am J Obstet Gynecol184: 350–35311228485|
|52.||Ferrandez L,, Martin M,, Fernandez M, (Year: 1990) Calcitonin, estradiol, and hydroxyproline as parameters in the early diagnosis of involutional osteoporosis. The importance of the “second calcitonin phenomenon”. Arch Orthop Trauma Surg109: 181–1852383444|
|53.||Schmolke B, (Year: 2001) Laboratory diagnosis of osteoporosis. Orthopade30: 425–43611515180|
Table 1 Body weight, BMD, Uterine index, and Serum estradiol from the OVX and the Sham rats.
|Parameters||The Sham group||The OVX group|
|Initial body weight (g)||252.62±11.47||250.72±16.35|
|Final body weight (g)||326.43±27.66||403.51±32.76**|
|The 12th week BMD(g/cm2)||0.224±0.019||0.201±0.010*|
|The 24th week BMD(g/cm2)||0.236±0.017||0.207±0.015*|
|Initial BMD(g/cm2)||lumbar vertebrae||0.202±0.011||0.201±0.015|
|The 12th week BMD(g/cm2)||0.213±0.016||0.188±0.017*|
|The 24th week BMD(g/cm2)||0.221±0.022||0.185±0.012**|
|Uterine index (mg/g)||2.312±0.252||0.532±0.124**|
|Serum estradiol (pk/mL)||35.53±6.45||13.53±5.47**|
Values with a superscript are significantly different from the Sham group (*P<0.05, **P<0.01).
Table 2 Fold changes and P value of endogenous metabolites detected by GC/TOF-MS.
|Metabolites||Retentionindex||Biologicalrole||12 weeksaftersurgery||24 weeks after surgery|
|Cholesterol||3178.2||a sterol and a lipid found in the cell membranes||2.02||0.047||1.53||0.006|
|Octadecadienoic acid||2216.4||an unsaturated n-6 fatty acid and essential fatty acid||1.19||0.533||1.85||0.005|
|Arachidonic acid||2375.2||a polyunsaturated, essential fatty acid and precursor ofprostaglandins, thromboxanes, and leukotrienes||2.25||0.005||3.29||0.002|
|Oleic acid||2215.3||An unsaturated fatty acid||1.09||0.793||1.72||0.013|
|9-(Z)-Hexadecenoic Acid||2029.7||an unsaturated fatty acid, a common constituentof the glycerides of adipose tissue||1.29||0.031||1.36||0.009|
|Hexadecanoic Acid||2049.2||a saturated fatty acid||−1.16||0.401||1.54||0.010|
|Octadecanoic Acid||2245.6||a saturated fatty acid||−1.10||0.340||1.18||0.047|
|Docosahexaenoic Acid||2566.9||an omega-3 fatty acid and associated with the heart disease||−1.22||0.350||−1.54||0.001|
|Glycerol||1282.4||an important component of triglycerides and of phospholipids, energy metabolism||1.18||0.362||1.43||0.013|
|Isoleucine||1305.8||Essential branched-chain aliphatic amino acid associated with energy metabolismand blood sugar regulation||1.24||0.012||1.78||0.018|
|homocysteine||1437.9||sulfur-containing amino acid, associated with cardiovasculardisease, Alzheimer’s disease and fracture||1.53||0.042||1.98||0.044|
|L-Hydroxyproline||1522.4||a common non-proteinogenic amino acid and Paget’s disease||2.20||0.006||2.43||0.005|
|Tryptophan||2244.0||an essential amino acid which is the precursor of serotonin||1.03||0.780||−1.44||0.001|
|Glutamine||1608.8||energy metabolism, nitrogen balances||1.01||0.969||−1.24||0.022|
|Cystine||1433.9||immune function and the synthesis of glutathione||1.61||0.030||1.44||0.226|
|Glycine||1438.4||a simple, nonessential amino acid||1.27||0.039||1.15||0.837|
|Alanine||1107.5||An important participant and regulator in glucose metabolism||−2.79||0.005||−2.77||0.017|
|Succinic Acid||1323.2||citric acid cycle, Glucose metabolism||−2.47||0.022||−1.14||0.814|
|Glucose||1935.6||a monosaccharide containing six carbon atoms||1.15||0.010||1.06||0.318|
|L-Lysine||1720.8||an essential amino acid||−1.89||0.011||−1.76||0.002|
|Ethanedioic Acid||1144.9||glyoxylic acid or ascorbic acid metabolism||1.79||0.004||1.25||0.302|
|2-Hydroxy-Butanoic Acid||1137.3||an organic acid that is involved in propanoate metabolism||−6.88||0.012||−2.26||0.007|
|Dodecanoic Acid||1656.1||a component of triglycerides||−1.54||0.023||1.24||0.026|
|2-Aminobutyric Acid||1185.5||a key intermediate in the biosynthesis of ophthalmic acid or ophthalmate.||−3.16||0.012||−3.67||0.003|
|3-Hydroxybutyric- Acid||1172.8||ketone body||−1.24||0.425||1.43||0.010|
|Glycerol-3-Phosphate||1787.1||a chemical intermediate in the glycolysis metabolic pathway||1.10||0.320||1.38||0.044|
Previous Document: In Vivo MR Imaging of Intraarterially Delivered Magnetically Labeled Mesenchymal Stem Cells in a Can...
Next Document: Analysis of the enzymatic properties of a broad family of alanine aminotransferases.