POLR2A blocks osteoclastic bone resorption and protects against osteoporosis by interacting with CREB1

Multi-modal molecular determinants of clinically relevant osteoporosis subtypes

Due to a rapidly aging global population, osteoporosis and the associated risk of bone fractures have become a wide-spread public health problem. However, osteoporosis is very heterogeneous, and the existing standard diagnostic measure is not sufficient to accurately identify all patients at risk of osteoporotic fractures and to guide therapy. Here, we constructed the first prospective multi-omics atlas of the largest osteoporosis cohort to date (longitudinal data from 366 participants at three time points), and also implemented an explainable data-intensive analysis framework (DLSF: Deep Latent Space Fusion) for an omnigenic model based on a multi-modal approach that can capture the multi-modal molecular signatures (M3S) as explicit functional representations of hidden genotypes. Accordingly, through DLSF, we identified two subtypes of the osteoporosis population in Chinese individuals with corresponding molecular phenotypes, i.e., clinical intervention relevant subtypes (CISs), in which bone mineral density benefits response to calcium supplements in 2-year follow-up samples. Many snpGenes associated with these molecular phenotypes reveal diverse candidate biological mechanisms underlying osteoporosis, with xQTL preferences of osteoporosis and its subtypes indicating an omnigenic effect on different biological domains. Finally, these two subtypes were found to have different relevance to prior fracture and different fracture risk according to 4-year follow-up data. Thus, in clinical application, M3S could help us further develop improved diagnostic and treatment strategies for osteoporosis and identify a new composite index for fracture prediction, which were remarkably validated in an independent cohort (166 participants).

A mutation in POLR2A gene associated with body size traits in Dezhou donkeys revealed with GWAS

The Dezhou donkey is a famous local donkey breed in China. The aim of the present study was to identify the genes associated with the body size traits of the Dezhou donkey and facilitate the breeding activities of the donkeys. A total of 349 donkeys from 2 generations (113 individuals in F0 and 236 in F1) were analyzed with restriction-site-associated DNA sequencing. A genome-wide association study revealed that the region between 13.7 and 15.6 Mb of chromosome 13 is significantly associated with body sizes. Candidate genes related to body size development, including POLR2A, CHRNB1, FGF11, and ZBTB4, were identified. The results of GO and KEGG analysis indicated that the genes involved in many GO terms were related to metabolic processes and developmental processes. Additionally, a T>C mutation (Chr13:14312485) was found at intron 10 of the POLR2A gene. The association analysis showed significant differences among genotypes for the size traits. The body size of the individuals with the TT genotype was significantly higher than that with the CC genotype. The results showed that the polymorphism of POLR2A has the potential to be used as a marker in the breeding programs of the Dezhou donkeys.

The CREB1 inhibitor 666-15 maintains cartilage homeostasis and mitigates osteoarthritis progression

Aims

cAMP response element binding protein (CREB1) is involved in the progression of osteoarthritis (OA). However, available findings about the role of CREB1 in OA are inconsistent. 666-15 is a potent and selective CREB1 inhibitor, but its role in OA is unclear. This study aimed to investigate the precise role of CREB1 in OA, and whether 666-15 exerts an anti-OA effect.

Methods

CREB1 activity and expression of a disintegrin and metalloproteinase with thrombospondin motifs 4 (ADAMTS4) in cells and tissues were measured by immunoblotting and immunohistochemical (IHC) staining. The effect of 666-15 on chondrocyte viability and apoptosis was examined by cell counting kit-8 (CCK-8) assay, JC-10, and terminal deoxynucleotidyl transferase-mediated dUTP nick end-labelling (TUNEL) staining. The effect of 666-15 on the microstructure of subchondral bone, and the synthesis and catabolism of cartilage, in anterior cruciate ligament transection mice were detected by micro-CT, safranin O and fast green (S/F), immunohistochemical staining, and enzyme-linked immunosorbent assay (ELISA).

Results

CREB1 was hyperactive in osteoarthritic articular cartilage, interleukin (IL)-1β-treated cartilage explants, and IL-1β- or carbonyl cyanide 3-chlorophenylhydrazone (CCCP)-treated chondrocytes. 666-15 enhanced cell viability of OA-like chondrocytes and alleviated IL-1β- or CCCP-induced chondrocyte injury through inhibition of mitochondrial dysfunction-associated apoptosis. Moreover, inhibition of CREB1 by 666-15 suppressed expression of ADAMTS4. Additionally, 666-15 alleviated joint degeneration in an ACLT mouse model.

Conclusion

Hyperactive CREB1 played a critical role in OA development, and 666-15 exerted anti-IL-1β or anti-CCCP effects in vitro as well as joint-protective effects in vivo. 666-15 may therefore be used as a promising anti-OA drug.

CREB1 alleviates the apoptosis and potentiates the osteogenic differentiation of zoledronic acid-treated human periodontal ligament stem cells via up-regulating VEGF

