Skip to main content

Cartilage-selective genes identified in genome-scale analysis of non-cartilage and cartilage gene expression



Cartilage plays a fundamental role in the development of the human skeleton. Early in embryogenesis, mesenchymal cells condense and differentiate into chondrocytes to shape the early skeleton. Subsequently, the cartilage anlagen differentiate to form the growth plates, which are responsible for linear bone growth, and the articular chondrocytes, which facilitate joint function. However, despite the multiplicity of roles of cartilage during human fetal life, surprisingly little is known about its transcriptome. To address this, a whole genome microarray expression profile was generated using RNA isolated from 18–22 week human distal femur fetal cartilage and compared with a database of control normal human tissues aggregated at UCLA, termed Celsius.


161 cartilage-selective genes were identified, defined as genes significantly expressed in cartilage with low expression and little variation across a panel of 34 non-cartilage tissues. Among these 161 genes were cartilage-specific genes such as cartilage collagen genes and 25 genes which have been associated with skeletal phenotypes in humans and/or mice. Many of the other cartilage-selective genes do not have established roles in cartilage or are novel, unannotated genes. Quantitative RT-PCR confirmed the unique pattern of gene expression observed by microarray analysis.


Defining the gene expression pattern for cartilage has identified new genes that may contribute to human skeletogenesis as well as provided further candidate genes for skeletal dysplasias. The data suggest that fetal cartilage is a complex and transcriptionally active tissue and demonstrate that the set of genes selectively expressed in the tissue has been greatly underestimated.


Skeletogenesis begins with condensation of mesenchymal chondroprogenitor cells to form the cartilage anlagen that pattern the early skeleton. Subsequently, for bones that grow by endochondral ossification, the chondrocytes differentiate further to establish the growth plates. At the joint surfaces, development of articular cartilage facilitates and maintains joint movement during fetal life. These multi-step processes require the coordinated expression of many genes, including genes encoding extracellular matrix proteins and morphogens, as well as proliferative, angiogenic, and apoptotic signals [1]. Most of our knowledge of the function of the genes involved has been derived from developmental studies in model systems and cell lines [2] as well as from the identification of disease genes in skeletal disorders.

Whole genome analysis of chondrocyte gene expression has the potential to reveal novel genes and gene expression programs which define the tissue. Although the complete set of genes expressed in human cartilage has not yet been described, analysis of human cartilage cDNA libraries has provided an initial in vivo picture of the cartilage transcriptome [36]. These investigations have also identified expression of both known and novel genes. Comparative microarray studies in rat cartilage [7] and several chondrocyte cell lines [8, 9] have provided a larger set of genes of potential importance in chondrocytes, including genes specific to the stages of chondrocyte differentiation. Wang et al. (2004) identified 92 genes with two-fold variation in expression between hypertrophic and proliferative growth plate chondrocytes. In this in vivo study, significant gene expression changes were principally associated with cell cycle, transcription, extracellular matrix structure, receptor and transporter functions. In microarray studies of mouse micromass cultures [8], 212 genes exhibited at least a ten-fold difference in gene expression as the cultures differentiated. Thus global characterization of gene expression is beginning to describe the identities of key regulatory molecules and their targets in chondrocytes.

Disrupting genes involved in the organization and maturation of the growth plate and/or the stability of articular cartilage results in inherited skeletal disorders that range from perinatal lethal phenotypes to mild disorders with early-onset osteoarthropathy as their major feature [10, 11]. Of the approximately 370 clinically distinguishable skeletal dysplasias [12], mutations in 115 genes have been associated with about 150 disorders. Many of these disease genes are expressed in a cartilage-selective pattern, and therefore identifying additional genes uniquely expressed in cartilage should yield new skeletal dysplasia candidate genes.

To identify a larger set of genes uniquely expressed in chondrocytes, this study describes a genome-scale gene expression profile for 18–22 week human fetal cartilage. There were 161 genes which appeared to be selectively expressed in fetal cartilage, comprising a variety of novel genes that may contribute to skeletal development. The data suggest a complex pattern of cartilage gene expression and indicate that the number of genes selectively expressed in cartilage has been greatly underestimated.


Identification of cartilage-selective genes

To define a set of genes preferentially or uniquely expressed in normal human fetal cartilage, cartilage probeset intensities were compared with probeset intensities across a variety of normal tissues. A two-step process was employed for gene identification, consisting of a training step and a validation step (see the additional data file 1, for a flow chart of an overview of the analysis). The tissue-selectivity of a representative sampling of the identified genes was confirmed by quantitative RT-PCR.

The training dataset consisted of five independent cartilage samples and 41 non-cartilage samples, all analyzed using Affymetrix U133 2.0 Plus arrays. The average correlation coefficient among the cartilage samples (R2) was 0.96. To identify unbiased relationships within the data, and to test the robustness of the normalization and tissue-specificity, an unsupervised approach [13], in which the genes and tissues were grouped based only on expression patterns, was employed. Probesets with the greatest variation across all tissues and whose expression in any two arrays differed by at least two standard deviations from their mean expression across the entire set of samples were selected. This selection yielded 9483 probesets.

Two-way hierarchical clustering based on similarity of expression of these 9483 probesets within the samples was performed (Figure 2). Samples from the same tissues clustered together, indicating that the normalization was sufficiently robust to allow tissue-selective expression patterns to be identified. Even with these relatively non-stringent selection criteria, the results showed a surprisingly large number of genes with a fetal cartilage-selective expression pattern. At least 89 probesets representing 64 genes with coordinately higher expression in cartilage relative to non-cartilage tissues appeared to drive the clustering of the two groups (Figure 2B). These probesets formed a gene expression node in the dendrogram which shared an overall expression correlation of 0.99. The genes represented by these probesets included some well established cartilage-selective genes, including aggrecan (AGC1), type × collagen (COL10A1), and matrilin 3 (MATN3), among others. Thus, a comparative approach with microarrays can identify genes whose expression is cartilage-selective.

