The effects of physiological and injurious hydrostatic pressure on murine ex vivo articular and growth plate cartilage explants: an RNAseq study

Abstract

Introduction:

Chondrocytes are continuously exposed to loads placed upon them. Physiological loads are pivotal to the maintenance of articular cartilage health, while abnormal loads contribute to pathological joint degradation. Similarly, the growth plate cartilage is subject to various loads during growth and development. Due to the high-water content of cartilage, hydrostatic pressure is considered one of the main biomechanical influencers on chondrocytes and has been shown to play an important role in the mechano-regulation of cartilage.

Methods:

Herein, we conducted RNAseq analysis of ex vivo hip cap (articular), and metatarsal (growth plate) cartilage cultures subjected to physiological (5 MPa) and injurious (50 MPa) hydrostatic pressure, using the Illumina platform (n = 4 replicates).

Results:

Several hundreds of genes were shown to be differentially modulated by hydrostatic pressure, with the majority of these changes evidenced in hip cap cartilage cultures (375 significantly upregulated and 322 downregulated in 5 MPa versus control; 1022 upregulated and 724 downregulated in 50 MPa versus control). Conversely, fewer genes were differentially affected by hydrostatic pressure in the metatarsal cultures (5 significantly upregulated and 23 downregulated in 5 MPa versus control; 7 significantly upregulated and 19 downregulated in 50 MPa versus control). Using Gene Ontology annotations for Biological Processes, in the hip cap data we identified a number of pathways that were modulated by both physiological and injurious hydrostatic pressure. Pathways upregulated in response to 50 MPa versus control, included those involved in the generation of precursor metabolites and cellular respiration. Biological processes that were downregulated in this tissue included ossification, connective tissue development, and chondrocyte differentiation.

Discussion:

Collectively our data highlights the divergent chondrocyte phenotypes in articular and growth plate cartilage. Further, we show that the magnitude of hydrostatic pressure application has distinct effects on gene expression and biological processes in hip cap cartilage explants. Finally, we identified differential expression of a number of genes that have previously been identified as osteoarthritis risk genes, including Ctsk, and Chadl. Together these data may provide potential genetic targets for future investigations in osteoarthritis research and novel therapeutics.

1 Introduction

Articular cartilage is a specialized connective tissue that covers the ends of bones in synovial joints and facilitates joint movement. It is load bearing and therefore protects underlying subchondral bone from excessive forces. The articular cartilage consists of chondrocytes which retain a stable phenotype to ensure the longevity of the tissue (1, 2). This is in contrast to the chondrocytes of the growth plate cartilage which undergo defined stages of maturation and differentiation to enable longitudinal bone growth (3).

Structurally, the articular cartilage can be divided into superficial, intermediate, and deep zones which are distinct in their organization of both the chondrocytes, surrounded by their individual pericellular matrix, and the collagen type-II and aggrecan-rich matrix (3). The articular cartilage functions to withstand physiological loading over the life-course. However, in the degenerative joint disease osteoarthritis, pathology is characterized by progressive articular cartilage degradation (4). Whilst osteoarthritis is well established to affect all tissues of the joint, the cellular and molecular mechanisms are incompletely understood (5–7). Various forms of mechanical stimuli are involved in the maintenance of the articular cartilage and thus the mechanoresponse of the chondrocyte plays an important role in the development of osteoarthritis (8–10). Compression, tensile and shear stress result in deformative loading, whereas osmotic and hydrostatic pressure induce stress without tissue or cellular deformation (8, 9, 11, 12). As a highly hydrated tissue, interstitial fluid pressurization within the articular cartilage is considered one of the main biomechanical influencers on chondrocytes (13–15). Throughout the cartilage zones, chondrocytes are subjected and respond to a hydrostatic pressure gradient, ranging from 0.1-10 MPa, to direct matrix remodeling, chondrogenesis and chondrocyte metabolism (13, 14). However, excessive hydrostatic pressure (≥20 MPa) outside the physiological range has been shown to induce apoptosis, alter cell morphology and metabolism, reduce extracellular matrix (ECM) synthesis, induce inflammatory cytokine production, and modulate oxidative stress (16–19).

In vitro, hydrostatic pressure can be applied experimentally to cells and tissues derived from both animals and humans to investigate mechanotransduction, for example in monolayer cultures (20–23), micromass or pellet cultures (24, 25), 3D cell scaffolds (26–29), and explant cultures (17, 22, 30, 31). The ability to provide either dynamic or continuous hydrostatic pressure, alter the magnitude and/or the duration of pressure provides an alternative approach to study the effects of mechanical stimulation (13, 32). Whilst there is little consensus within the field on the duration and pressure magnitudes in cultures, our previous meta-analysis has indicated that in human and animal- derived cells, low pressure (5 MPa) leads to anabolic responses, including elevated aggrecan expression and proteoglycan release, whereas a higher pressure (50 MPa) has a negative effect on proteoglycan production (33). Therefore, it is possible to investigate the effects of hydrostatic pressure at both physiological and pathophysiological levels.

To determine the effects of hydrostatic pressure on the molecular pathways involved in the regulation of chondrocyte physiology, transcriptomic analyses are often employed to identify responsive genes. Several studies in animal cells have utilized these approaches in the study of chondrocyte progenitor cells, immortalized chondrocytes, and primary chondrocytes within a hydrogel; however, transcriptome sequencing on ex vivo models has not yet been performed (21, 29, 34). Phenotypic changes are often observed in cells cultured in a monolayer, with cells de-differentiating or altering morphology, whereas ex vivo models allow examination of cells within their native environment (35). Herein, the aim of this study was to perform RNAseq analysis on two murine ex vivo cartilage models (hip cap and metatarsal) after exposure to physiological and injurious hydrostatic pressure, to examine the effects of hydrostatic pressure on gene expression in two different chondrocyte phenotypes.

2 Methods

2.1 Isolation and culture of ex vivo cartilage models

All mice utilized in these studies were kept in controlled conditions at the University of Brighton and all tissue isolation procedures were performed in accordance with the UK Animals (Scientific Procedures) Act of 1986 and regulations set by the UK Home Office and local institutional guidelines (PPL: PP3310437). Analyses were conducted blindly where possible to minimize the effects of subjective bias. Animal studies were conducted in line with the ARRIVE guidelines.

