Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome [email protected] Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass

Identification of known and novel long non-coding RNAs potentially responsible for the effects of BMD GWAS loci

Osteoporosis, characterized by low bone mineral density (BMD), is the most common complex disease affecting bone and constitutes a major societal health problem. Genome-wide association studies (GWASs) have identified over 1100 associations influencing BMD. It has been shown that perturbations to long non-coding RNAs (lncRNAs) influence BMD and the activities of bone cells; however, the extent to which lncRNAs are involved in the genetic regulation of BMD is unknown. Here, we combined the analysis of allelic imbalance (AI) in human acetabular bone fragments with a transcriptome-wide association study (TWAS) and expression quantitative trait loci (eQTL) colocalization analysis using data from the Genotype-Tissue Expression (GTEx) project to identify lncRNAs potentially responsible for GWAS associations. We identified 27 lncRNAs in bone that are located in proximity to a BMD GWAS association and harbor SNPs demonstrating AI. Using GTEx data we identified an additional 31 lncRNAs whose expression was associated (FDR correction0.1). The 58 lncRNAs are located in 43 BMD associations. To further support a causal role for the identified lncRNAs, we show that 23 of the 58 lncRNAs are differentially expressed as a function of osteoblast differentiation. Our approach identifies lncRNAs that are potentially responsible for BMD GWAS associations and suggest that lncRNAs play a role in the genetics of osteoporosis.First author draf

Transcriptome-wide association study and eQTL colocalization identify potentially causal genes responsible for bone mineral density GWAS associations

Genome-wide association studies (GWASs) for bone mineral density (BMD) have identified over 1,100 associations to date. However, identifying causal genes implicated by such studies has been challenging. Recent advances in the development of transcriptome reference datasets and computational approaches such as transcriptome-wide association studies (TWASs) and expression quantitative trait loci (eQTL) colocalization have proven to be informative in identifying putatively causal genes underlying GWAS associations. Here, we used TWAS/eQTL colocalization in conjunction with transcriptomic data from the Genotype-Tissue Expression (GTEx) project to identify potentially causal genes for the largest BMD GWAS performed to date. Using this approach, we identified 512 genes as significant (Bonferroni <= 0.05) using both TWAS and eQTL colocalization. This set of genes was enriched for regulators of BMD and members of bone relevant biological processes. To investigate the significance of our findings, we selected PPP6R3, the gene with the strongest support from our analysis which was not previously implicated in the regulation of BMD, for further investigation. We observed that Ppp6r3 deletion in mice decreased BMD. In this work, we provide an updated resource of putatively causal BMD genes and demonstrate that PPP6R3 is a putatively causal BMD GWAS gene. These data increase our understanding of the genetics of BMD and provide further evidence for the utility of combined TWAS/colocalization approaches in untangling the genetics of complex traits.First author draf

Single-cell transcriptomics of bone marrow stromal cells in Diversity Outbred mice: a model for population-level scRNA-seq studies

Genome-wide association studies (GWASs) have advanced our understanding of the genetics of osteoporosis; however, the challenge has been converting associations to causal genes. Studies have utilized transcriptomics data to link disease-associated variants to genes, but few population transcriptomics datasets have been generated on bone at the single-cell level. To address this challenge, we profiled the transcriptomes of bone marrow-derived stromal cells (BMSCs) cultured under osteogenic conditions from five Diversity Outbred (DO) mice using single-cell RNA-seq (scRNA-seq). The goal of the study was to determine if BMSCs could serve as a model to generate cell-type specific transcriptomic profiles of mesenchymal lineage cells from large populations of mice to inform genetic studies. By enriching for mesenchymal lineage cells in vitro, coupled with pooling of multiple samples and downstream genotype deconvolution, we demonstrate the scalability of this model for population-level studies. We demonstrate that dissociation of BMSCs from a heavily mineralized matrix had little effect on viability or their transcriptomic signatures. Furthermore, we show that BMSCs cultured under osteogenic conditions are diverse and consist of cells with characteristics of mesenchymal progenitors, marrow adipogenic lineage precursors (MALPs), osteoblasts, osteocyte-like cells, and immune cells. Importantly, all cells were similar from a transcriptomic perspective to cells isolated in vivo. We employed scRNA-seq analytical tools to confirm the biological identity of profiled cell-types. SCENIC was used to reconstruct gene regulatory networks (GRNs) and we showed that cell-types show GRNs expected of osteogenic and pre-adipogenic lineage cells. Further, CELLECT analysis showed that osteoblasts, osteocyte-like cells, and MALPs captured a significant component of bone mineral density (BMD) heritability. Together, these data suggest that BMSCs cultured under osteogenic conditions coupled with scRNA-seq can be used as a scalable and biologically-informative model to generate cell-type specific transcriptomic profiles of mesenchymal lineage cells in large populations. This article is protected by copyright. All rights reserved

Osteoblasts Generate Testosterone From DHEA and Activate Androgen Signaling in Prostate Cancer Cells

