Identification of Key Genes and Pathways Associated with Sex Differences in Osteoarthritis Based on Bioinformatics Analysis

Causal association of menstrual reproductive factors on the risk of osteoarthritis: A univariate and multivariate Mendelian randomization study

OBJECTIVE

Several observational studies have revealed a potential relationship between menstrual reproductive factors (MRF) and osteoarthritis (OA). However, the precise causal relationship remains elusive. This study performed Mendelian randomization (MR) to provide deeper insights into this relationship.

METHODS

Utilizing summary statistics of genome-wide association studies (GWAS), we conducted univariate MR to estimate 2 menstrual factors (Age at menarche, AAM; Age at menopause, AMP) and 5 reproductive factors (Age at first live birth, AFB; Age at last live birth, ALB; Number of live births, NLB; Age first had sexual intercourse, AFSI; Age started oral contraceptive pill, ASOC) on OA (overall OA, OOA; knee OA, KOA and hip OA, HOA). The sample size of MRF ranged from 123846 to 406457, and the OA sample size range from 393873 to 484598. Inverse variance weighted (IVW) method was used as the primary MR analysis methods, and MR Egger, weighted median was performed as supplements. Sensitivity analysis was employed to test for heterogeneity and horizontal pleiotropy. Finally, multivariable MR was utilized to adjust for the influence of BMI on OA.

RESULTS

After conducting multiple tests (P<0.0023) and adjusting for BMI, MR analysis indicated that a lower AFB will increase the risk of OOA (odds ratio [OR] = 0.97, 95% confidence interval [CI]: 0.95-0.99, P = 3.39×10-4) and KOA (OR = 0.60, 95% CI: 0.47-0.78, P = 1.07×10-4). ALB (OR = 0.61, 95% CI: 0.45-0.84, P = 2.06×10-3) and Age AFSI (OR = 0.66, 95% CI: 0.53-0.82, P = 2.42×10-4) were negatively associated with KOA. In addition, our results showed that earlier AMP adversely affected HOA (OR = 1.12, 95% CI: 1.01-1.23, P = 0.033), and earlier ASOC promote the development of OOA (OR = 0.97, 95% CI: 0.95-1.00, P = 0.032) and KOA (OR = 0.58, 95% CI: 0.40-0.84, P = 4.49×10-3). ALB (OR = 0.98, 95% CI: 0.96-1.00, P = 0.030) and AFSI (OR = 0.98, 95% CI: 0.97-0.99, P = 2.66×10-3) also showed a negative association with OOA but they all did not pass multiple tests. The effects of AAM and NLB on OA were insignificant after BMI correction.

CONCLUSION

This research Certificates that Early AFB promotes the development of OOA, meanwhile early AFB, ALB, and AFSI are also risk factors of KOA. Reproductive factors, especially those related to birth, may have the greatest impact on KOA. It provides guidance for promoting women's appropriate age fertility and strengthening perinatal care.

Identification of CaMK4 as a sex-difference-related gene in knee osteoarthritis

Background

Osteoarthritis (OA) is a common degenerative joint disease with a higher prevalence in females than in males. Sex may be a key factor affecting the progression of OA. This study aimed to investigate critical sex-difference-related genes in patients with OA and confirm their potential roles in OA regulation.

Methods

OA datasets GSE12021, GSE55457, and GSE36700 were downloaded from the Gene Expression Omnibus database to screen OA-causing genes that are differentially expressed in the two sexes. Cytoscape was used to construct a protein-protein interaction network and determine hub genes. Synovial tissues of patients (male and female) with OA and female non-OA healthy controls were obtained to confirm the expression of hub genes and screen the key genes among them. Destabilization of the medial meniscus (DMM)-induced OA mice model was established to verify the screened key genes. Hematoxylin and eosin (HE) staining and Safranin O-fast green dye staining were employed to observe synovial inflammation and pathological cartilage status.

Results

The abovementioned three datasets were intersected to obtain 99 overlapping differentially expressed genes, of which 77 were upregulated and 22 were downregulated in female patients with OA. The hub genes screened were EGF, AQP4, CDC42, NTRK3, ERBB2, STAT1, and CaMK4. Among them, Ca2+/calmodulin-dependent protein kinase-4 (CaMK4) was identified as a key sex-related gene for OA. It was significantly higher in female OA patients than in the cases of male patients. Moreover, CaMK4 was significantly increased in female patients with OA compared with the female non-OA group. These results suggest that CaMK4 plays an important role in the progression of OA. OA mouse models demonstrated that CaMK4 expression in the mice knee joint synovial tissue elevated after DMM, with aggravated synovial inflammation and significant cartilage damage. Cartilage damage improved after intraperitoneal administration of the CaMK4 inhibitor KN-93.

Conclusions

CaMK4 is a key sex-related gene influencing the progression and pathogenesis of OA and may be considered as a new target for OA treatment.

