Integrative analysis of genome-wide association studies and gene expression profiles identified candidate genes for osteoporosis in Kashin-Beck disease patients

Plasma metabolites and inflammatory proteins profiling predict outcome of Fufang Duzhong Jiangu granules treating Kashin-Beck disease

To investigate predictive biomarkers that could be used to identify patients' response to treatment, plasma metabolomics and proteomics analyses were performed in Kashin-Beck disease (KBD) patients treated with Fufang Duzhong Jiangu Granules (FDJG). Plasma was collected from 12 KBD patients before treatment and 1 month after FDJG treatment. LC-MS and olink proteomics were employed for obtaining plasma metabolomics profiling and inflammatory protein profiles. Patients were classified into responders and non-responders based on drug efficacy. Enrichment analyses of differential metabolites and proteins of the responders at baseline and after treatment were conducted to study the mechanism of drug action. Differential metabolites and proteins between the two groups were screened as biomarkers to predict the drug efficacy. The receiver operating characteristic curve was used to evaluate the prediction accuracy of biomarkers. The changes in metabolites and inflammatory proteins in responders after treatment reflected the mechanism of FDJG treatment for KBD, which may act on glycerophospholipid metabolism, d-glutamine and d-glutamate metabolism, nitrogen metabolism and NF-kappa B signaling pathway. Three metabolites were identified as potential predictors: N-undecanoylglycine, β-aminopropionitrile and PC [18:3(6Z,9Z,12Z)/20:4(8Z,11Z,14Z,17Z)]. For inflammatory protein, interleukin-8 was identified as a predictive biomarker to detect responders. Combined use of these four biomarkers had high predictive ability (area under the curve = 0.972).

A Metabolomics Study of Feces Revealed That a Disturbance of Selenium-Centered Metabolic Bioprocess Was Involved in Kashin-Beck Disease, an Osteoarthropathy Endemic to China

Background: Kashin-Beck disease (KBD) is a distinct osteoarthropathy in China with an unclear pathogenesis. This study aims to explore whether perturbations in the intestine metabolome could be linked to KBD individuals. Methods: An investigation was conducted in KBD endemic villages and fecal samples were collected. After applying inclusion and exclusion criteria, a total of 75 subjects were enrolled for this study, including 46 KBD (including 19 Grade I KBD and 27 Grade II KBD) and 29 controls. Untargeted metabolomics analysis was performed on the platform of UHPLC-MS. PLS-DA and OPLS-DA were conducted to compare the groups and identify the differential metabolites (DMs). Pathway analysis was conducted on MPaLA platform to explore the functional implication of the DMs. Results: Metabolomics analysis showed that compared with the control group, KBD individuals have a total of 584 differential metabolites with dysregulated levels such as adrenic acid (log2FC = -1.87, VIP = 4.84, p = 7.63 × 10-7), hydrogen phosphate (log2FC = -2.57, VIP = 1.27, p = 1.02 × 10-3), taurochenodeoxycholic acid (VIP = 1.16, log2FC = -3.24, p = 0.03), prostaglandin E3 (VIP = 1.17, log2FC = 2.67, p = 5.61 × 10-4), etc. Pathway analysis revealed several significantly perturbed pathways associated with KBD such as selenium micronutrient network (Q value = 3.11 × 10-3, Wikipathways), metabolism of lipids (Q value = 8.43 × 10-4, Reactome), free fatty acid receptors (Q value = 3.99 × 10-3, Reactome), and recycling of bile acids and salts (Q value = 2.98 × 10-3, Reactome). Subgroup comparisons found a total of 267 differential metabolites were shared by KBD vs. control, KBD II vs. control, and KBD I vs. control, while little difference was found between KBD II and KBD I (only one differential metabolite detected). Conclusions: KBD individuals showed distinct metabolic features characterized by perturbations in lipid metabolism and selenium-related bioprocesses. Our findings suggest that the loss of nutrients metabolism balance in intestine was involved in KBD pathogenesis. Linking the nutrients metabolism (especially selenium and lipid) to KBD cartilage damage should be a future direction of KBD study.

