Email updates

Keep up to date with the latest news and content from BMC Molecular Biology and BioMed Central.

Open Access Highly Accessed Research article

Selection of reliable reference genes for qPCR studies on chondroprotective action

Stefan Toegel1*, Wenwen Huang1, Claudia Piana1, Frank M Unger1, Michael Wirth1, Mary B Goldring2, Franz Gabor1 and Helmut Viernstein1

Author Affiliations

1 Department of Pharmaceutical Technology and Biopharmaceutics, University of Vienna, Althanstrasse 14, A-1090 Vienna, Austria

2 Beth Israel Deaconess Medical Center, New England Baptist Bone and Joint Institute, and Harvard Medical School, Boston, Massachusetts 02115, USA

For all author emails, please log on.

BMC Molecular Biology 2007, 8:13  doi:10.1186/1471-2199-8-13

The electronic version of this article is the complete one and can be found online at: http://www.biomedcentral.com/1471-2199/8/13


Received:14 November 2006
Accepted:26 February 2007
Published:26 February 2007

© 2007 Toegel et al; licensee BioMed Central Ltd.

This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract

Background

Chondroprotective agents (CPA) such as glucosamine, curcumin and diacerein represent potential remedies for the management of osteoarthritis and several studies have been performed on their effects in-vitro and in-vivo. For the investigation of chondroprotective action on chondrocyte gene expression, quantitative real-time RT-PCR is the method of choice. However, validation of applied normalization strategies represents a crucial and sometimes neglected step in the analysis process. Therefore, the present study aimed to determine the expression stability of common reference genes (ACTB, Beta actin; GAPDH, Glyceraldehyde-3-phosphate; B2M, Beta-2-microglobulin; HPRT1, Hypoxanthine phosphoribosyl-transferase I; SDHA, Succinate dehydrogenase complex, subunit A; YWHAZ, Tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein, zeta polypeptide) under the influence of glucosamine, curcumin and diacerein in the IL-1β-stimulated C-28/I2 chondrocyte model, using the geNorm software tool.

Results

CPA treatment of C-28/I2 chondrocytes significantly affected the expression level of many reference genes (p < 0.05). According to their expression stability, geNorm analysis revealed rankings of the 3 most stable genes (from most stable to least stable) as follows: GAPDH, B2M and SDHA in glucosamine treated samples and HPRT1, GAPDH and B2M in curcumin or diacerein treated samples. Interestingly, ACTB was one of the most variably expressed genes throughout all experiments.

Conclusion

Our study points out the problem of relying on commonly used reference genes without an accurate validation process. For normalization purposes in gene profiling studies on glucosamine action, the genes GAPDH, B2M and SDHA are recommended as single reference genes depending on the expression level of the target gene or more favourably in combination. For experiments with curcumin and diacerein the use of HPRT1, GAPDH and B2M should be considered.

Background

Osteoarthritis (OA) is a chronic, degenerative disorder of unknown cause characterized by gradual loss of articular cartilage. It is the most common of all joint diseases and represents a major social and economic burden since its prevalence increases with age [1]. Currently, classic therapeutic approaches are still limited to symptom relieving drugs or surgical intervention. Chondroprotection with drugs possessing disease-modifying qualities represents an alternative concept in the treatment of OA and its clinical potential has been the subject of numerous studies, both clinical and in-vitro [2-10]. Examples of chondroprotective agents (CPA) are glucosamine, chondroitin sulphate, curcumin, diacerein, rhein, and avocado/soybean unsaponifiables. However, medical opinion about the applicability and clinical efficacy of CPA in OA remains divided. In the case of glucosamine, for example, numerous clinical reports have indicated beneficial effects such as symptom relief in OA [2,3]. On the other hand, several reviews have provided criticism of study designs or processing [11,12]. Moreover, the recently published 'Glucosamine/chondroitin Arthritis Intervention Trial (GAIT)', which was intended to clarify the effectiveness of glucosamine and chondroitin in reducing the pain of osteoarthritis, has generated further controversy by revealing negative findings and high placebo response [13]. It should be further noted that to date limited information regarding synergistic effects, contraindications and drug interactions of CPA is available [14]. However, the rationale to intervene therapeutically in the course of OA with pharmacological agents able to improve biosynthetic processes of chondrocytes accounts for an ongoing interest in the evaluation of CPA. Moreover, these often natural-sourced agents are widely available, generally well tolerated and presumably possess more desirable safety profiles as compared to non-steroidal anti-inflammatory drugs (NSAIDS). A detailed clarification of the molecular effects of CPA is therefore regarded as a reasonable and promising approach leading to progress in understanding of OA modifying action.