Femoral heads were isolated from 4-week-old male C57/BL6J mice (Charles River), as previously described (Figure 1) (36). In brief, the hip joint was dislocated by applying slight pressure at the joint, and the femoral cap was avulsed using forceps. At this developmental stage, the predominant component of this tissue is the articular cartilage, therefore underlying subchondral bone was not included. Both hip caps were pooled from each individual mouse (n=4 mice/experimental group). Hip caps were cultured in Dulbecco’s Modified Eagle Medium with GlutaMAX, substituted with 100 U/ml penicillin, 100µg/ml streptomycin (Thermo Fisher Scientific) in a humidified atmosphere (37°C, 5% CO₂).

Figure 1

Embryonic metatarsal organ cultures provide a well‐established model of endochondral bone growth (Figure 1) (37). Metatarsals were isolated from E15 embryos of C57/BL6J (Charles River) mice as previously described (36). Six metatarsal bones were pooled per sample (n=4 samples/experimental group). Metatarsal bones were cultured in α-Minimum Essential Medium supplemented with 0.2% BSA Fraction V; 1 mmol/l β-glycerophosphate (βGP); 0.05 mg/ml L-ascorbic acid phosphate; 0.05 mg/ml gentamicin and 1.25 μg/ml fungizone (Thermo Fisher Scientific) in a humidified atmosphere (37°C, 5% CO₂).

2.2 Application of hydrostatic pressure

After 24 hours of culture, hips caps and metatarsals were placed into 5 ml sterile plastic syringes fitted with Luer lock end caps, taking care to eliminate all air bubbles (Suppl. Figure 1). Movement of the syringe plunger allowed for equilibration of pressure between syringe contents and the pressure vessel water (17). The syringes were placed in a water-filled pressure vessel at room temperature. Syringes were pressurized to 0 MPa (control), 5 MPa (physiological) or 50 MPa (injurious) hydrostatic pressure for 1 hour (Figure 1; Suppl. Figure 1). Following exposure to hydrostatic pressure, tissues were placed back into the incubator and cultured for a further 24 hours in the respective media, then flash frozen at -80°C until RNA extraction.

2.3 RNA extraction and sequencing

Tissue (<100 mg) were defrosted on ice and 1 ml Trizol (Qiagen) was added to each sample; tissues were homogenised using a mechanical disruptor, making sure to keep them cool by putting on ice every 15 seconds. Samples were incubated at room temperature for a minimum of 10 minutes to allow for cell lysis and centrifuged at 12,000 x g for 15 minutes at 4°C to pellet the excess tissue, whilst retaining RNA in solution. The supernatant was transferred to a clean tube and 200 μL of chloroform (Sigma) added. After vigorous shaking for 20 seconds, the samples were incubated at room temperature for 3 minutes and then centrifuged at 12,000 x g for 15 minutes at 4°C to enable phase separation. The upper, aqueous phase was transferred to a new tube, avoiding the interface. Following the addition of an equal volume of 70% ethanol, the samples were mixed thoroughly by vortexing and total RNA purified using RNeasy Mini spin columns (Qiagen), according to the manufacturer’s recommendations. Purified RNA was eluted in 30 μl of RNase-free water, repeating the elution twice by reapplying the elute. The concentration and purity of the RNA samples were assessed using a Nanodrop One C spectrophotometer (Labtech) and the quality of the RNA was assessed on a TapeStation 4200 (Agilent Technologies).

All samples passed purity quality control checks but exhibited RNA Integrity Number (RIN) equivalent values below the ideal minimum of 7 (average value 2.8). The low RIN values obtained are considered typical for these explant tissue samples and suggest some partial degradation of the total RNA. DV200 analysis using the Agilent TapeStation 4200 software showed a percentage of fragments between 200 and 10000 bp ranging between 53.51-87.4% in all RNA samples. Sequencing libraries were prepared using the Universal Plus™ Total RNASeq with NuQuant kit and a mouse rRNA depletion module (Tecan Genomics), required for partially degraded RNA samples. Library construction strategy was pair end and strand specific. Libraries were checked for quality using the TapeStation 4200, quantified, normalized and sequenced on the Illumina NextSeq500 sequencer using a high-output kit (17 libraries) and a mid-output kit (7 libraries).

2.4 Data analysis

Initial sequencing read quality control was conducted using fastqc (version 0.11.9) (38) and multiqc (version 1.8) (39). Trimming was performed using TrimGalore using a minimum quality threshold of 20, discarding any trimmed reads shorter than 20 nucleotides. Trimmed reads were quantified using kallisto quant and transcript quantifications were converted to gene level by tximport. The transcriptome mapping data for all samples was imported into R for data summarization at the gene level. The data was normalized and analyzed using the DESeq2 pipeline (40). Unsupervised clustering of the sample data was performed using the R packages pheatmap and pcaMethods. Significant genes were identified by analysis using a model design that considered the sequencing run and strandedness of the library as possible batch effects (design= ~SeqRun + Library + Condition) and applying a 5% significance threshold to p-values adjusted using the Benjamini and Hochberg procedure (a significance threshold referred to elsewhere in the text as padj<=0.05, or 5% FDR). For functional analysis of the groups of differentially expressed genes, clusterProfiler was utilized to identify significantly over-represented functional categories using a significance threshold of 5% on the Benjamini and Hochberg corrected p-values (41). Annotations for the Gene Ontology (GO) Biological Process (BP), from the R package org.Mm.eg.db (version 3.11.4) were used (42). Genes that were significantly differentially expressed between our samples were compared to recent genome-wide association studies of osteoarthritis that have identified a number of osteoarthritis risk genes (43, 44).

3 Results

Herein, we conducted RNAseq analysis of murine ex vivo hip cap (articular), and metatarsal (growth plate) cartilage cultures (n=4 replicates) subjected to physiological (5 MPa) and injurious (50 MPa) hydrostatic pressure. Unsupervised clustering of the gene expression data indicated a clear distinction between the hip cap and metatarsal sample data, but two of the hip cap cartilage samples (H502 [exposed to 50 MPa hydrostatic pressure] and HC4 [control, 0 MPa hydrostatic pressure]) appeared to be outliers, thus were excluded from all downstream statistical analyses (Figures 2A, B; Suppl. Figures 2, 3).

Figure 2

3.1 Gene expression profiles of articular and growth plate cartilage

