Mutations in genes encoding cartilage oligomeric matrix protein and matrilin-3 cause a spectrum of chondrodysplasias called multiple epiphyseal dysplasia (MED) and pseudoachondroplasia (PSACH). The majority of these diseases feature classical endoplasmic reticulum (ER) stress and activation of the unfolded protein response (UPR) as a result of misfolding of the mutant protein. However, the importance and the pathological contribution of ER stress in the disease pathogenesis are unknown. The aim of this study was to investigate the generic role of ER stress and the UPR in the pathogenesis of these diseases. A transgenic mouse line (ColIITgcog) was generated using the collagen II promoter to drive expression of an ER stress-inducing protein (Tgcog) in chondrocytes. The skeletal and histological phenotypes of these ColIITgcog mice were characterised. The expression and intracellular retention of Tgcog induced ER stress and activated the UPR as characterised by increased BiP expression, phosphorylation of eIF2α and spliced Xbp1. ColIITgcog mice exhibited decreased long bone growth and decreased chondrocyte proliferation rate. However, there was no disruption of chondrocyte morphology or growth plate architecture and perturbations in apoptosis were not apparent. Our data demonstrate that the targeted induction of ER stress in chondrocytes was sufficient to reduce the rate of bone growth, a key clinical feature associated with MED and PSACH, in the absence of any growth plate dysplasia. This study establishes that classical ER stress is a pathogenic factor that contributes to the disease mechanism of MED and PSACH. However, not all the pathological features of MED and PSACH were recapitulated, suggesting that a combination of intra- and extra-cellular factors are likely to be responsible for the disease pathology as a whole.

Pseudoachondroplasia and multiple epiphyseal dysplasia are genetic skeletal diseases resulting from mutations in cartilage structural proteins. Electron microscopy and immunohistochemistry previously showed that the appearance of the cartilage extracellular matrix (ECM) in targeted mouse models of these diseases is disrupted; however, the precise changes in ECM organization and the pathological consequences remain unknown. Our aim was to determine the effects of matrilin-3 and COMP mutations on the composition and extractability of ECM components to inform how these detrimental changes might influence cartilage organization and degeneration. Cartilage was sequentially extracted using increasing denaturants and the extraction profiles of specific proteins determined using SDS-PAGE/Western blotting. Furthermore, the relative composition of protein pools was determined using mass spectrometry for a non-biased semi-quantitative analysis. Western blotting revealed changes in the extraction of matrilins, COMP and collagen IX in mutant cartilage. Mass spectrometry confirmed quantitative changes in the extraction of structural and non-structural ECM proteins, including proteins with roles in cellular processes such as protein folding and trafficking. In particular, genotype-specific differences in the extraction of collagens XII and XIV and tenascins C and X were identified; interestingly, increased expression of several of these genes has recently been implicated in susceptibility and/or progression of murine osteoarthritis. We demonstrated that mutation of matrilin-3 and COMP caused changes in the extractability of other cartilage proteins and that proteomic analyses of Matn3 V194D, Comp T585M and Comp DelD469 mouse models revealed both common and discrete disease signatures that provide novel insight into skeletal disease mechanisms and cartilage degradation.

To characterize the circadian clock in murine cartilage tissue and identify tissue-specific clock target genes, and to investigate whether the circadian clock changes during aging or during cartilage degeneration using an experimental mouse model of osteoarthritis (OA).

Integrins are a functionally significant family of metazoan cell surface adhesion receptors. The receptors are dimers composed of an alpha and a beta chain. Vertebrate genomes encode an expanded set of integrin alpha and beta chains in comparison with protostomes such as drosophila or the nematode worm. The publication of the genome of a basal chordate, Ciona intestinalis, provides a unique opportunity to gain further insight into how and when the expanded integrin supergene family found in vertebrates evolved.