At the molecular level, osteoarthritic tissue destruction originates from an imbalance of catabolic cytokine and anabolic growth factor pathways in chondrocytes and synovial cells. Studies have shown that the sensitive homeostasis of cartilage extracellular matrix is effectively disturbed by the TNFα/IL-1β driven cytokine cascade [15-17]. Recent in-vitro studies suggest that CPA have direct influence on the cytokine cascade resulting in down-regulation of catabolic downstream mediators such as collagenase, aggrecanase, NO synthase and prostaglandin E2 [8-10]. For the investigation of chondroprotective action, IL-1β-stimulated chondrocyte models represent a helpful tool [7,8,10,18-20]. In these models, chondrocytes cultured in monolayer or alginate beads are incubated with an inflammatory mediator, such as IL-1β and the CPA of interest, alone and in combination, and selected targets are investigated at the protein or gene expression level.

During the analysis of the molecular biological consequences of chondroprotective action, real-time quantitative RT-PCR (qPCR) is the method of choice for monitoring alterations of gene expression patterns in chondrocytes [21-25]. qPCR is one of the most powerful analytical tools for the sensitive and accurate quantification of mRNA expression levels in cells of different populations or under different culture conditions. One of the main challenges encountered in qPCR data interpretation is the selection of an appropriate normalization strategy. Data normalization is required to control for the experimental error introduced during the multistage process of isolating and processing RNA [26]. Different options range from standardization of tissue weight or cell numbers to the use of artificial molecules incorporated into the samples [26]. However, the most popular method for internal control of qPCR data error is the use of reference genes with presumably stable expression levels under a variety of experimental conditions. Since several studies have demonstrated the potential regulation of those reference genes by the conditions of bioassay [27-31], leading qPCR investigators have declared that the validation of reference genes for each experimental situation is a crucial requirement for the acquisition of biologically meaningful data [26,29,32,33]. In this context, a recent study demonstrated that when a carefully validated reference gene is used to normalize qPCR data, the results can be significantly different from those obtained when an unvalidated reference gene is used [34]. Nevertheless, a survey of the recent literature makes it apparent that reference genes are often used for normalization of qPCR data without any mention of a validation process. This is also true for papers on chondrocyte function. Furthermore, only two studies investigating reference gene stabilities in chondrocytes under mechanical loading or between patient samples, respectively, have been published [35,36].

The present study was designed to obtain a methodological evaluation of the influence of chondroprotective agents on reference gene mRNA levels in the IL-1β stimulated chondrocyte model. We therefore cultured IL-1β-treated C-28/I2 chondrocytes in the presence of selected CPA, i.e. glucosamine, curcumin, or diacerein and assessed the mRNA levels of prominent reference genes by qPCR. Our results will provide information for future studies on CPA about appropriate reference genes for internal standardization of qPCR data in the IL-1β stimulated chondrocyte model.

Results

Optimization and quality assessment of the qPCR protocol

Since the efficiency of reverse transcription reactions depends on the total amount of RNA, the input RNA for the reverse transcription step should be identical if various groups of samples are compared [26]. This requires accurate quantification of the isolated RNA as well as assessment of the RNA integrity. In the present study, RiboGreen assays of RNA isolated from differently treated C-28/I2 chondrocytes revealed RNA amounts ranging from 6.96 μg to 36.3 μg. No correlation between the amount of total RNA and the different treatments was found. The integrity of the RNA preparations was verified by the inspection of the 28S and 18S ribosomal RNA bands using agarose gel electrophoresis. Efficiency of the DNase digest was additionally assessed by implementing minus RT controls for each RNA preparation in the qPCR experiments. Some of the minus RT control samples recorded Ct values of around 32 indicating minor genomic DNA contamination. However, the ΔCt – separating the minus RT controls from the unknown samples – of more than 13 Ct was regarded as high enough to be negligible [37]. In summary, all RNA samples were found suitable for qPCR reactions.

Optimization of primer concentrations is critical for increasing specificity and sensitivity of the qPCR assay [33]. The results of the described primer optimization experiments are presented in Table 1. For all qPCR reactions of the study, these respective primer concentrations were used. Afterwards, standard curves were prepared to examine the quality of the overall qPCR assay. Standard curves of all examined genes revealed high amplification efficiencies over 4 to 5 orders of magnitude (Table 1) covering the Ct range of all samples. The R2 values of the standard curves exceeded 0.996. Gene-specific amplification was confirmed by a single peak in the melting-curve analyses (data not shown). Only ACTB amplification resulted in the detection of a second peak, which is known to originate from the presence of regions with highly different GC content rather than from unspecific reactions [38]. Moreover, no primer-dimer was detectable in any of the melting-curves indicating optimal performance of the primers. The no-template controls (NTC), included in all qPCR runs, did not record any positive Ct value throughout the study.

Table 1. Function of reference genes, primer sequences, optimal primer concentrations and reaction efficiencies of qPCR experiments.

Transcriptional profiling of the reference genes