CTR9 drives osteochondral lineage differentiation of human mesenchymal stem cells via epigenetic regulation of BMP-2 signaling

Cell fate determination of human mesenchymal stem/stromal cells (hMSCs) is precisely regulated by lineage-specific transcription factors and epigenetic enzymes. We found that CTR9, a key scaffold subunit of polymerase-associated factor complex (PAFc), selectively regulates hMSC differentiation to osteoblasts and chondrocytes, but not to adipocytes. An in vivo ectopic osteogenesis assay confirmed the essentiality of CTR9 in hMSC-derived bone formation. CTR9 counteracts the activity of Enhancer Of Zeste 2 (EZH2), the epigenetic enzyme that deposits H3K27me3, in hMSCs. Accordingly, CTR9 knockdown (KD) hMSCs gain H3K27me3 mark, and the osteogenic differentiation defects of CTR9 KD hMSCs can be partially rescued by treatment with EZH2 inhibitors. Transcriptome analyses identified bone morphology protein-2 (BMP-2) as a downstream effector of CTR9. BMP-2 secretion, membrane anchorage, and the BMP-SMAD pathway were impaired in CTR9 KD MSCs, and the effects were rescued by BMP-2 supplementation. This study uncovers an epigenetic mechanism engaging the CTR9-H3K27me3-BMP-2 axis to regulate the osteochondral lineage differentiation of hMSCs.

An eQTL variant of ALDH1A2 is associated with Kashin-Beck disease in Chinese population

INTRODUCTION

The aims of the study were to investigate the relationship between aldehyde dehydrogenase 1 family member A2 (ALDH1A2) and Kashin-Beck disease (KBD), explore the effects of the rs3204689 polymorphism and methylation status on the expression levels of ALDH1A2, and further clarify the pathogenesis of KBD.

MATERIALS AND METHODS

The genotype of ALDH1A2 rs3204689 was detected by PCR-RFLP in 103 KBD patients and 109 healthy controls in the whole blood. The mRNA level of ALDH1A2 was measured by qRT-PCR, and the protein expression was detected using IHC staining and Western blotting. The MSP-PCR was used to identify the ALDH1A2 methylation level.

RESULTS

There were significant differences in G/G, G/C, and C/C frequencies of ALDH1A2 rs3204689 between the KBD and control groups (χ² = 7.113, P = 0.029); the minor allele G of ALDH1A2 was associated with the risk of KBD (χ² = 5.984, P = 0.014). The mRNA and protein levels of ALDH1A2 were increased in the whole blood and cartilage of KBD patients compared with the controls (P = 0.049, P < 0.0001, P = 0.019). Meanwhile, a statistically significant difference was observed between G/G, G/C and C/C genotype on mRNA expression (P = 0.039). The methylation level of the ALDH1A2 gene promoter region showed no significant difference between the KBD and control groups (χ² = 0.317, P = 0.573).

CONCLUSION

Our case-control study indicates that the common variant rs3204689 near ALDH1A2 is associated with KBD in Chinese population. The risk allele G of rs3204689 is statistically linked to the high expression of ALDH1A2, which is up-regulated in the cartilage and whole blood of KBD patients. Our findings suggest a potential role of ALDH1A2 in the pathogenesis of KBD.

A network pharmacology approach to explore and validate the potential targets of ginsenoside on osteoporosis

Background: Osteoporosis (OP) is determined as a chronic systemic bone disorder to increase the susceptibility to fracture. Ginsenosides have been found the anti-osteoporotic activity of in vivo and in vitro. However, its mechanism remains unknown.Methods: The potential mechanism of ginsenosides in anti-osteoporotic activity was identified by using network phamacology analysis. The active compounds of ginsenosides and their targets associated to OP were retrieved from Traditional Chinese Medicine Systems Pharmacology Database and Analysis Platform, Drug Bank, Pharmmapper, and Cytoscape. The Gene ontology (GO) enrichment analysis and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis target genes were performed in String, Phenopedia, DisGeNET database, and Metascape software. The protein to protein interaction were created by String database and Cytoscape software. The molecular docking was used to investigate the interactions between active coumpounds and potential targets by utilizing SwissDock tool, UCSF Chimera, and Pymol software. Results: A total of eight important active ingredients and 17 potential targets related to OP treatment were subjected to analyze. GO analysis showed the anti-osteoporosis targets of ginsenoside mainly play a role in the response to steroid hormone. KEGG enrichment analysis indicated that ginsenoside treats OP by osteoblast differentiation signal pathway. Lastly, the molecular docking outcomes indicated that ginsenoside rh2 had a good binding ability with four target proteins IL1B, TNF, IFNG, and NFKBIA. Conclusion: IL1B, TNF, IFNG, and NFKBIA are the most important targets and osteoblast differentiation is the most valuable signaling pathways in ginsenoside for the treatment of OP, which might be beneficial to elucidate the mechanism concerned to the action of ginsenoside and might supply a better understanding of its anti-OP effects.