Periodontitis represents a severe inflammatory illness in tooth supporting tissue. It has been supported that cAMP response element binding protein 1 (CREB1), a common transcription factor, extensively participates in osteogenic differentiation. Here, the current study was to look into the impacts of CREB1 on the process of periodontitis and its possible action mechanism. After human periodontal ligament stem cells (PDLSCs) were challenged with zoledronic acid (ZA), CREB1 expression was examined with RT-qPCR and western blotting. CCK-8 assay appraised cell activity. Following CREB1 elevation or/and vascular endothelial growth factor (VEGF) silencing in ZA-treated PDLSCs, CCK-8 and TUNEL assays separately estimated cell viability and apoptosis. Western blotting tested the expression of apoptosis- and osteogenic differentiation-associated proteins. ALP staining measured PDLSCs osteogenic ability and ARS staining estimated mineralized nodule formation. JASPAR predicted the potential binding of CREB1 with VEGF promoter, which was then testified by ChIP and luciferase reporter assays. RT-qPCR and western blotting tested VEGF expression. CREB1 expression was declined in ZA-exposed PDLSCs and CREB1 elevation exacerbated the viability and osteogenic differentiation while obstructed the apoptosis of PDLSCs. Additionally, CREB1 bond to VEGF promoter and transcriptionally activated VEGF expression. Further, VEGF absence partially stimulated the apoptosis while suppressed the osteogenic differentiation of CREB1-overexpressing PDLSCs treated by ZA. To be concluded, CREB1 might activate VEGF transcription to obstruct the apoptosis while contribute to the osteogenic differentiation of ZA-treated PDLSCs.

Immune-related gene IL17RA as a diagnostic marker in osteoporosis

Objectives: Bone immune disorders are major contributors to osteoporosis development. This study aims to identify potential diagnostic markers and molecular targets for osteoporosis treatment from an immunological perspective. Method: We downloaded dataset GSE56116 from the Gene Expression Omnibus database, and identified differentially expressed genes (DEGs) between normal and osteoporosis groups. Subsequently, differentially expressed immune-related genes (DEIRGs) were identified, and a functional enrichment analysis was performed. A protein-protein interaction network was also constructed based on data from STRING database to identify hub genes. Following external validation using an additional dataset (GSE35959), effective biomarkers were confirmed using RT-qPCR and immunohistochemical (IHC) staining. ROC curves were constructed to validate the diagnostic values of the identified biomarkers. Finally, a ceRNA and a transcription factor network was constructed, and a Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analysis was performed to explore the biological functions of these diagnostic markers. Results: In total, 307 and 31 DEGs and DEIRGs were identified, respectively. The enrichment analysis revealed that the DEIRGs are mainly associated with Gene Ontology terms of positive regulation of MAPK cascade, granulocyte chemotaxis, and cytokine receptor. protein-protein interaction network analysis revealed 10 hub genes: FGF8, KL, CCL3, FGF4, IL9, FGF9, BMP7, IL17RA, IL12RB2, CD40LG. The expression level of IL17RA was also found to be significantly high. RT-qPCR and immunohistochemical results showed that the expression of IL17RA was significantly higher in osteoporosis patients compared to the normal group, as evidenced by the area under the curve Area Under Curve of 0.802. Then, we constructed NEAT1-hsa-miR-128-3p-IL17RA, and SNHG1-hsa-miR-128-3p-IL17RA ceRNA networks in addition to ERF-IL17RA, IRF8-IL17RA, POLR2A-IL17RA and ERG-IL17RA transcriptional networks. Finally, functional enrichment analysis revealed that IL17RA was involved in the development and progression of osteoporosis by regulating local immune and inflammatory processes in bone tissue. Conclusion: This study identifies the immune-related gene IL17RA as a diagnostic marker of osteoporosis from an immunological perspective, and provides insight into its biological function.

Whole-exome sequencing in syndromic craniosynostosis increases diagnostic yield and identifies candidate genes in osteogenic signaling pathways

Craniosynostosis (CS) is a common congenital anomaly defined by premature fusion of one or more cranial sutures. Syndromic CS involves additional organ anomalies or neurocognitive deficits and accounts for 25%-30% of the cases. In a recent population-based study by our group, 84% of the syndromic CS cases had a genetically verified diagnosis after targeted analyses. A number of different genetic causes were detected, confirming that syndromic CS is highly heterogeneous. In this study, we performed whole-exome sequencing of 10 children and parents from the same cohort where previous genetic results were negative. We detected pathogenic, or likely pathogenic, variants in four additional genes (NFIA, EXTL3, POLR2A, and FOXP2) associated with rare conditions. In two of these (POLR2A and FOXP2), CS has not previously been reported. We further detected a rare predicted damaging variant in SH3BP4, which has not previously been related to human disease. All findings were clustered in genes involved in the pathways of osteogenesis and suture patency. We conclude that whole-exome sequencing expands the list of genes associated with syndromic CS, and provides new candidate genes in osteogenic signaling pathways.

A multiethnic whole genome sequencing study to identify novel loci for bone mineral density

At present, there have only been a few DNA sequencing-based studies to explore the genetic determinants of bone mineral density (BMD). We carried out the largest whole genome sequencing analysis to date for femoral neck and spine BMD (n = 4981), with one of the highest average sequencing depths implemented thus far at 22×, in a multiethnic sample (58% Caucasian and 42% African American) from the Louisiana Osteoporosis Study (LOS). The LOS samples were combined with summary statistics from the GEFOS consortium and several independent samples of various ethnicities to perform GWAS meta-analysis (n = 44 506). We identified 31 and 30 genomic risk loci for femoral neck and spine BMD, respectively. The findings substantiate many previously reported susceptibility loci (e.g. WNT16 and ESR1) and reveal several others that are either novel or have not been widely replicated in GWAS for BMD, including two for femoral neck (IGF2 and ZNF423) and one for spine (SIPA1). Although we were not able to uncover ethnicity specific differences in the genetic determinants of BMD, we did identify several loci which demonstrated sex-specific associations, including two for women (PDE4D and PIGN) and three for men (TRAF3IP2, NFIB and LYSMD4). Gene-based rare variant association testing detected MAML2, a regulator of the Notch signaling pathway, which has not previously been suggested, for association with spine BMD. The findings provide novel insights into the pathophysiological mechanisms of osteoporosis.