A real-time qPCR assay, based on SYBR Green detection, was designed for the transcriptional profiling of six frequently used reference genes (ACTB, GAPDH, B2M, HPRT1, SDHA, YWHAZ) in the cDNA samples. Figures 1A, B and 1C show the distribution of the Ct values of all samples under all culture conditions (Table 2) for glucosamine, curcumin and diacerein experiments, respectively. This 'whole error approach' is a simplified way to visualize the overall error introduced during the whole assay and gives an overview of the abundance of the genes in the samples. Since lower Ct values indicate higher numbers of starting cDNA molecules, the most abundant reference genes in C-28/I2 cells were ACTB, GAPDH and B2M followed by YWHAZ. The amplification of HPRT1 and SDHA occurred about 5 Ct values later. However, this 'whole error approach' does neither take into account the influence of qPCR efficiencies, nor does it provide satisfactory resolution for statistical evaluation of small differences between study groups. In a further step, the Ct values of each sample were related to the average Ct value of the respective untreated control group by the equation E-ΔCt. Tables 3, 4 and 5 show the relative reference gene expression levels of IL-1β treated chondrocytes under the influence of glucosamine, curcumin and diacerein, respectively. IL-1β did not provoke significant differences of any reference gene except for ACTB which was significantly up-regulated in the glucosamine (Table 3) and diacerein (Table 5) experiments (p < 0.05). Regarding the impact of the CPA used, a trend for dose-dependent down-regulation of all reference genes was observed for both the pre-treated and the simultaneously treated study groups. In contrast, time-dependent down-regulation seems to occur at the high concentrations of CPA used. We found that both glucosamine (Table 3) and curcumin (Table 4) treatment affected the expression level of all reference genes significantly (p < 0.05). Comparing two samples, a twofold change in gene expression level is commonly regarded as the benchmark of biological meaningfulness [29,35]. Therefore, a less than twofold change is considered as a requirement for suitability of a reference gene. In our study, a 48 h incubation of chondrocytes with high concentrations of glucosamine (Table 3) and curcumin (Table 4) (5 mM and 50 μM, respectively) induced most reference genes to exceed this benchmark of biological meaningfulness. The only exceptions were B2M in the case of glucosamine treatment (Table 3) and YWHAZ in the case of curcumin treatment (Table 4). The data further indicate differences in the reference gene regulation caused by different CPA. For diacerein-treated chondrocytes, GAPDH, B2M, HPRT1 and SDHA were not significantly affected by any culture condition used.

thumbnailFigure 1. qPCR cycle threshold values of 6 reference genes among 7 culture conditions (table 2). Distribution of expression levels of ACTB, GAPDH, B2M, HPRT1, SDHA, and YWHAZ in IL-1β treated C-28/I2 chondrocytes under the influence of glucosamine (A), curcumin (B) and diacerein (C) is shown as medians (lines), 25th percentile to the 75th percentile (boxes) and ranges (whiskers) (n = 63). Samples were normalized to the amount of total RNA.

Table 2. Experimental design

Table 3. Relative mRNA levels of C-28/I2 chondrocytes under the influence of IL-1β and glucosamine.

Table 4. Relative mRNA levels of C-28/I2 chondrocytes under the influence of IL-1β and curcumin.

Table 5. Relative mRNA levels of C-28/I2 chondrocytes under the influence of IL-1β and diacerein.

GeNorm analysis

In order to determine the most stable reference genes in our study, the data were analyzed using the geNorm software. Table 6 shows the ranking of the 6 candidate reference genes according to their expression stability under the influence of glucosamine, curcumin and diacerein. The optimal number of reference genes necessary for the calculation of reliable normalization factors was assessed calculating the pairwise variations Vn/Vn+1 according to Vandesompele et al. [32]. RT-PCR normalization factors were calculated on the basis of the 3 most stable genes and stepwise inclusion of the other genes in the order of their expression stability. The pairwise variation Vn/Vn+1 was calculated between the 2 sequential normalization factors (NFn and NFn+1). A large variation means that the added gene has a significant effect on the newly calculated normalization factor and should preferably be included for calculation. In all 3 experiments, the inclusion of the 4th reference gene did not contribute significantly to the variation of the normalization factor (V3/4 value < 0.15). Based on the proposed cut-off value of 0.15, below which the inclusion of an additional control gene is not required, the use of 3 control genes would therefore be sufficient for normalization in our study. The 3 most stable reference genes for glucosamine-treated samples were GAPDH, B2M and SDHA. In contrast, for curcumin and diacerein-treated samples, HPRT1, GAPDH and B2M provided the most stable expression levels with SDHA ranking the 6th and 5th positions, respectively. It is remarkable that ACTB, one of the most prominent and frequently used reference genes, proved to be one of the most variably expressed genes under the influence of selected CPA.

Table 6. Ranking of the reference genes as determined using the geNorm software tool.

Discussion

