Osteoporosis is a complex disease with a strong genetic contribution. A recently published genome-wide association study (GWAS) for estimated bone mineral density (eBMD) identified 1103 independent ...genome-wide significant association signals. Most of these variants are non-coding, suggesting that regulatory effects may drive many of the associations. To identify genes with a role in osteoporosis, we integrate the eBMD GWAS association results with those from our previous osteoclast expression quantitative trait locus (eQTL) dataset.
We identify sixty-nine significant cis-eQTL effects for eBMD GWAS variants after correction for multiple testing. We detect co-localisation of eBMD GWAS and osteoclast eQTL association signals for 21 of the 69 loci, implicating a number of genes including CCR5, ZBTB38, CPE, GNA12, RIPK3, IQGAP1 and FLCN. Summary-data-based Mendelian Randomisation analysis of the eBMD GWAS and osteoclast eQTL datasets identifies significant associations for 53 genes, with TULP4 presenting as a strong candidate for pleiotropic effects on eBMD and gene expression in osteoclasts. By performing analysis using the GARFIELD software, we demonstrate significant enrichment of osteoporosis risk variants among high-confidence osteoclast eQTL across multiple GWAS P value thresholds. Mice lacking one of the genes of interest, the apoptosis/necroptosis gene RIPK3, show disturbed bone micro-architecture and increased osteoclast number, highlighting a new biological pathway relevant to osteoporosis.
We utilise a unique osteoclast eQTL dataset to identify a number of potential effector genes for osteoporosis risk variants, which will help focus functional studies in this area.
Siglec‐15, a Siglec family member and type‐1 transmembrane protein, is expressed mainly in human macrophages and dendritic cells. It is comprised of a lysine‐containing transmembrane domain, two ...extracellular immunoglobulin (Ig)‐like domains and a short cytoplasmic domain. Siglec‐15 is highly conserved in vertebrates and acts as an immunoreceptor. It exerts diverse functions on osteoclast physiology as well as the tumor microenvironment. Siglec‐15 interacts with adapter protein DAP12 ‐ Syk signaling pathway to regulate the RANKL/RANK‐mediated PI3K, AKT, and ERK signaling pathways during osteoclast formation in vitro. Consistently, the lack of the Siglec‐15 gene in mice leads to impaired osteoclast activity and osteopetrosis in vivo. In addition, Siglec‐15 is expressed by tumor‐associated macrophages (TAMs) and regulates the tumor microenvironment by activating the SYK/MAPK signaling pathway. Interestingly, Siglec‐15 shares sequence homology to programmed death‐ligand 1 (PD‐L1) and has a potential immune‐regulatory role in cancer immunology. Thus, Siglec‐15 might also represent an alternative target for the treatment of cancers that do not respond to anti‐PD‐L1/PD‐1 immunotherapy. Understanding the role of Siglec‐15 in osteoclastogenesis and the tumor microenvironment will help us to develop new treatments for bone disorders and cancer.
Osteoporosis is a common bone disease that has a strong genetic component. Genome-wide linkage studies have identified the chromosomal region 3p14-p22 as a quantitative trait locus for bone mineral ...density (BMD). We have previously identified associations between variation in two related genes located in 3p14-p22, ARHGEF3 and RHOA, and BMD in women. In this study we performed knockdown of these genes using small interfering RNA (siRNA) in human osteoblast-like and osteoclast-like cells in culture, with subsequent microarray analysis to identify genes differentially regulated from a list of 264 candidate genes. Validation of selected findings was then carried out in additional human cell lines/cultures using quantitative real-time PCR (qRT-PCR). The qRT-PCR results showed significant down-regulation of the ACTA2 gene, encoding the cytoskeletal protein alpha 2 actin, in response to RHOA knockdown in both osteoblast-like (P<0.001) and osteoclast-like cells (P = 0.002). RHOA knockdown also caused up-regulation of the PTH1R gene, encoding the parathyroid hormone 1 receptor, in Saos-2 osteoblast-like cells (P<0.001). Other findings included down-regulation of the TNFRSF11B gene, encoding osteoprotegerin, in response to ARHGEF3 knockdown in the Saos-2 and hFOB 1.19 osteoblast-like cells (P = 0.003-0.02), and down-regulation of ARHGDIA, encoding the Rho GDP dissociation inhibitor alpha, in response to RHOA knockdown in osteoclast-like cells (P<0.001). These studies identify ARHGEF3 and RHOA as potential regulators of a number of genes in bone cells, including TNFRSF11B, ARHGDIA, PTH1R and ACTA2, with influences on the latter evident in both osteoblast-like and osteoclast-like cells. This adds further evidence to previous studies suggesting a role for the ARHGEF3 and RHOA genes in bone metabolism.
Osteoarthritis (OA) is a joint disease with a heritable component. Genetic loci identified via genome-wide association studies (GWAS) account for an estimated 26.3% of the disease trait variance in ...humans. Currently, there is no method for predicting the onset or progression of OA. We describe the first use of the Collaborative Cross (CC), a powerful genetic resource, to investigate knee OA in mice, with follow-up targeted multi-omics analysis of homologous regions of the human genome.
We histologically screened 275 mice for knee OA and conducted quantitative trait locus (QTL) mapping in the complete cohort (> 8 months) and the younger onset sub-cohort (8-12 months). Multi-omic analysis of human genetic datasets was conducted to investigate significant loci.
We observed a range of OA phenotypes. QTL mapping identified a genome-wide significant locus on mouse chromosome 19 containing Glis3, the human equivalent of which has been identified as associated with OA in recent GWAS. Mapping the younger onset sub-cohort identified a genome-wide significant locus on chromosome 17. Multi-omic analysis of the homologous region of the human genome (6p21.32) indicated the presence of pleiotropic effects on the expression of the HLA - DPB2 gene and knee OA development risk, potentially mediated through the effects on DNA methylation.
The significant associations at the 6p21.32 locus in human datasets highlight the value of the CC model of spontaneous OA that we have developed and lend support for an immune role in the disease. Our results in mice also add to the accumulating evidence of a role for Glis3 in OA.
Polycystic ovary syndrome (PCOS) is a complex multifactorial disorder with a substantial genetic component. However, the clinical manifestations of PCOS are heterogeneous with notable differences ...between lean and obese women, implying a different pathophysiology manifesting in differential body mass index (BMI). We performed a meta-analysis of genome-wide association study (GWAS) data from six well-characterised cohorts, using a case-control study design stratified by BMI, aiming to identify genetic variants associated with lean and overweight/obese PCOS subtypes.
The study comprised 254,588 women (5,937 cases and 248,651 controls) from individual studies performed in Australia, Estonia, Finland, the Netherlands and United States of America, and separated according to three BMI stratifications (lean, overweight and obese). Genome-wide association analyses were performed for each stratification within each cohort, with the data for each BMI group meta-analysed using METAL software. Almost half of the total study population (47%, n = 119,584) were of lean BMI (≤ 25 kg/m
). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10
) and rs2228260 within XBP1 (P = 3.68 × 10
). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10
). There were no significant loci observed for the overweight or obese sub-strata when analysed separately, however, when these strata were combined, an association signal led by rs569675099 within DENND1A reached genome-wide significance (P = 3.22 × 10
) and a gene-based association was identified with ERBB4 (P = 1.59 × 10
). Nineteen of 28 signals identified in previous GWAS, were replicated with consistent allelic effect in the lean stratum. There were less replicated signals in the overweight and obese groups, and only 4 SNPs were replicated in each of the three BMI strata.
Genetic variation at the XBP1, LINC02905 and ERBB4 loci were associated with PCOS within unique BMI strata, while DENND1A demonstrated associations across multiple strata, providing evidence of both distinct and shared genetic features between lean and overweight/obese PCOS-affected women. This study demonstrated that PCOS-affected women with contrasting body weight are not only phenotypically distinct but also show variation in genetic architecture; lean PCOS women typically display elevated gonadotrophin ratios, lower insulin resistance, higher androgen levels, including adrenal androgens, and more favourable lipid profiles. Overall, these findings add to the growing body of evidence supporting a genetic basis for PCOS as well as differences in genetic patterns relevant to PCOS BMI-subtype.
Osteoclasts are giant bone-digesting cells that harbor specialized lysosome-related organelles termed secretory lysosomes (SLs). SLs store cathepsin K and serve as a membrane precursor to the ruffled ...border, the osteoclast's 'resorptive apparatus'. Yet, the molecular composition and spatiotemporal organization of SLs remains incompletely understood. Here, using organelle-resolution proteomics, we identify member a2 of the solute carrier 37 family (Slc37a2) as a SL sugar transporter. We demonstrate in mice that Slc37a2 localizes to the SL limiting membrane and that these organelles adopt a hitherto unnoticed but dynamic tubular network in living osteoclasts that is required for bone digestion. Accordingly, mice lacking Slc37a2 accrue high bone mass owing to uncoupled bone metabolism and disturbances in SL export of monosaccharide sugars, a prerequisite for SL delivery to the bone-lining osteoclast plasma membrane. Thus, Slc37a2 is a physiological component of the osteoclast's unique secretory organelle and a potential therapeutic target for metabolic bone diseases.
Abstract
There has been a growing interest in the role of the subchondral bone and its resident osteoclasts in the progression of osteoarthritis (OA). A recent genome-wide association study (GWAS) ...identified 100 independent association signals for OA traits. Most of these signals are led by noncoding variants, suggesting that genetic regulatory effects may drive many of the associations. We have generated a unique human osteoclast-like cell-specific expression quantitative trait locus (eQTL) resource for studying the genetics of bone disease. Considering the potential role of osteoclasts in the pathogenesis of OA, we performed an integrative analysis of this dataset with the recently published OA GWAS results. Summary data-based Mendelian randomization (SMR) and colocalization analyses identified 38 genes with a potential role in OA, including some that have been implicated in Mendelian diseases with joint/skeletal abnormalities, such as BICRA, EIF6, CHST3, and FBN2. Several OA GWAS signals demonstrated colocalization with more than one eQTL peak, including at 19q13.32 (hip OA with BCAM, PRKD2, and BICRA eQTL). We also identified a number of eQTL signals colocalizing with more than one OA trait, including FAM53A, GCAT, HMGN1, MGAT4A, RRP7BP, and TRIOBP. An SMR analysis identified 3 loci with evidence of pleiotropic effects on OA-risk and gene expression: LINC01481, CPNE1, and EIF6. Both CPNE1 and EIF6 are located at 20q11.22, a locus harboring 2 other strong OA candidate genes, GDF5 and UQCC1, suggesting the presence of an OA-risk gene cluster. In summary, we have used our osteoclast-specific eQTL dataset to identify genes potentially involved with the pathogenesis of OA.
Osteoporosis is a disease that is characterised by reduced bone mineral density (BMD) and can be exacerbated by the excessive bone resorption of osteoclasts (OCs). Bioinformatic methods, including ...functional enrichment and network analysis, can provide information about the underlying molecular mechanisms that participate in the progression of osteoporosis. In this study, we harvested human OC-like cells differentiated in culture and their precursor peripheral blood mononuclear cells (PBMCs) and characterised the transcriptome of the two cell types using RNA-sequencing in order to identify differentially expressed genes. Differential gene expression analysis was performed in RStudio using the edgeR package. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed to identify enriched GO terms and signalling pathways, with inter-connected regions characterised using protein-protein interaction analysis. In this study, we identified 3201 differentially expressed genes using a 5% false discovery rate; 1834 genes were upregulated, whereas 1367 genes were downregulated. We confirmed a significant upregulation of several well-established OC genes including
,
,
,
,
, and
. The GO analysis suggested that upregulated genes are involved in cell division, cell migration, and cell adhesion, while the KEGG pathway analysis highlighted oxidative phosphorylation, glycolysis and gluconeogenesis, lysosome, and focal adhesion pathways. This study provides new information about changes in gene expression and highlights key biological pathways involved in osteoclastogenesis.
Osteoporosis is a common and debilitating bone disease that is characterised by low bone mineral density, typically assessed using dual-energy X-ray absorptiometry. Quantitative ultrasound (QUS), ...commonly utilising the two parameters velocity of sound (VOS) and broadband ultrasound attenuation (BUA), is an alternative technology used to assess bone properties at peripheral skeletal sites. The genetic influence on the bone qualities assessed by QUS remains an under-studied area. We performed a comprehensive genome-wide association study (GWAS) including low-frequency variants (minor allele frequency ≥0.005) for BUA and VOS using a discovery population of individuals with whole-genome sequence (WGS) data from the UK10K project (n = 1268). These results were then meta-analysed with those from two deeply imputed GWAS replication cohorts (n = 1610 and 13 749). In the gender-combined analysis, we identified eight loci associated with BUA and five with VOS at the genome-wide significance level, including three novel loci for BUA at 8p23.1 (PPP1R3B), 11q23.1 (LOC387810) and 22q11.21 (SEPT5) (P = 2.4 × 10-8 to 1.6 × 10-9). Gene-based association testing in the gender-combined dataset revealed eight loci associated with BUA and seven with VOS after correction for multiple testing, with one novel locus for BUA at FAM167A (8p23.1) (P = 1.4 × 10-6). An additional novel locus for BUA was seen in the male-specific analysis at DEFB103B (8p23.1) (P = 1.8 × 10-6). Fracture analysis revealed significant associations between variation at the WNT16 and RSPO3 loci and fracture risk (P = 0.004 and 4.0 × 10-4, respectively). In conclusion, by performing a large GWAS meta-analysis for QUS parameters of bone using a combination of WGS and deeply imputed genotype data, we have identified five novel genetic loci associated with BUA.