Zhiyao Ma1†
David Xinzheyang Li1,2†
Melanie Kunze1
Aillette Mulet-Sierra1
Lindsey Westover3
Adetola B. Adesida1*- 1Department of Surgery, Divisions of Orthopaedic Surgery, Surgical Research and Otolaryngology-Head and Neck Surgery, Faculty of Medicine and Dentistry, University of Alberta, Edmonton, AB, Canada
- 2Department of Civil and Environmental Engineering, University of Alberta, Edmonton, AB, Canada
- 3Department of Mechanical Engineering, University of Alberta, Edmonton, AB, Canada
Background: Osteoarthritis (OA) primarily affects mechanical load-bearing joints. The knee joint is the most impacted by OA. Knee OA (KOA) occurs in almost all demographic groups, but the prevalence and severity are disproportionately higher in females. The molecular mechanism underlying the pathogenesis and progression of KOA is unknown. The molecular basis of biological sex matters of KOA is not fully understood. Mechanical stimulation plays a vital role in modulating OA-related responses of load-bearing tissues. Mechanical unloading by simulated microgravity (SMG) induced OA-like gene expression in engineered cartilage, while mechanical loading by cyclic hydrostatic pressure (CHP), on the other hand, exerted a pro-chondrogenic effect. This study aimed to evaluate the effects of mechanical loading and unloading via CHP and SMG, respectively, on the OA-related profile changes of engineered meniscus tissues and explore biological sex-related differences.
Methods: Tissue-engineered menisci were made from female and male meniscus fibrochondrocytes (MFCs) under static conditions of normal gravity in chondrogenic media and subjected to SMG and CHP culture. Constructs were assayed via histology, immunofluorescence, GAG/DNA assays, RNA sequencing, and testing of mechanical properties.
Results: The mRNA expression of ACAN and COL2A1, was upregulated by CHP but downregulated by SMG. COL10A1, a marker for chondrocyte hypertrophy, was downregulated by CHP compared to SMG. Furthermore, CHP increased GAG/DNA levels and wet weight in both female and male donors, but only significantly in females. From the transcriptomics, CHP and SMG significantly modulated genes related to the ossification, regulation of ossification, extracellular matrix, and angiogenesis Gene Ontology (GO) terms. A clear difference in fold-change magnitude and direction was seen between the two treatments for many of the genes. Furthermore, differences in fold-change magnitudes were seen between male and female donors within each treatment. SMG and CHP also significantly modulated genes in OA-related KEGG pathways, such as mineral absorption, Wnt signalling pathway, and HIF-1 signalling pathway.
Conclusion: Engineered menisci responded to CHP and SMG in a sex-dependent manner. SMG may induce an OA-like profile, while CHP promotes chondrogenesis. The combination of SMG and CHP could serve as a model to study the early molecular events of KOA and potential drug-targetable pathways.
Introduction
Osteoarthritis (OA) is the most common form of degenerative disease and primarily affects loading-bearing joints, with the knee joint being the most prevalent (Nicolella et al., 2012; Boyan et al., 2013a). A hallmark feature of knee osteoarthritis (KOA) is the atrophy of articular cartilage. Surrounding joint tissues, including the menisci, will also undergo breakdown (Sun et al., 2010; Blackburn et al., 2016). These abnormal changes can lead to rapid loss of the knee joint’s function and mobility, making KOA a leading cause of physical disability (Murray et al., 2012). Although KOA occurs in almost all demographic groups, the prevalence and severity of KOA increases with age and is disproportionately higher in females than males (Badley and Kasman, 2004; Breedveld, 2004; O’Connor, 2007; Boyan et al., 2012; Pan et al., 2016). It has been reported by the World Health Organization that 9.6% of males and 18% of females above the age of 60 years have symptomatic OA (Osteoarthritis, 2014).
The molecular mechanisms and cellular events underlying the pathogenesis and progression of KOA are not well understood. As well, the molecular basis of biological sex matters has not been previously investigated. Currently, there is no consensus model to reflect the pathophysiology of KOA holistically. However, the molecular and cellular characteristics of KOA resemble the hypertrophic differentiation of chondrocytes as they progress to the bone during endochondral ossification (Dreier, 2010). Conveniently, this process includes the upregulation of hypertrophic markers COL10A1 and MMP13, which can be indicators for the KOA phenotype (Aigner et al., 1993; D’Angelo et al., 2000). Healthy chondrocytes, on the other hand, resist hypertrophic differentiation and lack expression of these hypertrophic markers (Leijten et al., 2012). Given the disproportionate incidence of KOA in females as compared to males, we reasonably expect the cellular and molecular characteristics of KOA to show sex-dependent differences.
Mechanical stimulation was reported by an abundance of studies to play a critical role in modulating OA-related responses of loading-bearing tissues. Applying mechanical loading to joints through regular exercise is essential to maintaining healthy cartilage and preventing breakdown from prolonged disuse (Bader et al., 20112011). Cyclic hydrostatic pressure (CHP) as a loading modality mimics physiological loading patterns and can be easily recreated in-vitro using specialized bioreactors (Elder and Athanasiou, 2009; Mellor et al., 2017). Studies have shown that CHP applied to engineered tissue constructs has induced a mostly pro-chondrogenic effect. For example, Zellner et al. applied dynamic hydrostatic pressure (cyclic at 1 Hz for 4 h per day, 0.55–5.03 MPa) for 7 days to cellular aggregates generated from inner and outer meniscus fibrochondrocytes (MFCs) (Zellner et al., 2015). After 14 additional days of static culture, aggregates loaded initially for the 7 days showed immunohistochemically enhanced chondrogenesis compared to unloaded controls (Zellner et al., 2015). Further, Gunja et al. applied dynamic hydrostatic pressure (0.1 Hz for 1 h every 3 days, 10 MPa) for 28 days to engineered tissue constructs using MFCs with added TGF-β1 growth factors (Gunja et al., 2009). Loaded tissue constructs with growth factors showed additive and synergistic increases in collagen deposition (approximately 2.5-fold), GAG deposition (2-fold), and enhanced compressive properties compared to unloaded controls without growth factors (Gunja et al., 2009).
Mechanical unloading of joints from long-term immobilization has been shown to induce cartilage atrophy that resembles characteristics of KOA. In a case study by Souza et al., joint immobilization of healthy individuals without prior history of OA resulted in magnetic resonance imaging (MRI) parameters of their knee articular cartilage that resemble KOA (Souza et al., 2012). When returning to standard weight-bearing, the MRI parameters for the joints were restored to baseline values consistent with healthy articular cartilage (Souza et al., 2012). Mechanical unloading has also been modelled by simulated microgravity (SMG) using rotating wall vessel bioreactors (Yu et al., 2011; Jin et al., 2013; Mayer-Wagner et al., 2014; Mellor et al., 2014; Mellor et al., 2017). Mayer-Wagner et al. applied simulated microgravity for 21 days to human mesenchymal stem cell (hMSC) pellets and found a decrease in histological staining of proteoglycans and collagen type-II compared to normal gravity controls (Mayer-Wagner et al., 2014). SMG pellets also showed a lower COL2A1/COL10A1 expression ratio suggesting that mechanical unloading via SMG reduced the chondrogenic differentiation of hMSCs (Mayer-Wagner et al., 2014). Finally, in a comparative study, CHP-loaded pellets from human adipose-derived stem cells showed increased expression of ACAN, SOX9, and COL2A1, and a 3-fold increase in GAG productions compared to unloaded SMG groups (Mellor et al., 2017). However, none of the above studies investigated the sex-dependent differences in the magnitude of differential modulation by mechanical loading and unloading via CHP and SMG, respectively.
A recent definition of OA from the Osteoarthritis Research Society International includes the menisci of the knee joint as a tissue undergoing breakdown and abnormal changes from the disease (Blackburn et al., 2016). Knee menisci undergoing OA has also been shown to have similar characteristics as knee articular cartilage undergoing OA, such as focal calcification and the increased expression of hypertrophic markers COL10A1 and MMP13 (Sun et al., 2010; Kiraly et al., 2017). This suggests that MFC may play an active role in the pathogenesis of KOA alongside knee articular chondrocytes, making MFC a reasonable cell option to model OA in-vitro.
Taken together, the goal of this study was to evaluate the effects of mechanical loading and unloading via CHP and SMG, respectively, and determine sex-dependent differences in the modulation of OA-related characteristics in bioengineered human meniscus tissues. This would serve as an OA disease model in-vitro to determine the cellular and molecular profiles responsible for the sex-dependent incidence of the disease.
Materials and Methods
The experiment is outlined in Figure 1. Most culture methods and assays were performed identically to those described in previous work (Liang et al., 2017; Szojka et al., 2021a; Szojka et al., 2021b).
FIGURE 1. Experiment outlines. Created with Biorender.com (2021).
Ethics Statement
Human non-osteoarthritic inner meniscus samples were collected from patients undergoing partial meniscectomies at the University of Alberta Hospital and Grey Nuns Community Hospital in Edmonton. The ethics of this study was approved by the Health Research Ethics Board of the University of Alberta. Non-identifying donor information is listed in Table 1.
TABLE 1. Non-identifying donor information.
Cell and Tissue Culture
Meniscus fibrochondrocytes (MFCs) were isolated from inner meniscus tissue samples by type II collagenase (0.15% w/v of 300 units/mg; Worthington) mediated digestion, followed by 48 h recovery. After recovery, cells were replated in tissue culture flasks at the density of 104 cells/cm2 and expanded for 1 week in high glucose Dulbecco’s modified Eagle’s medium (HG-DMEM) supplemented with 10% v/v heat-inactivated fetal bovine serum (FBS), 10 mM 4-(2-hydroxyethyl)-1-piperazineethanesulfonic acid (HEPES), 100 U/mL penicillin, 100 μg/ml streptomycin and 2 mM L-glutamine (PSG; Life Technologies, ON, Canada), 5 ng/ml of FGF-2 (Neuromics, MH, United States, catalogue #: PR80001) and 1 ng/ml of TGF-β1 (ProSpec, catalogue #: CYT-716) for 1 week
Expanded MFCs were resuspended in defined chondrogenic medium (HG-DMEM supplemented with HEPES, PSG, ITS +1 premix (Corning, Discovery Labware, Inc, MA, United States), 125 μg/ml of human serum albumin, 100 nM of dexamethasone, 365 μg/ml ascorbic acid 2-phosphate, 40 μg/ml of l-proline, and 10 ng/ml of TGF-β3) and seeded onto bovine type I collagen scaffolds (dimensions: 6 mm diameter, 3.5 mm height; Integra LifeSciences, NJ, United States) at the density of 5×106 cells/cm3. The cell-containing constructs were precultured statically in standard 24-well plates with 2.5 ml/construct of the chondrogenic medium described above for 2 weeks; media changes occurred once per week.
Mechanical Stimulation
After the 2-weeks preculture, tissue constructs were randomly assigned to a mechanical stimulation group. For the static control group, constructs were placed in a tissue culture tube (Sarstedt, Germany). For the mechanical unloading group, a commercially available bioreactor (RCCS-4; Synthecon Inc.) was used to culture tissue constructs in a simulated microgravity (SMG) environment. The rotation speed was adjusted over time to maintain constructs in suspension (30 rpm from day 1–2; 34 rpm from day 3–7; 37 rpm from day 8–13; 40 rpm from day 14–21). For the mechanical loading group, cyclic hydrostatic pressure (CHP) was applied to tissue constructs using a MechanoCulture TR (CellScale, Canada). Constructs were loaded 1 h per day and daily with 0.9 MPa cyclic hydrostatic pressure at the frequency of 1 Hz. When not loaded, tissue constructs were cultured in 6-well plates under static conditions. All experimental groups were cultured with chondrogenic medium with supplemented TGF-β3 growth factor, and the volume of medium per tissue construct (approximately 6.5 ml per tissue construct per week) was equivalent among the different groups. The mechanical stimulation was applied for 3 weeks, and medium change was performed once a week. At the end of the 3-weeks treatment, the CHP group was allowed 30 min of rest following the last loading event for gene expression changes to occur. Tissue replicates for RNA extraction were placed in TRIzol reagent and frozen at −80°C. Constructs from SMG and static control groups were harvested at approximately the same time.
Histology, Immunofluorescence, and Biochemical Analysis
The wet weight of tissue constructs (n = 5–8 replicates) intended for histology and biochemical analysis was recorded at the end of the experiment. Constructs (n = 2 replicates, only one replicate is presented) were then fixed in 1 ml of 10% v/v buffered formalin (Fisher Scientific, MA, United States) overnight at 4°C, dehydrated, embedded in paraffin wax, and sectioned at 5 µm thickness. Tissue sections approximately from the middle region of the constructs were stained with Safranin-O (Sigma-Aldrich, United States, #S2255-25G), Fast Green FGF (Sigma-Aldrich, United States, #F7258-25G), and Haematoxylin (Sigma-Aldrich, United States, #MHS32-1L) for histological examination of cell morphology and extracellular matrix deposition. Briefly, for immunofluorescence staining, tissue sections were labelled with primary antibody against human type I, type II, and type X collagen (1:200 dilution of anti-human rabbit type I collagen, Cedarlane, Canada, #CL50111AP-1; 1:200 dilution of mouse anti-human type II collagen, Developmental Studies Hybridoma Bank, United States, #II-II6B3; 1:100 dilution of rabbit anti-human type X collagen, Abcam, UK, #ab58632) and incubated overnight at 4°C (Anderson-Baron et al., 2021). On the next day, the secondary antibody (1:200 dilution of goat anti-rabbit, Abcam, UK, #ab150080; 1:200 dilution of goat anti-mouse, Abcam, UK, #ab150117) and DAPI (Cedarlane, Canada) was applied to visualize the stained components.
Biochemical assays quantified the total content of glycosaminoglycan (GAG) and DNA. Tissue constructs (n = 4 replicates for donors F1-3 and M1-3; n = 3 replicates for donors F4 and M4) were digested overnight with proteinase K (Sigma-Aldrich, United States, #P2308) at 56°C. GAG content was measured with a 1,9-dimethyl methylene blue assay (DMMB, Sigma-Aldrich, United States, #341088). Chondroitin sulphate (Sigma-Aldrich, United States, #C8529) was used to generate the standard curve. DNA content was measured with a CyQuant cell proliferation assay kit (ThermoFisher Scientific, United States, #C7026) with different dilutions of supplied bacteriophage λ DNA as the standard.
Mechanical Property Assessment
Detailed sequence of the mechanical testing protocol is included in Supplementary Material S3. A stepwise stress relaxation test (Supplementary Figure S1) was used to assess the mechanical properties of tissue constructs (n = 2 replicates for donor F1-3 and M1-3) with the BioDynamic 5210 system (TA Instruments, United States). The cross-section areas of tissue constructs were measured before mechanical tests. For the test, constructs were placed between two platens and the initial height was determined by bringing tissue to near contact with platens. Constructs were first preconditioned by 15 cycles of sine wave dynamic loading with the amplitude of 5% tissue height at the frequency of 1 Hz. The following stress relaxation test consisted of 3 incremental strain steps. In the first two steps, the constructs were subjected to a 10% strain ramp at the rate of 50% strain/sec followed by 5 min relaxation under constant strain. In the third step, the relaxation time was adjusted to 10 min. All tested constructs were able to reach equilibrium within the given relaxation period. Force was recorded as a function of time, and stress was calculated by normalizing force to construct’s cross-section area. The peak modulus was calculated by dividing the maximum stress measured immediately after each strain increment by the strain increment. The strain was applied in 10% increments up to a maximum of 30% strain.
RNA Extraction, qPCR, and Next-Generation Sequencing
Tissue constructs (n = 2 for donors F1-3 and M1-3; n = 3 replicates for donors F4 and M4) intended for transcriptome analysis were preserved in Trizol (Life Technologies, United States) immediately upon harvesting and stored at −80°C until RNA extraction. RNA was extracted and purified from ground constructs using PuroSPIN Total DNA Purification KIT (Luna Nanotech, Canada) following the manufacturer’s protocol. RNA was reversely transcribed into cDNA, and genes of interest were amplified by quantitative real-time polymerase chain reaction (RT-qPCR) using specific primers (Supplementary Table S1). The expression level of genes of interest was normalized to chosen housekeeping genes (i.e., B-actin, B2M, and YWHAZ) based on the coefficient of variation (CV) and M-value as measures of reference gene stability (Hellemans et al., 2007), and the data was presented using the 2-∆∆CT method (Livak and Schmittgen, 2001; Schmittgen and Livak, 2008). Next-generation RNA-sequencing was performed on the Illumina NextSeq 500 platform with paired-end 42 bp × 42 bp reads, and FastQ files were obtained for further bioinformatics analysis.
Bioinformatics
Next-generation sequencing data were analyzed with Partek® Flow® software (Version 10.0.21.0302, Copyright© 2021, Partek Inc, St. Louis, MO, United States). Raw input reads were first trimmed from the 3′ end to achieve a quality score beyond 20 and then aligned to the reference human genome hg38 using the STAR 2.7.3a aligner. Aligned data were quantified to a transcript model (hg38-RefSeq Transcripts 94–2020–05–01) using the Partek E/M algorithm. Genes with maximum read counts below 50 were filtered out to reduce noise. Quantified and filtered reads were normalized in sequential order using the Add: 1.0, TMM, and Log 2.0 methods. Statistical analysis was performed using analysis of variance (ANOVA) for biological sex and treatment. Within each sex, the donors were assigned as a random variable. Differentially expressed genes (DEGs) for each comparison were determined by p-values, adjusted p-values (q-values), and fold change (FC). Gene ontology enrichment, pathway enrichment, and the visualization of DEGs using Venn diagrams were all conducted in Partek.
Statistical analysis
Statistical analyses were performed in Prism 9 (GraphPad) and Partek® Flow® software. The statistical test used, and p-values and q-values are indicated in the respective figure legends. For analysis in Figure 3B, paired t-test were used between SMG and CHP groups within each sex and unpaired t-test were used between female and male groups within each mechanical treatment. For analysis in Figure 4A, a repeated measurement one-way ANOVA with Geisser-Greenhouse correction was used to compare mechanical treatment groups within each sex and a Tukey’s multiple comparison test was used to find the adjusted p-value between each comparison. Within each mechanical treatment, unpaired t-tests between female and male groups were used. For the male cohort, no statistical analysis was conducted for the SMG group in contraction due to limited data points. For analysis in Figure 4B, repeated measurement one-way ANOVA with Geisser-Greenhouse correction was used to compare different strain levels within each mechanical treatment group and a Tukey’s multiple comparison test was used to find the adjusted p-value between each comparison. Within each strain level, paired t-tests were used to compare CHP and SMG groups. For the male cohort, no statistical analysis was conducted for the SMG group due to limited data points.
Results
Dataset Overview
Transcriptome analysis included the expression profiles of 8 donors (4 females and 4 males), each individually exposed to static, mechanical loading (CHP), and mechanical unloading (SMG) conditions. After preprocessing as described in the methods, 13,361 genes were preserved for downstream analysis.
Transcriptome Profiles of the Engineered Meniscus to CHP and SMG
We first analyzed the overall effect of mechanical loading and unloading on the transcriptome profiles of all donors combined. Summarized results are shown in Figure 2. Fold-change of gene expression levels in the CHP and SMG groups were calculated by normalizing to its corresponding static group. Differentially expressed genes (DEGs) were defined as genes with expression fold-change over 2 and q-value less than 0.05. Mechanical loading from CHP significantly modulated 236 genes, while mechanical unloading from SMG significantly modulated 388 genes. The overlay of DEGs between the two mechanical stimulation groups showed only a small proportion of common DEGs (52 genes), whereas the majority of DEGs were uniquely modulated by CHP (184 genes) and SMG (336 genes) (Figure 2A). These results indicated that CHP and SMG distinctly modulate the transcriptome profile of donors in this study.
FIGURE 2. Effect of CHP and SMG on the transcriptome profile alteration of engineered meniscus tissues. (A) Overlap of differentially expressed genes (DEGs) by CHP and SMG. (B) Top non-redundant Gene Ontology (GO) terms enriched by DEGs of CHP and SMG. DEGs were identified based on all 8 donors. Top significant enriched GO terms were selected by p-value and plotted to the number of genes included in each term.
The most significant Gene Ontology (GO) terms enriched by the DEGs for CHP and SMG were examined next (Figure 2B). Although the top 3 most significantly enriched Gene Ontology (GO) in biological components were identical for CHP and SMG (“extracellular space,” “extracellular region,” and “extracellular matrix”), the included gene expression profiles were different between the treatment groups. The top 20 genes with the highest absolute fold change participating in the ECM relevant activities in CHP and SMG groups are listed in Table 2. For CHP, most of the genes are signalling molecules or proteins associated with ECM structure remodelling. For SMG, many of the strongly regulated genes played a more general role, such as various growth factors coding genes: IGFBP1, TGFA, and NGF. Among the top regulated genes, only NETO1 and OLFML2A were common between CHP and SMG; these two genes were upregulated in both treatment groups.
TABLE 2. Top 20 genes with the highest absolute fold change participating in the ECM relevant activities in CHP and SMG groups as compared to static controls.
Sex-Dependent Response of Engineered Meniscus to CHP and SMG
Next, we sought to explore the sex-dependent differences in the engineered meniscus responses to CHP and SMG. Therefore, we separated donors into female and male cohorts and evaluated several factors involved in normal cartilage physiology and OA-related alterations. The GAG and type II collagen content are often used to characterize the degree of chondrogenesis for cartilage tissues, and higher content is linked to higher chondrogenic capacity. Type X collagen, on the other hand, is a hypertrophic marker that suggests an OA-like phenotype. The histological staining for GAG by Safranin-O and immunofluorescence labelled type II collagen (Figure 3A) showed highly variable chondrogenic capacities within the female and male cohorts in the baseline static control group. Regardless, tissue constructs exposed to CHP and SMG showed an increase in type II collagen content, with CHP having a generally more pronounced effect. Type X collagen staining intensity was also modulated by mechanical stimulations compared to baseline. A clear increase in type X collagen content was observed in the SMG group, while the CHP group showed comparable type X collagen intensity as the baseline.
FIGURE 3. Effect of CHP and SMG on the chondrogenic and hypertrophic differentiation potential of female and male donors. (A) Histological and immunofluorescent staining analysis. (B) Regulation of selected genes expression level. Gene expression level of individual samples were measured by RT-qPCR. Fold change of expression level was calculated by normalizing expression level in the CHP or SMG group to its corresponding static group. Scale bar: 100 μm *Represents p < 0.05.
To further pursue sex matters in responses to CHP and SMG, the gene expression of selected markers was examined by RT-qPCR (Figure 3B). In addition to COL1A2, COL2A1, and COL10A1, the transcription factor SOX9 and cartilage-specific proteoglycan core protein ACAN were also quantified. The observed increase of type II collagen level in the CHP group was confirmed quantitatively by gene expression results. The average COL2A1 expression level was upregulated 215.9-fold in the CHP group for the female cohort compared to a 21.4-fold increase in CHP for the male cohort. For ACAN, SOX9, and COL1A2, the CHP group had a higher average fold-change than the SMG group, but the differences between the sexes were not significant. For COL10A1, the average fold-change of expression level was comparable between the CHP and SMG groups for both females and males. However, when taking individual female donors into account, only female donor 2 showed a significant increase in COL10A1 expression level in CHP (8.9-fold) compared to SMG (5.7-fold) group. Generally, CHP, as compared to SMG, reduced COL10A1 expression level for the female cohort. The RT-qPCR data were consistent with the histological observations, providing additional evidence on the cellular and molecular influence of CHP and SMG. The RT-qPCR data also agreed with the RNA-sequencing data, thereby providing additional validation for the observed trends (Supplementary Figure S2).
Sex-dependent differences were also assessed in other related aspects (Figure 4A). Quantitative GAG/DNA measurements showed that CHP increased while SMG decreased the GAG production per MFC. The tissue wet weights showed similar trends as the GAG/DNA ratios; CHP tissues weighed more, and SMG tissues weighed less on average than the static control groups. This difference in tissue wet weight was significant in only the female cohort between treatments. Further, at the end of the mechanical stimulation period, tissue constructs from all three experimental groups contracted to certain degrees. The percentage of contraction (percentage reduction in area) was quantified, and comparison among groups showed that SMG tissues had increased contraction compared to the CHP tissues. The contraction plot is presented as the percentage contracted in Supplementary Figure S4.
FIGURE 4. Effect of CHP and SMG on chondrogenesis related factors of female and male donors (A) Biochemical and morphological analysis; contraction is calculated as the % of area lost as compared to the original area. No statistical analysis was conducted for the contraction data with the male SMG group due to limited data points (B) Mechanical property analysis. No statistical analysis was conducted for the mechanical property data with the male SMG group due to limited data points (C) Pearson correlation heatmap of analyzed factors. Fold change value of characterized factors was calculated by normalizing the CHP or SMG group to its corresponding static group. Heatmap was generated by calculating the pairwise Pearson correlation coefficient of included factors. *Represents p < 0.05.
The differences among stimulation groups were significant within the female and male cohorts while also significant between female CHP and male CHP groups. For the mechanical properties of tissue constructs, SMG groups showed higher peak modulus at all tested strain levels for both sexes as expected based on contraction results, and the differences of the average fold-change decreased with increasing strain level. No significant differences were observed between female and male cohorts (Figure 4B).
To better understand the relationship among all factors of interest and the sex-dependent differences, a Pearson correlation network (Figure 4C) was generated for both the female and male cohorts with data from all three treatment groups. Some observations mentioned above were further confirmed with correlation coefficients, such as the positive correlation between GAG/DNA and wet weight, the negative correlation between GAG/DNA and contraction, and the positive correlation between SOX9 and COL2A1 expression levels. Overall, the correlation analysis yielded different patterns between the female and male cohorts, especially for two examined genes, ACTA2 and TAGLN, that are characteristic of the contractile phenotype in dedifferentiated articular chondrocytes compared to the other factors (Parreno et al., 2017). The level of contraction was positively correlated with ACTA2 and TAGLN expression levels in both female and male groups, confirming their indication for contraction levels. However, the female cohort showed differently a mostly strong negative correlation between the contractile genes and the rest of the factors compared to the male cohort. ACTA2 and TAGLN activity is associated with cytoskeletal composition and structure (Killion et al., 2017; Tang and Gerlach, 2017), and thus results may suggest sex-dependent differences of cytoskeletal activity in response to mechanical stimulation.
Comparison of Female and Male Transcriptome Response to CHP and SMG
Finally, we investigated the sex-dependent difference in the global transcriptome profile in response to mechanical loading and unloading. Consistent with the DEGs analysis for all donors combined, CHP and SMG uniquely regulated a large proportion of genes within each sex group (Figure 5A). We identified the top 25 enriched KEGG pathways for female and male tissues under CHP and SMG using the corresponding DEG sets (Figure 5B). For CHP, the most enriched KEGG pathway for both female and male groups was “mineral absorption.” Other relevant terms, such as “IL-17 signalling pathway”, “HIF-1 signalling pathway”, and “glycosaminoglycan biosynthesis,” were also observed for both sex groups. Although some of the enriched pathways were shared, a distinct gene profile with different magnitudes and sometimes the direction of modulation was observed for each sex group (Table 3). For example, the NOTUM gene in the Wnt-signalling pathway was significantly upregulated by 6.7-fold in CHP for the female cohort and only 1.8-fold in CHP for the male cohort. A large proportion of the top 25 KEGG pathways enriched by SMG overlapped with the pathways in CHP, but the corresponding gene profiles were different. For example, SMG did not significantly regulate the NOTUM, and it showed the opposite direction between female and male cohorts. The MMP3 gene was downregulated by SMG but upregulated by CHP. Interestingly, in addition to the OA-related KEGG pathways, the “fluid shear stress and atherosclerosis” pathway was enriched for both female and male cohorts, suggesting an effect of mechanical stimulation from CHP and SMG on engineered meniscus constructs.
FIGURE 5. Transcriptome response of female and male donors to CHP and SMG. (A) Overlap of DEGs of female and male donors exposed to CHP and SMG. (B) Top enriched KEGG pathways by identified DEGs. KEGG terms were selected by p-value and plotted to the enrichment score.
TABLE 3. Key enriched KEGG pathways and associated differentially expressed genes (DEGs) in CHP and SMG groups as compared to static controls within female and male donor cohorts. *q < 0.05, **q < 0.01, ***q < 0.001 represent statistical difference of the two groups in each fold change value.
In relation to mechanobiology of the meniscus, we investigated the transcriptome profile changes of various mechanosensitive molecules such as TRPV1/4, PIEZO1, TMEM63A/B/C, and RUNX2 (McNulty and Guilak, 2015; Servin-Vences et al., 2017; Agarwal et al., 2021; Hwang et al., 2021; Lee et al., 2021). Furthermore, the mechanotransduction function of FOSB in MFCs (Szojka et al., 2021a) was reported by Szojka et al., and its function in the IL-17 signalling pathway was well demonstrated (Benderdour et al., 2002). Finally, Vyhlidal et al. recognized the potential of caveolae molecules such as CAV1/2 in the mechanotransductive mechanism of the meniscus, but this needs to be further verified in future studies (Vyhlidal and Adesida, 2021). The transcriptome changes are summarized in Table 4. Within our dataset, only CAV1, CAV2, and FOSB are significant by q-value for loading regime comparison and sex-based comparison. CAV2 was regulated in the opposite direction by CHP and SMG. CHP upregulated FOSB by 49.5-fold compared to 3.4-fold by SMG. The difference was even larger in female tissues (137.8-fold by CHP and 2.8-fold by SMG), and in comparison, the fold change of FOSB expression was similar by CHP and SMG for male tissues (9.7-fold and 3.7-fold). Other mechanosensitive molecules were not significantly modulated from the treatments. The transcriptome profile changes suggest sex-dependent differences in the mechanotransduction mechanisms as well as varying capabilities to sense cytoskeletal structural changes.
TABLE 4. Relevant mechanosensitive molecules in CHP and SMG groups as compared to static controls within female and male donor cohorts. *q < 0.05, **q < 0.01, ***q < 0.001 represent statistical difference of the two groups in each fold change value.
Discussion
Currently, there is no existing model to reflect the physiological mechanism of OA holistically, but several in-vitro and in-vivo models have been developed to answer questions regarding the mechanisms of OA (Cope et al., 2019). This study aimed to evaluate the effects of mechanical loading and unloading via CHP and SMG, respectively, on the OA-related profile changes of engineered meniscus tissues and explore biological sex-related differences. This can serve as an in-vitro model to investigate the cellular and molecular profiles responsible for the sex-dependent incidence of OA disease.
Cartilage is a highly mechanosensitive tissue, and appropriate levels of mechanical stimulation are crucial for homeostasis and healthy cartilage development. Mechanical stimuli are transmitted by the pericellular matrix (PCM) (Melrose et al., 2006; Youn et al., 2006) to the chondrocyte surface and sensed by mechano-receptors, triggering a cascade of downstream activities (Zhao et al., 2020). The importance of mechanical loading under normal gravity environments has been demonstrated by studies examining the protective effects of moderate loading against tissue degradation as well as investigating the unwanted consequences from unloading. From several in-vitro and in-vivo models, mechanical loading has been shown to attenuate inflammatory cytokine-induced expression of matrix-degrading enzyme (Ohtsuki et al., 2019), upregulating the content of sulphated glycosaminoglycan (sGAG), aggrecan, cartilage oligomeric matrix protein, type II collagen, and lubricin (PRG4) (Antunes et al., 2020), and modulating relevant pathways such as the HIF-1 (Huang et al., 2017) and IL-4 (He et al., 2019) signalling pathways. On the contrary, prolonged mechanical unloading by space flight has been shown to accelerate cartilage degeneration (Herranz et al., 2013; Fitzgerald, 2017; Fitzgerald et al., 2019). Although several studies suggested that SMG was beneficial for preserving a chondrogenic phenotype through induction of 3D aggregates from monolayer cultures (Aleshcheva et al., 2015; Wuest et al., 2018; Wehland et al., 2020) and promoted cartilaginous components deposition for scaffold cultures (Freed and Vunjak-Novakovic, 1997), the degree of hypertrophic differentiation was not investigated in these studies. Through regulation of key genes and molecular pathways, SMG was also reported to increase activities associated with cartilage catabolism (Stamenković et al., 2010; Wang et al., 2010; Laws et al., 2016) and promoted hypertrophic differentiation of chondrocytes (Jin et al., 2013; Weiss et al., 2017).
In addition to cartilage, the meniscus also plays a critical role in the biomechanics of the knee joint. Evidence has shown that the role of meniscus fibrochondrocytes (MFCs) in response to mechanical signals affects the physiological, pathological, and repair response of the meniscus (McNulty and Guilak, 2015). The PCM of MFCs is also involved in mechanotransduction, although it may play a protective role against larger stresses and strains (Gupta and Haut Donahue, 2006). Various in-vivo studies have documented that mechanical stimulation can drive both anti- and pro-inflammatory responses in MFCs (Ferretti et al., 2005; Killian et al., 2014), as well as the detrimental effects of mechanical unloading from joint immobilization on meniscus development, function, and repair (Videman et al., 1979; Klein et al., 1989; Anderson et al., 1993; Djurasovic et al., 1998; Mikic et al., 2000; Bray et al., 2001). Additionally, several in-vitro studies have shown the anabolic effect of mechanical loading on MFCs with enhanced meniscus ECM components like ACAN expression and collagens and mechanical properties (Aufderheide and Athanasiou, 2006; Puetzer et al., 2012). Finally, there have been various in-vitro meniscus repair models that utilized mechanical loading to suppress IL-1 mediated increases in MMP activity, enhance sGAG production, and increase integrative strength of the engineered tissue constructs (McNulty et al., 2007; McNulty and Guilak, 2008; McNulty et al., 2010; Riera et al., 2011; McNulty and Guilak, 2015).
Consistent with these previous reports, our results showed that mechanical loading via CHP increased the deposition of type II collagen and aggrecan, supported by immunofluorescence staining and gene expression analysis. In addition, the Safranin-O staining and biochemical quantification of GAG production per cell confirmed that CHP increased the chondrogenic potential of meniscus MFCs. As expected, the wet weight of engineered constructs in the CHP group was the highest among the three treatment groups since the GAG is a major water-binding component of cartilage (KIANI et al., 2002). The ratio of COL2A1 to COL1A2 has also been used to evaluate chondrogenic capacity (Wuest et al., 2018; Wehland et al., 2020), and in our results, CHP substantially increased the COL2A1/COL1A2 ratio compared to SMG groups. Although upregulation of type II collagen deposition and COL2A1 expression was observed in the SMG groups compared to static controls, the magnitude of upregulation was not less than in the CHP groups. The GAG/DNA ratio of the SMG group was also lower than static controls in addition to a stronger intensity of type X collagen staining. Taken together, our data highlights the detrimental effect of SMG on the chondrogenic capacity of the engineered meniscus and suggests that mechanical unloading increases the hypertrophic differentiation of MFCs, driving them to display OA-like characteristics.
To the best of our knowledge, this is one of the first studies to investigate the effect of mechanical loading and unloading on cartilage models by examining global transcriptome profile alterations. The ECM of chondrocytes is crucial for regulating key functions through receptor-mediated matrix-cell interactions. The composition of the ECM and the bound signalling molecules largely influences the chondrogenic capacity of embedded chondrocytes (Gao et al., 20142014). By overlaying the DEGs, we found that CHP and SMG regulated MFC functions through largely different mechanisms. As listed in Table 2, several Wnt-signaling pathway-related genes were strongly regulated by CHP, in addition to various matrix remodelling enzymes. But for SMG, the genes with the highest fold-change functioned in more general ways, such as the several growth factors encoding genes identified to regulate general development processes. The KEGG pathway analysis also showed that several chondrogenesis- and OA-related pathways such as mineral absorption, Wnt-signalling pathway, HIF-1 signalling pathways, and IL-17 signalling pathway were enriched by both CHP and SMG. However, the expression profile of related genes in terms of magnitude and direction was quite different (Table 3). It is worth mentioning that no sex hormone-related pathways were present in the top enriched pathways. Thus, the observed sex-dependent differences are independent of differences between sex-related hormones. Taken together, the comparative transcriptome analysis suggested a distinct effect of CHP and SMG on regulating chondrogenesis, but further investigations are needed to determine the specific underlying mechanisms.
The cellular and molecular mechanisms behind the well-documented sex discrepancy in OA incident rates are poorly understood. Many factors are believed to contribute to the higher incidence and severity of OA in females, such as age, psychosocial status, metabolic variables, hormonal differences, anatomical variations, and inflammatory disease (Ferre et al., 2019). While the focus of many previous studies was mainly on bone shape (Yang et al., 2014; Wise et al., 2016; Frysz et al., 2020), gait kinematics (Phinyomark et al., 2016; Ro et al., 2017; Allison et al., 2018), and sex hormones (Zazulak et al., 2006; Brennan et al., 2010; Boyan et al., 2013b; Wang et al., 2013; Antony et al., 2016), little effort was invested into determining differences in the global transcriptome profile. Two studies reported that female OA patients have higher levels of inflammatory cytokines in the synovial fluid than males (Hooshmand et al., 2015; Kosek et al., 2018), but the underlying signalling mechanism was not investigated. To explore the sex-dependent differences in OA pathogenesis, we separated our donor cohort based on sex and compared the OA-related characteristics.
Although the female and male donor cohorts included in this study showed highly variable trends in terms of histological and immunofluorescence staining for cartilage markers like aggrecan and type II collagen, an expected similar trend was observed for COL2A1, SOX9, ACAN, and COL2A1/COL1A2 average fold-change levels across donor and treatment groups. The average fold change level for COL10A1 was comparable between the CHP and SMG groups for both sexes. However, one female donor showed significantly higher fold-changes for COL10A1 while the remaining female cohort showed a generally reduced expression of COL10A1 in CHP compared to SMG, but not significant. The observed deposition of type X collagen in immunofluorescence staining shows high variability, and sex-dependent differences are difficult to elucidate. Our results suggest overall that at the protein level, the effect of mechanical stimulation is more dominant than sex differences for the deposition of cartilaginous components.
A significant sex-dependent difference was observed for tissue contraction between female and male CHP groups. Additionally, the Pearson correlation network generated for the female and male donor cohorts shows clear sex-dependent differences in expression trends between contractile genes ACTA2 and TAGLN in relation to the other factors investigated in this study. Firstly, the measured contraction based on the percentage of reduction in tissue area positively correlated with ACTA2 and TAGLN expression levels in both female and male cohorts. This suggests that these contractile genes are indeed correlated with physical contraction levels in our study. Interestingly, ACTA2 and TAGLN expression levels in relation to other factors show distinct sex-dependent differences. Among the female cohort, the contractile genes showed a mostly strong negative correlation with the other factors, while the genes in the male cohort showed a generally positive correlation. Mechanical loading has been suggested to mediate the function of chondrocytes by stimulating the reorganization of cytoskeleton (Killion et al., 2017; Tang and Gerlach, 2017). And simulated microgravity was also reported to alter the structure of cytoskeleton components (Wehland et al., 2020) and regulate the expression of several cytoskeletal genes (Aleshcheva et al., 2015). Since ACTA2 and TAGLN activity is associated with cytoskeletal composition and structure, the opposite trend in correlation with other factors between the female and male cohort may suggest sex-dependent differences from cytoskeletal activity in response to varying levels of mechanical stimulation. Furthermore, the mechano-sensitive gene FOSB in the IL-17 signalling pathway was identified in our dataset to show a 137.8-fold upregulation from CHP in the female cohort while only observing a 9.7-fold increase in the male cohort. This may suggest sex-dependent differences in the mechanotransduction mechanisms as well as varying capabilities to sense cytoskeletal structural changes arising from mechanical stimulation.
There are several limitations to our model explored in this study. Firstly, we did not have any tissue samples from OA patients to confirm the OA-phenotype observed in our model. However, the modulation of several OA markers (Kiraly et al., 2017) from SMG such as the upregulation of MMP13 and the increased staining of collagen type X, reasonably suggest that the SMG treatment is pushing the engineered tissues constructs towards an OA-phenotype, and this effect is likely to increase with longer treatment periods. Another limitation is that even prolonged joint mobility under normal gravity can induce an OA-phenotype, and thus the static group in this study may also serve as an alternate condition for simulating OA. However, the static group here is meant as a baseline control between the two mechanical treatments and serves to reduce donor-to-donor variability by normalizing measurements to the static control within each donor. Finally, a limitation in our model is that although the engineered meniscus constructs contain the necessary component of fibrocartilage, the tissue microenvironment experienced by the cells is likely different from that of the native meniscus. In particular, the cells in the native meniscus may respond differently to receptor-mediated cell-ECM interactions from mechanical load than in the engineered meniscus models, especially due to differences in the matrix stiffness. We hope to address this limitation in our future studies by evaluating cytoskeleton factors and measurements and using substrates with tunable stiffnesses such as hydrogels.
Conclusion
Taken together, our data suggest that engineered meniscus tissues responded to mechanical loading and unloading via CHP and SMG in a sex-dependent manner. Mechanical unloading via SMG was shown to induce an OA-like profile, while mechanical loading via CHP promotes elements of chondrogenesis. Within each mechanical stimulation group, female and male donor cohorts show sex-dependent differences in the magnitude and direction of many differentially expressed genes, as well as tissue contraction and correlation of contractile genes with the other factors investigated in this study. The combination of CHP and SMG can feasibly serve as an in-vitro model to study the cellular and molecular mechanisms of KOA and provide a platform for exploring potential drug-targetable pathways as therapeutics.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found below: GEO and GSE192982.
Author Contributions
All listed authors meet the Frontier in Bioengineering and Biotechnology criteria for authorship. ZM, DL, and AA designed the study. ZM, DL, and AM-S performed tissue culture. ZM, DL, and MK performed RT-qPCR expression analysis. ZM, DL, and AA were responsible for RNA-seq data analysis with Partek Flow software. ZM, DL, and LW were responsible for mechanical testing and analysis. ZM and DL performed the statistical analysis and prepared tables and figures. ZM, DL, and AA wrote the manuscript with input from all co-authors. LW and AA were responsible for acquiring financial support and supervision of the study.
Funding
ZM: NSERC (NSERC RGPIN-2018-06290 Adesida) DL: NSERC Undergraduate Student Research Awards MK: Alberta Cancer Foundation-Mickleborough Interfacial Biosciences Research Program (ACF-MIBRP 27128 Adesida) AM-S: Canadian Institutes of Health Research (CIHR MOP 125921 Adesida) LW: University of Alberta Pilot Seed Grant Program (UOFAB PSGP); University of Alberta Women and Children’s Health Research Institute Innovation Grant (UOFAB WCHRIIG 3126) AA: Natural Sciences and Engineering Research Council (NSERC RGPIN-2018-06290); Canadian Institutes of Health Research (CIHR MOP 125921 Adesida); Canada Foundation for Innovation (CFI 33786); University Hospital of Alberta Foundation (UHF; RES0028185 Adesida); Edmonton Orthopaedic Research Committee, Cliff Lede Family Charitable Foundation (RES00045921 Adesida); University of Alberta Pilot Seed Grant Program (UOFAB PSGP); University of Alberta Women and Children’s Health Research Institute Innovation Grant (UOFAB WCHRIIG 3126); the Alberta Cancer Foundation-Mickleborough Interfacial Biosciences Research Program (ACF-MIBRP 27128 Adesida). Research grant funding for the work was provided by Natural Sciences and Engineering Research Council (NSERC RGPIN-2018-06290 Adesida), NSERC RTI-2019-00310 Adesida; Canada Foundation for Innovation (CFI 33786); University Hospital of Alberta Foundation (UHF; RES0028185 Adesida); Edmonton Orthopaedic Research Committee, Cliff Lede Family Charitable Foundation (RES00045921 Adesida); University of Alberta Pilot Seed Grant Program (UOFAB PSGP); University of Alberta Women and Children’s Health Research Institute Innovation Grant (UOFAB WCHRIIG 3126); the Alberta Cancer Foundation-Mickleborough Interfacial Biosciences Research Program (ACF-MIBRP 27128 Adesida).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We thank Tara Stach and the team at the Biomedical Research Centre from the University of British Columbia for their help with RNA sequencing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2022.823679/full#supplementary-material
References
Agarwal, P., Lee, H.-p., Smeriglio, P., Grandi, F., Goodman, S., Chaudhuri, O., et al. (2021). A Dysfunctional TRPV4-Gsk3β Pathway Prevents Osteoarthritic Chondrocytes from Sensing Changes in Extracellular Matrix Viscoelasticity. Nat. Biomed. Eng. 5, 1472–1484. doi:10.1038/s41551-021-00691-3
Aigner, T., Reichenberger, E., Bertling, W., Kirsch, T., Stöß, H., and Mark, K. (1993). Type X Collagen Expression in Osteoarthritic and Rheumatoid Articular Cartilage. Virchows Archiv B Cell Pathol 63, 205–211. doi:10.1007/bf02899263
Aleshcheva, G., Wehland, M., Sahana, J., Bauer, J., Corydon, T. J., Hemmersbach, R., et al. (2015). Moderate Alterations of the Cytoskeleton in Human Chondrocytes after Short‐term Microgravity Produced by Parabolic Flight Maneuvers Could Be Prevented by Up‐regulation of BMP‐2 and SOX‐9. FASEB j. 29, 2303–2314. doi:10.1096/fj.14-268151
Allison, K., Hall, M., Wrigley, T. V., Pua, Y.-H., Metcalf, B., and Bennell, K. L. (2018). Sex-specific Walking Kinematics and Kinetics in Individuals with Unilateral, Symptomatic Hip Osteoarthritis: A Cross Sectional Study. Gait & Posture 65, 234–239. doi:10.1016/j.gaitpost.2018.07.179
Anderson, D. R., Gershuni, D. H., Nakhostine, M., and Danzig, L. A. (1993). The Effects of Non-weight-bearing and Limited Motion on the Tensile Properties of the Meniscus. Arthrosc. J. Arthroscopic Relat. Surg. 9, 440–445. doi:10.1016/s0749-8063(05)80319-6
Anderson-Baron, M., Liang, Y., Kunze, M., Mulet-Sierra, A., Osswald, M., Ansari, K., et al. (2021). Suppression of Hypertrophy during in vitro Chondrogenesis of Cocultures of Human Mesenchymal Stem Cells and Nasal Chondrocytes Correlates with Lack of in vivo Calcification and Vascular Invasion. Front. Bioeng. Biotechnol. 8. doi:10.3389/fbioe.2020.572356
Antony, B., Venn, A., Cicuttini, F., March, L., Blizzard, L., Dwyer, T., et al. (2016). Association of Body Composition and Hormonal and Inflammatory Factors with Tibial Cartilage Volume and Sex Difference in Cartilage Volume in Young Adults. Arthritis Care Res. 68, 517–525. doi:10.1002/acr.22715
Antunes, B. P., Vainieri, M. L., Alini, M., Monsonego-Ornan, E., Grad, S., and Yayon, A. (2020). Enhanced Chondrogenic Phenotype of Primary Bovine Articular Chondrocytes in Fibrin-Hyaluronan Hydrogel by Multi-Axial Mechanical Loading and FGF18. Acta Biomater. 105, 170–179. doi:10.1016/j.actbio.2020.01.032
Aufderheide, A. C., and Athanasiou, K. A. (2006). A Direct Compression Stimulator for Articular Cartilage and Meniscal Explants. Ann. Biomed. Eng. 34, 1463–1474. doi:10.1007/s10439-006-9157-x
Bader, D. L., Salter, D. M., and Chowdhury, T. T. (20112011). Biomechanical Influence of Cartilage Homeostasis in Health and Disease. Arthritis 2011, 1–16. doi:10.1155/2011/979032
Badley, E. M., and Kasman, N. M. (2004). The Impact of Arthritis on Canadian Women. BMC Women's Health 4, S18. doi:10.1186/1472-6874-4-s1-s18
Benderdour, M., Tardif, G., Pelletier, J. P., Di Battista, J. A., Reboul, P., Ranger, P., et al. (2002). Interleukin 17 (IL-17) Induces Collagenase-3 Production in Human Osteoarthritic Chondrocytes via AP-1 Dependent Activation: Differential Activation of AP-1 Members by IL-17 and IL-1beta. J. Rheumatol. 29, 1262–1272.
Blackburn, S., Research User Group, C., Rhodes, C., Higginbottom, A., and Dziedzic, K. (2016). The OARSI Standardised Definition of Osteoarthritis: A Lay Version. Osteoarthritis and Cartilage 24, S192. doi:10.1016/j.joca.2016.01.379
Boyan, B. D., Hart, D. A., Enoka, R. M., Nicolella, D. P., Resnick, E., Berkley, K. J., et al. (2013). Hormonal Modulation of Connective Tissue Homeostasis and Sex Differences in Risk for Osteoarthritis of the Knee. Biol. Sex Dif 4, 3. doi:10.1186/2042-6410-4-3
Boyan, B. D., Tosi, L., Coutts, R., Enoka, R., Hart, D. A., Nicolella, D. P., et al. (2012). Sex Differences in Osteoarthritis of the Knee. J. Am. Acad. Orthopaedic Surgeons 20, 668–669. doi:10.5435/jaaos-20-10-668
Boyan, B. D., Tosi, L. L., Coutts, R. D., Enoka, R. M., Hart, D. A., Nicolella, D. P., et al. (2013). Addressing the Gaps: Sex Differences in Osteoarthritis of the Knee. Biol. Sex Dif 4, 4. doi:10.1186/2042-6410-4-4
Bray, R. C., Smith, J. A., Eng, M. K., Leonard, C. A., Sutherland, C. A., and Salo, P. T. (2001). Vascular Response of the Meniscus to Injury: Effects of Immobilization. J. Orthop. Res. 19, 384–390. doi:10.1016/s0736-0266(00)00037-1
Breedveld, F. C. (2004). Osteoarthritis--the Impact of a Serious Disease. Rheumatology 43, 4i–8. doi:10.1093/rheumatology/keh102
Brennan, S. L., Cicuttini, F. M., Shortreed, S., Forbes, A., Jones, G., Stuckey, S. L., et al. (2010). Women Lose Patella Cartilage at a Faster Rate Than Men: A 4.5-year Cohort Study of Subjects with Knee OA. Maturitas 67, 270–274. doi:10.1016/j.maturitas.2010.07.008
Cope, P. J., Ourradi, K., Li, Y., and Sharif, M. (2019). Models of Osteoarthritis: the Good, the Bad and the Promising. Osteoarthritis and Cartilage 27, 230–239. doi:10.1016/j.joca.2018.09.016
D’Angelo, M., Yan, Z, Nooreyazdan, M., Pacifici, M, Sarment, D. S., Billings, P. C., et al. (2000). MMP-13 Is Induced during Chondrocyte Hypertrophy. J. Cell Biochem. 77, 678–693.
Djurasovic, M., Aldridge, J. W., Grumbles, R., Rosenwasser, M. P., Howell, D., and Ratcliffe, A. (1998). Knee Joint Immobilization Decreases Aggrecan Gene Expression in the Meniscus. Am. J. Sports Med. 26, 460–466. doi:10.1177/03635465980260032101
Dreier, R. (2010). Hypertrophic Differentiation of Chondrocytes in Osteoarthritis: the Developmental Aspect of Degenerative Joint Disorders. Arthritis Res. Ther. 12, 216. doi:10.1186/ar3117
Elder, B. D., and Athanasiou, K. A. (2009). Hydrostatic Pressure in Articular Cartilage Tissue Engineering: From Chondrocytes to Tissue Regeneration. Tissue Eng. B: Rev. 15, 43–53. doi:10.1089/ten.teb.2008.0435
Ferre, I. M., Roof, M. A., Anoushiravani, A. A., Wasterlain, A. S., and Lajam, C. M. (2019). Understanding the Observed Sex Discrepancy in the Prevalence of Osteoarthritis. JBJS Rev. 7, e8. doi:10.2106/jbjs.rvw.18.00182
Ferretti, M., Srinivasan, A., Deschner, J., Gassner, R., Baliko, F., Piesco, N., et al. (2005). Anti-inflammatory Effects of Continuous Passive Motion on Meniscal Fibrocartilage. J. Orthop. Res. 23, 1165–1171. doi:10.1016/j.orthres.2005.01.025
Fitzgerald, J. (2017). Cartilage Breakdown in Microgravity-A Problem for Long-Term Spaceflight? Npj Regen. Med. 2, 10. doi:10.1038/s41536-017-0016-1
Fitzgerald, J., Endicott, J., Hansen, U., and Janowitz, C. (2019). Articular Cartilage and Sternal Fibrocartilage Respond Differently to Extended Microgravity. npj Microgravity 5, 3. doi:10.1038/s41526-019-0063-6
Freed, L. E., and Vunjak-Novakovic, G. (1997). Microgravity Tissue Engineering. In Vitro Cell.Dev.Biol.-Animal 33, 381–385. doi:10.1007/s11626-997-0009-2
Frysz, M., Gregory, J., Aspden, R. M., Paternoster, L., and Tobias, J. H. (2020). Sex Differences in Proximal Femur Shape: Findings from a Population-Based Study in Adolescents. Sci. Rep. 10, 4612. doi:10.1038/s41598-020-61653-4
Gao, Y., Liu, S, Huang, J., Guo, W, Chen, J., Zhang, L., et al. (20142014). The ECM-Cell Interaction of Cartilage Extracellular Matrix on Chondrocytes. Biomed. Res. Int. 1–8. doi:10.1155/2014/648459
Gunja, N. J., Uthamanthil, R. K., and Athanasiou, K. A. (2009). Effects of TGF-Β1 and Hydrostatic Pressure on Meniscus Cell-Seeded Scaffolds. Biomaterials 30, 565–573. doi:10.1016/j.biomaterials.2008.10.007
Gupta, T., and Haut Donahue, T. L. (2006). Role of Cell Location and Morphology in the Mechanical Environment Around Meniscal Cells. Acta Biomater. 2, 483–492. doi:10.1016/j.actbio.2006.05.009
He, Z., Leong, D. J., Xu, L., Hardin, J. A., Majeska, R. J., Schaffler, M. B., et al. (2019). CITED2 Mediates the Cross‐talk between Mechanical Loading and IL‐4 to Promote Chondroprotection. Ann. N.Y. Acad. Sci. 1442, 128–137. doi:10.1111/nyas.14021
Hellemans, J., Mortier, G., de Paepe, A., Speleman, F., and Vandesompele, J. (2007). qBase Relative Quantification Framework and Software for Management and Automated Analysis of Real-Time Quantitative PCR Data. Genome Biol. 8, R19. doi:10.1186/gb-2007-8-2-r19
Herranz, R., Anken, R., Boonstra, J., Braun, M., Christianen, P. C. M., de Geest, M., et al. (2013). Ground-Based Facilities for Simulation of Microgravity: Organism-specific Recommendations for Their Use, and Recommended Terminology. Astrobiology 13, 1–17. doi:10.1089/ast.2012.0876
Hooshmand, S., Juma, S., Khalil, D. A., Shamloufard, P., and Arjmandi, B. H. (2015). Women with Osteoarthritis Have Elevated Synovial Fluid Levels of Insulin-like Growth Factor (IGF)-1 and IGF-Binding Protein-3. J. Immunoassay Immunochemistry 36, 284–294. doi:10.1080/15321819.2014.947431
Huang, Z., Zhou, M., Wang, Q., Zhu, M., Chen, S., and Li, H. (2017). Mechanical and Hypoxia Stress Can Cause Chondrocytes Apoptosis through Over-activation of Endoplasmic Reticulum Stress. Arch. Oral Biol. 84, 125–132. doi:10.1016/j.archoralbio.2017.09.021
Hwang, H. S., Park, I. Y., Hong, J. I., Kim, J. R., and Kim, H. A. (2021). Comparison of Joint Degeneration and Pain in Male and Female Mice in DMM Model of Osteoarthritis. Osteoarthritis and Cartilage 29, 728–738. doi:10.1016/j.joca.2021.02.007
Jin, L., Feng, G., Reames, D. L., Shimer, A. L., Shen, F. H., and Li, X. (2013). The Effects of Simulated Microgravity on Intervertebral Disc Degeneration. Spine J. 13, 235–242. doi:10.1016/j.spinee.2012.01.022
Kiani, C., Chen, L., Wu, Y. J., Yee, A. J., and Yang, B. B. (2002). Structure and Function of Aggrecan. Cell Res 12, 19–32. doi:10.1038/sj.cr.7290106
Killian, M. L., Haut, R. C., and Haut Donahue, T. L. (2014). Acute Cell Viability and Nitric Oxide Release in Lateral Menisci Following Closed-Joint Knee Injury in a Lapine Model of post-traumatic Osteoarthritis. BMC Musculoskelet. Disord. 15, 297. doi:10.1186/1471-2474-15-297
Killion, C. H., Mitchell, E. H., Duke, C. G., and Serra, R. (2017). Mechanical Loading Regulates Organization of the Actin Cytoskeleton and Column Formation in Postnatal Growth Plate. MBoC 28, 1862–1870. doi:10.1091/mbc.e17-02-0084
Kiraly, A. J., Roberts, A., Cox, M., Mauerhan, D., Hanley, E., and Sun, Y. (2017). Comparison of Meniscal Cell-Mediated and Chondrocyte-Mediated Calcification. Toorthj 11, 225–233. doi:10.2174/1874325001711010225
Klein, L., Heiple, K. G., Torzilli, P. A., Goldberg, V. M., and Burstein, A. H. (1989). Prevention of Ligament and Meniscus Atrophy by Active Joint Motion in a Non-weight-bearing Model. J. Orthop. Res. 7, 80–85. doi:10.1002/jor.1100070111
Kosek, E., Finn, A., Ultenius, C., Hugo, A., Svensson, C., and Ahmed, A. S. (2018). Differences in Neuroimmune Signalling between Male and Female Patients Suffering from Knee Osteoarthritis. J. Neuroimmunology 321, 49–60. doi:10.1016/j.jneuroim.2018.05.009
Laws, C. J., Berg-Johansen, B., Hargens, A. R., and Lotz, J. C. (2016). The Effect of Simulated Microgravity on Lumbar Spine Biomechanics: an in vitro Study. Eur. Spine J. 25, 2889–2897. doi:10.1007/s00586-015-4221-6
Lee, W., Nims, R. J., Savadipour, A., Zhang, Q., Leddy, H. A., Liu, F., et al. (2021). Inflammatory Signaling Sensitizes Piezo1 Mechanotransduction in Articular Chondrocytes as a Pathogenic Feed-Forward Mechanism in Osteoarthritis. Proc. Natl. Acad. Sci. USA 118, e2001611118. doi:10.1073/pnas.2001611118
Leijten, J. C. H., Emons, J., Sticht, C., van Gool, S., Decker, E., Uitterlinden, A., et al. (2012). Gremlin 1, Frizzled-Related Protein, and Dkk-1 Are Key Regulators of Human Articular Cartilage Homeostasis. Arthritis Rheum. 64, 3302–3312. doi:10.1002/art.34535
Liang, Y., Idrees, E., Andrews, S. H. J., Labib, K., Szojka, A., Kunze, M., et al. (2017). Plasticity of Human Meniscus Fibrochondrocytes: A Study on Effects of Mitotic Divisions and Oxygen Tension. Sci. Rep. 7, 12148. doi:10.1038/s41598-017-12096-x
Livak, K. J., and Schmittgen, T. D. (2001). Analysis of Relative Gene Expression Data Using Real-Time Quantitative PCR and the 2−ΔΔCT Method. Methods 25, 402–408. doi:10.1006/meth.2001.1262
Mayer-Wagner, S., Hammerschmid, F., Redeker, J. I., Schmitt, B., Holzapfel, B. M., Jansson, V., et al. (2014). Simulated Microgravity Affects Chondrogenesis and Hypertrophy of Human Mesenchymal Stem Cells. Int. Orthopaedics (Sicot) 38, 2615–2621. doi:10.1007/s00264-014-2454-3
McNulty, A. L., Estes, B. T., Wilusz, R. E., Weinberg, J. B., and Guilak, F. (2010). Dynamic Loading Enhances Integrative Meniscal Repair in the Presence of Interleukin-1. Osteoarthritis and Cartilage 18, 830–838. doi:10.1016/j.joca.2010.02.009
McNulty, A. L., and Guilak, F. (2008). Integrative Repair of the Meniscus: Lessons from in vitro Studies. Biorheology 45, 487–500. doi:10.3233/bir-2008-0489
McNulty, A. L., and Guilak, F. (2015). Mechanobiology of the Meniscus. J. Biomech. 48, 1469–1478. doi:10.1016/j.jbiomech.2015.02.008
McNulty, A. L., Moutos, F. T., Weinberg, J. B., and Guilak, F. (2007). Enhanced Integrative Repair of the Porcine Meniscus in vitro by Inhibition of Interleukin-1 or Tumor Necrosis Factor α. Arthritis Rheum. 56, 3033–3043. doi:10.1002/art.22839
Mellor, L. F., Baker, T. L., Brown, R. J., Catlin, L. W., and Oxford, J. T. (2014). Optimal 3D Culture of Primary Articular Chondrocytes for Use in the Rotating Wall Vessel Bioreactor. Aviat Space Environ. Med. 85, 798–804. doi:10.3357/asem.3905.2014
Mellor, L. F., Steward, A. J., Nordberg, R. C., Taylor, M. A., and Loboa, E. G. (2017). Comparison of Simulated Microgravity and Hydrostatic Pressure for Chondrogenesis of hASC. Aerospace Med. Hum. Perform. 88, 377–384. doi:10.3357/amhp.4743.2017
Melrose, J., Roughley, P., Knox, S., Smith, S., Lord, M., and Whitelock, J. (2006). The Structure, Location, and Function of Perlecan, a Prominent Pericellular Proteoglycan of Fetal, Postnatal, and Mature Hyaline Cartilages. J. Biol. Chem. 281, 36905–36914. doi:10.1074/jbc.m608462200
Mikic, B., Johnson, T. L., Chhabra, A. B., Schalet, B. J., Wong, M., and Hunziker, E. B. (2000). Differential Effects of Embryonic Immobilization on the Development of Fibrocartilaginous Skeletal Elements. J. Rehabil. Res. Dev. 37, 127–133.
Murray, C. J., Vos, T., Lozano, R., Naghavi, M., Flaxman, A. D., Michaud, C., et al. (2012). Disability-adjusted Life Years (DALYs) for 291 Diseases and Injuries in 21 Regions, 1990-2010: a Systematic Analysis for the Global Burden of Disease Study 2010. Lancet 380, 2197–2223. doi:10.1016/S0140-6736(12)61689-4
Nicolella, D. P., O’Connor, M. I., Enoka, R. M., Boyan, B. D., Hart, D. A., Resnick, E., et al. (2012). Mechanical Contributors to Sex Differences in Idiopathic Knee Osteoarthritis. Biol. Sex Dif 3, 28. doi:10.1186/2042-6410-3-28
O’Connor, M. I. (2007). Sex Differences in Osteoarthritis of the Hip and Knee. JAAOS - J. Am. Acad. Orthopaedic Surgeons 15.
Ohtsuki, T., Shinaoka, A., Kumagishi-Shinaoka, K., Asano, K., Hatipoglu, O. F., Inagaki, J., et al. (2019). Mechanical Strain Attenuates Cytokine-Induced ADAMTS9 Expression via Transient Receptor Potential Vanilloid Type 1. Exp. Cell Res. 383, 111556. doi:10.1016/j.yexcr.2019.111556
Osteoarthritis (2014). Bone and Joint Canada. http://boneandjointcanada.com/osteoarthritis/.
Pan, Q., O’Connor, M. I., Coutts, R. D., Hyzy, S. L., Olivares-Navarrete, R., Schwartz, Z., et al. (2016). Characterization of Osteoarthritic Human Knees Indicates Potential Sex Differences. Biol. Sex. Differ. 7, 27. doi:10.1186/s13293-016-0080-z
Parreno, J., Raju, S., Wu, P.-h., and Kandel, R. A. (2017). MRTF-A Signaling Regulates the Acquisition of the Contractile Phenotype in Dedifferentiated Chondrocytes. Matrix Biol. 62, 3–14. doi:10.1016/j.matbio.2016.10.004
Phinyomark, A., Osis, S. T., Hettinga, B. A., Kobsar, D., and Ferber, R. (2016). Gender Differences in Gait Kinematics for Patients with Knee Osteoarthritis. BMC Musculoskelet. Disord. 17, 157. doi:10.1186/s12891-016-1013-z
Puetzer, J. L., Ballyns, J. J., and Bonassar, L. J. (2012). The Effect of the Duration of Mechanical Stimulation and Post-Stimulation Culture on the Structure and Properties of Dynamically Compressed Tissue-Engineered Menisci. Tissue Eng. A 18, 1365–1375. doi:10.1089/ten.tea.2011.0589
Riera, K. M., Rothfusz, N. E., Wilusz, R. E., Weinberg, J. B., Guilak, F., and McNulty, A. L. (2011). Interleukin-1, Tumor Necrosis Factor-Alpha, and Transforming Growth Factor-Beta 1 and Integrative Meniscal Repair: Influences on Meniscal Cell Proliferation and Migration. Arthritis Res. Ther. 13, R187. doi:10.1186/ar3515
Ro, D. H., Lee, D. Y., Moon, G., Lee, S., Seo, S. G., Kim, S. H., et al. (2017). Sex Differences in Knee Joint Loading: Cross-Sectional Study in Geriatric Population. J. Orthop. Res. 35, 1283–1289. doi:10.1002/jor.23374
Schmittgen, T. D., and Livak, K. J. (2008). Analyzing Real-Time PCR Data by the Comparative CT Method. Nat. Protoc. 3, 1101–1108. doi:10.1038/nprot.2008.73
Servin-Vences, M. R., Moroni, M., Lewin, G. R., and Poole, K. (2017). Direct Measurement of TRPV4 and PIEZO1 Activity Reveals Multiple Mechanotransduction Pathways in Chondrocytes. eLife 6. doi:10.7554/eLife.21074
Souza, R. B., Baum, T., Wu, S., Feeley, B. T., Kadel, N., Li, X., et al. (2012). Effects of Unloading on Knee Articular Cartilage T1rho and T2 Magnetic Resonance Imaging Relaxation Times: A Case Series. J. Orthop. Sports Phys. Ther. 42, 511–520. doi:10.2519/jospt.2012.3975
Stamenković, V., Keller, G., Nesic, D., Cogoli, A., and Grogan, S. P. (2010). Neocartilage Formation in 1 G, Simulated, and Microgravity Environments: Implications for Tissue Engineering. Tissue Eng. Part. A. 16, 1729–1736. doi:10.1089/ten.tea.2008.0624
Sun, Y., Mauerhan, D. R., Honeycutt, P. R., Kneisl, J. S., Norton, H. J., Zinchenko, N., et al. (2010). Calcium Deposition in Osteoarthritic Meniscus and Meniscal Cell Culture. Arthritis Res. Ther. 12, R56. doi:10.1186/ar2968
Szojka, A. R. A., Li, D. X., Sopcak, M. E. J., Ma, Z, Kunze, M, Mulet-Sierra, A., et al. (2021). Mechano-Hypoxia Conditioning of Engineered Human Meniscus. Front. Bioeng. Biotechnol. 9. doi:10.3389/fbioe.2021.739438
Szojka, A. R., Marqueti, R. d. C., Li, D. X., Molter, C. W., Liang, Y., Kunze, M., et al. (2021). Human Engineered Meniscus Transcriptome after Short-Term Combined Hypoxia and Dynamic Compression. J. Tissue Eng. 12, 204173142199084. doi:10.1177/2041731421990842
Tang, D. D., and Gerlach, B. D. (2017). The Roles and Regulation of the Actin Cytoskeleton, Intermediate Filaments and Microtubules in Smooth Muscle Cell Migration. Respir. Res. 18, 54. doi:10.1186/s12931-017-0544-7
Videman, T., Eronen, I., Friman, C., and Langenskiöld, A. (1979). Glycosaminoglycan Metabolism of the Medial Meniscus, the Medial Collateral Ligament and the Hip Joint Capsule in Experimental Osteoarthritis Caused by Immobilization of the Rabbit Knee. Acta Orthopaedica Scand. 50, 465–470. doi:10.3109/17453677908989791
Vyhlidal, M. J., and Adesida, A. B. (2021). Mechanotransduction in Meniscus Fibrochondrocytes: What about Caveolae? J. Cell Physiol. doi:10.1002/jcp.30616
Wang, Q., Zheng, Y.-P., Wang, X.-Y., Huang, Y.-P., Liu, M.-Q., Wang, S.-Z., et al. (2010). Ultrasound Evaluation of Site-specific Effect of Simulated Microgravity on Articular Cartilage. Ultrasound Med. Biol. 36, 1089–1097. doi:10.1016/j.ultrasmedbio.2010.04.018
Wang, X. D., Kou, X. X., Meng, Z., Bi, R. Y., Liu, Y., Zhang, J. N., et al. (2013). Estrogen Aggravates Iodoacetate-Induced Temporomandibular Joint Osteoarthritis. J. Dent Res. 92, 918–924. doi:10.1177/0022034513501323
Wehland, M., Steinwerth, P., Aleshcheva, G., Sahana, J., Hemmersbach, R., Lützenberg, R., et al. (2020). Tissue Engineering of Cartilage Using a Random Positioning Machine. Ijms 21, 9596. doi:10.3390/ijms21249596
Weiss, W. M., Mulet-Sierra, A., Kunze, M., Jomha, N. M., and Adesida, A. B. (2017). Coculture of Meniscus Cells and Mesenchymal Stem Cells in Simulated Microgravity. npj Microgravity 3, 28. doi:10.1038/s41526-017-0032-x
Wise, B. L., Liu, F., Kritikos, L., Lynch, J. A., Parimi, N., Zhang, Y., et al. (2016). The Association of Distal Femur and Proximal Tibia Shape with Sex: The Osteoarthritis Initiative. Semin. Arthritis Rheum. 46, 20–26. doi:10.1016/j.semarthrit.2016.02.006
Wuest, S., Caliò, M., Wernas, T., Tanner, S., Giger-Lange, C., Wyss, F., et al. (2018). Influence of Mechanical Unloading on Articular Chondrocyte Dedifferentiation. Ijms 19, 1289. doi:10.3390/ijms19051289
Yang, B., Yu, J.-K., Zheng, Z.-Z., Lu, Z.-H., and Zhang, J.-Y. (2014). Comparative Study of Sex Differences in Distal Femur Morphology in Osteoarthritic Knees in a Chinese Population. PLoS ONE 9, e89394. doi:10.1371/journal.pone.0089394
Youn, I., Choi, J. B., Cao, L., Setton, L. A., and Guilak, F. (2006). Zonal Variations in the Three-Dimensional Morphology of the Chondron Measured in situ Using Confocal Microscopy. Osteoarthritis and Cartilage 14, 889–897. doi:10.1016/j.joca.2006.02.017
Yu, B., Yu, D., Cao, L., Zhao, X., Long, T., Liu, G., et al. (2011). Simulated Microgravity Using a Rotary Cell Culture System Promotes Chondrogenesis of Human Adipose-Derived Mesenchymal Stem Cells via the P38 MAPK Pathway. Biochem. Biophysical Res. Commun. 414, 412–418. doi:10.1016/j.bbrc.2011.09.103
Zazulak, B. T., Paterno, M., Myer, G. D., Romani, W. A., and Hewett, T. E. (2006). The Effects of the Menstrual Cycle on Anterior Knee Laxity. Sports Med. 36, 847–862. doi:10.2165/00007256-200636100-00004
Zellner, J., Mueller, M., Xin, Y., Krutsch, W., Brandl, A., Kujat, R., et al. (2015). Dynamic Hydrostatic Pressure Enhances Differentially the Chondrogenesis of Meniscal Cells from the Inner and Outer Zone. J. Biomech. 48, 1479–1484. doi:10.1016/j.jbiomech.2015.02.003
Keywords: disease modelling, human-engineered meniscus, mechanical stimulation, simulated microgravity, cyclic hydrostatic pressure
Citation: Ma Z, Li DX, Kunze M, Mulet-Sierra A, Westover L and Adesida AB (2022) Engineered Human Meniscus in Modeling Sex Differences of Knee Osteoarthritis in Vitro. Front. Bioeng. Biotechnol. 10:823679. doi: 10.3389/fbioe.2022.823679
Received: 29 November 2021; Accepted: 21 January 2022;
Published: 15 February 2022.
Edited by:
Mark Buckley, University of Rochester, United StatesReviewed by:
Csaba Matta, University of Debrecen, HungaryRobert Nims, Washington University in St. Louis, United States
Jennifer Puetzer, Virginia Commonwealth University, United States
Copyright © 2022 Ma, Li, Kunze, Mulet-Sierra, Westover and Adesida. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Adetola B. Adesida, YWRlc2lkYUB1YWxiZXJ0YS5jYQ==
†These authors have contributed equally to this work and share first authorship