The claim of real-time qPCR to assess changes of mRNA expression levels in a quantitative way requires the choice of an adequate normalization strategy [26]. Normalization must be performed to compensate sample-to-sample and run-to-run variations. Reference genes theoretically represent a convenient way to normalize qPCR data since their use controls for all steps of gene expression analysis. However, the quality of normalized quantitative expression data cannot be better than the quality of the normalizer itself. Any variation in the normalizer will obscure real changes and produce additional errors. Therefore, the stability of the reference gene of choice represents the benchmark for the resolution of the qPCR assay and directly influences the significance of the results. By geometric averaging of multiple control genes, the geNorm software tool represents an alternative for accurate normalization of qPCR data, in particular, when the significance of subtle gene expression differences is evaluated [32]. The determination of the most stable reference genes in a given experimental setup and the subsequent calculation of normalization factors thus open up the possibility of reliably studying small expression differences even in low-copy number genes.

In the present study, we determined the expression stability of common reference genes under the influence of glucosamine, curcumin and diacerein in C-28/I2 chondrocytes. Immortalized C-28/I2 cells represent an alternative to primary human chondrocytes circumventing certain shortcomings of primary cells such as high inter-individual variability between different donors and dedifferentiation during in-vitro culture [39,40]. Moreover, C-28/I2 cells are known to be reactive to IL-1β [41] suggesting that they represent a suitable model for investigating the effects of CPA on the IL-1β-driven catabolic cascade. The present study design was developed to mirror experimental protocols that were used in recent studies [18-25] investigating the activity of CPA in chondrocytes in the presence of IL-1β. In order to create high variability, we used 3 different CPA for the treatment of chondrocyte cultures, i.e. glucosamine, curcumin and diacerein. The range of concentrations used in the presented experiments were in accordance with recently published in-vitro studies [20,23,42]. To mimic different approaches used previously, the CPA were added to the culture medium either simultaneously with IL-1β or prior to the addition of IL-1β.

The presented data report insignificant impact of IL-1β on the variability of most reference gene expression levels (p > 0.05). Only ACTB levels were significantly changed in the glucosamine and diacerein experiments. In this regard, the C-28/I2 model is in accordance with results of Chan et al. [25] who showed that GAPDH was not regulated by IL-1β in primary bovine chondrocytes. Moreover, our results demonstrated significantly variable expression levels of all reference genes in C-28/I2 chondrocytes under most culture conditions containing CPA. It is suggested that this variability was due to both dose and time dependent effects induced by the CPA. In an attempt to provide information on the most uniformly expressed reference genes for studies analyzing chondroprotective activity, we determined the ranking of 6 candidate genes according to their expression stability depending on the CPA under investigation. This ranking was not equivalent between the 3 different CPA. The high variability of ACTB in our study further points out the problem of relying on commonly used reference genes without an accurate validation process. Therefore, using the geNorm software, we defined the sets of reference genes suitable for studies on the presented CPA. GAPDH, B2M and SDHA are proposed for studies on glucosamine, while HPRT1, GAPDH and B2M can be used to investigate curcumin or diacerein activity. If glucosamine, curcumin and diacerein are to be compared, the same set as proposed for curcumin or diacerein can be used. It is generally recommended that the gene of interest and the selected reference gene should roughly be expressed at the same level. Therefore, the choice of a single reference gene (e.g. GAPDH or SDHA in glucosamine studies) may be determined by the abundance of the gene of interest. When a study aims at the fine measurement of small gene expression differences even the geometric mean of 2 or 3 reference genes could be of advantage.

Conclusion

In conclusion, our report underlines the importance of accurate reference gene validation for each qPCR study. However, this validation process represents an evident disadvantage of the reference gene normalization strategy, since remarkable effort and costs are associated with its realization. In order to provide a reference for future studies on chondrocyte responses to different agents, we have defined the sets of suitable reference genes for the molecular biological investigation of glucosamine, curcumin and diacerein activities. Our findings will facilitate the use of the reference gene method for the normalization and data interpretation of CPA effects in the IL-1β stimulated chondrocyte model, in particular when subtle gene expression differences are to be investigated.

Methods

Materials

Glucosamine hydrochloride, curcumin, diacerein and 0.25 % Trypsin-EDTA were purchased from Sigma (St. Louis, MO, USA). Recombinant human IL-1β was obtained from Strathmann Biotec (Hamburg, D). Primers for real-time PCR experiments and the mi-Gel Extraction Kit were from Metabion (Martinsried, D). Dulbecco's modified Eagle's medium and 100x ITS solution were purchased from Gibco (Lofer, A). The StrataScript First-Strand Synthesis System and the Brilliant SYBR Green QPCR Master Mix were from Stratagene (La Jolla, CA, USA). All other chemicals were of analytical grade and obtained from Merck (Darmstadt, D) unless otherwise specified.

Cell culture and experimental design

Immortalized C-28/I2 chondrocytes were used for all experiments in this study. In order to determine the influence of glucosamine, curcumin and diacerein on reference gene expression levels in chondrocytes, three sets of individual experiments – one for each compound – were performed.