Shared genetic liability between major depressive disorder and osteoarthritis

Aims

Deciphering the genetic relationships between major depressive disorder (MDD) and osteoarthritis (OA) may facilitate an understanding of their biological mechanisms, as well as inform more effective treatment regimens. We aim to investigate the mechanisms underlying relationships between MDD and OA in the context of common genetic variations.

Methods

Linkage disequilibrium score regression was used to test the genetic correlation between MDD and OA. Polygenic analysis was performed to estimate shared genetic variations between the two diseases. Two-sample bidirectional Mendelian randomization analysis was used to investigate causal relationships between MDD and OA. Genomic loci shared between MDD and OA were identified using cross-trait meta-analysis. Fine-mapping of transcriptome-wide associations was used to prioritize putatively causal genes for the two diseases.

Results

MDD has a significant genetic correlation with OA (r_g = 0.29) and the two diseases share a considerable proportion of causal variants. Mendelian randomization analysis indicates that genetic liability to MDD has a causal effect on OA (b_xy = 0.24) and genetic liability to OA conferred a causal effect on MDD (b_xy = 0.20). Cross-trait meta-analyses identified 29 shared genomic loci between MDD and OA. Together with fine-mapping of transcriptome-wide association signals, our results suggest that Estrogen Receptor 1 (ESR1), SRY-Box Transcription Factor 5 (SOX5), and Glutathione Peroxidase 1 (GPX1) may have therapeutic implications for both MDD and OA.

Conclusion

The study reveals substantial shared genetic liability between MDD and OA, which may confer risk for one another. Our findings provide a novel insight into phenotypic relationships between MDD and OA.

Cite this article: Bone Joint Res 2022;11(1):12–22.

Integrative Analysis of Genome-Wide Association Studies and DNA Methylation Profile Identified Genetic Control Genes of DNA Methylation for Kashin-Beck Disease

OBJECTIVE

Epigenetic modifications of DNA are regarded as a crucial factor for understanding the molecular basis of complex phenotypes. This study aims to uncover insight into the epigenetic modifications for Kashin-Beck disease (KBD) by integrating genome-wide association studies (GWAS), methylation quantitative trait loci (meQTLs), and DNA methylation profiles data.

DESIGN

The knee articular cartilages of 5 KBD patients and 5 healthy controls were collected for DNA methylation profiling, using Illumina Infinium HumanMethylation450 BeadChip. Mass spectrograph validation of identified differently methylated genes was conducted using independent samples of 4 KBD patients and 3 healthy controls, together with a previous sample of 2743 Han Chinese individuals of GWAS study for KBD and a study of 697 normal subjects for meQTLs annotation datasets. KBD GWAS single nucleotide polymorphisms (SNPs) and normal meQTLs SNPs were integrated with DNA methylation profiles of KBD articular cartilage to identify genetic control (GC) genes of DNA methylation for KBD. Quantitative polymerase chain reaction (qPCR) was performed to validate the mRNA expression of several identified candidate genes.

RESULTS

A total of 162 CpG sites, 253 SNPs, and 123 GC genes for KBD were identified. Enrichment analysis detected 642 marked GO terms and 19 KEGG pathways (P < 0.05). Six potential key GC genes were conducted for qPCR experiment (ERG, MN1, MITF, WISP1, TRIO, and NOSTRIN).