Extracellular matrix (ECM) is a key metazoan characteristic. In addition to providing structure and orientation to tissues, it is involved in many cellular processes such as adhesion, migration, proliferation and differentiation. Here we provide a comprehensive analysis of ECM molecules focussing on when vertebrate specific matrices evolved. We identify 60 ECM genes and 20 associated processing enzymes in the genome of the urochordate Ciona intestinalis. A comparison with vertebrate and protostome genomes has permitted the identification of both a core set of metazoan matrix genes and vertebrate-specific innovations in the ECM. We have identified a few potential cases of de novo vertebrate ECM gene innovation, but the majority of ECM genes have resulted from duplication of pre-existing genes present in the ancestral vertebrate. In conclusion, the modern complexity we see in vertebrate ECM has come about largely by duplication and modification of pre-existing matrix molecules. Extracellular matrix genes and their processing enzymes appear to be over-represented in the vertebrate genome suggesting that these genes played an active role enabling and underpinning the evolution of vertebrates.

Osteoarthritis has been associated with a plethora of pathological factors and one which has recently emerged is chondrocyte endoplasmic reticulum (ER) stress. ER stress is sensed by key ER-resident stress sensors, one of which is activating transcription factor 6 (ATF6). The purpose of this study is to determine whether increased ER stress plays a role in OA.

The unfolded protein response (UPR) is a conserved cellular response to the accumulation of proteinaceous material in endoplasmic reticulum (ER), active both in health and disease to alleviate cellular stress and improve protein folding. Multiple epiphyseal dysplasia (EDM5) is a genetic skeletal condition and a classic example of an intracellular protein aggregation disease, whereby mutant matrilin-3 forms large insoluble aggregates in the ER lumen, resulting in a specific 'disease signature' of increased expression of chaperones and foldases, and alternative splicing of the UPR effector XBP1. Matrilin-3 is expressed exclusively by chondrocytes thereby making EDM5 a perfect model system to study the role of protein aggregation in disease. In order to dissect the role of XBP1 signalling in aggregation-related conditions we crossed a p.V194D Matn3 knock-in mouse model of EDM5 with a mouse line carrying a cartilage specific deletion of XBP1 and analysed the resulting phenotype. Interestingly, the growth of mice carrying the Matn3 p.V194D mutation compounded with the cartilage specific deletion of XBP1 was severely retarded. Further phenotyping revealed increased intracellular retention of amyloid-like aggregates of mutant matrilin-3 coupled with dramatically decreased cell proliferation and increased apoptosis, suggesting a role of XBP1 signalling in protein accumulation and/or degradation. Transcriptomic analysis of chondrocytes extracted from wild type, EDM5, Xbp1-null and compound mutant lines revealed that the alternative splicing of Xbp1 is crucial in modulating levels of protein aggregation. Moreover, through detailed transcriptomic comparison with a model of metaphyseal chondrodysplasia type Schmid (MCDS), an UPR-related skeletal condition in which XBP1 was removed without overt consequences, we show for the first time that the differentiation-state of cells within the cartilage growth plate influences the UPR resulting from retention of a misfolded mutant protein and postulate that modulation of XBP1 signalling pathway presents a therapeutic target for aggregation related conditions in cells undergoing proliferation.

Metaphyseal chondrodysplasia, Schmid type (MCDS) is characterized by mild short stature and growth plate hypertrophic zone expansion, and caused by collagen X mutations. We recently demonstrated the central importance of ER stress in the pathology of MCDS by recapitulating the disease phenotype by expressing misfolding forms of collagen X (Schmid) or thyroglobulin (Cog) in the hypertrophic zone. Here we characterize the Schmid and Cog ER stress signaling networks by transcriptional profiling of microdissected mutant and wildtype hypertrophic zones. Both models displayed similar unfolded protein responses (UPRs), involving activation of canonical ER stress sensors and upregulation of their downstream targets, including molecular chaperones, foldases, and ER-associated degradation machinery. Also upregulated were the emerging UPR regulators Wfs1 and Syvn1, recently identified UPR components including Armet and Creld2, and genes not previously implicated in ER stress such as Steap1 and Fgf21. Despite upregulation of the Chop/Cebpb pathway, apoptosis was not increased in mutant hypertrophic zones. Ultrastructural analysis of mutant growth plates revealed ER stress and disrupted chondrocyte maturation throughout mutant hypertrophic zones. This disruption was defined by profiling the expression of wildtype growth plate zone gene signatures in the mutant hypertrophic zones. Hypertrophic zone gene upregulation and proliferative zone gene downregulation were both inhibited in Schmid hypertrophic zones, resulting in the persistence of a proliferative chondrocyte-like expression profile in ER-stressed Schmid chondrocytes. Our findings provide a transcriptional map of two chondrocyte UPR gene networks in vivo, and define the consequences of UPR activation for the adaptation, differentiation, and survival of chondrocytes experiencing ER stress during hypertrophy. Thus they provide important insights into ER stress signaling and its impact on cartilage pathophysiology.