2.3 × 105 chondrocytes were seeded into 25 cm2 culture flasks (Corning, NY, USA) and cultured for 5 days in DMEM supplemented with 10 % fetal calf serum (Biochrom AG, Berlin, D), 2 μl/ml gentamycin and 50 μg/ml ascorbate. On days 5 and 6 the culture medium was changed to serum-free DMEM supplemented with 1 % (v/v) 100× ITS (10 μg/ml insulin + 5.5 μg/ml transferrin + 6.7 ng/ml sodium selenite) and 2 μl/ml gentamycin. Additionally, the particular CPA and/or IL-1β (10 ng/ml) were added as indicated in Table 2. Thereby, 7 different conditions for each CPA experiment were created: one control group without additional supplements, one IL-1β treated group, two groups simultaneously treated with both IL-1β and CPA at two different concentrations, two groups pre-incubated with CPA at two different concentrations prior to IL-1β addition, and one group treated with the respective CPA at the higher concentration without the addition of IL-1β. Each of these study groups consisted of three identically treated 'biological' replicates. Throughout cell culturing, cells were maintained in a humidified 5 % CO2/95 % air atmosphere at 37°C.

RNA isolation and cDNA synthesis

On day 7 of the culture period, chondrocytes were harvested by trypsination using 0.05 % Trypsin-EDTA. The cells were washed with PBS and total RNA was extracted immediately using the NucleoSpin RNA II Kit (Macherey-Nagel, Dueren, D) according to the manufacturer's instructions. On-column DNase digestion was performed for all samples during the isolation process. To control for the integrity of the isolated RNA preparations, all samples were qualitatively assessed according to standard methods prior to use. Briefly, a 1.5% denaturing agarose gel was loaded with 2 μg of total RNA samples and stained with ethidium bromide after electrophoresis. The quantity of extracted RNA was determined using the Quant-It RiboGreen reagent (Molecular Probes, Eugene, USA) and the quantitative plate-read function of the Stratagene MxPro QPCR software. An R2 > 0.990 of the standard curve plots was considered acceptable for accurate RNA quantification. Reverse transcription (RT) was performed using the StrataScript First-Strand Synthesis System, following the manufacturer's instructions. Briefly, 2 μg of total RNA were reverse transcribed in a 20 μl reaction volume using the manufacturer-supplied oligo(dT) primers. Minus RT controls were prepared for each sample using the identical procedure except for the omission of reverse transcriptase. All cDNA preparations and minus RT controls were diluted at a ratio of 1:10 with RNase free water prior to real-time PCR.

Reference gene selection and primer optimization

Six reference genes (ACTB, GAPDH, B2M, HPRT1, SDHA, YWHAZ) belonging to different functional classes were selected to reduce the chance that these genes might be co-regulated (Table 1). Primers for GAPDH were designed using AlleleID 2.01 software (Premier Biosoft, CA, USA) and the Genbank sequence NM_002046. In order to avoid non-specific product formation, the primers were designed across exon junctions and with minimized self and cross dimer ΔG values. Additionally, BLAST analyses were performed to verify specificity. Primer sequences for all other reference genes were used as described by Vandesompele et al. [32]. Amplicon sizes were 140 bp for ACTB, 165 bp for GAPDH, 86 bp for B2M, 94 bp for HPRT1, 86 bp for SDHA and 94 bp for YWHAZ. Gene specific primers were used to amplify respective PCR products in order to provide sufficient template material for primer concentration optimization. PCR amplifications were performed using GoTaq DNA Polymerase (Promega, Wisconsin, USA) in a 50 μl reaction volume (10 μl of 5x Colourless Go Taq Reaction Buffer, 1 μl of PCR Nucleotide Mix, 5 μl of gene specific forward and reverse primer, 0.25 μl of GoTaq DNA polymerase, 2 μl of template DNA, and nuclease-free water to 50 μl). The PCR program consisted of an initial denaturation step at 95°C for 2 min, followed by 35 cycles at 95°C for 45 s, at 55°C for 45 s and at 72°C for 1 min, and a final extension step at 72°C for 5 min. Amplified PCR products were diluted with bromophenol blue solution at a ratio of 1:4 and purified by gel electrophoresis using a 1% agarose gel. After staining with ethidium bromide, the designated bands were cut out of the agarose gel and the products were extracted from the gel using the mi-Gel Extraction Kit according to the manufacturer's instructions. Optimization of primer concentrations for the qPCR experiments was performed by qPCR assays as described below using different primer concentrations and a 1:103 dilution of purified PCR products as templates. Optimal primer pairs for each reference gene were identified selecting the combination that gave the lowest Ct value and the highest ΔRn value.

Real-time qPCR