Prior to differential gene expression analyses focusing on the effects of hydrostatic pressure, the gene expression profiles of the two different cartilage explants were investigated to assess the genes and pathways that may be differentially expressed between a transient (growth plate) and an inherently stable (articular) cartilage phenotype (Suppl. Data 1). There were 2775 genes upregulated and 3368 genes downregulated in hip cap cartilage in comparison to metatarsal cartilage (Figure 3A). Upregulated genes with the greatest log₂ fold change included ribosomal protein L9 (Rpl9-ps4, 39.4-fold), collagen type X (Col10a1, 7.7-fold), and frizzled-related protein (Frzb, 7.2-fold) (Table 1). Downregulated genes with the greatest log₂ fold change included microfibrillar-associated protein 4 (Mfap4, 8.6-fold), insulin-like growth factor binding protein 2 (Igfbp2, 7.3-fold) and fibroblastic growth factor 10 (Fgf10, 7.2-fold) (Table 1).

Figure 3

Next, we sought to examine whether these differentially expressed genes were enriched in particular biological processes. Using annotations for GO BP, the data revealed a number of significantly enriched processes, which include ossification (GO:0001503; 124 genes), bone development (GO:0060348; 81 genes), cartilage development (GO:0051216; 83 genes), connective tissue development (GO:0061448; 98 genes), and extracellular matrix organization (GO:0030198; 84 genes), in hip cap cultures in comparison to metatarsals (Suppl. Table 1). Conversely, those that were downregulated included muscle tissue development (GO:0060537; 140 genes) and muscle cell differentiation (GO:0042692; 127 genes), as well as synapse organization (GO:0050808; 142 genes) (Suppl. Table 1).

When comparing the two datasets, the hip cap data yielded many more significant changes than the metatarsal data, and the greater spread of log₂ fold changes taking place in the hip cap samples suggests that the hip cap cartilage explants are more responsive to changes in pressure than the metatarsal explants (Figure 3). Therefore, subsequent analyses focused on the data from the hip cap explants, with the highest up- and down-regulated genes, either commonly or uniquely expressed between each group, in the metatarsal data sets detailed in Suppl. Tables 2, 3.

3.2 Effects of physiological and injurious hydrostatic pressure on gene expression in hip cap cartilage explants

Compared to control, there were 375 genes significantly upregulated with 5 MPa hydrostatic pressure and 322 significantly downregulated in hip cap cultures (Figure 3A; Suppl. Data 1). With injurious hydrostatic pressure (50 MPa), there were 1022 significantly upregulated and 724 significantly downregulated genes (Figure 3A; Suppl. Data 1). Whilst some of the genes were consistently up- or down-regulated across the two hydrostatic pressures in both the hip cap and the metatarsal datasets, most of these significant genes were uniquely expressed by the hip cap datasets at 5 or 50 MPa (Figure 3B). Genes commonly modulated by hydrostatic pressure in these cultures with the greatest log₂ fold changes are detailed in Suppl. Table 4. The genes uniquely expressed in response to hydrostatic pressure magnitudes include Car2 (upregulated in 5 MPa versus control, 1.5-fold), Mlip (upregulated in 50 MPa vs control, 2.6-fold), Tg (downregulated in 5 MPa versus control, 2.2-fold) and Ryr3 (downregulated in 50 MPa versus control, 2.6-fold) (Table 2).

3.3 GO BP enrichment analysis of differentially expressed genes

Using annotations for GO BP, the data revealed significantly enriched processes including regulation of cytokine production (GO:001819; 21 genes), Ras protein signal transduction (GO:007265; 20 genes) and ATP metabolic processes (GO0046034; 18 genes) with 5 MPa hydrostatic pressure application (Table 3; Suppl. Data 2). Conversely, process including cellular component disassembly (GO:0022411; 16 genes) and nuclear transport (GO0051169; 15 genes) were downregulated (Table 3; Suppl. Data 2). With injurious hydrostatic pressure (50 MPa), enriched pathways included generation of precursor metabolites and energy (GO:0006091; 39 genes), and cellular respiration (GO:0045333, 31 genes) (Figure 4; Table 3; Suppl. Data 2). Other upregulated GO BP relevant to the known functions of chondrocytes included regulation of developmental growth (GO0048638; 33 genes), and regulation of cell size (GO0008361; 21 genes) (Suppl. Data 2). Whereas those downregulated included ossification (GO:0001503; 25 genes), cartilage development (GO:0051216; 18 genes), connective tissue development (GO:0061448; 21 genes), and chondrocyte differentiation (GO:0002062; 17 genes) (Table 3; Suppl. Data 2). Further analysis of these enriched pathways in injurious hydrostatic pressure highlighted differential expression of several genes known to be involved in osteoarthritis, such as Fgf2, Ep300, Ngf, Adam9, Igfbp3, Sox9, Comp, Col6a1, Col6a2 and Col11a1.

Figure 4

3.4 Differential expression of previously identified osteoarthritis risk genes

Recent genome-wide association studies of osteoarthritis have identified a number of osteoarthritis risk genes (43, 44). We therefore sought to compare whether these genes were differentially expressed in response to hydrostatic pressure in our datasets (Table 4). Only one of these genes was differentially expressed in our 5 MPa versus control datasets (Wscd2, 0.6-fold downregulation; data not shown). However, with injurious (50 MPa) hydrostatic pressure application, there were 12 genes differentially expressed (Table 4). These included cathepsin K (Ctsk, 0.9-fold upregulation), and chondroadherin-like (Chadl, 0.9-fold downregulation) (Table 4; Suppl. Figure 4).

4 Discussion

In this study we conducted RNAseq analysis of two different ex vivo cartilage explants (metatarsal and hip cap), to examine the effects of two magnitudes of hydrostatic pressure on gene expression. We observed clear differences between the cartilage types, including the upregulation of key genes such as Frzb and Col10a1 in the hip cap explants. Extensive changes in gene expression were observed with hydrostatic pressure in the hip cap cartilage groups, however this was to a weaker extent in the metatarsal explants. Within the hip cap data set, enriched GO BP in the genes that were significantly downregulated in response to injurious hydrostatic pressure (50 MPa) versus control, included those involved in cartilage, bone and connective tissue development. Interestingly, these pathways were also increased when comparing the hip cap to the metatarsal data, suggesting that injurious hydrostatic pressure may promote a more transient-like phenotype in the hip cap cultures. This is further supported by our observed enrichment of the GO BPs for developmental growth and cell size in hip caps exposed to 50 MPa hydrostatic pressure. Indeed, it is well established that in osteoarthritis, the inherently stable articular cartilage undergoes changes that reflect a more developmental cartilage phenotype, such as that in the growth plate (3, 7). Therefore, lessons can be learnt from a better understanding of these two phenotypes, and their similarities and differences in our pursuit of maintaining articular cartilage health in ageing. This is of particular importance given the lack of regenerative capability of the articular cartilage, thus meaning therapies for osteoarthritis remain limited (4, 7).

