|Relationship of age and body mass index to the expression of obesity and osteoarthritis-related genes in human meniscus.|
|Jump to Full Text|
|PMID: 23318714 Owner: NLM Status: MEDLINE|
|OBJECTIVE: Aging and obesity contribute to the initiation and progression of osteoarthritis with little information on their relation to gene expression in joint tissues, particularly the meniscus. Here, we test the hypothesis that patient age and body mass index (BMI) correlate with the expression of osteoarthritis- and obesity-related gene signatures in the meniscus.
DESIGN: Meniscus was obtained from patients (N=68) undergoing arthroscopic partial meniscectomy. The mRNA expression of 24 osteoarthritis-related and 4 obesity-related genes in meniscus was assessed by quantitative real-time PCR. The relationship between gene expression and patient age and BMI was analyzed using Spearman's rank-order correlation. Hierarchical cluster dendrogram and heat map were generated to study inter-gene associations.
RESULTS: Age was negatively correlated (P<0.05) with the expression of MMP-1 (r=-0.447), NFκB2 (r=-0.361), NFκBIA (r=-0.312), IκBA (r=-0.308), IL-8 (r=-0.305), ADAMTS-4 (r=-0.294), APLN (apelin) (r=-0.250) and IL-6 (r=-0.244). Similarly, BMI was negatively correlated with the expression of APLN (r=-0.328), ACAN (r=-0.268) and MMP-1 (r=-0.261). After adjusting for the correlation between age and BMI (r=0.310; P=0.008), the only independent effect of BMI on gene expression was for APLN (r=-0.272). However, age had an independent effect on the expression on ADAMTS-4 (r=-0.253), MMP-1 (r=-0.399), IL-8 (r=-0.327), COL1A1 (r=-0.287), NFκBIA (r=-0.278), NFκB2 (r=-0.312) and IκBA (r=-0.299). The gene correlation analysis identified four clusters of potentially relevant genes: transcription factors, matrix-degrading enzymes, cytokines and chemokines, and obesity genes.
CONCLUSION: Age and BMI were negatively correlated with several osteoarthritis- and obesity-related genes. Although the bulk of these changes appeared to be driven by age, expression of APLN was related to BMI. Inter-gene correlation analysis implicated a common role for strongly correlated genes. Although age-related variations in gene expression appear to be more relevant than obesity-related differences for the role of the meniscus in osteoarthritis development, further investigation into the role of APLN in meniscus and joint health is warranted.
|M F Rai; L J Sandell; J M Cheverud; R H Brophy|
Related Documents :
|9479504 - Chromosomal localization of the human diadenosine 5',5"'-p1,p4-tetraphosphate pyrophosp...
20513434 - Molecular maps of the reorganization of genome-nuclear lamina interactions during diffe...
10585774 - Cloning and mapping of a novel human serum/glucocorticoid regulated kinase-like gene, s...
10944454 - Mapping and determination of the cdna sequence of the erc gene preferentially expressed...
18420194 - The specific contribution of hypoxia-inducible factor-2alpha to hypoxic gene expression...
17553164 - Androgen-regulated genes differentially modulated by the androgen receptor coactivator ...
|Type: Journal Article; Research Support, N.I.H., Extramural; Research Support, Non-U.S. Gov't Date: 2013-01-15|
|Title: International journal of obesity (2005) Volume: 37 ISSN: 1476-5497 ISO Abbreviation: Int J Obes (Lond) Publication Date: 2013 Sep|
|Created Date: 2013-09-10 Completed Date: 2014-05-05 Revised Date: 2014-07-10|
Medline Journal Info:
|Nlm Unique ID: 101256108 Medline TA: Int J Obes (Lond) Country: England|
|Languages: eng Pagination: 1238-46 Citation Subset: IM|
|APA/MLA Format Download EndNote Download BibTex|
Aging* / metabolism
Body Mass Index*
Cartilage, Articular / metabolism*
Gene Expression Profiling
Gene Expression Regulation
I-kappa B Proteins / metabolism
Intercellular Signaling Peptides and Proteins / metabolism*
Interleukin-8 / metabolism
Matrix Metalloproteinase 1 / metabolism
Menisci, Tibial / injuries, metabolism*
NF-kappa B p52 Subunit / metabolism
Obesity / epidemiology, genetics, metabolism*
Osteoarthritis / epidemiology, genetics, metabolism*
Procollagen N-Endopeptidase / metabolism
Protein Array Analysis
Real-Time Polymerase Chain Reaction
United States / epidemiology
|P30 AR057235/AR/NIAMS NIH HHS; P30-AR057235/AR/NIAMS NIH HHS; R01 AR036994/AR/NIAMS NIH HHS; R01 AR063757/AR/NIAMS NIH HHS; R01-AR036994/AR/NIAMS NIH HHS; T32 AR060719/AR/NIAMS NIH HHS; T32-AR060719/AR/NIAMS NIH HHS|
|0/APLN protein, human; 0/I-kappa B Proteins; 0/Intercellular Signaling Peptides and Proteins; 0/Interleukin-8; 0/NF-kappa B p52 Subunit; 0/NFKB2 protein, human; 139874-52-5/NF-kappaB inhibitor alpha; EC 3.4.24.-/ADAM Proteins; EC 3.4.24.-/aggrecanase-1; EC 126.96.36.199/Procollagen N-Endopeptidase; EC 188.8.131.52/Matrix Metalloproteinase 1|
Journal ID (nlm-journal-id): 101256108
Journal ID (pubmed-jr-id): 32579
Journal ID (nlm-ta): Int J Obes (Lond)
Journal ID (iso-abbrev): Int J Obes (Lond)
nihms-submitted publication date: Day: 22 Month: 12 Year: 2012
Electronic publication date: Day: 15 Month: 1 Year: 2013
Print publication date: Month: 9 Year: 2013
pmc-release publication date: Day: 01 Month: 3 Year: 2014
Volume: 37 Issue: 9
First Page: 1238 Last Page: 1246
PubMed Id: 23318714
|Relationship of age and body mass index to the expression of obesity and osteoarthritis related genes in human meniscus|
|Muhammad Farooq Rai, Ph.D.1|
|Linda J. Sandell, Ph.D.12|
|James M. Cheverud, Ph.D.3|
|Robert H. Brophy, M.D.4*|
1Department of Orthopaedic Surgery, Washington University School of Medicine at Barnes-Jewish Hospital, 660 S. Euclid Ave. MS 8233, St. Louis MO 63110 United States
2Cell Biology and Physiology, Washington University School of Medicine at Barnes-Jewish Hospital, 660 S. Euclid Ave. MS 8233, St. Louis MO 63110 United States
3Anatomy and Neurobiology, Washington University School of Medicine at Barnes-Jewish Hospital, 660 S. Euclid Ave. MS 8233, St. Louis MO 63110 United States
4Department of Orthopaedic Surgery, Washington University School of Medicine at Barnes-Jewish Hospital, 14532 South Outer Forty Drive, Chesterfield, MO 63017 United States
|*Corresponding author, Robert H. Brophy, M.D., Department of Orthopaedic Surgery, Washington University School of Medicine at Barnes-Jewish Hospital, 14532 South Outer Forty Drive, Chesterfield, MO 63017 United States Ph: 314-514-3564; Fax: 314-514-3689; email@example.com
Osteoarthritis is a degenerative joint disorder with severe consequences in terms of both disability and the financial burden in society1, 2. Today, osteoarthritis is not considered to be merely a disease of cartilage but a disease of the whole joint and its related tissues3. The pathologic changes in the arthritic knee are not confined only to articular cartilage (degradation, calcification) but also affect subchondral bone (thickening, osteophyte formation), synovium (synovitis), joint capsule (hypertrophy), ligaments (degeneration) and meniscus (extrusion, degeneration, calcification)4, 5.
The meniscus (pl. menisci) is a semi-lunar biconcave fibrocartilaginous disk that resides within the medial and lateral tibiofemoral articulations (i.e. interposed between the femoral condyles and tibial plateaux) and contributes to tibiofemoral congruence, joint lubrication, stability and load distribution6–10. In each knee joint, there are two menisci: a medial and a lateral and both are attached to bone via insertional ligaments and provide protection to underlying articular cartilage. The menisci increase the contact area in the femorotibial joint and thus decrease the stress on the articular cartilage6, 11, 12. Dispersal of synovial fluid over and through the menisci decreases friction across articular surfaces and contributes to joint lubrication10. Normal joint movement causes compression of the menisci and the resulting extrusion of synovial fluid bathes the articular cartilage to provide nutrition13.
A meniscus tear may occur as a consequence of trauma, twisting of the knee, deceleration and landing from a jump or getting up from a squatting or crouching position14, 15. Younger patients are more likely to have acute, traumatic meniscal tears whereas older patients are more likely to have degenerative tears16. It is believed that meniscal tears are associated with osteoarthritis development8, 17–21, but it is unclear how degenerative changes in the menisci affect cartilage homeostasis22–25.
A tear is usually treated by arthroscopic repair or partial meniscectomy. While meniscal repair is preferred when possible26, partial meniscectomy is the most commonly performed procedure 16 and is currently preferred over total meniscectomy27, 28, particularly since total meniscectomy has been associated with more cartilage degeneration and osteoarthritis than partial meniscectomy29, 30.
Several factors contribute to the initiation and progression of primary knee osteoarthritis including aging and obesity31–33. Aging undoubtedly contributes to osteoarthritis development through cartilage degeneration 34 but the explicit mechanisms for this contribution are not fully understood. Aging may contribute to the development of osteoarthritis through age-related changes at cellular and tissue levels making the joint more susceptible to damage with a decreased ability to maintain homeostasis5. Similarly, obesity, although a modifiable risk factor, represents another predictor of osteoarthritis35. Although the association between obesity and osteoarthritis is complex, multiple biochemical mediators (cytokines, chemokines and adipokines) that crosstalk with obesity and osteoarthritis are reported to contribute to this relationship36–38.
Our previous findings demonstrated that the expression of osteoarthritis-related genes in a torn meniscus decreases with age without sex-based differences39. In the literature, several studies have identified significant sexual dimorphism in the development and progression of knee osteoarthritis40, 41. However, there is limited information available on sex-based differences in gene expression in meniscus other than our recent finding of no substantial differences between women and men in the expression of osteoarthritis-related genes.
Little is known how gene expression signatures of a torn meniscus relate to osteoarthritis in the context of aging and obesity. Therefore, in this study we aim to gain insight into the interaction of osteoarthritis- and obesity-related gene signatures in the meniscus. Here, we test the hypothesis that gene expression patterns in menisci correlate with the age and body mass index (BMI) of patients. These findings have the potential to shed further light on how meniscal injury, age and obesity relate to the initiation and progression of osteoarthritis in the knee.
All procedures were approved by the Washington University School of Medicine human subjects institutional review board. Informed consent was obtained about the use of tissues for research purpose as approved by Human Research Protection Office.
A total of 68 meniscal samples were collected at the Department of Orthopaedic Surgery (Barnes-Jewish Hospital, Chesterfield, MO), between June 2010 and January 2012, at the time of arthroscopic partial meniscectomy for a symptomatic tear in their meniscus. Demographic data for age, sex, height and body weight was recorded and BMI was calculated (Table-1). None of the patients had advanced osteoarthritis or any ligament injury at the time of the meniscal surgery. Eligible patients indicated for partial meniscectomy voluntarily consented to participate in this study. Appropriate resection of non-repairable meniscus tears was performed and samples of the torn tissue were removed from the knee and transported to the laboratory in normal saline. We recorded whether the sample came from the inner half of the meniscus, the outer half of the meniscus or both.
The tissues, upon arrival in the lab, were washed with phosphate-buffered saline (PBS; HyClone; Thermo Fisher Scientific, Rockford, IL) and weighed. TRIzol reagent (Invitrogen, Carlsbad, CA) was added to each specimen (1mL/50–100 mg), which were then stored at −80°C prior to total RNA extraction.
The thawed tissues were homogenized with use of a Polytron homogenizer (Brinkmann Instruments, Westbury, NY). Total RNA was isolated using a combination of phenol chloroform method and RNeasy mini spin columns (Qiagen Inc., Valencia, CA). Briefly, aliquots of the homogenized suspension were transferred to microfuge tubes and incubated at room-temperature for 5 min. Chloroform (200 µl) was added to each tube, vortexed vigorously, incubated for 5 min and centrifuged at 12000g for 15 min at 4°C. The upper aqueous phase, that contained RNA, was transferred to a clean tube and an equal volume of 70% RNase-free ethanol was added and the contents were mixed. The aqueous phase/ethanol mixture was applied to spin columns and spun for 15 sec at 8000 rpm at room-temperature. After decanting the flow-through, the column was washed twice with RW1 and RPE buffers. Finally, RNA was eluted in 30 µl of RNase-free water, quantified by a NanoDrop (Thermo Fisher Scientific, Pittsburg, PA) and stored at −80°C until reverse transcribed for complementary DNA (cDNA).
Total RNA was first subjected to DNase I (Life Technologies, Grand Island, NY) digestion to avoid contamination from genomic DNA. In 0.2-mL tubes, 150–300 ng of RNA, 2 units of amplification grade DNase I (1 U/µl), and 2 µl of 10X DNase I reaction buffer were added to a total 20-µl volume in RNase-free water and incubated at room-temperature for 15 min. DNase I was inactivated with 2 µl of 25 mM EDTA at 65°C for 10 min.
First strand cDNA synthesis was carried out by use of random hexamers and the SuperScript II First-Strand Synthesis System (Invitrogen, Carlsbad CA). In a 40-µl reaction volume, 2 µl of the random primer (50 ng/µl) and 2 µl of dNTP mix (10 mM) were added to DNase-I-treated RNA and incubated at 65°C for 5 min. The reaction mixture was chilled on ice and reverse-transcription was continued by adding 4 µl of 0.1M DTT (dithiothreitol), 400 units of SuperScript II (200 units/µl) and 5X first-strand buffer in the reaction mixture. Finally, the tube was incubated at 42°C for 50 min followed by deactivation at 70°C for 15 min.
Several genes representative if pathways involved in osteoarthritis pathophysiology were selected including proinflammatory cytokines (IL-1α, IL-1β, IL-6, TNFα), chemokines (IL-8, CCL3, CCL3L1, CCL20, CXCL1, CXCL3, CXCL6), matrix degrading enzymes [aggrecanases (ADAMTS-4, ADAMTS-5), matrix metalloproteinases (MMP-1, MMP-3, MMP-9, MMP-13)], matrix genes (COL1A1, COL2A1, ACAN, BMP-2) and transcription factors (NFκBIA, NFκB2, IκBA). The sequence information for these primers has been previously published39. We also studied selected markers for obesity including adiponectin (ADIPOQ), apelin (APLN), leptin (LEP) 42 and resistin (RETN)43 (see respective references for primer sequences). All primers were obtained from Invitrogen, Carlsbad, CA.
Quantitative real-time PCR (qPCR) was performed on a 7500 Fast Real-Time PCR System (Applied Biosystems, Carlsbad, CA) on cDNA samples to assess expression levels of the genes listed above. The housekeeping gene glyceraldehyde 3-phosphate dehydrogenase (GAPDH) acted as an endogenous reference gene for normalization of fluorescence threshold (Ct) values of target genes. In a 20-µl reaction volume, 10 µl of SYBR Green PCR Master Mix (Applied Biosystems, Carlsbad, CA), 1.5 µl of cDNA, and 200 nM of primers were added. Samples were amplified with an initial activation step at 95°C for 10 min, followed by forty cycles of denaturation at 95°C for 15 sec and annealing at 60°C for 60 sec. The Ct values for house-keeping gene and the genes of interest were measured and normalized to GAPDH for each sample (ΔCt). Rather than employing the traditional ΔΔCt method to identify increased or decreased expression of samples relative to a single reference sample and then comparing candidate genes to the multiple control genes, a regression strategy was used to correct all samples for interfering factors at once. This strategy eliminates amplification of error that would be created by introducing ratios when comparing candidate gene expression to control gene expression. Ratios of raw Ct values for the candidate gene relative to the Ct value for GAPDH removes the effect of total cellular material from inter-sample comparisons only if the slope of the regression between the scores is 1.0; otherwise, the ratio remains correlated with its denominator.
In order to determine the associations among gene expression levels, the values were subjected to a hierarchical cluster analysis using Ward Minimum Variance Method to produce a dendrogram. A manual heat map was then generated to display the pairwise correlations of candidate genes across all gene categories. To achieve this, a color code was assigned to the correlation strengths among various genes and finally the heat map was aligned with the hierarchical cluster dendrogram.
All statistical analyses were performed using Systat-12 software (Systat, Chicago, IL.). Since the distribution of the gene expression values deviated from a normal distribution and contained several extreme values that would dominate any ordinary least squares linear regression analysis, Spearman's rank-order correlation was used for examining the associations between BMI and age and gene expression values on an ordinal scale.
The correlation between BMI and age is significantly positive, so we used rank-order partial correlation analysis to test for correlations between gene expression and BMI independent of age and between gene expression and age independent of BMI. The first-order partial correlation coefficient is given by
Probability values for the Spearman’s partial correlations were computed by treating
The demographic information of the study population, from which meniscal samples were obtained, is presented in Table-1. The tissues were collected from 68 patients (20 females and 48 males) at the time of arthroscopic surgery. The mean age of patients was 45.3 years with a range of 14–72 years with BMI 27.70±5.17 kg/m2 (mean ± S.D.).
A one-way analysis of variance was first used to determine sex-based differences between 20 female and 48 male patients. The analysis showed no sexual dimorphism in gene expression (data not shown), therefore our results for age and BMI are based on combined data set (N=68) from both female and male subjects. The samples came predominantly from the inner rim (N=46) while only one sample came from the outer meniscus in isolation and the rest included a combination of the inner and outer zones of the meniscus. This distribution precluded any analysis based on zonal differences in gene expression.
The present study confirms our previous 39 that ostearthritis- and obesity-related gene expression is significantly negatively correlated with age. The results from each gene category detailed below are depicted in Fig. 2 and Supplementary Table-1.
Although, the mRNA expression of various proinflammatory cytokines including IL-1α (r=−0.219), IL-1β (r=−0.119), IL-6 (r=−0.243) and TNFα (r=−0.151) were negatively correlated with age, only IL-6 expression was statistically significant (r=-0.244; P=0.045).
Chemokine expression levels are negatively correlated with age but the only statistically significant negative correlation was found to be between age and IL-8 expression (r=-0.31; P=0.011).
Both ADAMTS-4 (r=−0.294) and MMP-1 (r=−0.447) were strongly significantly negatively associated with age (P<0.05) while ADAMTS-5 and other MMPs levels were negatively correlated with age without being statistically significant.
NFκB2 (r=−0.361), NFκBIA (r=−0.312) and IκBA (r=−0.308) were all significantly negatively correlated with age (P<0.05).
The quantification of extracellular matrix genes showed that COL1A1 (r=−0.234), COL2A1 (r=−0.030), ACAN (r=−0.139) and BMP-2 (r=−0.144) were negatively correlated with age at a sub-significant level. COL1A1 had a borderline significant level of correlation (P=0.054).
The quantification of genes related to obesity showed that their expression was also negatively correlated with age. Gene expression for APLN (r=−0.250) was significantly negatively correlated with age (P=0.040) while that of ADIPOQ (r=−0.250) was at a borderline significance level (P=0.052) and LEP (r=−0.160) and RETN (r=−0.202) showed no statistical significance.
We found that a majority of the selected genes were also negatively correlated with BMI (Figure 3, Supplementary Table 2). There were no significant correlations for proinflammatory cytokines, chemokines, or aggrecanases. The expression of all MMPs except for MMP-13 was also found to be negatively correlated with BMI with only MMP-1 (r=−0.261) reaching statistical significant (P=0.032). For transcription factors, only NFκB2 (r=−0.183) was of borderline significance (P=0.053). The expression of ACAN (r=−0.267) was negatively correlated with BMI (P=0.027). Among obesity related genes studied here only APLN (r=−0.328) was found to be significantly negatively correlated (P=0.006) with BMI.
Age and BMI were significantly positively correlated (r=0.310; P=0.008) in this cohort. This correlation indicated that majority of the older patients were obese (Supplementary Fig.1) and vice versa.
Rank-order partial correlation of gene expression level with age, controlling for BMI, was often significant and negative (Table 2). Among these are ADAMTS-4 (r=−0.253; P=0.038), MMP-1 (r=−0.399; P=0.0008), IL-8 (r=−0.326; P=0.007), COL1A1 (r=−0.287; P=0.018), NFκBIA (r=−0.278; P=0.023), NFκB2 (r=−0.312; P=0.01), IκBA (r=−0.299; P=0.014) and ADIPOQ (r=−0.224; P = 0.068) at borderline significance. However, when BMI is considered independently from age, only APLN (r=−0.224) remained significantly negatively correlated with BMI (P=0.027). Other gene expression levels that were originally negatively correlated with BMI are no longer significant after controlling for age. These results indicate that age effects on gene expression in the meniscus are stronger than the effects of BMI on gene expression. Most of the apparent negative correlation between gene expression levels and BMI can be accounted for by age differences.
Gene expression levels were correlated pairwise for all candidate genes to identify potentially common regulatory mechanisms based on concordant or discordant expression patterns (Fig. 4). We observed strong inter-gene correlations among a majority of cytokine/chemokine genes. There was also a strong inter-correlation among various transcription factors involved in most if not all cytokine/chemokine pathways in osteoarthritis pathophysiology. In addition, weak inter-gene correlations were found between matrix genes (collagens) and matrix degrading enzymes (metalloproteinases and aggrecanases).
In this study, we sought to determine the gene expression signatures in torn menisci (obtained from patients undergoing arthroscopic partial meniscectomy) and to study their correlation with patients’ age and BMI. Gene expression was generally negatively correlated with both age and BMI. Although most of the apparent negative correlation between gene expression levels and BMI can be accounted for by age differences, this is the first study we are aware of to identify differences in gene expression in the injured meniscus which relate to patient BMI.
The published evidence suggests a relationship between obesity and osteoarthritis36–38. Dysregulation of lipid homeostasis is one of the mechanisms potentially leading to osteoarthritis development 37. Among the twenty-eight genes analyzed in this study, six genes (IL-8, MMP-13, CXCL1, CCL20, COL1A1 and COL2A1) had an insignificant positive correlation with BMI while expression of the remaining twenty-two genes was negatively correlated with BMI, with only three genes (APLN, ACAN and MMP-1) significantly so. Only the expression of APLN was exclusively due to obesity status independent of age. As stated above, APLN is more closely related to obesity than osteoarthritis, but its precise role remains uncertain44. However, it is reported that APLN expression is regulated by insulin45 and might be involved in regulation of glucose homeostasis and obesity46. Since APLN is known to play a catabolic role in articular cartilage44 and bone47 metabolisms, APLN may have a catabolic role in meniscal homeostasis, particularly in obese patients, and may constitute a risk factor for the development of osteoarthritis. The expression of APLN in the injured menisci of normal weight individuals is presumably appropriate, although this cannot be established without comparison to uninjured menisci.
What is less clear is how APLN interacts with meniscal and other joint tissue, and how the magnitude of APLN expression in the torn meniscus relates to joint health. Perhaps paradoxically, the torn menisci of obese patients express less APLN than torn menisci of patients with a lower BMI. At least one study has shown worse long term clinical outcomes for patients with elevated BMI after partial lateral meniscectomy in otherwise normal knees48. The decreased expression of APLN in obese patients suggests a biological pathway by which obesity affects the knee joint and may influence the risk of developing osteoarthritis. Future investigation should look at how the expression of APLN affects meniscal and other joint tissue. Another important question is whether expression of APLN is plastic and could be modified by future weight loss (or gain).
Aging is known to affect the basal-level gene expression pattern in joint tissues including the meniscus49. We have also reported that gene expression in meniscus tears varies by age39, perhaps because of lower cellularity in torn meniscus in patients younger than forty years50. Here, we validate our results with increased sample size and with more statistical power to show that for most genes, expression levels are negatively associated with a patient’s age. A moderate significant negative correlation of MMP-1, NFκB2, NFκBIA, IκBA, IL-8, ADAMTS-4, IL-6 and APLN (apelin) expression with age was observed. Several of these genes were previously found to show less expression in older patients. Except for APLN, which is an obesity-related gene, all other genes are involved in osteoarthritis pathogenesis and their expression pattern is an important consideration in the molecular characterization of osteoarthritis51–53. The precise role of APLN in osteoarthritis is not well-studied, but it is suggested that it may be involved in cartilage metabolism and constitutes a risk factor in the pathophysiology of osteoarthritis44. In parallel, although not statistically significant, the expression of other genes (Fig. 2; Supplementary Table-1) was also negatively correlated with age. It was interesting to note that none of the osteoarthritis- or obesity-related genes were significantly positively correlated with age. Taken together, these findings demonstrate that an elevated intrinsic response exists in torn menisci from younger individuals compared to older individuals. These results also support our previous findings where it was reported that osteoarthritis-related proinflammatory and catabolic genes are downregulated in older patients39. In addition, it has been shown that there is a greater prevalence of traumatic tears in younger patients18 that are associated with an increased risk for initiation and perpetuation of osteoarthritis8, 17–19. Thus, our findings provide a molecular rationale for the increased risk of developing osteoarthritis in young individuals with a meniscus tear.
The gene expression pattern as it relates to aging and obesity has clinical relevance in several ways. For example, negative correlation of gene expression with age may relate to the reduced long-term clinical outcomes seen in older patients after meniscal repair54. Although short-term clinical success rates are good in older individuals55, these patients have reportedly lower rates of meniscal healing and higher long-term failure rates post-repair. As a lack of viable cells in the meniscus is associated with degeneration and repeat tearing56, one plausible physiologic reason for higher late failure rates in older patients may be the decreased intrinsic and perimeniscal cellularity found in the torn menisci of this patient group50. Differences in gene expression and intrinsic inflammation may contribute to this difference in healing potential as well. The fact that some of the genes are negatively correlated with obesity independent of age (Table-2) indicates that factors other than obesity and age likely play a role in differential gene expression patterns. Among these factors, physical activity level could be the most significant because it has been reported that a high percentage of obese individuals with knee osteoarthritis are sedentary suggesting that a lack of physical activity may increase the susceptibility to cellular inflammation57, 58.
Osteoarthritis is a sexually dimorphic trait more common among females than males40, 41, thus giving rise to the concept that sex is a contributing factor for osteoarthritis development. This might lead one to expect higher expression of osteoarthritis-related genes in women than men following a meniscal tear. Our previous findings provided weak evidence in support of this hypothesis where we detected that only CCL3L1 was significantly upregulated in women compared to men39. In the present study, however, we failed to appreciate any significant difference in the gene expression profile between the two sexes even with greater sample size and more statistical power. It is unclear what the clinical implications of these findings may be in this context but these differences are unlikely to contribute to the greater prevalence of knee osteoarthritis in women.
The co-expression analysis and its ability to define distinct physiological sets of genes are portrayed by hierarchical cluster dendrogram and heat map (Fig. 4). The genes with similar physiological role are clustered together with varying degree of correlation strengths. These findings implicate that common regulatory mechanism for genes or group of genes may exist. For instance, we found that cytokine and chemokines were strongly inter-correlated similar to transcription factors. In contrast, obesity-related genes were weakly correlated with matrix genes and matrix-degrading enzymes.
One limitation of this study is the inability to test for differences in gene expression between the inner and outer zones of the meniscus. The meniscus contains several different cell types that exhibit different gene expression patterns and are highly localized to the inner vs outer regions59, 60. While the majority of our samples are from the inner meniscus, potential regional variation could influence our results.
Another major limitation of this study is the lack of data on gene expression in normal menisci and its potential relationship to age and BMI. While it would be interesting to have data on expression from normal meniscal tissue, such information is not essential for a couple of reasons. First and foremost, this study focuses on the injured meniscus. The intention is not to compare the injured meniscus to the uninjured meniscus but rather to identify BMI and age related differences in gene expression of injured menisci. Patients with meniscus tears are the population of interest and the impact of BMI and age on gene expression is potentially very clinically relevant in terms of the prognosis for these patients. Given the heightened risk of developing osteoarthritis in patients with meniscus tears, this information may help stratify this risk independently of the potential relationship of gene expression in torn menisci to gene expression in normal menisci. Secondly, age and BMI may have similar or dissimilar effects on gene expression in the injured meniscus compared to the uninjured and that question is not the focus of our present investigation. Since sampling normal meniscal tissue is not entirely benign, we do not consider it an ethical approach in patients with intact menisci. Considering that measuring gene expression in normal meniscal tissue is not paramount or practical, the relationship, if any, between gene expression of normal meniscus tissue and torn meniscus tissue is interesting but not integral to the current investigation.
Since aging and obesity both contribute to development of knee osteoarthritis, we predicted that the expression pattern of genes of interest will help understand tissue-specific responses to age and obesity. We found that patients’ age and BMI were negatively correlated with the expression of several osteoarthritis- and obesity-related genes in the meniscus. While a bulk of these changes appear to be driven by age, expression of APLN in the meniscus is related to BMI. Further investigation into the role of APLN in meniscus and joint health is warranted. These findings also suggest that age-related changes are more relevant than obesity-related differences for the role of the meniscus in the pathogenesis of osteoarthritis. The potential effect of other factors such as physical activity, ethnicity, smoking, diabetes and inflammation on meniscal gene expression should also be investigated.
2 Click here for additional data file (NIHMS425600-supplement-2.docx)
3 Supplementary Fig. 1. Association between age and BMI
Spearman’s rank-order correlation was calculated to plot the relationship between age and BMI. Solid line show the best fitted correlation curve, while dotted lines show 95% confidence interval based on the curve. In this study, age and BMI (obesity) were found to be strongly positively correlated.
Click here for additional data file (NIHMS425600-supplement-3.docx)
FN1CONFLICT OF INTEREST
Supplementary information is available at IJO’s website.
Funding for this study was provided by an Orthopaedic Research and Education Foundation (OREF) grant (R.H.B.), National Institute of Arthritis and Musculoskeletal and Skin Diseases grant RO1-AR036994 (L.J.S.) and Musculoskeletal Research Center, grant P30-AR057235. M.F.R. is supported by Ruth L. Kirschstein National Research Service Award Fellowship (T32-AR060719) from National Institutes of Health. The content is solely the responsibility of the authors and does not necessarily represent the official views of the National Institute of Arthritis, Musculoskeletal and Skin Diseases or the National Institutes of Health.
|1.||Woolf AD,Pfleger B. Burden of major musculoskeletal conditionsBulletin of the World Health OrganizationYear: 200381964665614710506|
|2.||van Baar ME,Dekker J,Lemmens JA,Oostendorp RA,Bijlsma JW. Pain and disability in patients with osteoarthritis of hip or knee: the relationship with articular, kinesiological, and psychological characteristicsJ RheumatolYear: 19982511251339458215|
|3.||Brandt KD,Radin EL,Dieppe PA,van de Putte L. Yet more evidence that osteoarthritis is not a cartilage diseaseAnn Rheum DisYear: 200665101261126416973787|
|4.||Loeser RF,Goldring SR,Scanzello CR,Goldring MB. Osteoarthritis: A disease of the joint as an organArthritis RheumYear: 20126461697170722392533|
|5.||Loeser RF. Age-related changes in the musculoskeletal system and the development of osteoarthritisClin Geriatr MedYear: 201026337138620699160|
|6.||Walker PS,Erkman MJ. The role of the menisci in force transmission across the kneeClin Orthop Relat ResYear: 19751091841921173360|
|7.||Masouros SD,McDermott ID,Amis AA,Bull AM. Biomechanics of the meniscus-meniscal ligament construct of the kneeKnee Surg Sports Traumatol ArthroscYear: 200816121121113218802689|
|8.||Englund M,Guermazi A,Lohmander LS. The meniscus in knee osteoarthritisRheum Dis Clin North AmYear: 200935357959019931804|
|9.||Chevrier A,Nelea M,Hurtig MB,Hoemann CD,Buschmann MD. Meniscus structure in human, sheep, and rabbit for animal models of meniscus repairJ Orthop ResYear: 20092791197120319242978|
|10.||Macconaill MA. The Function of Intra-Articular Fibrocartilages, with Special Reference to the Knee and Inferior Radio-Ulnar JointsJ AnatYear: 193266Pt 221022717104369|
|11.||Messner K,Gao J. The menisci of the knee joint. Anatomical and functional characteristics, and a rationale for clinical treatmentJ AnatYear: 1998193Pt 21611789827632|
|12.||Renstrom P,Johnson RJ. Anatomy and biomechanics of the menisciClinics in sports medicineYear: 1990935235382199066|
|13.||Mow VC,Holmes MH,Lai WM. Fluid transport and mechanical properties of articular cartilage: a reviewJ BiomechYear: 19841753773946376512|
|14.||Lim SY,Peh WC. Magnetic resonance imaging of sports injuries of the kneeAnnals of the Academy of Medicine, SingaporeYear: 2008374354361|
|15.||Hayes CW,Coggins CA. Sports-related injuries of the knee: an approach to MRI interpretationClinics in sports medicineYear: 200625465967916962421|
|16.||McDermott ID,Amis AA. The consequences of meniscectomyJ Bone Joint Surg BrYear: 200688121549155617159163|
|17.||Englund M. The role of the meniscus in osteoarthritis genesisRheum Dis Clin North AmYear: 200834357357918687273|
|18.||Englund M,Guermazi A,Lohmander SL. The role of the meniscus in knee osteoarthritis: a cause or consequence?Radiologic clinics of North AmericaYear: 200947470371219631077|
|19.||Englund M. Meniscal tear--a feature of osteoarthritisActa orthopaedica Scandinavica. SupplementumYear: 200475312145 backcover.. 15188684|
|20.||Ding C,Martel-Pelletier J,Pelletier JP,Abram F,Raynauld JP,Cicuttini F,et al. Meniscal tear as an osteoarthritis risk factor in a largely non-osteoarthritic cohort: a cross-sectional studyJ RheumatolYear: 200734477678417361984|
|21.||Zeichen J,Hankemeier S,Knobloch K,Jagodzinski M. [Arthroscopic partial meniscectomy]Oper Orthop TraumatolYear: 2006185–638039217171326|
|22.||Noble J,Hamblen DL. The pathology of the degenerate meniscus lesionJ Bone Joint Surg BrYear: 19755721801861173585|
|23.||Casscells SW. The torn or degenerated meniscus and its relationship to degeneration of the weight-bearing areas of the femur and tibiaClin Orthop Relat ResYear: 1978132196200679540|
|24.||Lewandrowski KU,Muller J,Schollmeier G. Concomitant meniscal and articular cartilage lesions in the femorotibial jointAm J Sports MedYear: 19972544864949240982|
|25.||Bennett LD,Buckland-Wright JC. Meniscal and articular cartilage changes in knee osteoarthritis: a cross-sectional double-contrast macroradiographic studyRheumatology (Oxford)Year: 200241891792312154209|
|26.||Brophy RH,Matava MJ. Surgical options for meniscal replacementJ Am Acad Orthop SurgYear: 201220526527222553098|
|27.||Howell JR,Handoll HH. Surgical treatment for meniscal injuries of the knee in adultsCochrane Database Syst RevYear: 20002CD00135310796436|
|28.||Paxton ES,Stock MV,Brophy RH. Meniscal repair versus partial meniscectomy: a systematic review comparing reoperation rates and clinical outcomesArthroscopyYear: 20112791275128821820843|
|29.||Andersson-Molina H,Karlsson H,Rockborn P. Arthroscopic partial and total meniscectomy: A long-term follow-up study with matched controlsArthroscopyYear: 200218218318911830813|
|30.||Hede A,Larsen E,Sandberg H. Partial versus total meniscectomy. A prospective, randomised study with long-term follow-upJ Bone Joint Surg BrYear: 19927411181211732238|
|31.||Aigner T,Richter W. OA in 2011: Age-related OA--a concept emerging from infancy?Nat Rev RheumatolYear: 201282707222231238|
|32.||Lohmander LS,Gerhardsson de Verdier M,Rollof J,Nilsson PM,Engstrom G. Incidence of severe knee and hip osteoarthritis in relation to different measures of body mass: a population-based prospective cohort studyAnn Rheum DisYear: 200968449049618467514|
|33.||Lementowski PW,Zelicof SB. Obesity and osteoarthritisAm J Orthop (Belle Mead NJ)Year: 200837314815118438470|
|34.||Lotz M,Loeser RF. Effects of aging on articular cartilage homeostasisBoneYear: 2012|
|35.||Jiang L,Rong J,Wang Y,Hu F,Bao C,Li X,et al. The relationship between body mass index and hip osteoarthritis: a systematic review and meta-analysisJoint Bone SpineYear: 201178215015520580591|
|36.||Rai MF,Sandell LJ. Inflammatory mediators: tracing links between obesity and osteoarthritisCritical reviews in eukaryotic gene expressionYear: 201121213114222077152|
|37.||Pottie P,Presle N,Terlain B,Netter P,Mainard D,Berenbaum F. Obesity and osteoarthritis: more complex than predicted!Ann Rheum DisYear: 200665111403140517038451|
|38.||Sharma L,Chang A. Overweight: advancing our understanding of its impact on the knee and the hipAnn Rheum DisYear: 200766214114217242017|
|39.||Brophy RH,Rai MF,Zhang Z,Torgomyan A,Sandell LJ. Molecular analysis of age and sex-related gene expression in meniscal tears with and without a concomitant anterior cruciate ligament tearJ Bone Joint Surg AmYear: 201294538539322362494|
|40.||Rosner IA,Goldberg VM,Moskowitz RW. Estrogens and osteoarthritisClin Orthop Relat ResYear: 198621377832430748|
|41.||Verbrugge LM. Women, men, and osteoarthritisArthritis Care ResYear: 1995842122208605259|
|42.||Geiger K,Muendlein A,Stark N,Saely CH,Wabitsch M,Fraunberger P,et al. Hypoxia induces apelin expression in human adipocytesHorm Metab ResYear: 201143638038521448846|
|43.||Nohira T,Nagao K,Kameyama K,Nakai H,Fukumine N,Okabe K,et al. Identification of an alternative splicing transcript for the resistin gene and distribution of its mRNA in human tissueEur J EndocrinolYear: 2004151115115415248836|
|44.||Hu PF,Chen WP,Tang JL,Bao JP,Wu LD. Apelin plays a catabolic role on articular cartilage: in vivo and in vitro studiesInt J Mol MedYear: 201026335736320664951|
|45.||Kleinz MJ,Davenport AP. Emerging roles of apelin in biology and medicinePharmacol TherYear: 2005107219821115907343|
|46.||Rayalam S,Della-Fera MA,Krieg PA,Cox CM,Robins A,Baile CA. A putative role for apelin in the etiology of obesityBiochem Biophys Res CommunYear: 2008368381581918275845|
|47.||Xie H,Yuan LQ,Luo XH,Huang J,Cui RR,Guo LJ,et al. Apelin suppresses apoptosis of human osteoblastsApoptosis : an international journal on programmed cell deathYear: 200712124725417136493|
|48.||Scheller G,Sobau C,Bulow JU. Arthroscopic partial lateral meniscectomy in an otherwise normal knee: Clinical, functional, and radiographic results of a long-term follow-up studyArthroscopyYear: 200117994695211694926|
|49.||Loeser RF,Olex AL,McNulty MA,Carlson CS,Callahan MF,Ferguson CM,et al. Microarray analysis reveals age-related differences in gene expression during the development of osteoarthritis in miceArthritis RheumYear: 201264370571721972019|
|50.||Mesiha M,Zurakowski D,Soriano J,Nielson JH,Zarins B,Murray MM. Pathologic characteristics of the torn human meniscusAm J Sports MedYear: 200735110311217092929|
|51.||Sun Y,Mauerhan DR,Honeycutt PR,Kneisl JS,Norton JH,Hanley EN Jr,et al. Analysis of meniscal degeneration and meniscal gene expressionBMC Musculoskelet DisordYear: 2010111920109188|
|52.||Sandell LJ,Xing X,Franz C,Davies S,Chang LW,Patra D. Exuberant expression of chemokine genes by adult human articular chondrocytes in response to IL-1ßOsteoarthritis CartilageYear: 2008|
|53.||Hellio Le Graverand MP,Vignon E,Otterness IG,Hart DA. Early changes in lapine menisci during osteoarthritis development: Part II: molecular alterationsOsteoarthritis CartilageYear: 200191657211178949|
|54.||Eggli S,Wegmuller H,Kosina J,Huckell C,Jakob RP. Long-term results of arthroscopic meniscal repair. An analysis of isolated tearsAm J Sports MedYear: 19952367157208600740|
|55.||Kotsovolos ES,Hantes ME,Mastrokalos DS,Lorbach O,Paessler HH. Results of all-inside meniscal repair with the FasT-Fix meniscal repair systemArthroscopyYear: 20062213916399454|
|56.||Brindle T,Nyland J,Johnson DL. The meniscus: review of basic principles with application to surgery and rehabilitationJournal of athletic trainingYear: 200136216016916558666|
|57.||Issa RI,Griffin TM. Pathobiology of obesity and osteoarthritis: integrating biomechanics and inflammationPathobiology of aging & age related diseasesYear: 201222012|
|58.||Mork PJ,Holtermann A,Nilsen TI. Effect of body mass index and physical exercise on risk of knee and hip osteoarthritis: longitudinal data from the Norwegian HUNT StudyJ Epidemiol Community HealthYear: 2012|
|59.||Spindler KP,Miller RR,Andrish JT,McDevitt CA. Comparison of collagen synthesis in the peripheral and central region of the canine meniscusClin Orthop Relat ResYear: 19943032562638194243|
|60.||Fuller ES,Smith MM,Little CB,Melrose J. Zonal differences in meniscus matrix turnover and cytokine responseOsteoarthritis CartilageYear: 2012201495922062355|
Keywords: Meniscus tear, gene expression, osteoarthritis, age, body mass index, obesity.
Previous Document: Molecular machines directly observed by high-speed atomic force microscopy.
Next Document: Growth in early life and the development of obesity by age 9 years: are there critical periods and a...