Mutant matrilin-3 (V194D) forms non-native disulphide bonded aggregates in the rER of chondrocytes from cell and mouse models of multiple epiphyseal dysplasia (MED). Intracellular retention of mutant matrilin-3 causes endoplasmic reticulum (ER) stress and induces an unfolded protein response (UPR) including the upregulation of two genes recently implicated in ER stress: Armet and Creld2. Nothing is known about the role of Armet and Creld2 in human genetic diseases. In this study, we used a variety of cell and mouse models of chondrodysplasia to determine the genotype-specific expression profiles of Armet and Creld2. We also studied their interactions with various mutant proteins and investigated their potential roles as protein disulphide isomerases (PDIs). Armet and Creld2 were up-regulated in cell and/or mouse models of chondrodysplasias caused by mutations in Matn3 and Col10a1, but not Comp. Intriguingly, both Armet and Creld2 were also secreted into the ECM of these disease models following ER stress. Armet and Creld2 interacted with mutant matrilin-3, but not with COMP, thereby validating the genotype-specific expression. Substrate-trapping experiments confirmed Creld2 processed PDI-like activity, thus identifying a putative functional role. Finally, alanine substitution of the two terminal cysteine residues from the A-domain of V194D matrilin-3 prevented aggregation, promoted mutant protein secretion and reduced the levels of Armet and Creld2 in a cell culture model. We demonstrate that Armet and Creld2 are genotype-specific ER stress response proteins with substrate specificities, and that aggregation of mutant matrilin-3 is a key disease trigger in MED that could be exploited as a potential therapeutic target.

Disease mechanisms leading to different forms of chondrodysplasia include extracellular matrix (ECM) alterations and intracellular stress resulting in abnormal changes to chondrocyte proliferation and survival. Delineating the relative contribution of these two disease mechanisms is a major challenge in understanding disease pathophysiology in genetic skeletal diseases and a prerequisite for developing effective therapies. To determine the influence of intracellular stress and changes in chondrocyte phenotype to the development of chondrodysplasia, we targeted the expression of the G2320R mutant form of thyroglobulin to the endoplasmic reticulum (ER) of resting and proliferating chondrocytes. Previous studies on this mutant protein have shown that it induces intracellular aggregates and causes cell stress and death in the thyroid gland. The expression and retention of this exogenous mutant protein in resting and proliferating chondrocytes resulted in a chronic cell stress response, growth plate dysplasia and reduced bone growth, without inducing any alterations to the architecture and organization of the cartilage ECM. More significantly, the decreased bone growth seemed to be the direct result of reduced chondrocyte proliferation in the proliferative zone of growth plates in transgenic mice, without transcriptional activation of a classical unfolded protein response (UPR) or apoptosis. Overall, these data show that mutant protein retention in the ER of resting and proliferative zone chondrocytes is sufficient to cause disrupted bone growth. The specific disease pathways triggered by mutant protein retention do not necessarily involve a prototypic UPR, but all pathways impact upon chondrocyte proliferation in the cartilage growth plate.