CONCLUSIONS

The results suggest that GC genes of DNA methylation may lead to the erosion of cartilage in KBD, which may help us in understanding the epigenetic alteration of KBD.

Exome Sequencing Reveals a Phenotype Modifying Variant in ZNF528 in Primary Osteoporosis With a COL1A2 Deletion

We studied a family with severe primary osteoporosis carrying a heterozygous p.Arg8Phefs*14 deletion in COL1A2, leading to haploinsufficiency. Three affected individuals carried the mutation and presented nearly identical spinal fractures but lacked other typical features of either osteogenesis imperfecta or Ehlers-Danlos syndrome. Although mutations leading to haploinsufficiency in COL1A2 are rare, mutations in COL1A1 that lead to less protein typically result in a milder phenotype. We hypothesized that other genetic factors may contribute to the severe phenotype in this family. We performed whole-exome sequencing in five family members and identified in all three affected individuals a rare nonsense variant (c.1282C > T/p.Arg428*, rs150257846) in ZNF528. We studied the effect of the variant using qPCR and Western blot and its subcellular localization with immunofluorescence. Our results indicate production of a truncated ZNF528 protein that locates in the cell nucleus as per the wild-type protein. ChIP and RNA sequencing analyses on ZNF528 and ZNF528-c.1282C > T indicated that ZNF528 binding sites are linked to pathways and genes regulating bone morphology. Compared with the wild type, ZNF528-c.1282C > T showed a global shift in genomic binding profile and pathway enrichment, possibly contributing to the pathophysiology of primary osteoporosis. We identified five putative target genes for ZNF528 and showed that the expression of these genes is altered in patient cells. In conclusion, the variant leads to expression of truncated ZNF528 and a global change of its genomic occupancy, which in turn may lead to altered expression of target genes. ZNF528 is a novel candidate gene for bone disorders and may function as a transcriptional regulator in pathways affecting bone morphology and contribute to the phenotype of primary osteoporosis in this family together with the COL1A2 deletion. © 2020 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals LLC on behalf of American Society for Bone and Mineral Research (ASBMR).

The integrative analysis of DNA methylation and mRNA expression profiles confirmed the role of selenocompound metabolism pathway in Kashin-Beck disease

Kashin-Beck disease (KBD) is an endemic chronic osteochondropathy. The etiology of KBD remains unknown. In this study, we conducted an integrative analysis of genome-wide DNA methylation and mRNA expression profiles between KBD and normal controls to identify novel candidate genes and pathways for KBD. Articular cartilage samples from 17 grade III KBD patients and 17 healthy controls were used in this study. DNA methylation profiling of knee cartilage and mRNA expression profile data were obtained from our previous studies. InCroMAP was performed to integrative analysis of genome-wide DNA methylation profiles and mRNA expression profiles. Gene ontology (GO) enrichment analysis was conducted by online DAVID 6.7. The quantitative real-time polymerase chain reaction (qPCR), Western blot, immunohistochemistry (IHC), and lentiviral vector transfection were used to validate one of the identified pathways. We identified 298 common genes (such as COL4A1, HOXA13, TNFAIP6 and TGFBI), 36 GO terms (including collagen function, skeletal system development, growth factor), and 32 KEGG pathways associated with KBD (including Selenocompound metabolism pathway, PI3K-Akt signaling pathway, and TGF-beta signaling pathway). Our results suggest the dysfunction of many genes and pathways implicated in the pathogenesis of KBD, most importantly, both the integrative analysis and in vitro study in KBD cartilage highlighted the importance of selenocompound metabolism pathway in the pathogenesis of KBD for the first time.

A trans-ethnic two-stage polygenetic scoring analysis detects genetic correlation between osteoporosis and schizophrenia

BACKGROUNDS

To explore the genetic correlation between schizophrenia (SCZ) and osteoporosis (OP).

DESIGN, SETTING, PARTICIPANTS, MEASUREMENTS