Figure 1

Human fetal cartilage section dissected for RNA expression profiling. (A) Normal distal femur cartilage from an 18–22 week fetus. Brackets define the cartilage that was dissected and used for RNA profiling. (B) Toluidine blue stained longitudinal section of the distal head of the femur magnified 8.5 fold. The dissected portion included chondrocytes from the articular, reserve, proliferating, and hypertrophic zones.

Figure 2

Unsupervised analysis of all tissues using U133 2.0 microarrays. (A) Two-way hierarchical clustering of 46 normal tissues (including five fetal cartilage and 41 non-cartilage tissues) and 9483 probesets which vary more than two standard deviations from the mean expression of the probeset. As detailed in the Methods, to generate a cartilage-centric figure, the median cartilage signals were normalized to zero prior to the clustering. This resulted in cartilage gene expression being displayed to the left and expression among the non-cartilage tissues to the right and produced deep branching among the cartilage samples. Mean-centered clustering reflected the high correlation among the cartilage samples (average regression correlation of 0.92) and did not qualitatively alter the correlations among the non-cartilage samples (data not shown). (B) Zoomed image of gene expression from 89 probesets representing 64 genes containing many known cartilage-specific genes which distinguish fetal cartilage expression.

To define a ranked list of genes significantly expressed in cartilage, a supervised analysis [13], comparing cartilage versus non-cartilage gene expression, was employed. This consisted of a two-class analysis with a modified t-test (SAM) (See additional data file 2, for the complete results of this analysis). There were 2634 probesets representing 1720 genes with at least three-fold differential expression when comparing cartilage and non-cartilage tissues, with a false discovery rate of zero. Of these, 2446 of the 2634 probesets demonstrated higher expression in cartilage with respect to non-cartilage tissues, while the remaining 188 probesets were expressed at significantly higher levels in the other tissues. As observed for the hierarchical clustering, probesets representing well-known cartilage markers, including COL2A1, AGC1, COMP, COL9A3, and MMP3 were among the top genes listed. In addition, lubricin (PRG4), an articular cartilage-specific marker, was also identified, confirming the ability to identify genes specific to fetal articular cartilage. Indeed, among the top 35 probesets more highly expressed in cartilage, only four probesets, representing unannotated genes, were derived from genes not previously known to be expressed in cartilage.

In silico validation

Three array platforms were used to validate the 2446 probesets identified in the supervised analysis and generate a robust list of cartilage-selective genes (Table 1). A majority of these probesets (2245 probesets (> 92%)) were identified in 124 U133A and 74 U133B arrays using the Celsius database (see Materials and Methods), and represented expression from 34 normal tissues. A small proportion of the probesets (201/2446) are not found on the Affymetrix™ Human Genome U133A/B Arrays, so these probesets were identified in the analysis of 26 U133 Plus 2.0 arrays, representing eight non-cartilage tissues. A summary of the validation and the tissue distribution are available as additional files.

Table 1 Cartilage-selective genes validated in silico on U133A (left), U133B (center), and U133 Plus 2.0 (right) platforms.

Of the three platforms, the U133A dataset was the most robust with regard to the number of arrays, biological replicates, diversity of tissues, probes identified, and gene annotation. From this platform, 1363 of the 2446 probesets identified in the SAM analysis as expressed at a higher level in cartilage were obtained. Two hundred seventy-four of the 1363 probesets (274/1363), representing 237 genes, exhibited at least five-fold higher expression when compared to non-cartilage tissues and were ranked by cartilage-specificity using an analog of coefficient of variation (CV) (see Methods). Of these, 56 probesets, representing 49 genes, were identified with a CV < 50% in non-cartilage samples, constituting the cartilage-selective gene set from this platform (Table 1, left). Twenty of these genes have mutations that have been associated with skeletal phenotypes in humans and/or mice, representing 44% of the probes selected from this platform.

Eight hundred eighty-two of the 2446 probesets were identified from the U133B validation set. Two hundred fourteen of these probesets, representing 158 genes, were well measured in cartilage with at least five-fold higher expression in cartilage relative to non-cartilage tissues. Of these, 77 probesets had a CV less than 50% in non-cartilage samples, representing 71 cartilage-selective genes (Table 1, center), including 3 genes also identified using the U133A platform (COL11A1, EDIL3, and PDPN).

A subset of the cartilage-selective genes was represented only on the Human Genome U133 Plus 2.0 arrays and were selected from the analysis of 28 non-cartilage samples. In total 201/2446 probesets were not represented in the U133A/B array set. Of these 201 probesets, 96 probesets, representing 85 genes, had a five-fold higher expression in cartilage than non-cartilage samples. By including the CV selection criterion, 52 probesets, representing 50 cartilage-selective genes were identified and added to the complete tally (Table 1, right), including 6 genes also identified using the U133A and U133B arrays (IRAK2, NRP2, WTAP, PITPNC1, AKR1C2, and PTK2).

In summary, 480 genes demonstrating enriched or specific expression in cartilage were selected from the comparison of cartilage and non-cartilage tissues with data derived from the U133A (n = 237), U133B (n = 158) and U133 Plus 2.0 (n = 85) platforms. Of these, a non-redundant set of 161 genes (Table 2), including 11 uncharacterized genes and 16 genes represented by unannotated probesets, were classified as cartilage-selective. These data greatly expand the number of genes known to be selectively expressed in cartilage and emphasize the unique pattern of gene expression that determines its properties.

