Document Detail

Relationship of age and body mass index to the expression of obesity and osteoarthritis-related genes in human meniscus.
Jump to Full Text
MedLine Citation:
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 :
10405324 - A novel human gene encoding an f-box/wd40 containing protein maps in the shfm3 critical...
23374714 - Replication initiator 1 in adipose tissue function and human obesity.
16568094 - Map1d, a novel methionine aminopeptidase family member is overexpressed in colon cancer.
8812434 - Transcription elongation factor sii (tcea) maps to human chromosome 3p22 --> p21.3.
14706644 - Study of drug effects of calcium channel blockers on retinal degeneration of rd mouse.
23831514 - Tfii-i regulates target genes in the pi-3k and tgf-β signaling pathways through a novel...
Publication Detail:
Type:  Journal Article; Research Support, N.I.H., Extramural; Research Support, Non-U.S. Gov't     Date:  2013-01-15
Journal Detail:
Title:  International journal of obesity (2005)     Volume:  37     ISSN:  1476-5497     ISO Abbreviation:  Int J Obes (Lond)     Publication Date:  2013 Sep 
Date Detail:
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    
Other Details:
Languages:  eng     Pagination:  1238-46     Citation Subset:  IM    
Export Citation:
APA/MLA Format     Download EndNote     Download BibTex
MeSH Terms
ADAM Proteins / metabolism
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*
Middle Aged
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
Grant Support
Reg. No./Substance:
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 N-Endopeptidase; EC Metalloproteinase 1

From MEDLINE®/PubMed®, a database of the U.S. National Library of Medicine

Full Text
Journal Information
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)
ISSN: 0307-0565
ISSN: 1476-5497
Article Information
Download PDF

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
ID: 3751987
DOI: 10.1038/ijo.2012.221
ID: NIHMS425600

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;


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 distribution610. 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, 1721, but it is unclear how degenerative changes in the menisci affect cartilage homeostasis2225.

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 obesity3133. 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 relationship3638.

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.

Ethics statements

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.

Patients and clinical samples

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.

Tissue processing

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.

Total RNA isolation

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).

First strand cDNA synthesis

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.

Selection of candidate genes and oligonucleotide information

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.

Quantification of gene expression levels

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.

Inter-relations of gene expression levels

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.

Statistical Analysis

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

[Formula ID: FD1]
where X refers to the gene expression trait, Y to the factor of interest (either BMI or age), and Z to the factor controlled for in the analysis (either BMI or age).

Probability values for the Spearman’s partial correlations were computed by treating

[Formula ID: FD2]
as coming from a t-distribution with (n−k−2) degrees of freedom, where r is the partial correlation and k is the number of variables being partialled out.

Characteristics of the study population

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.).

Gene expression is not correlated with sex

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.

Gene expression is negatively correlated with age

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.

Proinflammatory cytokines

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).

Matrix degrading enzymes

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.

Transcription factors

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).

Cartilage/meniscus matrix genes

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).

Obesity-related genes

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.

Gene expression is negatively associated with BMI

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.

Age has more impact on gene expression pattern than does BMI

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.

Genes of similar categories cluster together

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 osteoarthritis3638. 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 osteoarthritis5153. 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, 1719. 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.

Conclusion and future directions

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.

Supplementary Material 1 Click here for additional data file (NIHMS425600-supplement-1.docx)

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)





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

Article Categories:
  • Article

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...