We conducted a trans-ethnic two-stage genetic correlation analysis of OP and SCZ, totally invoking 2286 Caucasia subjects in discovery stage and 4124 Chinese subjects in replication stage. The bone mineral density (BMD) and bone area values of ulna & radius, hip and spine were measured using Hologic 4500W dual energy X-ray absorptiometry machine. SCZ was diagnosed according to DSM-IV criteria. For the genome-wide association study (GWAS) of Caucasian OP, Chinese OP and Chinese SCZ, SNP genotyping was performed using Affymetrix SNP 6.0 array. For the GWAS of Caucasian SCZ, SNP genotyping was conducted using the Affymetrix 5.0 array, Affymetrix 6.0 array and Illumina 550 K array. Polygenetic risk scoring (PRS) analysis was conducted by PRSice software. Also, Linkage disequilibrium score regression (LD Score regression) analysis was performed to evaluate the genetic correlation between OP and SCZ. Multi-trait analysis of GWAS (MTAG) was performed to detect novel candidate genes for osteoporosis and SCZ.

RESULTS

In the Caucasia discovery samples, significant genetic correlations were observed for ulna & radius BMD vs. SCZ (P value = 0.010), ulna & radius area vs. SCZ (P value = 0.031). In the Chinese replication samples, we observed significant correlation for ulna & radius area vs. SCZ (P value = 0.019). In addition, LD Score regression also identified significant genetic correlations between SCZ and bone phenotypes in Caucasian and Chinese sample respectively. MTAG analysis identified several novel candidate genes, such as CTNNA2 (MTAG P value = 2.24 × 10-6) for SCZ and FADS2 (MTAG P value = 2.66 × 10-7) for osteoporosis.

CONCLUSIONS

Our study results support the overlapped genetic basis for osteoporosis and SCZ, and provide novel clues for elucidating the biological mechanism of increased osteoporosis risk in SCZ patients.

Chemokines and Bone

Chemokines are a family of small proteins, subdivided by their conserved cysteine residues and common structural features. Chemokines interact with their cognate G-protein-coupled receptors to elicit downstream signals that result in cell migration, proliferation, and survival. This review presents evidence for how the various CXC and CC subfamily chemokines influence bone hemostasis by acting on osteoclasts, osteoblasts, and progenitor cells. Also discussed are the ways in which chemokines contribute to bone loss as a result of inflammatory diseases such as rheumatoid arthritis, HIV infection, and periodontal infection. Both positive and negative effects of chemokines on bone formation and bone loss are presented. In addition, the role of chemokines in altering the bone microenvironment through effects on angiogenesis and tumor invasion is discussed. Very few therapeutic agents that influence bone formation by targeting chemokines or chemokine receptors are available, although a few are currently being evaluated.

MicroRNA-214-5p/TGF-β/Smad2 signaling alters adipogenic differentiation of bone marrow stem cells in postmenopausal osteoporosis

Postmenopausal osteoporosis (OPM) is a common type of osteoporosis in females. It is a systemic, chronic bone disease that presents as microstructure degradation of osseous tissue, decreased bone mineral density and increased osteopsathyrosis caused by hypoovarianism and reduced estrogen levels in the body following menopause. In the present study, the role of microRNA (miR)-214-5p in the regulation of the expression of bone marrow stem cells (BMSCs) was investigated, and its molecular mechanism of osteogenic induction in vitro was assessed. When dexamethasone-induced adipogenic differentiation was performed, miR-214-5p expression was increased compared with the control group, as determined by RT-qPCR. Furthermore, oil red O staining, RT-qPCR and western blot analysis demonstrated that overexpression of miR-214-5p promoted adipogenic differentiation, inhibited alkaline phosphatase (ALP), runt-related transcription factor 2 (Runx2), osteocalcin (OC) and collagen α-1 (I) chain (COL1A1) mRNA expression, and suppressed transforming growth factor (TGF)-β, phosphorylated (p)-Smad2 and collagen type IV α1 chain (COL4A1) protein expression in BMSCs. Additionally, downregulation of miR-214-5p increased the ALP, Runx2, OC and COL1 mRNA expression and increased TGF-β, Smad2 and COL4A1 protein expression in BMSCs. Furthermore, a TGF-β inhibitor was employed to inhibit TGF-β expression in BMSCs following miR-214-5p downregulation, which led to reduced Smad2, TGF-β and COL4A1 protein expression, and ALP, Runx2, OC and COL1 mRNA expression was also reduced, compared with the miR-214-5p downregulation only group. It was demonstrated that miR-214-5p may weaken osteogenic differentiation of BMSCs through regulating COL4A1. In conclusion, the results of the present study indicated that miR-214-5p may promote the adipogenic differentiation of BMSCs through regulation of the TGF-β/Smad2/COL4A1 signaling pathway, and potentially may be used to develop a novel drug for postmenopausal osteoporosis.