Table 2 Non-redundant set of 161 cartilage-selective genes organized by chromosomal location.

qRT-PCR validation

Quantitative RT-PCR was used to independently assess the tissue-selectivity of the genes identified in the microarray analysis. For each of the three microarray platforms, the probesets with a CV less than 50% were divided into 10% intervals (0–10% CV, 10–20% etc.) (Table 1), and one gene from the middle of each interval was selected for analysis by qRT-PCR.

All of the thirteen of genes analyzed demonstrated higher expression in cartilage than in the seven non-cartilage tissues studied (Table 3). Also, with one exception (OSMR), the selection threshold of at least five-fold higher expression in cartilage tissues as compared with the average expression among all non-cartilage tissues imposed for the microarray analysis, was observed. For most of the genes studied by qRT-PCR, there was little expression in the seven non-cartilage tissues (median Ct = 33.2), indicating that including the coefficient of variation in the ranking algorithm preferentially identifies genes selectively expressed in cartilage. Also, there was an inverse correlation between the gene rank and the standard deviation in expression level among non-cartilage tissues, indicating that genes with a higher rank were more selectively expressed in cartilage. Finally, there was a trend of decreasing cartilage selectivity moving from the U133A to U133B to U133 2.0 qRT-PCR validations, likely reflecting the decreasing robustness of comparison datasets in the respective platforms. Overall the qRT-PCR experiments replicated and validated findings derived from the comparative microarray data.

Table 3 Summary of qRT-PCR amplification of cartilage-selective genes in fetal cartilage and seven non-cartilage tissues.


Using genome-scale microarrays, gene expression in human fetal cartilage was compared with a robust set of other normal tissues. Hierarchical clustering showed remarkable similarity among the 18–22 week fetal cartilage expression profiles and demonstrated that a subset of the cartilage transcriptome is composed of a unique gene set not generally expressed in the other tissues studied. Using SAM, 2446 probesets measured preferential expression of 1712 genes with at least three-fold higher expression in cartilage as compared with other tissues. 1028 (42%) of these probesets matched genes identified in a cartilage growth plate cDNA library [4] validating their expression in cartilage via an independent dataset. The identification of genes known to have restricted patterns of expression in cartilage confirmed the presence of RNA derived from the reserve (GREM1), hypertrophic (BMP6, COL10A1), and terminally differentiated (MMP13) chondrocytes, in addition to genes expressed throughout all zones of the growth plate. This analysis suggested that there is differential transcriptional regulation of many genes in fetal cartilage and that the data could be used to identify genes selectively expressed in the tissue.

Tissue-selective genes have been previously defined as genes with enriched expression in a particular tissue [14] and characterized with algorithms dependent on the degree of differential expression relative to other tissues, including t-test [15], SAM [16], fold change [14, 17, 18], and enrichment scores [14]. While these approaches successfully identify tissue-selective genes, the reliance on fold change reduces the significance of many selectively expressed genes with low fold change. To compensate for this and identify cartilage-selective genes expressed at lower levels, the approach presented here placed increased significance on the preferentially expressed genes that showed the least variation of expression in non-cartilage tissues. This was made possible by the use of publicly-released reference gene expression data performed on the same platform and led to the reliable identification of genes with lower fold changes, but high cartilage selectivity. The impact of the use of coefficient of variation on the ranked gene list is apparent in Tables 1 and 2. In the U133A dataset, nine of the top 25 genes were ranked higher than 100 in significance in the SAM ranking. The average fold change of the probes for these nine genes was 10.7, while the average fold change of the probes for the other 14 of the top 25 genes was 42.2. One of these probes, COL10A1 was among the top four cartilage-selective genes using the CV algorithm but ranked at 284 by SAM (Table 1). In the U133B dataset, which contains a higher percentage of unannotated genes, 4 of the top 50 probes had a SAM ranking below 100, and the average SAM ranking was 576. Overall, to identify only the most cartilage-selective genes, a threshold of 50% coefficient of variation was used across all three platforms, yielding 161 cartilage-selective genes. A subset of 13 of the 161 cartilage-selective genes was studied by quantitative RT-PCR in cartilage and eight non-cartilage tissues to independently assess tissue selectivity. The data confirmed the cartilage-selectivity of genes with less than 50% CV, validating the selection procedure and suggesting that the gene expression patterns determined by microarray analysis are representative.

The coefficient of variation selection approach could, in theory, equally select for three different patterns of expression: cartilage-specific genes; genes with a consistent level of baseline expression in non-cartilage tissues; and genes with significant but equal expression in all tissues (e.g. housekeeping genes). In this data analysis, however, the most highly ranked genes consistently demonstrated little or no expression in non-cartilage tissues. The data thus demonstrate that incorporating coefficient of variation preferentially selected for genes not significantly expressed in non-cartilage tissues, yielding genes likely to have important and perhaps unique roles in cartilage.

Regardless of expression level, a cartilage-selective expression pattern suggests that the product of each identified gene may have a functional role in the development of the skeleton. Concordant with this hypothesis, mutations in 25 of the 161 selected genes have been associated with skeletal phenotypes in humans and/or mice. Included among them were the products of the well characterized genes encoding aggrecan and the cartilage-specific collagens, gene products known to have a prominent role in skeletal development and endochondral ossification. By this measure, the remaining genes may be candidate genes for skeletal dysplasias in which the disease gene has yet to be identified. As new skeletal dysplasia loci are defined, coincidence between a locus and a cartilage-selective gene may promote rapid identification of the disease gene. Knockout of the orthologous genes in mice would also facilitate exploring the role of each gene in skeletal development.