Articular cartilage covers the ends of the bones in synovial joints, and the chondrocytes within maintain a stable phenotype to ensure joint health and longevity. This is in contrast to the growth plate cartilage, which is more transient in nature, with chondrocytes undergoing differentiation processes which drive endochondral ossification and longitudinal bone growth (3). The chondrocytes of these two cartilaginous structures express different programs, further defined by our RNAseq analysis in hip cap (articular) and metatarsal (growth plate) cartilage. Amongst the most differentially expressed genes in our studies were Frzb, and Col10a1 (both upregulated) and Igfbp2, and Fgf10 (both downregulated). Col10a1 is a key determinant of chondrocyte hypertrophy, with mutant or abnormal human Col10a1 expression associated with abnormalities in this process (45–47). The increase in Col10a1 in our hip cap explants therefore suggests a greater degree of hypertrophy than in our metatarsal explants. Abnormal Col10a1 expression is a well-established feature in osteoarthritis (48–50). Similarly, two SNPs in Frzb, an antagonist of the canonical WNT pathway, have been associated with osteoarthritis (51–53). Further, in pre-clinical models, osteoarthritis severity scores are significantly higher in the joints with deletion of Frzb compared to littermates (54). Together, our data are consistent with previous studies considering the different phenotypes of these cells, thus suggesting diverging phenotypes of these cell populations (55–57).

The high-water content of cartilage (approx. 70-80% water per wet mass) is maintained by an abundance of proteoglycans in the matrix. Chondrocytes in both the growth plate and the articular cartilage are subjected to a number of mechanical forces, including compressive and shear stresses, during loading (9, 13). These mechanical signals then modulate biochemical activity and changes in chondrocyte behavior (22). The majority of research to date has focused on understanding compressive forces on the health of the articular cartilage, however most of this force transforms to hydrostatic pressure due to the interstitial fluid content of joints (14, 58). As such, it can be assumed that hydrostatic pressure is the more prevalent stress to which chondrocytes are exposed. Chondrocytes demonstrate an improved cartilaginous physiology when exposed to hydrostatic pressure, as indicated by their increased ECM production (13). This therefore suggests that understanding the complexities of hydrostatic pressure could be a potential avenue for tissue regeneration in osteoarthritis.

Despite the application of hydrostatic pressure being experimentally controllable, studies have varied in their magnitude, style and duration of hydrostatic pressure application. Our previous meta-analysis informed these factors in the experimental set up for our RNAseq study herein (33). In articular cartilage during normal movement, typical hydrostatic pressure loading of 0.5–10 MPa have been measured (13, 59). Our meta-analysis in 3D cultured chondrocytes confirmed that, based on aggrecan gene expression data, 4–5 MPa can significantly enhance proteoglycan production (33). Conversely, our meta-analysis detailed that the hydrostatic pressure magnitude of 50 MPa had a negative effect on proteoglycans (33). As such, we deemed the magnitudes of physiological (5 MPa) and injurious (50 MPa) hydrostatic pressure to be applicable in our pursuit of understanding gene changes in our explants.

In an RNAseq study performed on monolayer cultures, Zhu et al. used human articular chondrocytes to compare hydrostatic pressure (0.1 MPa) and perfusion methods on the chondrocyte phenotype, with the aim of understanding methods for reducing chondrocyte dedifferentiation in culture (60). Their RNAseq analysis revealed upregulation of well-known chondrocyte genes with hydrostatic pressure and conclude that a low hydrostatic pressure can be beneficial to chondrocytes (60). Further, a previous microarray study examined the effects of continuous hydrostatic pressure (25 MPa) on the chondrogenic ATDC5 cell line, again cultured in monolayer (21). Similarities can be observed between the genes they observe to be modulated by hydrostatic pressure and ours described herein, including differential expression of apoptosis-related and cartilage matrix genes (21). However, Montagne et al. applied a continuous hydrostatic pressure for 24 hours, which is in comparison to our study whereby we applied a single load for 1 hour and is akin to a single injurious event. Further, our examination of two different magnitudes of hydrostatic pressure and in physiologically-relevant cartilage explants adds further strength to our study. In addition, several genes known to play a key role in progression of osteoarthritis (e.g., Fgf2, Ep300, Ngf, Adam9, Igfbp3, Sox9, Comp, Col6a1, Col6a2 and Col11a1) were modulated in our injurious hydrostatic pressure hip cap datasets, thereby validating this approach.

Overall, our results seem to indicate osteoarthritic-like effects of injurious hydrostatic pressure on our hip cap cartilage explants. Among the modulated genes identified in our study, several genes which have been identified as osteoarthritis risk genes from recent GWAS studies were differentially expressed, however verification of these by in situ hybridization or RT-qPCR would be beneficial (43, 44). There was only one gene (Wscd2, WSC Domain-Containing Protein 2) modulated in the 5 MPa versus control dataset, with the majority being in the 50 MPa comparison. Interestingly, Wscd2 has previously been identified as an osteocyte transcriptome signature gene and downregulated in murine bone with ageing, although its role in cartilage has, to our knowledge, not yet fully been defined (61, 62).

Of these risk genes modulated by 50 MPa hydrostatic pressure, the gene that underwent the highest fold upregulation was cathepsin K (Ctsk), a protein expressed by osteoclasts used for collagen degradation (63). This finding is consistent with the previous microarray study by Montague et al. in which Ctsk was found to be strongly induced following the exposure of hydrostatic pressure for 4 hours (21). Indeed, Ctsk has been shown to be overexpressed in the articular cartilage and subchondral bone in osteoarthritis (64, 65). Further, Ctsk deletion in a murine surgical osteoarthritis model (destabilization of the medial meniscus) protected against disease progression (66), as did pharmacological treatment with a cathepsin K inhibitor (SB-553484) in a canine model (67). Pre-clinical findings have been translated to clinical trials with the selective cathepsin K inhibitor MIV-711 reducing bone and cartilage disease progression in individuals with symptomatic, radiographic knee osteoarthritis (68).