Integrating genome-wide DNA methylation and mRNA expression profiles identified different molecular features between Kashin-Beck disease and primary osteoarthritis

Background

Kashin-Beck disease (KBD) is an endemic osteochondropathy of unknown etiology. Osteoarthritis (OA) is a form of degenerative joint disease sharing similar clinical manifestations and pathological changes to articular cartilage with KBD.

Methods

A genome-wide DNA methylation profile of articular cartilage from five KBD patients and five OA patients was first performed using the Illumina Infinium HumanMethylation450 BeadChip. Together with a previous gene expression profiling dataset comparing KBD cartilage with OA cartilage, an integrative pathway enrichment analysis of the genome-wide DNA methylation and the mRNA expression profiles conducted in articular cartilage was performed by InCroMAP software.

Results

We identified 241 common genes altered in both the DNA methylation profile and the mRNA expression profile of articular cartilage of KBD versus OA, including CHST13 (NM_152889, fold-change = 0.5979, P_methy = 0.0430), TGFBR1 (NM_004612, fold-change = 2.077, P_methy = 0.0430), TGFBR2 (NM_001024847, fold-change = 1.543, P_methy = 0.037), TGFBR3 (NM_001276, fold-change = 0.4515, P_methy = 6.04 × 10⁻⁴), and ADAM12 (NM_021641, fold-change = 1.9768, P_methy = 0.0178). Integrative pathway enrichment analysis identified 19 significant KEGG pathways, including mTOR signaling (P = 0.0301), glycosaminoglycan biosynthesis-chondroitin sulfate/dermatan sulfate (P = 0.0391), glycosaminoglycan biosynthesis-keratan sulfate (P = 0.0278), and PI3K-Akt signaling (P = 0.0243).

Conclusion

This study identified different molecular features between Kashin-Beck disease and primary osteoarthritis and provided novel clues for clarifying the pathogenetic differences between KBD and OA.

Electronic supplementary material

The online version of this article (10.1186/s13075-018-1531-1) contains supplementary material, which is available to authorized users.

P4 medicine and osteoporosis: a systematic review

BACKGROUND

Osteoporosis is the most frequent bone metabolic disease. In order to improve early detection, prediction, prevention, diagnosis, and treatment of the disease, a new model of P4 medicine (personalized, predictive, preventive, and participatory medicine) could be applied. The aim of this work was to systematically review the publications of four different types of "omics" studies related to osteoporosis, in order to discover novel predictive, preventive, diagnostic, and therapeutic targets for better management of the geriatric population.

METHODS

To systematically search the PubMed database, we created specific groups of criteria for four different types of "omics" information on osteoporosis: genomic, transcriptomic, proteomic, and metabolomic. We then analyzed the intersections between them in order to find correlations and common pathways or molecules with important roles in osteoporosis, and with a potential application in disease prediction, prevention, diagnosis, or treatment.

RESULTS

Altogether, 180 publications of "omics" studies in the field of osteoporosis were found and reviewed at first selection. After introducing the inclusion and exclusion criteria (the secondary selection), 46 papers were included in the systematic review.

CONCLUSIONS

The intersection of reviewed papers identified five genes (ESR1, IBSP, CTNNB1, SOX4, and IDUA) and processes like the Wnt pathway, JAK/STAT signaling, and ERK/MAPK, which should be further validated for their predictive, diagnostic, or other clinical value in osteoporosis. Such molecular insights will enable us to fit osteoporosis into the P4 strategy and could increase the effectiveness of disease prediction and prevention, with a decrease in morbidity in the geriatric population.