Classification of the biological roles of the products of the cartilage-selective gene set reveals genes with diverse functions including structural proteins of the cartilage extracellular matrix, enzymes that modify them, and 41 gene products with unannotated function. There were 65 genes that are components of signaling pathways, and only 43% of these were identified by sequence analysis of a comparable fetal cartilage cDNA library [4]. Among the genes were elements of the nitric oxide, VEGF, TNF/RANK, and gp130 pathways, all of which have known roles in the growth plate [1923]. Mutations in the genes encoding some of the molecules in these pathways, including RPS6SKA3, LIFR, TNFRSF11A and IKBKG, have been associated with human skeletal dysplasias [12], again suggesting that the remaining genes may also serve critical roles in endochondral ossification.

Multiple genes encoding members of the LIF/gp130 signaling pathway met the definition of cartilage-selective genes. LIF is a cytokine that is expressed in terminally-differentiated growth plate chondrocytes [24] and signals through the gp130/LIFR complex. Homozygosity for loss of function mutations in the LIF receptor produces the recessively inherited skeletal dysplasia, Stuve-Wiedemann syndrome [25]. In addition to their skeletal features, these patients have cardiovascular, pulmonary, gastrointestinal, neurologic and metabolic abnormalities, likely attributable to the role that LIFR plays in embryonic or fetal development. Genes on the cartilage-selective gene list upstream of the receptor include RELA and RELB, NF-KB survival transcription factors that increase transcription of LIF [26], as well as the LIF gene itself. Through the LIFR/gp130 complex, LIF can regulate both the JAK/STAT and ERK MAP kinase pathways. Pathway components downstream of the receptor include ATF1, part of the ATF1/CREB transcription factor complex that participates in ERK MAP kinase signaling [27, 28]. The ATF1/CREB complex is also regulated by phosphorylation by the product of the RPS6SKA3 gene [29, 30], another gene in the MAPK/ERK pathway that is associated with a skeletal phenotype. The gene encoding RPS, a phosphorylation target of RPS6SKA3 [30], was also cartilage-selective, but the role of this protein in growth plate differentiation has yet to be determined. Finally, the gene encoding FOSL1, a FOS-like transcription factor activated by the ERK/MAPK pathway which binds cJUN to form a transcription complex [31, 32], was among the cartilage-selective genes identified. Thus comparative microarray analysis has identified multiple components of a regulatory pathway that can be explored to further evaluate their importance in growth plate differentiation and endochondral ossification.

While this study has provided a deep set of genes that exhibit a cartilage-selective expression pattern, there are some limitations to the analysis. First, the study focused on total cartilage RNA, including all types of growth plate chondrocytes, as well as articular cartilage. As a result, it cannot be determined if the selected genes are expressed in all types of chondrocytes or only a subset of cells. In this context, nine of the cartilage selective genes have been shown to be more highly expressed in hypertrophic cells relative to proliferating chondrocytes in the rat and/or mouse [7, 8]. Second, the cartilage samples were derived from a single anatomic site and a narrow window of fetal development, so it is unclear to what extent the observed gene expression pattern can be generalized. Third, neither all possible non-cartilage tissues nor each type of cell within each tissue were studied, so cartilage-selectivity could be affected if additional fetal and/or adult tissues that express the identified genes were found. This may be particularly important for other connective tissues such as bone, tendon and ligament which contain cells known to express some of the cartilage-selective genes identified here.

Not all genes selectively expressed in developing cartilage will necessarily be identified using this approach. For instance, the COL2A1 gene fell just below the rigorous 50% CV standard set to define cartilage-selectivity. The underlying reasons for this are complex. Probe performance as well as the known expression of COL2A1 in fetal liver and heart, are likely to have had an effect, as both factors could have contributed to the variation in expression in non-cartilage tissues. In addition, the approach presented here treated the three expression platforms, U133A, U133B and U133 2.0 equally from the viewpoint of the threshold for cartilage-selectivity. Because the comparative dataset of normal tissues was both broader and deeper for the 133A platform, additional genes from this platform, albeit with greater than 50% CV, could be considered to be tissue-selective (e.g. COL2A1). Thus, a platform independent threshold would likely yield additional genes of interest within the U133A dataset. Finally, tissue-specific genes were identified using only microarrays and a single generalized algorithm. Additional genes selectively expressed in cartilage could be identified by less stringent criteria or other methods.


Genome-scale comparative expression analysis using human fetal cartilage and a broad set of normal human tissues has identified 161 cartilage-selective genes, including 27 uncharacterized genes. The data identify novel gene products that may provide essential roles in normal skeletogenesis and suggest new candidates for the over 100 inherited skeletal disorders in which the disease gene has not been identified. The results demonstrate that fetal cartilage is a complex and transcriptionally active tissue, and that the set of genes selectively expressed in cartilage has been greatly underestimated.


A flow chart outlining methods and results as well as other supplemental information is provided in additional data file 1.

Cartilage specimen collection and processing

Seven independent 18–22 week normal human fetal cartilage samples were studied under an Institutional Review Board approved protocol. Cartilage from the distal femur was dissected to remove bone and any adherent non-cartilage tissue (Figure 1). RNA was isolated and purified as previously described [4] and the quality and quantity of RNA were confirmed using an Agilent 2100 bioanalyzer and a Nanodrop ND-1000 spectrophotometer, respectively. Probe labeling, microarray hybridization, washing and scanning were carried out as detailed in Affymetrix protocols [33]. Five samples were used to probe Affymetrix™ U133 Plus 2.0 microarrays; and two samples were used to probe the Affymetrix™ Human Genome U133A/B set. Annotations were from version 11/15/06. The data are publicly available in the GEO database series [GEO:GSE6565]. An additional sample was fixed in formalin, sectioned and stained with toluidine blue.