Chadl, which encodes for chondroadherin-like protein, plays a role in collagen binding and in the negative regulation of chondrocyte (69). In our studies, its expression underwent the highest fold downregulation with 50 MPa hydrostatic pressure. This is consistent with a previous RNAseq study which examined the subchondral bone of patients who underwent total joint replacement due to osteoarthritis (70). In this study both Chadl and Il11, also identified in our studies, were identified as the most consistently differentially expressed genes and thus have the potential to be targeted for clinical therapies.

Whilst several ion channels known to be involved in chondrocyte mechanotransduction (e.g., Piezo1, Trpv4, Trpv5) (9) were unchanged in our datasets, upregulation of Piezo2 and downregulation of Trpm4 was observed in hip caps exposed to both magnitudes of hydrostatic pressure (Suppl. Data 1). Interestingly, reliable detection of Piezo2 transcripts in primary murine chondrocytes appears to be conflicting in the literature (71, 72). Trpm4 has been identified in cartilage samples from osteoarthritic patients (73), however its role in cartilage mechanotransduction is unclear. Downregulation of Trpm5 and P2rx7 was only observed in hip caps exposed to 5 MPa compared to control (Suppl. Data 1). This suggests that whilst our ex vivo models are sensitive to some changes in ion channel expression with hydrostatic pressure, other mechanisms may exist.

Our study is unique in using two different cartilage explants, both of which offer a physiological model system. We have also applied hydrostatic pressure at magnitudes based on findings from our previous meta-analysis to ensure these are representative of both physiological and injurious load (33). However, we do recognize the limitation in our sample size presented herein. Therefore, the biological interpretation of our findings should be considered appropriately, with the need for a more detailed consideration of the differences observed. For example, it would be pertinent to use a temporal approach to the application of hydrostatic pressure as in this study we applied a single load for 1 hour and is akin to a single injurious event, rather than the continual degradation seen in osteoarthritis. It would also be of further interest to utilize cartilage from an osteoarthritis model (e.g., STR/ort mouse), or ultimately from human samples, to both validate our results here, and also examine the effects of hydrostatic pressure on gene expression in disease pathology. Despite these limitations, the current study was able to statistically differentiate the effects of hydrostatic pressure on chondrocytes.

In conclusion, we identified distinct differential gene expression signatures in hip cap and metatarsal cartilage explants, indicative of the divergent phenotypes of their residing chondrocytes. Our RNAseq studies examining the cartilage response to hydrostatic pressure provided evidence for injurious hydrostatic pressure to be associated with decreases in processes including cartilage development and chondrocyte differentiation. Together this informs on the potential benefits of hydrostatic pressure in cartilage tissue engineering strategies, which need to carefully consider the magnitude of application and the effects on gene expression. Further, we identified the differential expression of a number of genes that have previously been identified as osteoarthritis risk genes, including Ctsk and Chadl, further highlighting their potential as therapeutic targets. These data will therefore contribute to a better understanding of the role of hydrostatic pressure and the chondrocyte phenotype in health and osteoarthritis.

References

• 1

  ArcherCWFrancis-WestP. The chondrocyte. Int J Biochem Cell Biol (2003) 35(4):401–4. doi: 10.1016/S1357-2725(02)00301-1

  • CrossRef
  • Google Scholar

• 2

  PapTKorb-PapA. Cartilage damage in osteoarthritis and rheumatoid arthritis–two unequal siblings. Nat Rev Rheumatol (2015) 11(10):606–15. doi: 10.1038/nrrheum.2015.95

  • CrossRef
  • Google Scholar

• 3

  StainesKAPollardASMcGonnellIMFarquharsonCPitsillidesAA. Cartilage to bone transitions in health and disease. J Endocrinol (2013) 219(1):R1–R12. doi: 10.1530/JOE-13-0276

  • CrossRef
  • Google Scholar

• 4

  GoldringMB. Articular cartilage degradation in osteoarthritis. HSS J (2012) 8(1):7–9. doi: 10.1007/s11420-011-9250-z

  • CrossRef
  • Google Scholar

• 5

  PooleAR. Osteoarthritis as a whole joint disease. HSS J (2012) 8(1):4–6. doi: 10.1007/s11420-011-9248-6

  • CrossRef
  • Google Scholar

• 6

  SandellLJAignerT. Articular cartilage and changes in arthritis. An introduction: cell biology of osteoarthritis. Arthritis Res (2001) 3(2):107–13. doi: 10.1186/ar148

  • CrossRef
  • Google Scholar

• 7

  PitsillidesAABeierF. Cartilage biology in osteoarthritis–lessons from developmental biology. Nat Rev Rheumatol (2011) 7(11):654–63. doi: 10.1038/nrrheum.2011.129

  • CrossRef
  • Google Scholar

• 8

  HodgkinsonTKellyDCCurtinCMO'BrienFJ. Mechanosignalling in cartilage: an emerging target for the treatment of osteoarthritis. Nat Rev Rheumatol (2022) 18(2):67–84. doi: 10.1038/s41584-021-00724-w

  • CrossRef
  • Google Scholar

• 9

  PettenuzzoSArduinoABelluzziEPozzuoliAFontanellaCGRuggieriPet al. Biomechanics of chondrocytes and chondrons in healthy conditions and osteoarthritis: A review of the mechanical characterisations at the microscale. Biomedicines (2023) 11(7):1942. doi: 10.3390/biomedicines11071942

  • CrossRef
  • Google Scholar

• 10

  BelluzziETodrosSPozzuoliARuggieriPCarnielELBerardoA. Human cartilage biomechanics: experimental and theoretical approaches towards the identification of mechanical properties in healthy and osteoarthritic conditions. Processes (2023) 11(4):1014. doi: 10.3390/pr11041014

  • CrossRef
  • Google Scholar

• 11

  BachrachNMMowVCGuilakF. Incompressibility of the solid matrix of articular cartilage under high hydrostatic pressures. J Biomech (1998) 31(5):445–51. doi: 10.1016/S0021-9290(98)00035-9

  • CrossRef
  • Google Scholar

• 12

  O'ConorCJCaseNGuilakF. Mechanical regulation of chondrogenesis. Stem Cell Res Ther (2013) 4(4):61. doi: 10.1186/scrt211

  • CrossRef
  • Google Scholar

• 13

  ElderBDAthanasiouKA. Hydrostatic pressure in articular cartilage tissue engineering: from chondrocytes to tissue regeneration. Tissue Eng Part B Rev (2009) 15(1):43–53. doi: 10.1089/ten.teb.2008.0435

  • CrossRef
  • Google Scholar