Bone metastasis is a complication of prostate cancer in up to 90% of men afflicted with advanced disease. Therapies that reduce androgen exposure remain at the forefront of treatment. However, most prostate cancers transition to a state whereby reducing testicular androgen action becomes ineffective. A common mechanism of this transition is intratumoral production of testosterone (T) using the adrenal androgen precursor dehydroepiandrosterone (DHEA) through enzymatic conversion by 3β- and 17β- hydroxysteroid dehydrogenases (3βHSD and 17βHSD). Given the ability of prostate cancer to form blastic metastases in bone, we hypothesized that osteoblasts might be a source of androgen synthesis. RNA expression analyses of murine osteoblasts and human bone confirmed that at least one 3βHSD and 17βHSD enzyme isoform was expressed, suggesting that osteoblasts are capable of generating androgens from adrenal DHEA. Murine osteoblasts were treated with 100- nM and 1- μM DHEA or vehicle control. Conditioned media from these osteoblasts were assayed for intermediate and active androgens by liquid chromatography- tandem mass spectrometry. As DHEA was consumed, the androgen intermediates androstenediol and androstenedione were generated and subsequently converted to T. Conditioned media of DHEA- treated osteoblasts increased androgen receptor (AR) signaling, prostate- specific antigen (PSA) production, and cell numbers of the androgen- sensitive prostate cancer cell lines C4- 2B and LNCaP. DHEA did not induce AR signaling in osteoblasts despite AR expression in this cell type. We describe an unreported function of osteoblasts as a source of T that is especially relevant during androgen- responsive metastatic prostate cancer invasion into bone. © 2021 American Society for Bone and Mineral Research (ASBMR). This article has been contributed to by US Government employees and their work is in the public domain in the USA.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/169292/1/jbmr4313_am.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/169292/2/jbmr4313.pd

Osteoblasts Generate Testosterone From DHEA

Bone metastasis is a complication of prostate cancer in up to 90% of men afflicted with advanced disease. Therapies that reduce androgen exposure remain at the forefront of treatment. However, most prostate cancers transition to a state whereby reducing testicular androgen action becomes ineffective. A common mechanism of this transition is intratumoral production of testosterone (T) using the adrenal androgen precursor dehydroepiandrosterone (DHEA) through enzymatic conversion by 3β- and 17β- hydroxysteroid dehydrogenases (3βHSD and 17βHSD). Given the ability of prostate cancer to form blastic metastases in bone, we hypothesized that osteoblasts might be a source of androgen synthesis. RNA expression analyses of murine osteoblasts and human bone confirmed that at least one 3βHSD and 17βHSD enzyme isoform was expressed, suggesting that osteoblasts are capable of generating androgens from adrenal DHEA. Murine osteoblasts were treated with 100- nM and 1- μM DHEA or vehicle control. Conditioned media from these osteoblasts were assayed for intermediate and active androgens by liquid chromatography- tandem mass spectrometry. As DHEA was consumed, the androgen intermediates androstenediol and androstenedione were generated and subsequently converted to T. Conditioned media of DHEA- treated osteoblasts increased androgen receptor (AR) signaling, prostate- specific antigen (PSA) production, and cell numbers of the androgen- sensitive prostate cancer cell lines C4- 2B and LNCaP. DHEA did not induce AR signaling in osteoblasts despite AR expression in this cell type. We describe an unreported function of osteoblasts as a source of T that is especially relevant during androgen- responsive metastatic prostate cancer invasion into bone. © 2021 American Society for Bone and Mineral Research (ASBMR). This article has been contributed to by US Government employees and their work is in the public domain in the USA.Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/169292/1/jbmr4313_am.pdfhttp://deepblue.lib.umich.edu/bitstream/2027.42/169292/2/jbmr4313.pd

Mouse genome-wide association and systems genetics identifies Lhfp as a regulator of bone mass.

Bone mineral density (BMD) is a strong predictor of osteoporotic fracture. It is also one of the most heritable disease-associated quantitative traits. As a result, there has been considerable effort focused on dissecting its genetic basis. Here, we performed a genome-wide association study (GWAS) in a panel of inbred strains to identify associations influencing BMD. This analysis identified a significant (P = 3.1 x 10-12) BMD locus on Chromosome [email protected] Mbp that replicated in two separate inbred strain panels and overlapped a BMD quantitative trait locus (QTL) previously identified in a F2 intercross. The association mapped to a 300 Kbp region containing four genes; Gm2447, Gm20750, Cog6, and Lhfp. Further analysis found that Lipoma HMGIC Fusion Partner (Lhfp) was highly expressed in bone and osteoblasts. Furthermore, its expression was regulated by a local expression QTL (eQTL), which overlapped the BMD association. A co-expression network analysis revealed that Lhfp was strongly connected to genes involved in osteoblast differentiation. To directly evaluate its role in bone, Lhfp deficient mice (Lhfp-/-) were created using CRISPR/Cas9. Consistent with genetic and network predictions, bone marrow stromal cells (BMSCs) from Lhfp-/- mice displayed increased osteogenic differentiation. Lhfp-/- mice also had elevated BMD due to increased cortical bone mass. Lastly, we identified SNPs in human LHFP that were associated (P = 1.2 x 10-5) with heel BMD. In conclusion, we used GWAS and systems genetics to identify Lhfp as a regulator of osteoblast activity and bone mass