Non-cartilage microarray data

This project made use of the Celsius database [34, 35], which is a database of publicly available microarray datasets from Gene Expression Omnibus, Array Express, and individual databases. Only CEL files are entered into the database, permitting reprocessing using identical algorithms to enable experimental comparisons. Only data from Affymetrix™ Human Genome U133A/B and Plus 2.0 platforms that contained clear annotation that they were derived from normal human tissues were selected for this analysis.

Microarray analysis

Data normalization and transformation

Raw data were normalized using the RMA algorithm with default parameters, available as part of the Bioconductor R library [36, 37]. In brief, each CEL file was processed separately with an invariant pool of 50 arrays from a matching platform. Higher signal intensities observed in a subset of U133A non-cartilage samples from one provider [38] were additionally normalized by subtracting the median from other non-cartilage samples. After normalization, the training dataset was log2 transformed prior to analysis.

Median derived analog of CV for cartilage-selectivity ranking

In highly expressed cartilage genes, the degree of cartilage selectivity was defined as a median derived analog of CV (average deviation/median) applied to expression of these genes in non-cartilage tissues. The analog of CV was used to allow for greater tolerance of gene expression in some cartilage containing tissues without affecting the cartilage-selective assessment. A CV of 50% was empirically determined as a mathematically acceptable threshold for cartilage selectivity for probes across all three validation datasets (U133A, U133B, and U133 Plus 2.0).


For the unsupervised analysis, probe intensities in all tissues were subtracted by the median probe intensity in cartilage, so all expression was defined relative to cartilage (i.e. the median of cartilage expression was set at zero) for each selected probe. Probesets with the greatest variation across all tissues and whose expression in at least two samples differed by two standard deviations from their mean expression across the entire set of samples were selected. Two-way hierarchical clustering was performed using Pearson's correlation to group genes and arrays based on the similarity of their expression patterns [39]. For the supervised analysis, the significance analysis of microarrays (SAM) two class method [40] was applied. 100 hundred permutations were used and at least three-fold variation between cartilage and non-cartilage expression was required.


The probeset expression profiles for the 2,446 probesets identified in the training set (see Results) were acquired from 224 arrays representing 34 different tissues on three different platforms: Affymetrix U133A, U133B and U133 Plus 2.0. The data from the first platform, U133A, consisted of 1363 probeset profiles from 124 arrays, representing two normal fetal cartilage and 122 normal non-cartilage (32 tissues) samples (see additional data file 3, for the tissue distribution of samples used in this analysis). Arrays represented in this dataset were mostly from two large normal tissue expression profiling projects [38]. The second dataset, U133B, consisted of 882 expression profiles identified on 72 U133B microarrays from two normal fetal cartilage and 74 normal non-cartilage samples. These samples were primarily from the UCLA normal tissue microarray project (Chen, Day and Nelson, unpublished). The U133 2.0 dataset consisted of 201 expression profiles obtained from 26 Human Genome U133 Plus 2.0 arrays representing expression from five normal fetal cartilage and 21 non-cartilage (eight different tissue types) samples. The five cartilage arrays for the U133 2.0 platform were technical replicates of the arrays used in the training dataset. The non-cartilage arrays were a subset of the training dataset set aside for this validation only. From the 2446 probesets selected, probesets that exhibited at least a five-fold difference between the average cartilage intensity and the median signal intensity of all non-cartilage tissues, were selected. Cartilage-specificity was then determined using a median-derived analog of coefficient of variation (CV) as described above. Probesets with less than 50% CV were defined as reflecting cartilage-selective expression.


One microgram of RNA from seven tissues (brain, prostate, kidney, liver, heart, thyroid, and testis) in the FirstChoice® Human Total RNA survey panel (Ambion) was reverse transcribed using a high-capacity cDNA archive kit (ABI) and random primers. For cartilage, RNA from three independent cartilage samples was pooled and reverse transcribed. Amplification reactions were performed in triplicate using 100 ng of each cDNA. Thirty-five cycles of amplification were carried out in an ABI 7300 using the validated QuantiTect Gene Expression Assays and SYBR Green PCR kit (Qiagen). To assess specificity, amplification products were subjected to melting curve analysis and gel electrophoresis. The 2- [delta] [delta]Ct method was employed to calculate relative amplification. This was performed using an average of endogenous references (18S, GAPDH, and HPRT1) to improve normalization across the panel of tissues used [41]. For genes where no amplification was detected in a tissue, a Ct value of 35 was assigned, reflecting the maximum number of cycles carried out.



coefficient of variation