• 14

  PattappaGZellnerJJohnstoneBDochevaDAngeleP. Cells under pressure - the relationship between hydrostatic pressure and mesenchymal stem cell chondrogenesis. Eur Cell Mater (2019) 37:360–81. doi: 10.22203/eCM.v037a22

  • CrossRef
  • Google Scholar

• 15

  MizunoSOgawaR. Using changes in hydrostatic and osmotic pressure to manipulate metabolic function in chondrocytes. Am J Physiol Cell Physiol (2011) 300(6):C1234–45. doi: 10.1152/ajpcell.00309.2010

  • CrossRef
  • Google Scholar

• 16

  CheleschiSTentiSBarbarinoMGiannottiSBellisaiFFratiEet al. Exploring the crosstalk between hydrostatic pressure and adipokines: an in vitro study on human osteoarthritic chondrocytes. Int J Mol Sci (2021) 22(5):2745. doi: 10.3390/ijms22052745

  • CrossRef
  • Google Scholar

• 17

  HallACUrbanJPGehlKA. The effects of hydrostatic pressure on matrix synthesis in articular cartilage. J Orthop Res (1991) 9(1):1–10. doi: 10.1002/jor.1100090102

  • CrossRef
  • Google Scholar

• 18

  NakamuraSAraiYTakahashiKATerauchiROhashiSMazdaOet al. Hydrostatic pressure induces apoptosis of chondrocytes cultured in alginate beads. J Orthop Res (2006) 24(4):733–9. doi: 10.1002/jor.20077

  • CrossRef
  • Google Scholar

• 19

  TakahashiKKuboTAraiYKitajimaITakigawaMImanishiJet al. Hydrostatic pressure induces expression of interleukin 6 and tumour necrosis factor alpha mRNAs in a chondrocyte-like cell line. Ann Rheum Dis (1998) 57(4):231–6. doi: 10.1136/ard.57.4.231

  • CrossRef
  • Google Scholar

• 20

  InoueHAraiYKishidaTTerauchiRHonjoKNakagawaSet al. Hydrostatic pressure influences HIF-2 alpha expression in chondrocytes. Int J Mol Sci (2015) 16(1):1043–50. doi: 10.3390/ijms16011043

  • CrossRef
  • Google Scholar

• 21

  MontagneKOnumaYItoYAikiYFurukawaKSUshidaT. High hydrostatic pressure induces pro-osteoarthritic changes in cartilage precursor cells: A transcriptome analysis. PloS One (2017) 12(8):e0183226. doi: 10.1371/journal.pone.0183226

  • CrossRef
  • Google Scholar

• 22

  ParkkinenJJIkonenJLammiMJLaakkonenJTammiMHelminenHJ. Effects of cyclic hydrostatic pressure on proteoglycan synthesis in cultured chondrocytes and articular cartilage explants. Arch Biochem Biophys (1993) 300(1):458–65. doi: 10.1006/abbi.1993.1062

  • CrossRef
  • Google Scholar

• 23

  MontagneKFurukawaKSTaninakaYNgaoBUshidaT. Modulation of the long non-coding RNA Mir155hg by high, but not moderate, hydrostatic pressure in cartilage precursor cells. PloS One (2022) 17(12):e0275682. doi: 10.1371/journal.pone.0275682

  • CrossRef
  • Google Scholar

• 24

  KawanishiMOuraAFurukawaKFukubayashiTNakamuraKTateishiTet al. Redifferentiation of dedifferentiated bovine articular chondrocytes enhanced by cyclic hydrostatic pressure under a gas-controlled system. Tissue Eng (2007) 13(5):957–64. doi: 10.1089/ten.2006.0176

  • CrossRef
  • Google Scholar

• 25

  RiederBWeihsAMWeidingerASzwarcDNurnbergerSRedlHet al. Hydrostatic pressure-generated reactive oxygen species induce osteoarthritic conditions in cartilage pellet cultures. Sci Rep (2018) 8(1):17010. doi: 10.1038/s41598-018-34718-8

  • CrossRef
  • Google Scholar

• 26

  KunitomoTTakahashiKAAraiYSakaoKHonjoKSaitoMet al. Influence of extracellular matrix on the expression of inflammatory cytokines, proteases, and apoptosis-related genes induced by hydrostatic pressure in three-dimensionally cultured chondrocytes. J Orthop Sci (2009) 14(6):776–83. doi: 10.1007/s00776-009-1393-0

  • CrossRef
  • Google Scholar

• 27

  ToyodaTSeedhomBBYaoJQKirkhamJBrookesSBonassWA. Hydrostatic pressure modulates proteoglycan metabolism in chondrocytes seeded in agarose. Arthritis Rheumatol (2003) 48(10):2865–72. doi: 10.1002/art.11250

  • CrossRef
  • Google Scholar

• 28

  WengerRHansMGWelterJFSolchagaLASheuYRMalemudCJ. Hydrostatic pressure increases apoptosis in cartilage-constructs produced from human osteoarthritic chondrocytes. Front Biosci (2006) 11:1690–5. doi: 10.2741/1914

  • CrossRef
  • Google Scholar

• 29

  ZhaoXHuaYWangTCiZZhangYWangXet al. In vitro cartilage regeneration regulated by a hydrostatic pressure bioreactor based on hybrid photocrosslinkable hydrogels. Front Bioeng Biotechnol (2022) 10:916146. doi: 10.3389/fbioe.2022.916146

  • CrossRef
  • Google Scholar

• 30

  ChangSHMoriDKobayashiHMoriYNakamotoHOkadaKet al. Excessive mechanical loading promotes osteoarthritis through the gremlin-1-NF-kappaB pathway. Nat Commun (2019) 10(1):1442. doi: 10.1038/s41467-019-09491-5

  • CrossRef
  • Google Scholar

• 31

  HenstockJRRotherhamMRoseJBEl HajAJ. Cyclic hydrostatic pressure stimulates enhanced bone development in the foetal chick femur in vitro. Bone (2013) 53(2):468–77. doi: 10.1016/j.bone.2013.01.010

  • CrossRef
  • Google Scholar

• 32

  Dehghan-BanianiDMehrjouBChuPKLeeWYWWuH. Recent advances in "Functional engineering of articular cartilage zones by polymeric biomaterials mediated with physical, mechanical, and biological/chemical cues". Adv Healthc Mater (2023) 12(10):e2202581. doi: 10.1002/adhm.202202581

  • CrossRef
  • Google Scholar