Schmid metaphyseal chondrodysplasia (MCDS) involves dwarfism and growth plate cartilage hypertrophic zone expansion resulting from dominant mutations in the hypertrophic zone collagen, Col10a1. Mouse models phenocopying MCDS through the expression of an exogenous misfolding protein in the endoplasmic reticulum (ER) in hypertrophic chondrocytes have demonstrated the central importance of ER stress in the pathology of MCDS. The resultant unfolded protein response (UPR) in affected chondrocytes involved activation of canonical ER stress sensors, IRE1, ATF6, and PERK with the downstream effect of disrupted chondrocyte differentiation. Here, we investigated the role of the highly conserved IRE1/XBP1 pathway in the pathology of MCDS. Mice with a MCDS collagen X p.N617K knock-in mutation (ColXN617K) were crossed with mice in which Xbp1 was inactivated specifically in cartilage (Xbp1CartΔEx2), generating the compound mutant, C/X. The severity of dwarfism and hypertrophic zone expansion in C/X did not differ significantly from ColXN617K, revealing surprising redundancy for the IRE1/XBP1 UPR pathway in the pathology of MCDS. Transcriptomic analyses of hypertrophic zone cartilage identified differentially expressed gene cohorts in MCDS that are pathologically relevant (XBP1-independent) or pathologically redundant (XBP1-dependent). XBP1-independent gene expression changes included large-scale transcriptional attenuation of genes encoding secreted proteins and disrupted differentiation from proliferative to hypertrophic chondrocytes. Moreover, these changes were consistent with disruption of C/EBP-β, a master regulator of chondrocyte differentiation, by CHOP, a transcription factor downstream of PERK that inhibits C/EBP proteins, and down-regulation of C/EBP-β transcriptional co-factors, GADD45-β and RUNX2. Thus we propose that the pathology of MCDS is underpinned by XBP1 independent UPR-induced dysregulation of C/EBP-β-mediated chondrocyte differentiation. Our data suggest that modulation of C/EBP-β activity in MCDS chondrocytes may offer therapeutic opportunities.

The small GTPase RhoA is a major regulator of actin reorganization during the formation of stress fibers; thus identifying molecules that regulate Rho activity is necessary for a complete understanding of the mechanisms that determine cell contractility. Here, we have identified Arhgap28 as a Rho GTPase activating protein (RhoGAP) that switches RhoA to its inactive form. We generated an Arhgap28-LacZ reporter mouse that revealed gene expression in soft tissues at E12.5, pre-bone structures of the limb at E15.5, and prominent expression restricted mostly to ribs and limb long bones at E18.5 days of development. Expression of recombinant Arhgap28-V5 in human osteosarcoma SaOS-2 cells caused a reduction in the basal level of RhoA activation and disruption of actin stress fibers. Extracellular matrix assembly studies using a 3-dimensional cell culture system showed that Arhgap28 was upregulated during Rho-dependent assembly of the ECM. Taken together, these observations led to the hypothesis that an Arhgap28 knockout mouse model would show a connective tissue phenotype, perhaps affecting bone. Arhgap28-null mice were viable and appeared normal, suggesting that there could be compensation from other RhoGAPs. Indeed, we showed that expression of Arhgap6 (a closely related RhoGAP) was upregulated in Arhgap28-null bone tissue. An upregulation in RhoA expression was also detected suggesting that Arhgap28 may be able to additionally regulate Rho signaling at a transcriptional level. Microarray analyses revealed that Col2a1, Col9a1, Matn3, and Comp that encode extracellular matrix proteins were downregulated in Arhgap28-null bone. Although mutations in these genes cause bone dysplasias no bone phenotype was detected in the Arhgap-28 null mice. Together, these data suggest that the regulation of Rho by RhoGAPs, including Arhgap28, during the assembly and development of mechanically strong tissues is complex and may involve multiple RhoGAPs.

We describe a new method, PhenomeExpress, for the analysis of transcriptomic datasets to identify pathogenic disease mechanisms. Our analysis method includes input from both protein-protein interaction and phenotype similarity networks. This introduces valuable information from disease relevant phenotypes, which aids the identification of sub-networks that are significantly enriched in differentially expressed genes and are related to the disease relevant phenotypes. This contrasts with many active sub-network detection methods, which rely solely on protein-protein interaction networks derived from compounded data of many unrelated biological conditions and which are therefore not specific to the context of the experiment. PhenomeExpress thus exploits readily available animal model and human disease phenotype information. It combines this prior evidence of disease phenotypes with the experimentally derived disease data sets to provide a more targeted analysis. Two case studies, in subchondral bone in osteoarthritis and in Pax5 in acute lymphoblastic leukaemia, demonstrate that PhenomeExpress identifies core disease pathways in both mouse and human disease expression datasets derived from different technologies. We also validate the approach by comparison to state-of-the-art active sub-network detection methods, which reveals how it may enhance the detection of molecular phenotypes and provide a more detailed context to those previously identified as possible candidates.