significance analysis of microarrays


  1. 1.

    Goldring MB, Tsuchimochi K, Ijiri K: The control of chondrogenesis. J Cell Biochem. 2006, 97 (1): 33-44. 10.1002/jcb.20652.

    CAS  PubMed  Article  Google Scholar 

  2. 2.

    Pacifici M, Koyama E, Iwamoto M, Gentili C: Development of articular cartilage: what do we know about it and how may it occur?. Connect Tissue Res. 2000, 41 (3): 175-184.

    CAS  PubMed  Article  Google Scholar 

  3. 3.

    Pogue R, Sebald E, King L, Kronstadt E, Krakow D, Cohn DH: A transcriptional profile of human fetal cartilage. Matrix Biol. 2004, 23 (5): 299-307. 10.1016/j.matbio.2004.07.003.

    CAS  PubMed  Article  Google Scholar 

  4. 4.

    Krakow D, Sebald ET, Pogue R, Rimoin LP, King L, Cohn DH: Analysis of clones from a human cartilage cDNA library provides insight into chondrocyte gene expression and identifies novel candidate genes for the osteochondrodysplasias. Mol Genet Metab. 2003, 79 (1): 34-42. 10.1016/S1096-7192(03)00064-7.

    CAS  PubMed  Article  Google Scholar 

  5. 5.

    Zhang H, Marshall KW, Tang H, Hwang DM, Lee M, Liew CC: Profiling genes expressed in human fetal cartilage using 13,155 expressed sequence tags. Osteoarthritis Cartilage. 2003, 11 (5): 309-319. 10.1016/S1063-4584(03)00032-3.

    CAS  PubMed  Article  Google Scholar 

  6. 6.

    Tagariello A, Schlaubitz S, Hankeln T, Mohrmann G, Stelzer C, Schweizer A, Hermanns P, Lee B, Schmidt ER, Winterpacht A, Zabel B: Expression profiling of human fetal growth plate cartilage by EST sequencing. Matrix Biol. 2005, 24 (8): 530-538. 10.1016/j.matbio.2005.08.002.

    CAS  PubMed  Article  Google Scholar 

  7. 7.

    Wang Y, Middleton F, Horton JA, Reichel L, Farnum CE, Damron TA: Microarray analysis of proliferative and hypertrophic growth plate zones identifies differentiation markers and signal pathways. Bone. 2004, 35 (6): 1273-1293. 10.1016/j.bone.2004.09.009.

    CAS  PubMed  Article  Google Scholar 

  8. 8.

    James CG, Appleton CT, Ulici V, Underhill TM, Beier F: Microarray analyses of gene expression during chondrocyte differentiation identifies novel regulators of hypertrophy. Mol Biol Cell. 2005, 16 (11): 5316-5333. 10.1091/mbc.E05-01-0084.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  9. 9.

    Stokes DG, Liu G, Coimbra IB, Piera-Velazquez S, Crowl RM, Jimenez SA: Assessment of the gene expression profile of differentiated and dedifferentiated human fetal chondrocytes by microarray analysis. Arthritis Rheum. 2002, 46 (2): 404-419. 10.1002/art.10106.

    CAS  PubMed  Article  Google Scholar 

  10. 10.

    Lohmander LS: What can we do about osteoarthritis?. Arthritis Res. 2000, 2 (2): 95-100. 10.1186/ar74.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  11. 11.

    Reginato AM, Olsen BR: The role of structural genes in the pathogenesis of osteoarthritic disorders. Arthritis Res. 2002, 4 (6): 337-345. 10.1186/ar595.

    PubMed Central  PubMed  Article  Google Scholar 

  12. 12.

    Superti-Furga A, Unger S: Nosology and classification of genetic skeletal disorders: 2006 revision. Am J Med Genet A. 2006

    Google Scholar 

  13. 13.

    Brazma A, Vilo J: Gene expression data analysis. FEBS Lett. 2000, 480 (1): 17-24. 10.1016/S0014-5793(00)01772-5.

    CAS  PubMed  Article  Google Scholar 

  14. 14.

    Liang S, Li Y, Be X, Howes S, Liu W: Detecting and profiling tissue-selective genes. Physiol Genomics. 2006, 26 (2): 158-162. 10.1152/physiolgenomics.00313.2005.

    CAS  PubMed  Article  Google Scholar 

  15. 15.

    Hsiao LL, Dangond F, Yoshida T, Hong R, Jensen RV, Misra J, Dillon W, Lee KF, Clark KE, Haverty P, Weng Z, Mutter GL, Frosch MP, Macdonald ME, Milford EL, Crum CP, Bueno R, Pratt RE, Mahadevappa M, Warrington JA, Stephanopoulos G, Stephanopoulos G, Gullans SR: A compendium of gene expression in normal human tissues. Physiol Genomics. 2001, 7 (2): 97-104.

    CAS  PubMed  Article  Google Scholar 

  16. 16.

    Shyamsundar R, Kim YH, Higgins JP, Montgomery K, Jorden M, Sethuraman A, van de Rijn M, Botstein D, Brown PO, Pollack JR: A DNA microarray survey of gene expression in normal human tissues. Genome Biol. 2005, 6 (3): R22-10.1186/gb-2005-6-3-r22.

    PubMed Central  PubMed  Article  Google Scholar 

  17. 17.

    Zhao SH, Recknor J, Lunney JK, Nettleton D, Kuhar D, Orley S, Tuggle CK: Validation of a first-generation long-oligonucleotide microarray for transcriptional profiling in the pig. Genomics. 2005, 86 (5): 618-625. 10.1016/j.ygeno.2005.08.001.

    CAS  PubMed  Article  Google Scholar 

  18. 18.

    Saito-Hisaminato A, Katagiri T, Kakiuchi S, Nakamura T, Tsunoda T, Nakamura Y: Genome-wide profiling of gene expression in 29 normal human tissues with a cDNA microarray. DNA Res. 2002, 9 (2): 35-45. 10.1093/dnares/9.2.35.

    CAS  PubMed  Article  Google Scholar 

  19. 19.

    Teixeira CC, Ischiropoulos H, Leboy PS, Adams SL, Shapiro IM: Nitric oxide-nitric oxide synthase regulates key maturational events during chondrocyte terminal differentiation. Bone. 2005, 37 (1): 37-45. 10.1016/j.bone.2005.03.010.

    CAS  PubMed  Article  Google Scholar 

  20. 20.

    Sims NA, Jenkins BJ, Quinn JM, Nakamura A, Glatt M, Gillespie MT, Ernst M, Martin TJ: Glycoprotein 130 regulates bone turnover and bone size by distinct downstream signaling pathways. J Clin Invest. 2004, 113 (3): 379-389. 10.1172/JCI200419872.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  21. 21.

    Zelzer E, Mamluk R, Ferrara N, Johnson RS, Schipani E, Olsen BR: VEGFA is necessary for chondrocyte survival during bone development. Development. 2004, 131 (9): 2161-2171. 10.1242/dev.01053.

    CAS  PubMed  Article  Google Scholar 

  22. 22.

    Kanegae Y, Tavares AT, Izpisua Belmonte JC, Verma IM: Role of Rel/NF-kappaB transcription factors during the outgrowth of the vertebrate limb. Nature. 1998, 392 (6676): 611-614. 10.1038/33429.

    CAS  PubMed  Article  Google Scholar 

  23. 23.

    Wu S, De Luca F: Inhibition of the Proteasomal Function in Chondrocytes Down-Regulates Growth Plate Chondrogenesis and Longitudinal Bone Growth. Endocrinology. 2006, 147 (8): 3761-3768. 10.1210/en.2005-1672.

    CAS  PubMed  Article  Google Scholar 

  24. 24.

    Grimaud E, Blanchard F, Charrier C, Gouin F, Redini F, Heymann D: Leukaemia inhibitory factor (lif) is expressed in hypertrophic chondrocytes and vascular sprouts during osteogenesis. Cytokine. 2002, 20 (5): 224-230. 10.1006/cyto.2002.2002.

    CAS  PubMed  Article  Google Scholar 

  25. 25.

    Dagoneau N, Scheffer D, Huber C, Al-Gazali LI, Di Rocco M, Godard A, Martinovic J, Raas-Rothschild A, Sigaudy S, Unger S, Nicole S, Fontaine B, Taupin JL, Moreau JF, Superti-Furga A, Le Merrer M, Bonaventure J, Munnich A, Legeai-Mallet L, Cormier-Daire V: Null leukemia inhibitory factor receptor (LIFR) mutations in Stuve-Wiedemann/Schwartz-Jampel type 2 syndrome. Am J Hum Genet. 2004, 74 (2): 298-305. 10.1086/381715.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  26. 26.

    Leung TH, Hoffmann A, Baltimore D: One nucleotide in a kappaB site can determine cofactor specificity for NF-kappaB dimers. Cell. 2004, 118 (4): 453-464. 10.1016/j.cell.2004.08.007.

    CAS  PubMed  Article  Google Scholar 

  27. 27.

    Belmonte N, Phillips BW, Massiera F, Villageois P, Wdziekonski B, Saint-Marc P, Nichols J, Aubert J, Saeki K, Yuo A, Narumiya S, Ailhaud G, Dani C: Activation of extracellular signal-regulated kinases and CREB/ATF-1 mediate the expression of CCAAT/enhancer binding proteins beta and -delta in preadipocytes. Mol Endocrinol. 2001, 15 (11): 2037-2049. 10.1210/me.15.11.2037.

    CAS  PubMed  Google Scholar 

  28. 28.

    Gliki G, Abu-Ghazaleh R, Jezequel S, Wheeler-Jones C, Zachary I: Vascular endothelial growth factor-induced prostacyclin production is mediated by a protein kinase C (PKC)-dependent activation of extracellular signal-regulated protein kinases 1 and 2 involving PKC-delta and by mobilization of intracellular Ca2+. Biochem J. 2001, 353 (Pt 3): 503-512. 10.1042/0264-6021:3530503.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  29. 29.

    Gupta P, Prywes R: ATF1 phosphorylation by the ERK MAPK pathway is required for epidermal growth factor-induced c-jun expression. J Biol Chem. 2002, 277 (52): 50550-50556. 10.1074/jbc.M209799200.

    CAS  PubMed  Article  Google Scholar 

  30. 30.

    Wang Y, Prywes R: Activation of the c-fos enhancer by the erk MAP kinase pathway through two sequence elements: the c-fos AP-1 and p62TCF sites. Oncogene. 2000, 19 (11): 1379-1385. 10.1038/sj.onc.1203443.

    CAS  PubMed  Article  Google Scholar 

  31. 31.

    Owens JM, Matsuo K, Nicholson GC, Wagner EF, Chambers TJ: Fra-1 potentiates osteoclastic differentiation in osteoclast-macrophage precursor cell lines. J Cell Physiol. 1999, 179 (2): 170-178. 10.1002/(SICI)1097-4652(199905)179:2<170::AID-JCP7>3.0.CO;2-K.

    CAS  PubMed  Article  Google Scholar 

  32. 32.

    Eferl R, Hoebertz A, Schilling AF, Rath M, Karreth F, Kenner L, Amling M, Wagner EF: The Fos-related antigen Fra-1 is an activator of bone matrix formation. Embo J. 2004, 23 (14): 2789-2799. 10.1038/sj.emboj.7600282.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  33. 33.

    Affymetrix Technical Support Protocols. []

  34. 34.

    Day A, Carlson MRJ, Dong J, O'Connor BD, Nelson SF: Celsius: A community resource for Affymetrix microarry data . Genome Biology.

  35. 35.

    Celsius Microarray Database. []

  36. 36.

    Bolstad BM, Irizarry RA, Astrand M, Speed TP: A comparison of normalization methods for high density oligonucleotide array data based on variance and bias. Bioinformatics. 2003, 19 (2): 185-193. 10.1093/bioinformatics/19.2.185.

    CAS  PubMed  Article  Google Scholar 

  37. 37.

    Gentleman RC, Carey VJ, Bates DM, Bolstad B, Dettling M, Dudoit S, Ellis B, Gautier L, Ge Y, Gentry J, Hornik K, Hothorn T, Huber W, Iacus S, Irizarry R, Leisch F, Li C, Maechler M, Rossini AJ, Sawitzki G, Smith C, Smyth G, Tierney L, Yang JY, Zhang J: Bioconductor: open software development for computational biology and bioinformatics. Genome Biol. 2004, 5 (10): R80-10.1186/gb-2004-5-10-r80.

    PubMed Central  PubMed  Article  Google Scholar 

  38. 38.

    Su AI, Wiltshire T, Batalov S, Lapp H, Ching KA, Block D, Zhang J, Soden R, Hayakawa M, Kreiman G, Cooke MP, Walker JR, Hogenesch JB: A gene atlas of the mouse and human protein-encoding transcriptomes. Proc Natl Acad Sci U S A. 2004, 101 (16): 6062-6067. 10.1073/pnas.0400782101.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  39. 39.

    Spellman PT, Sherlock G, Zhang MQ, Iyer VR, Anders K, Eisen MB, Brown PO, Botstein D, Futcher B: Comprehensive identification of cell cycle-regulated genes of the yeast Saccharomyces cerevisiae by microarray hybridization. Mol Biol Cell. 1998, 9 (12): 3273-3297.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  40. 40.

    Tusher VG, Tibshirani R, Chu G: Significance analysis of microarrays applied to the ionizing radiation response. Proc Natl Acad Sci U S A. 2001, 98 (9): 5116-5121. 10.1073/pnas.091062498.

    CAS  PubMed Central  PubMed  Article  Google Scholar 

  41. 41.

    Vandesompele J, De Preter K, Pattyn F, Poppe B, Van Roy N, De Paepe A, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3 (7): RESEARCH0034-10.1186/gb-2002-3-7-research0034.

    PubMed Central  PubMed  Article  Google Scholar 