• 33

  HodderEGuppyFCovillDBushP. The effect of hydrostatic pressure on proteoglycan production in articular cartilage in vitro: a meta-analysis. Osteoarthritis Cartilage (2020) 28(8):1007–19. doi: 10.1016/j.joca.2020.03.021

  • CrossRef
  • Google Scholar

• 34

  SironenREloMKaarnirantaKHelminenHJLammiMJ. Transcriptional activation in chondrocytes submitted to hydrostatic pressure. Biorheology (2000) 37(1-2):85–93.

  • Google Scholar

• 35

  SamvelyanHJHughesDStevensCStainesKA. Models of osteoarthritis: relevance and new insights. Calcif Tissue Int (2021) 109(3):243–56. doi: 10.1007/s00223-020-00670-x

  • CrossRef
  • Google Scholar

• 36

  StainesKABrownGFarquharsonC. The ex vivo organ culture of bone. Methods Mol Biol (2019) 1914:199–215. doi: 10.1007/978-1-4939-8997-3_10

  • CrossRef
  • Google Scholar

• 37

  HoustonDAStainesKAMacRaeVEFarquharsonC. Culture of murine embryonic metatarsals: A physiological model of endochondral ossification. J Vis Exp (2016) 118. doi: 10.3791/54978-v

  • CrossRef
  • Google Scholar

• 38

  AndrewsSK,FSegonds-PichonABigginsLKruegerCWingettS. FastQC. Babraham, UK: Babraham Institute (2012).

  • Google Scholar

• 39

  EwelsPMagnussonMLundinSKallerM. MultiQC: summarize analysis results for multiple tools and samples in a single report. Bioinformatics (2016) 32(19):3047–8. doi: 10.1093/bioinformatics/btw354

  • CrossRef
  • Google Scholar

• 40

  LoveMIHuberWAndersS. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2. Genome Biol (2014) 15(12):550. doi: 10.1186/s13059-014-0550-8

  • CrossRef
  • Google Scholar

• 41

  YuGWangLGHanYHeQY. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS (2012) 16(5):284–7. doi: 10.1089/omi.2011.0118

  • CrossRef
  • Google Scholar

• 42

  CarlsonM. org.Mm.eg.db: Genome wide annotation for Mouse. R package version 3.8.2. Bioconductor (2019). doi: 10.18129/B9.bioc.org.Mm.eg.db

  • CrossRef
  • Google Scholar

• 43

  BoerCGHatzikotoulasKSouthamLStefansdottirLZhangYCoutinho de AlmeidaRet al. Deciphering osteoarthritis genetics across 826,690 individuals from 9 populations. Cell (2021) 184(18):4784–818 e17. doi: 10.1016/j.cell.2021.07.038

  • CrossRef
  • Google Scholar

• 44

  TachmazidouIHatzikotoulasKSouthamLEsparza-GordilloJHaberlandVZhengJet al. Identification of new therapeutic targets for osteoarthritis through genome-wide analyses of UK Biobank data. Nat Genet (2019) 51(2):230–6. doi: 10.1038/s41588-018-0327-1

  • CrossRef
  • Google Scholar

• 45

  AinNUMakitieONazS. Autosomal recessive chondrodysplasia with severe short stature caused by a biallelic COL10A1 variant. J Med Genet (2018) 55(6):403–7. doi: 10.1136/jmedgenet-2017-104885

  • CrossRef
  • Google Scholar

• 46

  IkegawaSNishimuraGNagaiTHasegawaTOhashiHNakamuraY. Mutation of the type X collagen gene (COL10A1) causes spondylometaphyseal dysplasia. Am J Hum Genet (1998) 63(6):1659–62. doi: 10.1086/302158

  • CrossRef
  • Google Scholar

• 47

  WarmanMLAbbottMApteSSHefferonTMcIntoshICohnDHet al. A type X collagen mutation causes Schmid metaphyseal chondrodysplasia. Nat Genet (1993) 5(1):79–82. doi: 10.1038/ng0993-79

  • CrossRef
  • Google Scholar

• 48

  GratalPMedieroASanchez-PernauteOPrieto-PotinILamuedraAHerrero-BeaumontGet al. Chondrocyte enlargement is a marker of osteoarthritis severity. Osteoarthritis Cartilage (2019) 27(8):1229–34. doi: 10.1016/j.joca.2019.04.013

  • CrossRef
  • Google Scholar

• 49

  HeYManon-JensenTArendt-NielsenLPetersenKKChristiansenTSamuelsJet al. Potential diagnostic value of a type X collagen neo-epitope biomarker for knee osteoarthritis. Osteoarthritis Cartilage (2019) 27(4):611–20. doi: 10.1016/j.joca.2019.01.001

  • CrossRef
  • Google Scholar

• 50

  von der MarkKFrischholzSAignerTBeierFBelkeJErdmannSet al. Upregulation of type X collagen expression in osteoarthritic cartilage. Acta Orthop Scand Suppl (1995) 266:125–9. doi: 10.3109/17453679509157667

  • CrossRef
  • Google Scholar

• 51

  LoriesRJBoonenSPeetersJde VlamKLuytenFP. Evidence for a differential association of the Arg200Trp single-nucleotide polymorphism in FRZB with hip osteoarthritis and osteoporosis. Rheumatol (Oxford) (2006) 45(1):113–4. doi: 10.1093/rheumatology/kei148

  • CrossRef
  • Google Scholar

• 52

  LoughlinJDowlingBChapmanKMarcellineLMustafaZSouthamLet al. Functional variants within the secreted frizzled-related protein 3 gene are associated with hip osteoarthritis in females. Proc Natl Acad Sci USA (2004) 101(26):9757–62. doi: 10.1073/pnas.0403456101

  • CrossRef
  • Google Scholar

• 53

  MinJLMeulenbeltIRiyaziNKloppenburgMHouwing-DuistermaatJJSeymourABet al. Association of the Frizzled-related protein gene with symptomatic osteoarthritis at multiple sites. Arthritis Rheumatol (2005) 52(4):1077–80. doi: 10.1002/art.20993

  • CrossRef
  • Google Scholar

• 54

  ThysenSLuytenFPLoriesRJ. Loss of Frzb and Sfrp1 differentially affects joint homeostasis in instability-induced osteoarthritis. Osteoarthritis Cartilage (2015) 23(2):275–9. doi: 10.1016/j.joca.2014.10.010

  • CrossRef
  • Google Scholar