All qPCR reactions were performed in 25 μl reaction mixtures containing 1 μl cDNA, 12.5 μl Brilliant SYBR Green QPCR Master Mix, primer pairs as listed in Table 1, and nuclease-free water to 25 μl. Each biological replicate was run in triplicate on a Mx3000P QPCR system (Stratagene). Thermocycling conditions consisted of an initial polymerase activation step at 95°C for 10 min, followed by 35 cycles at 95°C for 30 s, at 55°C for 1 min, and at 72°C for 1 min. Afterwards, melting curves were generated to confirm a single gene-specific peak and to detect primer-dimer formation by heating the samples stepwise from 55°C to 95°C while continuously monitoring the fluorescence. NTC were included in each run to control for contaminations. Minus RT controls were run to verify the efficiency of the on-column DNase digestion. The results were analyzed using the MxPro real-time QPCR software (Stratagene). Baselines and thresholds were automatically set by the software and used after manual inspection. The crossing point of the amplification curve with the threshold represented the cycle threshold (Ct). PCR efficiencies for each primer pair were derived from standard curves using five-fold dilution series covering a 4 to 5 log dynamic range starting from one randomly selected undiluted cDNA sample of the untreated control group. The reactions were run in duplicate and no-template controls were included.

Data presentation and calculations

Results were exported to Microsoft Excel and GraphPad Prism for further analyses. For each test compound, the distribution of the expression levels (Ct values) for each reference gene under the 7 experimental conditions was displayed as Box and Whiskers plots. In a next step, data within each experiment were normalized to the average level of mRNA expression (Ct value) of the respective control group using the relationship E-ΔCt (E, efficiency; ΔCt, difference between Ct values). Statistics were performed using one-way ANOVA with post-hoc Dunnett tests comparing each study group with the untreated control group (significance p < 0.05). The geNorm application software for Microsoft Excel was additionally used as described by Vandesompele et al. [32] to identify the most stable reference gene under the described conditions, and to determine the optimal number of reference genes required for reliable normalization of qPCR data.

Authors' contributions

ST conceived of the study, participated in the qPCR experiments, carried out the statistical analyses and drafted the manuscript. WH and CP carried out the cell culture and gene expression analyses. MBG contributed to writing the manuscript. MW, FMU, FG and HV supervised the study design and the experimental processes. All authors read and approved the manuscript.

Acknowledgements

Stefan Toegel would like to thank Stratagene for the Stratagene Research Award 2005. Part of the work and Claudia Piana were supported by the Integrated Project CellProm (NMP4-CT-2004-500039) granted from the sixth framework program of the European Community. Dr. Goldring's research is supported by NIH grant AG-22021.