Download references


The authors thank Louis Fridkis and Brian O'Connor for their contributions to the CELSIUS database. The authors also thank the UCLA microarray core for assistance with the generation and analysis of the microarray data. This work was supported in part by grants from the NIH (HD22657 and RR00425 to DHC and DK) and (HL072367 and U24NS052108 to SFN) and DK was supported by the Joseph Drown Foundation. AD was supported by a grant from the NSF UCLA-IGERT (DGE-9987641). DHC and DK are recipients of Winnick Family Clinical Scholars Awards.

Author information



Corresponding author

Correspondence to Daniel H Cohn.

Additional information

Authors' contributions

VF designed and carried out the bioinformatics analysis and qRT-PCR, and drafted the manuscript. AD participated in data acquisition and normalization. DK obtained, dissected and isolated the cartilage RNA and critically revised the manuscript, ZAC participated in the bioinformatics analysis. ZC participated in generation of microarray data. SN participated in the design and coordination of the study and critically revised the manuscript. DC conceived the study, participated in its design and coordination, performed microarray analysis and critically revised the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

Additional File 1: A flow chart illustrating a summary of the analysis. An unsupervised (black arrows) and supervised analysis (blue arrows) were performed with gene expression from 46 U133 2.0 Affymetrix arrays. An independent validation set comprised of 224 Affymetrix arrays (dashed arrows) was also used to test the 1713 genes for the most robust fetal cartilage selective genes. (PDF 14 KB)