• 55

  ChauMLuiJCLandmanEBSpathSSVortkampABaronJet al. Gene expression profiling reveals similarities between the spatial architectures of postnatal articular and growth plate cartilage. PloS One (2014) 9(7):e103061. doi: 10.1371/journal.pone.0103061

  • CrossRef
  • Google Scholar

• 56

  HaseebAKcRAngelozziMde CharleroyCRuxDTowerRJet al. SOX9 keeps growth plates and articular cartilage healthy by inhibiting chondrocyte dedifferentiation/osteoblastic redifferentiation. Proc Natl Acad Sci USA (2021) 118(8):e2019152118. doi: 10.1073/pnas.2019152118

  • CrossRef
  • Google Scholar

• 57

  LuiJCChauMChenWCheungCSHansonJRodriguez-CanalesJet al. Spatial regulation of gene expression during growth of articular cartilage in juvenile mice. Pediatr Res (2015) 77(3):406–15. doi: 10.1038/pr.2014.208

  • CrossRef
  • Google Scholar

• 58

  MowVCWangCCHungCT. The extracellular matrix, interstitial fluid and ions as a mechanical signal transducer in articular cartilage. Osteoarthritis Cartilage (1999) 7(1):41–58. doi: 10.1053/joca.1998.0161

  • CrossRef
  • Google Scholar

• 59

  HodgeWAFijanRSCarlsonKLBurgessRGHarrisWHMannRW. Contact pressures in the human hip joint measured in vivo. Proc Natl Acad Sci USA (1986) 83(9):2879–83. doi: 10.1073/pnas.83.9.2879

  • CrossRef
  • Google Scholar

• 60

  ZhuGMayer-WagnerSSchroderCWoiczinskiMBlumHLavagiIet al. Comparing effects of perfusion and hydrostatic pressure on gene profiles of human chondrocyte. J Biotechnol (2015) 210:59–65. doi: 10.1016/j.jbiotec.2015.06.409

  • CrossRef
  • Google Scholar

• 61

  KayaSBaileyKNSchurmanCAEvansDSAllistonT. Bone-cartilage crosstalk informed by aging mouse bone transcriptomics and human osteoarthritis genome-wide association studies. Bone Rep (2023) 18:101647. doi: 10.1016/j.bonr.2022.101647

  • CrossRef
  • Google Scholar

• 62

  YoultenSEKempJPLoganJGGhirardelloEJSergioCMDackMRGet al. Osteocyte transcriptome mapping identifies a molecular landscape controlling skeletal homeostasis and susceptibility to skeletal disease. Nat Commun (2021) 12(1):2444. doi: 10.1038/s41467-021-22517-1

  • CrossRef
  • Google Scholar

• 63

  CostaAGCusanoNESilvaBCCremersSBilezikianJP. Cathepsin K: its skeletal actions and role as a therapeutic target in osteoporosis. Nat Rev Rheumatol (2011) 7(8):447–56. doi: 10.1038/nrrheum.2011.77

  • CrossRef
  • Google Scholar

• 64

  KonttinenYTMandelinJLiTFSaloJLassusJLiljestromMet al. Acidic cysteine endoproteinase cathepsin K in the degeneration of the superficial articular hyaline cartilage in osteoarthritis. Arthritis Rheumatol (2002) 46(4):953–60. doi: 10.1002/art.10185

  • CrossRef
  • Google Scholar

• 65

  MorkoJPSoderstromMSaamanenAMSalminenHJVuorioEI. Up regulation of cathepsin K expression in articular chondrocytes in a transgenic mouse model for osteoarthritis. Ann Rheum Dis (2004) 63(6):649–55. doi: 10.1136/ard.2002.004671

  • CrossRef
  • Google Scholar

• 66

  SokiFNYoshidaRPagliaDNDuongLTHansenMFDrissiH. Articular cartilage protection in Ctsk(-/-) mice is associated with cellular and molecular changes in subchondral bone and cartilage matrix. J Cell Physiol (2018) 233(11):8666–76. doi: 10.1002/jcp.26745

  • CrossRef
  • Google Scholar

• 67

  ConnorJRLePageCSwiftBAYamashitaDBendeleAMMaulDet al. Protective effects of a cathepsin K inhibitor, SB-553484, in the canine partial medial meniscectomy model of osteoarthritis. Osteoarthritis Cartilage (2009) 17(9):1236–43. doi: 10.1016/j.joca.2009.03.015

  • CrossRef
  • Google Scholar

• 68

  ConaghanPGBowesMAKingsburySRBrettAGuillardGRizoskaBet al. Disease-modifying effects of a novel cathepsin K inhibitor in osteoarthritis: A randomized controlled trial. Ann Intern Med (2020) 172(2):86–95. doi: 10.7326/M19-0675

  • CrossRef
  • Google Scholar

• 69

  TillgrenVHoJCOnnerfjordPKalamajskiS. The novel small leucine-rich protein chondroadherin-like (CHADL) is expressed in cartilage and modulates chondrocyte differentiation. J Biol Chem (2015) 290(2):918–25. doi: 10.1074/jbc.M114.593541

  • CrossRef
  • Google Scholar

• 70

  TuerlingsMvan HoolwerffMHoutmanESuchimanELakenbergNMeiHet al. RNA sequencing reveals interacting key determinants of osteoarthritis acting in subchondral bone and articular cartilage: identification of IL11 and CHADL as attractive treatment targets. Arthritis Rheumatol (2021) 73(5):789–99. doi: 10.1002/art.41600

  • CrossRef
  • Google Scholar

• 71

  LeeWLeddyHAChenYLeeSHZelenskiNAMcNultyALet al. Synergy between Piezo1 and Piezo2 channels confers high-strain mechanosensitivity to articular cartilage. Proc Natl Acad Sci USA (2014) 111(47):E5114–22. doi: 10.1073/pnas.1414298111

  • CrossRef
  • Google Scholar

• 72

  Servin-VencesMRMoroniMLewinGRPooleK. Direct measurement of TRPV4 and PIEZO1 activity reveals multiple mechanotransduction pathways in chondrocytes. Elife (2017) 6:e21074. doi: 10.7554/eLife.21074

  • CrossRef
  • Google Scholar

• 73

  HalonenLPemmariANummenmaaEHamalainenMMoilanenTVuolteenahoKet al. Human osteoarthritic chondrocytes express nineteen different TRP-genes-TRPA1 and TRPM8 as potential drug targets. Int J Mol Sci (2023) 24(12):10057. doi: 10.3390/ijms241210057

  • CrossRef
  • Google Scholar