Abstract
Background:
Although osteoarthritis (OA) is the commonest joint disorder and has a rising prevalence as the population ages, no drugs are available that prevent or delay the onset and progression of disease. Recent studies identified the DIO2 gene encoding type 2 deiodinase (D2) as a susceptibility locus for OA, and further data suggest deiodinase-regulated local availability of triiodothyronine (T3) in the joint plays an important role in cartilage maintenance and repair. To investigate the hypothesis that reduced tissue T3 availability protects joints from development of OA, the joint phenotypes of adult mice lacking D2 (D2KO) or lacking both D1 and D2 (D1D2KO), the only enzymes that catalyze conversion of the prohormone thyroxine to active T3, were determined.
Methods:
Knee joints were prepared from male 16-week-old adult wild type (WT; n=9), D2KO (n=5), and D1D2KO (n=3) mice. Articular cartilage pathology was scored using the Osteoarthritis Research Society International (OARSI) histopathology scale for murine OA to determine the severity and extent of disease. Digital X-ray microradiography was used to determine the area and mineral content of subchondral bone immediately beneath the articular cartilage surface.
Results:
There were no differences in maximum and standardized OA scores, cartilage erosion indices, or articular cartilage cellularity among WT, D2KO, and D1D2KO mice. Subchondral bone area did not differ among genotypes, but mineral content was markedly increased in both D2KO and D1D2KO mice compared to WT.
Conclusions:
Although adult D2KO mice have normal articular cartilage and no other features of spontaneous joint damage, they exhibit increased subchondral bone mineral content.
Introduction
O
Recent data suggest that altered thyroid hormone availability in articular cartilage or joint tissues may play a key role in predisposition to OA (5 –9). Normal euthyroid status is maintained by the hypothalamic-pituitary-thyroid axis via a negative feedback loop in which thyroid hormones produced by the thyroid gland inhibit synthesis and secretion of thyrotropin (TSH) from the anterior pituitary. The main hormone secreted by the thyroid gland is the prohormone thyroxine (T4), which is converted to the physiologically active metabolite 3,5,3′-L-triiodothyronine (T3) by removal of an outer ring iodine atom. Thyroid hormones enter peripheral tissues, including bone and cartilage (10), via specific membrane transporters. In target cells, the type 2 iodothyronine deiodinase (D2) converts T4 to T3, whereas the type 3 enzyme (D3) irreversibly inactivates T3 and prevents activation of T4 by removal of an inner ring iodine. Thus, the balance of D2 and D3 activities determines local T3 availability and tissue thyroid hormone responsiveness (11 –13).
The D2 enzyme is labile due to substrate-induced ubiquitination and proteasome-mediated degradation (12,13), and its specific activity is difficult to measure, especially at low levels. Nevertheless, in-depth studies of iodothyronine metabolism in mouse skeletal tissues, cell lysates, and intact primary cell cultures have shown that specific D2 enzyme activity is restricted to mature bone-forming osteoblasts and not detectable in bone-resorbing osteoclasts or growth plate chondrocytes (10). A recent study using quantitative reverse transcription polymerase chain reaction (RT-qPCR) identified Dio2 mRNA expression in rat articular cartilage at levels comparable to brown adipose tissue, although D2 enzyme activity was not measured (14). By contrast, functional D3 enzyme, a much more stable protein with readily detectable activity, was identified in osteoblasts, osteoclasts, and growth plate chondrocytes, with activity highest in cells derived from developing and younger animals (10). It is unknown whether D3 is expressed in articular cartilage or whether D2 expression in articular cartilage is regulated during joint development. The type 1 deiodinase (D1), which converts T4 to T3 mainly in the liver and kidney and generates a large proportion of the circulating T3, is not expressed in bone or cartilage and appears to have no physiological role in the skeleton (10).
T3 regulates cartilage maturation and is essential for endochondral ossification, bone modeling, and longitudinal growth (15). A role for thyroid hormones in the pathogenesis of OA has been proposed. Although chondrocytes resist terminal differentiation in healthy articular cartilage, a process resembling endochondral ossification occurs during the progression of OA in which articular chondrocytes undergo hypertrophic differentiation with acceleration of cartilage mineralization (16,17). Increased activities of ADAMTS5 (a disintegrin and metalloproteinase with thrombospondin motifs-5) (18,19) and MMP13 (matrix metalloproteinase-13) (20) in OA cause degradation of cartilage matrix, and these enzymes are regulated in the growth plate by T3 (21 –24), which also induces hypertrophic chondrocyte differentiation at this location (14,21,25).
A genome-wide linkage analysis of siblings with generalized OA found a nonsynonymous coding variant (rs225014; Thr92Ala) that identified DIO2 as a disease-susceptibility locus (6). Replication studies in Asian and European cohorts confirmed this by identifying a DIO2 haplotype containing the minor allele of rs225014 and major allele of rs12885300 that was associated with symptomatic OA (6). The relationship between this haplotype and OA, however, was not replicated in the Rotterdam Study population (26), and DIO2 was not identified as a susceptibility locus in a more recent meta-analysis (27). Nevertheless, further data suggest the risk allele of rs225014 may be expressed at higher levels than the protective allele in OA cartilage from heterozygous patients (28), possibly because of epigenetic modifications (29). Increased DIO2 mRNA and protein expression has also been documented by RT-qPCR and immunohistochemical staining in cartilage obtained from joints affected by end-stage OA (14,28). Additionally, the rs12885300 DIO2 polymorphism has been suggested to influence the association between hip shape and OA susceptibility by possibly increasing vulnerability of articular cartilage to abnormal hip morphology (30). Furthermore, studies in transgenic rats that overexpress human D2 in chondrocytes indicate transgenic animals have an increased susceptibility to OA following surgical destabilization of the knee (14). Finally, a recent meta-analysis studying factors that regulate thyroid hormone metabolism and action identified a possible role for DIO3 as a disease-modifying locus in OA (5).
Together, these data suggest that deiodinase-regulated local availability of T3 in joint tissues may play an important role in articular cartilage renewal and repair, particularly as the polymorphisms in DIO2 and DIO3 that are associated with OA do not influence systemic thyroid status (31). Overall, current findings suggest that increased T3 availability in OA may be detrimental to joint maintenance and articular cartilage homeostasis. Thus, it was hypothesized that reduced tissue T3 availability protects joints from damage and the development of OA. To investigate, the joint phenotypes of adult mice lacking D2 (D2KO) or lacking both D1 and D2 (D1D2KO) were determined.
Materials and Methods
Mice
D1KO or D2KO mice were generated after targeted disruption of Dio1 or Dio2 in 129/SvJ embryonic stem cells. Mice were bred and maintained as homozygous colonies in a C57BL/6 background at Dartmouth Medical School animal facility using protocols approved by the Institutional Review Board (32,33). D1KO and D2KO mice were crossed to obtain double heterozygous F1 progeny, which were used to generate a colony of D1D2KO homozygous double mutants (12,34). Mice were housed under conditions of controlled lighting and temperature and sacrificed at 16 weeks of age. Bone samples were shipped to Imperial College London for phenotyping studies in which nine wild type (WT), five D2KO, and three D1D2KO male mice were used for detailed histological analysis of the joints. Histological sections from four wild type, four D2KO, and three D1D2KO mice were prepared at directly comparable levels for subchondral bone analysis.
Preparation of histological sections
Lower limbs were fixed in 10% neutral buffered formalin for 24 h and decalcified in 10% EDTA, pH7.4, for 7–10 days. Knee joints were embedded in paraffin and 5 μm thick midline sections were cut at 80 μm intervals. Slides were stained for cartilage proteoglycans using Safranin-O Fast-Green (35) and imaged using a Leica MZ75 binocular microscope (Leica AG, Heerbrugg, Switzerland), KL1500 light source, and DFC320 digital camera.
Analysis of articular cartilage
Articular cartilage pathology was scored by a blinded observer using the Osteoarthritis Research Society International (OARSI) histopathology scale for murine OA (35). A score of 0 represents normal cartilage; 0.5 represents loss of Safranin-O staining without structural changes; 1 represents small fibrillations without loss of cartilage; and 2 represents vertical clefts down to the layer immediately below the superficial layer and some loss of surface lamina. Vertical clefts or erosions to the calcified cartilage extending to <25%, 25–50%, 50–75%, or >75% of the articular surface are assigned scores of 3, 4, 5, or 6 respectively. Using this system, the severity and extent of osteoarthritis is expressed as a maximum score for each joint and as a cumulative score over the whole joint (35). Because the numbers of histological sections obtained differed between individual mice, a standardized OA score was determined. Thus, a reference histological atlas of a 16-week-old male C57BL/6 WT mouse knee was constructed using 5 μm thick coronal sections cut at 80 μm intervals through the intact joint to yield 16 levels of sectioning, of which nine passed directly through the joint cavity (Supplementary Fig. S1; Supplementary Data are available online at

Articular cartilage and subchondral bone analysis in wild type (WT) C57BL/6 mice. (
Histological analysis of subchondral bone area
Histological sections were imaged in 8-bit Tiff format. The subchondral bone region of interest (ROI) between the articular cartilage and epiphyseal growth plate of the tibia was contoured using ImageJ (
Subchondral bone mineral content
Quantitative digital X-ray microradiography was used to determine relative bone mineral content in tibial subchondral bone as described (38,39). Increasing gradations of mineralization density were represented in 16 equal intervals by a pseudocolor scheme (Fig. 2). The bone compartment proximal to the tibial growth plate was selected as the subchondral bone ROI. Relative and cumulative frequency distributions of gray levels in ROIs were determined and compared among genotypes.

Analysis of subchondral bone mineral content in WT C57BL/6 mice. (
Statistics
Individual comparisons between WT and mutant mice for the maximum OA score and range were analyzed using the Mann–Whitney U-test. Comparisons between WT and mutant mice for mean standardized OA score, erosion index, and subchondral bone area were analyzed by unpaired Student's t-test. p-Values<0.05 were considered significant. The Kolmogorov–Smirnov test was used to compare gray-level cumulative frequency distributions in digital X-ray micrographs between WT and mutant mice. p-Values for the D statistic were D>6.01, p<0.05; D≥7.20, p<0.01; and D≥8.62, p<0.001 (39).
Results
Length of articular cartilage
The articular cartilage surface in the reference atlas was determined in sections from the nine levels that passed directly through the joint cavity (Supplementary Fig. S1 and Fig. 1). The measured complete articular cartilage surface in the reference joint was 27.4 mm, and this length was used to calculate standardized OA scores. In sections from experimental samples, the measured articular cartilage surface was: WT 5.0±1.2, n=9; D2KO 6.0±1.1, n=5; D1D2KO 10.99±2.99, n=3 (mm±SEM; ANOVA; not significant).
Articular cartilage is not altered in D2KO mice
There were no differences in maximum and standardized OA scores or erosion indices between WT and D2KO mice (Fig. 3). The maximum OA scores were: WT 3 (0–3), n=9; D2KO 4 (0.5–4), n=5 (maximum score (range); Mann–Whitney U-test; not significant). The standardized OA scores were: WT 7.8±2.1, n=9; D2KO 10.9±3.6, n=5 (mean standardized score±SEM; Student's t-test; not significant). The erosion indices were WT 2.7±1.4, n=9; D2KO 6.1±2.7, n=5 (mean %±SEM; Student's t-test; not significant). There were also no differences in these parameters in D1D2KO mice: maximum OA score (0–2), n=3; standardized OA score 9.6±5.4, n=3; erosion index 5.5±3.3, n=3 (Mann–Whitney U-test and Student's t-tests vs. WT or D2KO mice; not significant).

Osteoarthritis (OA) pathology and scores in WT and D2KO mice. (
Articular cartilage cellularity also did not differ between WT and D2KO mice: WT 1412±82, n=9; D2KO 1705±161, n=5 (mean number of chondrocytes per mm2±SEM; Student's t-test; not significant). Articular cartilage cellularity also did not differ in D1D2KO mice: 1493±251, n=3 (Student's t-test vs. WT or D2KO mice; not significant).
Subchondral bone area is unchanged but bone mineral content is increased in D2KO mice
The subchondral bone area did not differ between WT and D2KO mice: WT 64±7.6%, n=4; D2KO 56±10.3%, n=4 (mean BA/TA; Student's t-test; not significant; Fig. 3) and was also not different in D1D2KO mice: 52±3.8%, n=3 (Student's t-test vs. WT or D2KO mice; not significant). Despite this, subchondral bone mineral content was markedly increased in D2KO mice (n=4) compared to WT (n=5), Kolmogorov–Smirnov test (p<0.001; Fig. 4), and a similar increase was evident in D1D2KO mice (n=2; p<0.001).

Subchondral bone mineral content in WT and D2KO mice. (
Discussion
It is demonstrated that 16-week-old adult male D2KO mice, which have a normal circulating T3 concentration (33) but lack the D2 enzyme that regulates intracellular T3 availability, have increased subchondral bone mineral content but normal articular cartilage and no other features of OA. Additional deletion of D1 in D1D2KO mice, which also have a normal circulating T3 concentration (12), did not affect the phenotype observed in D2KO mice, supporting the notion that D1 has no physiological role in the skeleton (10).
D2KO mice have a generalized increase in bone mass and mineralization resulting in bones with reduced toughness and increased susceptibility to fracture (34). The abnormality results from disruption of a critical function for D2 in osteoblasts, in which a higher T3 concentration relative to other skeletal cells is maintained in order to optimize bone mineralization (34). In this study, it is demonstrated that adult D2KO mice also have increased mineral content in subchondral bone but no abnormalities of articular cartilage or other pathological features of spontaneous OA.
It is important to consider whether the subchondral bone abnormality in D2KO mice reflects a potential role for D2 in the pathogenesis of OA. This is a complex question because OA may result from primary abnormalities in articular chondrocytes, changes in adjacent subchondral bone that lead to secondary cartilage damage and degeneration, or from defects in both tissues. Although it is well established that both subchondral bone remodeling and articular cartilage damage occur in OA, it is not clear whether these alterations occur simultaneously or whether subchondral bone changes result in cartilage damage or vice versa (3,9,40,41). For example, cathepsin K knockout mice have increased subchondral bone mineralization and are protected from OA (42), whereas rabbits treated with parathyroid hormone have increased subchondral bone mineral density but worsened articular cartilage integrity (43). The current findings demonstrate abnormal subchondral bone in the absence of articular cartilage damage, suggesting that any effect of D2 deletion on OA pathogenesis would likely be mediated via primary effects in subchondral bone osteoblasts.
In this respect, it is essential to consider the cells in which D2 enzyme activity is actually present; chondrocytes, osteoblasts, or both cell types. In detailed studies demonstrating that functional D2 enzyme is restricted to bone-forming osteoblasts, it was shown that measurable D2 mRNA expression does not necessarily correlate with detectable enzyme activity (10). An important reason for this discrepancy is that D2 protein is labile, undergoing rapid degradation following exposure to increasing concentrations of T4 in a local negative feedback loop (13,44). A further difficulty is that available antibodies for detection of D2 by immunohistochemistry or Western blotting are insufficiently sensitive or unreliable. It has not been established, therefore, whether functional D2 is genuinely expressed in cartilage. Thus, Dio2 mRNA expression was demonstrated in growth plate cartilage but measureable enzyme activity using a high sensitivity assay could not be detected (10), whereas others have shown Dio2 mRNA expression in rat articular cartilage (14) and increased levels of Dio2 mRNA (14) and protein (28) in human articular cartilage from joints resected for end-stage OA, although in each case enzyme activity was not determined. The significance of apparently increased D2 expression in end-stage OA cartilage is therefore uncertain; it is also not known whether increased D2 expression would represent a secondary response to joint destruction or whether it might precede cartilage damage and be a causative factor in disease progression.
In this context, transgenic rats overexpressing Dio2 in chondrocytes had increased susceptibility to OA following provocation by surgical destabilization of the knee (14). Nevertheless, a causal relationship between increased D2 expression and OA susceptibility was not established because D2 activity was not determined in transgenic rats overexpressing the enzyme (14). This is an important point because in a previous study, transgenic mice overexpressing D2 in the heart surprisingly displayed only mild thyrotoxic changes in cardiomyocytes (45), presumably because the phenotype was mitigated by T4-induced degradation of the overexpressed enzyme (44). Furthermore, in a recent study, siRNA-mediated inhibition of DIO2 mRNA expression in primary human chondrocytes resulted in decreased expression of liver X receptor α (LXRα) and increased expression of interleukin-1β (IL1β) and IL1β-induced cyclooxygenase-2 (COX2) expression (46). At odds with the notion that deletion or inhibition of D2 might protect against OA, these findings indicate that suppression of D2 results in a pro-inflammatory response in cartilage. The authors suggest that D2 has an anti-inflammatory role in chondrocytes resulting in inactivation of LXRα and suppression of IL-1β and COX2 (46). Thus, the apparent increase in D2 expression observed in end-stage human OA (14,28) could represent a consequence of disease progression resulting from activation of pro-inflammatory pathways such as the nuclear factor κB (NFκB) pathway, which is known to stimulate D2 expression (47,48).
In light of these conflicting findings, it is prudent to reconsider the original hypothesis that reduced tissue T3 availability protects joints from damage and the development of OA. Genetic studies proposed DIO2 as a susceptibility gene for OA based on identification of the rs225014 (Thr92Ala) polymorphism. The functional effect of this polymorphism on enzyme activity is, however, uncertain. Initial in vitro transfection studies demonstrated the Thr92Ala substitution in D2 did not affect protein half-life, enzyme activity, or its sensitivity to substrate-mediated inhibition (49). Nevertheless, analysis of thyroid and skeletal muscle samples from subjects homozygous for the Ala risk allele revealed D2 enzyme velocity was reduced by 50% compared to heterozygous individuals or subjects homozygous for the protective Thr allele (49), although it was subsequently suggested this finding may have resulted from a methodological artifact (50,51). The initial genetic studies suggesting that reduced D2 activity is associated with increased susceptibility to OA (6) were based on the assumption that the Thr92Ala polymorphism results in reduced enzyme activity (49). However, studies showing that transgenic rats overexpressing D2 in articular cartilage had increased susceptibility to OA (14) were not consistent with this conclusion, and recent findings of allelic imbalance and increased expression of D2 protein in OA cartilage (28,29) further suggest that increased D2 activity may in fact be detrimental for articular cartilage maintenance. On the other hand, a large genome-wide association study and a recent meta-analysis both failed to identify DIO2 as a disease-susceptibility locus for OA (26,27), while a recent study in primary human chondrocytes actually suggests an anti-inflammatory role for D2 in cartilage (46). Currently, it remains unclear, therefore, whether variation in DIO2 plays a role in osteoarthritis pathogenesis, although the independent identification of DIO3 as a disease susceptibility locus (5) strongly suggests a role for control of local tissue T3 availability in the regulation of joint homeostasis.
Overall, it is not certain whether the original hypothesis that reduced tissue T3 availability protects joints from damage and the development of OA is likely to be correct or whether the reverse is more probable. Alternatively, it is possible that tissue T3 deficiency or excess are each detrimental for optimal articular cartilage homeostasis, thus accounting for the apparently contradictory literature. Indeed, OA pathophysiology is complex, with examples of increased or decreased signaling in both the Wnt/β-catenin (52,53) and transforming growth factor-β pathways (54,55) having deleterious effects on joint integrity. The finding of increased subchondral bone mineral content in D2KO mice will thus require careful study to determine its pathophysiological significance; that is, whether loss of D2 protects joints from OA or whether it results in increased disease susceptibility. Future studies will require the use of an OA provocation model, such as surgical destabilization of the medial meniscus (56), or study of the consequences of ageing to determine whether deletion of D2 accelerates or delays the onset and progression of predictable joint destruction. Additional studies to investigate the temporo-spatial specific role of D2 in cartilage could also be performed by crossing mice harboring floxed Dio2 alleles (57) with Agc1-CreERT2 mice (58) that conditionally express Cre-recombinase in cartilage following treatment with tamoxifen.
Footnotes
Acknowledgments
The authors thank Professors Valerie Galton (Dartmouth Medical School, Lebanon, NH) and Donald St. Germain (Maine Medical Center Research Institute, Scarborough, ME) for kindly providing samples from D2KO and D1D2KO mice, and Mahrokh Nohadani, Ann Sandison and Philippe Donatien (Imperial College London) for their generous help with histology. J.A.W. received funding for this work form an Arthritis Research UK Clinical Research Fellowship.
Author Disclosure Statement
No competing financial interests exist.
