Current research efforts, although limited, clearly indicate that the gut microbiota may be implicated in bone metabolism, and therefore, further exploration of this relationship is a promising area of focus in bone health and osteoporosis research. Although most existing studies investigate this relationship using animal models, human studies are both needed and on the horizon.

Human osteoblasts are multifunctional bone cells, which play essential roles in bone formation, angiogenesis regulation, as well as maintenance of hematopoiesis. However, the categorization of primary osteoblast subtypes in vivo in humans has not yet been achieved. Here, we used single-cell RNA sequencing (scRNA-seq) to perform a systematic cellular taxonomy dissection of freshly isolated human osteoblasts from one 31-year-old male with osteoarthritis and osteopenia after hip replacement. Based on the gene expression patterns and cell lineage reconstruction, we identified three distinct cell clusters including preosteoblasts, mature osteoblasts, and an undetermined rare osteoblast subpopulation. This novel subtype was found to be the major source of the nuclear receptor subfamily 4 group A member 1 and 2 (NR4A1 and NR4A2) in primary osteoblasts, and the expression of NR4A1 was confirmed by immunofluorescence staining on mouse osteoblasts in vivo . Trajectory inference analysis suggested that the undetermined cluster, together with the preosteoblasts, are involved in the regulation of osteoblastogenesis and also give rise to mature osteoblasts. Investigation of the biological processes and signaling pathways enriched in each subpopulation revealed that in addition to bone formation, preosteoblasts and undetermined osteoblasts may also regulate both angiogenesis and hemopoiesis. Finally, we demonstrated that there are systematic differences between the transcriptional profiles of human and mouse osteoblasts, highlighting the necessity for studying bone physiological processes in humans rather than solely relying on mouse models. Our findings provide novel insights into the cellular heterogeneity and potential biological functions of human primary osteoblasts at the single-cell level.

Although GWAS have been successful in identifying some osteoporosis associated loci, the findings explain only a small fraction of the total genetic variance. In this study we use a recently developed novel pleiotropic conditional false discovery rate (cFDR) method to identify novel genetic loci associated with two risk traits for osteoporotic fracture (the clinical outcome and end result of osteoporosis), Height (HT) and Femoral Neck (FNK) BMD. The cFDR method allows us to improve the detection of associated variants by incorporating any potentially shared genetic mechanisms between the two associated traits. We analyzed the summary statistics from two GWAS meta-analyses for single nucleotide polymorphisms (SNPs) that are associated with HT and FNK BMD. Using the cFDR method, we show enrichment in the identification of SNPs associated with each trait conditioned on their strength of association with the second trait. The findings revealed 18 SNPs that are associated with both HT and FNK BMD, 4 of which had not previously been reported to play a role in bone health. The novel SNPs located at KIF1B and the intergenic region between FERD3L and TWISTNB are noteworthy as these genes may be associated with processes that are functionally important in bone metabolism. By leveraging GWAS results from related phenotypes we identified several novel loci that may contribute to the proportion of variability explained for each trait, although we cannot speculate about these potential contributions to heritability based on this analysis alone.

Bone marrow-derived mesenchymal stem cells (BM-MSCs) are multipotent stromal cells that have a critical role in the maintenance of skeletal tissues such as bone, cartilage, and the fat in bone marrow. In addition to providing microenvironmental support for hematopoietic processes, BM-MSCs can differentiate into various mesodermal lineages including osteoblast/osteocyte, chondrocyte, and adipocyte that are crucial for bone metabolism. While BM-MSCs have high cell-to-cell heterogeneity in gene expression, the cell subtypes that contribute to this heterogeneity in vivo in humans have not been characterized. To investigate the transcriptional diversity of BM-MSCs, we applied single-cell RNA sequencing (scRNA-seq) on freshly isolated CD271 + BM-derived mononuclear cells (BM-MNCs) from two human subjects. We successfully identified LEPR hi CD45 low BM-MSCs within the CD271 + BM-MNC population, and further codified the BM-MSCs into distinct subpopulations corresponding to the osteogenic, chondrogenic, and adipogenic differentiation trajectories, as well as terminal-stage quiescent cells. Biological functional annotations of the transcriptomes suggest that osteoblast precursors induce angiogenesis coupled with osteogenesis, and chondrocyte precursors have the potential to differentiate into myocytes. We also discovered transcripts for several clusters of differentiation (CD) markers that were either highly expressed (e.g., CD167b, CD91, CD130 and CD118) or absent (e.g., CD74, CD217, CD148 and CD68) in BM-MSCs, representing potential novel markers for human BM-MSC purification. This study is the first systematic in vivo dissection of human BM-MSCs cell subtypes at the single-cell resolution, revealing an insight into the extent of their cellular heterogeneity and roles in maintaining bone homeostasis.