ADAMTS, constituting a recently discovered family of secreted zinc-dependent metalloproteases, have been shown to have critical physiological roles through identification of a number of natural animal and human gene mutations. The identification of six ADAMTS genes in the basal chordate Ciona intestinalis provides new insight into how, when and in what order the vertebrate orthologues have evolved. The phylogenetic assignments, based on sequences conserved across all genes, are supported by conserved domain structures within defined sub-families. The phylogeny and the frequent localisation of ADAMTS genes in paralogous regions of the genome are consistent with the vertebrate lineages having arisen by large scale or genome duplication. The high level of conservation in the protease active site of vertebrate orthologues within some sub-families suggests subfunctionalisation, whereas the greater divergence in others would favour the evolution of novel substrate specificities and these observations are borne-out where substrate-specificity is known. The expansion and sub-specialization of the ADAMTS family is a component of the increased complexity of extracellular matrix that is associated with the evolution of vertebrates.

Pseudoachondroplasia (PSACH) results from mutations in cartilage oligomeric matrix protein (COMP) and the p.D469del mutation within the type III repeats of COMP accounts for approximately 30% of PSACH. To determine disease mechanisms of PSACH in vivo, we introduced the Comp D469del mutation into the mouse genome. Mutant animals were normal at birth but grew slower than their wild-type littermates and developed short-limb dwarfism. In the growth plates of mutant mice chondrocyte columns were reduced in number and poorly organized, while mutant COMP was retained within the endoplasmic reticulum (ER) of cells. Chondrocyte proliferation was reduced and apoptosis was both increased and spatially dysregulated. Previous studies on COMP mutations have shown mutant COMP is co-localized with chaperone proteins, and we have reported an unfolded protein response (UPR) in mouse models of PSACH-MED (multiple epiphyseal dysplasia) harboring mutations in Comp (T585M) and Matn3, Comp etc (V194D). However, we found no evidence of UPR in this mouse model of PSACH. In contrast, microarray analysis identified expression changes in groups of genes implicated in oxidative stress, cell cycle regulation, and apoptosis, which is consistent with the chondrocyte pathology. Overall, these data suggest that a novel form of chondrocyte stress triggered by the expression of mutant COMP is central to the pathogenesis of PSACH.

Mutations causing metaphyseal chondrodysplasia type Schmid (MCDS) (e.g., Col10a1p.N617K) induce the pathology by a mechanism involving increased endoplasmic reticulum (ER) stress triggering an unfolded protein response (UPR) in hypertrophic chondrocytes (Rajpar et al. 2009). Here we correlate the expression of mutant protein with the onset of the UPR and disease pathology (hypertrophic zone [HZ] expansion) in MCDS and ColXTg(cog) mouse lines from E14.5 to E17.5. Embryos homozygous for the Col10a1p.N617K mutation displayed a delayed secretion of mutant collagen X accompanied by a UPR at E14.5, delayed ossification of the primary center at E15.5, and an expanded HZ at E17.5. Heterozygote embryos expressed mutant collagen X from E14.5 but exhibited no evidence of a UPR or an HZ expansion until after E17.5. Embryos positive for the ER stress-inducing ColXTg(cog) allele expressed Tg(cog) at E14.5, but the onset of the UPR was not apparent until E15.5 in homozygous and E17.5 in hemizygous embryos. Only homozygous embryos exhibited an HZ expansion at E17.5. The differential onset of the UPR and pathology, dependent on mutation type and gene dosage, indicates that hypertrophic chondrocytes have a latent capacity to deal with ER stress, which must be exceeded to trigger the UPR and HZ expansion.