How miR-31-5p and miR-33a-5p Regulates SP1/CX43 Expression in Osteoarthritis Disease: Preliminary Insights

Osteoarthritis (OA) is a degenerative bone disease that involved micro and macro-environment of joints. To date, there are no radical curative treatments for OA and novel therapies are mandatory. Recent evidence suggests the role of miRNAs in OA progression. In our previous studies, we demonstrated the role of miR-31-5p and miR-33a families in different bone regeneration signaling. Here, we investigated the role of miR-31-5p and miR-33a-5p in OA progression. A different expression of miR-31-5p and miR-33a-5p into osteoblasts and chondrocytes isolated from joint tissues of OA patients classified in based on different Kellgren and Lawrence (KL) grading was highlighted; and through a bioinformatic approach the common miRNAs target Specificity proteins (Sp1) were identified. Sp1 regulates the expression of gap junction protein Connexin43 (Cx43), which in OA drives the modification of i) osteoblasts and chondrocytes genes expression, ii) joint inflammation cytokines releases and iii) cell functions. Concerning this, thanks to gain and loss of function studies, the possible role of Sp1 as a modulator of CX43 expression through miR-31-5p and miR-33a-5p action was also evaluated. Finally, we hypothesize that both miRNAs cooperate to modulate the expression of SP1 in osteoblasts and chondrocytes and interfering, consequently, with CX43 expression, and they might be further investigated as new possible biomarkers for OA.

Women With Osteoarthritis Are at Increased Risk of Ischemic Stroke: A Population-Based Cohort Study

Background

Osteoarthritis (OA) is more prevalent in women with age. Comorbidities are prevalent in OA patients. In this study, we conducted a follow-up study to evaluate whether women with OA are at an increased risk of ischemic stroke using insurance claims data of Taiwan.

Methods

We identified 13,520 women with OA aged 20–99 newly diagnosed in 2000–2006 and 27,033 women without OA for comparison, frequency matched by age and diagnosis date. Women with baseline history of hypertension and other disorders associated with stroke were excluded for this study. Incident ischemic stroke was assessed by the end of 2013. A nested case-control analysis was used to identify factors associated with the stroke in the OA cohort.

Results

The incidence rate of ischemic stroke in the OA cohort was 1.5-fold greater than that in comparisons (1.93 versus 1.26 per 1,000 person-years), with an adjusted hazard ratio of 1.34 (95% confidence interval [CI], 1.09–1.66). The nested case-control analysis showed that stroke cases were twice as likely to develop hypertension during the follow-up period than controls without stroke. The ischemic stroke risk was significantly associated with hypertension (odds ratio [OR] 1.84; 95% CI, 1.37–2.46) and atrial fibrillation (OR 2.25; 95% CI, 1.24–4.09). Ischemic stroke was not associated with the use of non-steroidal anti-inflammatory drugs or aspirin.

Conclusion

Women with OA are at an elevated risk of ischemic stroke. A close monitoring of hypertension, atrial fibrillation, and other stroke related comorbidities is required for stroke prevention for OA patients.

Identification of Key Genes Associated with Changes in the Host Response to Severe Burn Shock: A Bioinformatics Analysis with Data from the Gene Expression Omnibus (GEO) Database

Background

Patients with severe burns continue to display a high mortality rate during the initial shock period. The precise molecular mechanism underlying the change in host response during severe burn shock remains unknown. This study aimed to identify key genes leading to the change in host response during burn shock.

Methods

The GSE77791 dataset, which was utilized in a previous study that compared hydrocortisone administration to placebo (NaCl 0.9%) in the inflammatory reaction of severe burn shock, was downloaded from the Gene Expression Omnibus (GEO) database and analyzed to identify the differentially expressed genes (DEGs). Functional enrichment analyses of Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways were performed. The protein–protein interaction (PPI) network of DEGs was constructed using the Search Tool for the Retrieval of Interacting Genes (STRING) database and then visualized in Cytoscape. In addition, important modules in this network were selected using the Molecular Complex Detection (MCODE) algorithm, and hub genes were identified in cytoHubba, a Cytoscape plugin.

Results

A total of 1059 DEGs (508 downregulated genes and 551 upregulated genes) were identified from the dataset. The DEGs enriched in GO terms and KEGG pathways were related to immune response. The PPI network contained 439 nodes and 2430 protein pairs. Finally, important modules and hub genes were identified using the different Cytoscape plugins. The key genes in burn shock were identified as arginase 1 (ARG1), cytoskeleton-associated protein (CKAP4), complement C3a receptor (C3AR1), neutrophil elastase (ELANE), gamma-glutamyl hydrolase (GGH), orosomucoid (ORM1), and quiescin sulfhydryl (QSOX1).

Conclusion

The DEGs, functional terms and pathways, and hub genes identified in the present study can help shed light on the molecular mechanism underlying the changes in host response during burn shock and provide potential targets for early detection and treatment of burn shock.