No abstract

Background Clinical and epidemiological findings point to an association between type 2 diabetes (T2D) and low birth weight. However, the nature of the relationship is largely unknown. The aim of this study was to identify novel single nucleotide polymorphisms (SNPs) in T2D and birth weight, and their pleiotropic loci. Methods A pleiotropy-informed conditional false discovery rate (cFDR) method was applied to two independent genome-wide association studies (GWAS) summary statistics of T2D (n = 149 821) and birth weight (n = 26 836). Results A conditional Q–Q plot showed strong enrichment of genetic variants in T2D conditioned on different levels of association with birth weight. 133 T2D-associated SNPs, including 120 novel SNPs, were identified with a significance threshold of cFDR < 0.05; 13 significant birth weight-associated SNPs, including 12 novel SNPs (cFDR < 0.05) were identified. Conjunctional cFDR (ccFDR) analysis identified nine pleiotropic loci, including seven novel loci, shared by both T2D and birth weight (ccFDR < 0.05). Two novel SNPs located at the CDK5 regulatory subunit-associated protein 1-like 1 (CDKAL1; rs1012635; cFDR < 0.05) and adenylate cyclase 5 (ADCY5; rs4677887; cFDR < 0.05) genes are of note. These two genes increase the risk of T2D and low birth weight through the pathway of the “fetal insulin hypothesis.” Conclusion Several pleiotropic loci were identified between T2D and birth weight by leveraging GWAS results. The results make it possible to explain a greater proportion of trait heritability and improve our understanding of the shared pathophysiology between T2D and birth weight.

Our study provided evidence for shared genetic loci between T2D and AD in European subjects by using cFDR and ccFDR analyses. These results may provide novel insight into the etiology and potential therapeutic targets of T2D and/or AD.

Bone mineral density (BMD) is a complex trait with high missing heritability. Numerous evidences have shown that BMD variation has a relationship with coronary artery disease (CAD). This relationship may come from a common genetic basis called pleiotropy. By leveraging the pleiotropy with CAD, we may be able to improve the detection power of genetic variants associated with BMD. Using a recently developed conditional false discovery rate (cFDR) method, we jointly analyzed summary statistics from two large independent genome wide association studies (GWAS) of lumbar spine (LS) BMD and CAD. Strong pleiotropic enrichment and 7 pleiotropic SNPs were found for the two traits. We identified 41 SNPs for LS BMD (cFDR<0.05), of which 20 were replications of previous GWASs and 21 were potential novel SNPs that were not reported before. Four genes encompassed by 9 cFDR-significant SNPs were partially validated in the gene expression assay. Further functional enrichment analysis showed that genes corresponding to the cFDR-significant LS BMD SNPs were enriched in GO terms and KEGG pathways that played crucial roles in bone metabolism (adjP < 0.05). In protein-protein interaction analysis, strong interactions were found between the proteins produced by the corresponding genes. Our study demonstrated the reliability and high-efficiency of the cFDR method on the detection of trait-associated genetic variants, the present findings shed novel insights into the genetic variability of BMD as well as the shared genetic basis underlying osteoporosis and CAD.