Additional File 2: Supervised analysis and summary of in silico validation. Ranked list of cartilage genes more highly expressed by at least three fold in cartilage than non-cartilage tissues in the training dataset (five fetal cartilage samples compared to 41 normal non-cartilage samples). Ranked order is based on expression profiles obtained from 46 U133 Plus 2.0 arrays analyzed with SAM 2 class analysis with 100 permutations and with a False Discovery Rate (FDR) of 0. (B) Each probeset was independently evaluated in the validation datasets five fold higher expression using independent samples and three independent platforms as outlined in methods. Present indicates probe is identified in validation platform; Enriched indicates gene is expressed five fold higher in cartilage than non-cartilage tissues; Cartilage selectivity indicates at least five fold higher expression in cartilage than non-cartilage with a CV score of < 50% in non-cartilage samples. (C) "X" denotes gene was identified in a fetal cartilage cDNA library. (PDF 98 KB)

Additional File 3: Tissue distribution training and validation sets. (A) 31 Non-cartilage tissues and 124 arrays were used for the validation of cartilage selective genes identified on the U133A chip. Two fetal cartilage samples were compared against 122 non-cartilage samples. (B) 27 non-cartilage Tissues and 74 arrays were used for the validation of cartilage selective genes identified on the U133B chip. Two fetal cartilage samples were compared against 72 non-cartilage arrays (C) Eight non-cartilage tissues and 26 arrays used for the validation of cartilage selective genes identified on the U133B chip. Five fetal cartilage samples were compared against 72 non-cartilage arrays. (PDF 8 KB)

Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

Authors’ original file for figure 1

Authors’ original file for figure 2

Rights and permissions

Reprints and Permissions

About this article

Cite this article

Funari, V.A., Day, A., Krakow, D. et al. Cartilage-selective genes identified in genome-scale analysis of non-cartilage and cartilage gene expression. BMC Genomics 8, 165 (2007).

Download citation


  • Growth Plate
  • Additional Data File
  • Endochondral Ossification
  • Skeletal Dysplasia
  • Growth Plate Chondrocytes