References

  1. Haq I, Murphy E, Dacre J: Osteoarthritis.

    Postgrad Med J 2003, 79:377-383. PubMed Abstract | Publisher Full Text OpenURL

  2. DeLosReyes BG, Koda RT, Lien EJ: Glucosamine and chondroitin sulfates in the treatment of osteoarthitis: a survey.

    Prog Drug Res 2000, 55:81-103. PubMed Abstract OpenURL

  3. Setnikar I: Antireactive properties of "chondroprotective" drugs.

    Int J Tissue React 1992, 14:253-261. PubMed Abstract OpenURL

  4. Dougados M, Nguyen M, Berdah L, Mazieres B, Vignon E, Lequesne M: Evaluation of the structure modifying effects of diacerein in hip osteoarthritis.

    Arthritis Rheum 2001, 44:2539-2547. PubMed Abstract | Publisher Full Text OpenURL

  5. Rein E, Kharazmi A, Winther K: A herbal remedy, Hyben Vital (stand. powder of a subspecies of Rosa canina fruits), reduces pain and improves general wellbeing in patients with osteoarthritis – a double-blind, placebo-controlled, randomised trial.

    Phytomedicine 2004, 11:383-391. PubMed Abstract | Publisher Full Text OpenURL

  6. Joe B, Lokesh BR: Prophylactic and therapeutic effects of n-3 polyunsaturated fatty acids, capsaicin, and curcumin on adjuvant induced arthritis in rats.

    J Nutr Biochem 1997, 8:397-407. Publisher Full Text OpenURL

  7. Martin G, Bogdanowicz P, Domagala F, Ficheux H, Pujol JP: Rhein inhibits interleukin-1β-induced activation of MEK/ERK pathway and DNA binding of NF-kB and AP-1 in chondrocytes cultured in hypoxia: a potential mechanism for its disease-modifying effect in osteoarthritis.

    Inflammation 2003, 27:233-246. PubMed Abstract | Publisher Full Text OpenURL

  8. Dodge GR, Jimenez SA: Glucosamine sulfate modulates the levels of aggrecan and matrix metalloproteinase-3 synthesized by cultured human osteoarthritis articular chondrocytes.

    Osteoarthr Cartilage 2003, 11:424-432. Publisher Full Text OpenURL

  9. Bassleer C, Henrotin Y, Franchimont P: In-vitro evaluation of drugs proposed as chondroprotective agents.

    Int J Tissue React 1992, 14:231-241. PubMed Abstract OpenURL

  10. Liacini A, Sylvester J, Zafarullah M: Triptolide suppresses proinflammatory cytokine-induced matrix metalloproteinase and aggrecanase-1 gene expression in chondrocytes.

    Biochem Bioph Res Co 2005, 327:320-327. Publisher Full Text OpenURL

  11. Felson DT, McAlindon TE: Glucosamine and chondroitin for osteoarthritis: to recommend or not to recommend?

    Arthritis Care Res 2000, 13:179-182. PubMed Abstract | Publisher Full Text OpenURL

  12. Richette P, Bardin T: Structure-modifying agents for osteoarthritis: an update.

    Joint Bone Spine 2004, 71:18-23. PubMed Abstract | Publisher Full Text OpenURL

  13. Clegg DO, Reda DJ, Harris CL, Klein MA, O'Dell JR, Hooper MM, Bradley JD, Bingham CO 3rd, Weisman MH, Jackson CG, Lane NE, Cush JJ, Moreland LW, Schumacher HR Jr, Oddis CV, Wolfe F, Molitor JA, Yocum DE, Schnitzer TJ, Furst DE, Sawitzke AD, Shi H, Brandt KD, Moskowitz RW, Williams HJ: Glucosamine, chondroitin sulfate, and the two in combination for painful knee osteoarthritis.

    N Engl J Med 2006, 354:795-808. PubMed Abstract | Publisher Full Text OpenURL

  14. Walker-Bone K: 'Natural remedies' in the treatment of osteoarthritis.

    Drug Aging 2003, 20:517-526. Publisher Full Text OpenURL

  15. Wang J, Verdonk P, Elewaut D, Veys EM, Verbruggen G: Homeostasis of the extracellular matrix of normal and osteoarthritic human articular cartilage chondrocytes in vitro.

    Osteoarthr Cartilage 2003, 11:801-809. Publisher Full Text OpenURL

  16. Garnero P, Rousseau JC, Delmas PD: Molecular basis and clinical use of biochemical markers of bone, cartilage and synovium in joint diseases.

    Arthritis Rheum 2000, 43:953-968. PubMed Abstract | Publisher Full Text OpenURL

  17. Sandell LJ, Aigner T: Articular cartilage and changes in arthritis. An introduction: cell biology of osteoarthritis.

    Arthritis Res 2001, 3:107-113. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  18. Richardson DW, Dodge GR: Dose-dependent effects of corticosteroids on the expression of matrix-related genes in normal and cytokine-treated articular chondrocytes.

    Inflamm Res 2003, 52:39-49. PubMed Abstract | Publisher Full Text OpenURL

  19. Curtis CL, Rees SG, Cramp J, Flannery CR, Hughes CE, Little CB, Williams R, Wilson C, Dent CM, Harwood JL, Caterson B: Effects of n-3 fatty acids on cartilage metabolism.

    P Nutr Soc 2002, 61:381-389. Publisher Full Text OpenURL

  20. Schulze-Tanzil G, Mobasheri A, Sendzik J, John T, Shakibaei M: Effects of curcumin (Diferuloylmethane) on nuclear factor kB signaling in Interleukin-1β-stimulated chondrocytes.

    Ann NY Acad Sci 2004, 1030:578-586. PubMed Abstract | Publisher Full Text OpenURL

  21. Pelletier JP, Boileau C, Boily M, Brunet J, Mineau F, Geng C, Reboul P, Laufer S, Lajeunesse D, Martel-Pelletier J: The protective effect of licofelone on experimental osteoarthritis is correlated with the downregulation of gene expression and protein synthesis of several major cartilage catabolic factors: MMP-13, cathepsin K and aggrecanases.

    Arthritis Res Ther 2005, 7:R1091-1102. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  22. Uitterlinden EJ, Jahr H, Koevoet JL, Jenniskens YM, Bierma-Zeinstra SM, Degroot J: Glucosamine decreases expression of anabolic and catabolic genes in human osteoarthritic cartilage explants.

    Osteoarthritis Cartilage 2006, 14:250-257. PubMed Abstract | Publisher Full Text OpenURL

  23. Sanchez C, Mathy-Hartert M, Deberg MA, Ficheux H, Reginster JYL, Henrotin YE: Effects of rhein on human articular chondrocytes in alginate beads.

    Biochem Pharmacol 2003, 65:377-388. PubMed Abstract | Publisher Full Text OpenURL

  24. Miller MJS, Ahmed S, Bobrowski P, Haqqi TM: The chondroprotective actions of a natural product are associated with the activation of IGF-1 production by human chondrocytes despite the presence of IL-1β.

    BMC Complementary and Alternative Medicine 2006, 6:13. BioMed Central Full Text OpenURL

  25. Chan PS, Caron JP, Orth MW: Short-term gene expression changes in cartilage explants stimulated with IL-1β plus glucosamine and chondroitin sulfate.

    J Rheumatol 2006, 33:1329-1340. PubMed Abstract | Publisher Full Text OpenURL

  26. Huggett J, Dheda K, Bustin S, Zumla A: Real-time RT-PCR normalisation; strategies and considerations.

    Genes Immun 2005, 6:279-284. PubMed Abstract | Publisher Full Text OpenURL

  27. Schmittgen TD, Zakrajsek BA: Effect of experimental treatment on housekeeping gene expression: validation by real-time, quantitative RT-PCR.

    J Biochem Bioph Meth 2000, 46:69-81. Publisher Full Text OpenURL

  28. Zhong H, Simons JW: Direct comparison of GAPDH, β-actin, cyclophilin, and 28S rRNA as internal standards for quantifying RNA levels under hypoxia.

    Biochem Bioph Res Co 1999, 259:523-526. Publisher Full Text OpenURL

  29. Dheda K, Huggett JF, Bustin SA, Johnson MA, Rook G, Zumla A: Validation of housekeeping genes for normalizing RNA expression in real-time PCR.

    Biotechniques 2004, 37:112-119. PubMed Abstract OpenURL

  30. Gilsbach R, Kouta M, Boenisch H, Bruess M: Comparison of in vitro and in vivo reference genes for internal standardization of real-time PCR data.

    Biotechniques 2006, 40(2):173-7. PubMed Abstract OpenURL

  31. Bogaert L, VanPoucke M, DeBaere , Peelman L, Gasthuys F, Martens A: Selection of a set of reliable reference genes for quantificative real-time PCR in normal equine skin and equine sarcoids.

    BMC Biotechnology 2006, 6:24. PubMed Abstract | BioMed Central Full Text | PubMed Central Full Text OpenURL

  32. Vandesompele J, DePreter K, Pattyn F, Poppe B, VanRoy N, DePaepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes.

    Genome Biology 2002, 3:34. BioMed Central Full Text OpenURL

  33. Bustin SA: A-Z of Quantitative PCR. San Diego: IUL Press; 2004. OpenURL

  34. Dheda K, Huggett JF, Chang JS, Kim LU, Bustin SA, Johnson MA, Rook GAW, Zumla A: The implications of using an inappropriate reference gene for real-time reverse transcription PCR data normalization.

    Anal Biochem 2005, 344:141-143. PubMed Abstract | Publisher Full Text OpenURL

  35. Lee CR, Grad S, MacLean JJ, Iatridis JC, Alini M: Effect of mechanical loading on mRNA levels of common endogenous controls in articular chondrocytes and intervertebral disk.

    Anal Biochem 2005, 341:372-375. PubMed Abstract | Publisher Full Text OpenURL

  36. Asp J, Brantsing C, Lovstedt K, Benassi MS, Inerot S, Gamberi G, Picci P, Lindahl A: Evaluation of p16 and Id1 status and endogenous reference genes in human chondrosarcoma by real-time PCR.

    Int J Oncol 2005, 27:1577-1582. PubMed Abstract OpenURL

  37. Bustin SA, Nolan T: Pitfalls of quantitative real-time reverse-transcription polymerase chain reaction.

    J Biomol Tech 2004, 15(3):155-6. PubMed Abstract | Publisher Full Text OpenURL

  38. Pattyn F, Speleman F, DePaepe A, Vandesompele J: RTPrimerDB: the Real-Time PCR primer and probe database. [http://medgen.ugent.be/rtprimerdb/] webcite

    Nucleic Acids Res 2003, 31:122-123. PubMed Abstract | Publisher Full Text | PubMed Central Full Text OpenURL

  39. Pufe T, Harde V, Petersen W, Goldring MB, Tillmann B, Mentlein R: Vascular endothelial growth factor (VEGF) induces matrix metalloproteinase expression in immortalized chondrocytes.

    J Pathol 2004, 202:367-374. PubMed Abstract | Publisher Full Text OpenURL

  40. Finger F, Schoerle C, Zien A, Gebhard , Goldring MB, Aigner T: Molecular phenotyping of human chondrocyte cell lines T/C-28a2, T/C-28a4 and C-28/I2.

    Arthritis Rheum 2003, 48:3395-3403. PubMed Abstract | Publisher Full Text OpenURL

  41. Goldring MB, Birkhead JR, Suen LF, Yamin R, Mizuno S, Glowacki J, Arbiser JL, Apperley JF: Interleukin-1β-modulated gene expression in immortalized human chondrocytes.

    J Clin Invest 1994, 94:2307-2316. PubMed Abstract | PubMed Central Full Text OpenURL

  42. Byron CR, Orth MW, Venta PJ, Lloyd JW, Caron JP: Influence of glucosamine on matrix metalloproteinase expression and activity in lipopolysaccharide-stimulated equine chondrocytes.

    Am J Vet Res 2003, 64:666-671. PubMed Abstract | Publisher Full Text OpenURL