The metzincins are a large gene superfamily of proteases characterized by the presence of a zinc protease domain, and include the ADAM, ADAMTS, BMP1/TLL, meprin and MMP genes. Metzincins are involved in the proteolysis of a wide variety of proteins, including those of the extracellular matrix. The metzincin gene superfamily comprises eighty proteins in the human genome and ninety-three in the mouse. When and how the level of complexity apparent in the vertebrate metzincin gene superfamily arose has not been determined in detail. Here we present a comprehensive analysis of vertebrate metzincins using genes from both Ciona intestinalis and Danio rerio to provide new insights into the complex evolution of this gene superfamily.

Mutations, mostly in the region of the COL10A1 gene encoding the C-terminal non-collagenous domain, cause the dwarfism metaphyseal chondrodysplasia type Schmid (MCDS). In most cases, the disease mechanism involves the misfolding of the mutant protein causing increased endoplasmic reticulum (ER) stress and an unfolded protein response (UPR). However, in an iliac crest biopsy, the COL10A1 p.Y632X mutation was found to produce instability of the mutant mRNA such that little mutant protein may be produced. To investigate the disease mechanism further, a gene-targeted mouse model of the Col10a1 p.Y632X mutation was generated. In this model, the mutant mRNA showed no instability, and in mice heterozygous for the mutation, mutant and wild-type mRNAs were present at equal concentrations. The protein was translated from the mutant allele and retained within the cell, triggering increased ER stress and a UPR. The mutation produced a relatively severe form of MCDS. Nevertheless, treatment of the mice with carbamazepine (CBZ), a drug which stimulates intracellular proteolysis and alleviates ER stress, effectively reduced the disease severity in this model of MCDS caused by a premature stop codon in the Col10a1 gene. Specifically, the drug reduced ER stress in the growth plate, restored growth plate architecture toward the wild-type state, significantly increased bone growth and within 2 weeks of treatment corrected the MCDS-induced hip distortion. These results indicate that CBZ is likely to be effective in ongoing clinical trials against all forms of MCDS whether caused by premature stop codons or substitutions.

Pseudoachondroplasia (PSACH) is an autosomal dominant skeletal dysplasia caused by mutations in cartilage oligomeric matrix protein (COMP) and characterised by short limbed dwarfism and early onset osteoarthritis. Mouse models of PSACH show variable retention of mutant COMP in the ER of chondrocytes, however, in each case a different stress pathway is activated and the underlying disease mechanisms remain largely unknown. T585M COMP mutant mice are a model of moderate PSACH and demonstrate a mild ER stress response. Although mutant COMP is not retained in significant quantities within the ER of chondrocytes, both BiP and the pro-apoptotic ER stress-related transcription factor CHOP are mildly elevated, whilst bcl-2 levels are decreased, resulting in increased and spatially dysregulated chondrocyte apoptosis. To determine whether the abnormal chondrocyte apoptosis observed in the growth plate of mutant mice is CHOP-mediated, we bred T585M COMP mutant mice with CHOP-null mice to homozygosity, and analysed the resulting phenotype. Although abnormal apoptosis was alleviated in the resting zone following CHOP deletion, the mutant growth plates were generally more disorganised. Furthermore, the bone lengths of COMP mutant CHOP null mice were significantly shorter at 9 weeks of age when compared to the COMP mutant mice, including a significant difference in the skull length. Overall, these data demonstrate that CHOP-mediated apoptosis is an early event in the pathobiology of PSACH and suggest that the lack of CHOP, in conjunction with a COMP mutation, may lead to aggravation of the skeletal phenotype via a potentially synergistic effect on endochondral ossification.

Integrins comprise a large family of alpha,beta heterodimeric, transmembrane cell adhesion receptors that mediate diverse essential biological functions. Higher vertebrates possess a single beta1 gene, and the beta1 subunit associates with a large number of alpha subunits to form the major class of extracellular matrix (ECM) receptors. Despite the fact that the zebrafish (Danio rerio) is a rapidly emerging model organism of choice for developmental biology and for models of human disease, little is currently known about beta1 integrin sequences and functions in this organism.