Search

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 CTSK, DCSTAMP, ACP5, MMP9, ITGB3, and ATP6V0D2. 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.

Lung function impairment underlies chronic obstructive pulmonary disease (COPD) and predicts mortality. In the largest multi-ancestry GWAS meta-analysis of lung function to date, comprising 580,869 participants, 1020 independent association signals identified 559 genes supported by ≥2 criteria from a systematic variant-to-gene mapping framework. These genes were enriched in 29 pathways. Individual variants showed heterogeneity across ancestries, age and smoking groups, and collectively as a genetic risk score (GRS) showed strong association with COPD across ancestry groups. We undertook phenome-wide association studies (PheWAS) for selected associated variants, and trait and pathway-specific GRS to infer possible consequences of intervening in pathways underlying lung function. We highlight new putative causal variants, genes, proteins and pathways, including those targeted by existing drugs. These findings bring us closer to understanding the mechanisms underlying lung function and COPD, and should inform functional genomics experiments and potentially future COPD therapies.

Osteoclasts are giant bone-digesting cells that harbour 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 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 ruffled border. Thus, Slc37a2 is a physiological component of the osteoclast’s unique secretory organelle and a potential therapeutic target for metabolic bone diseases.

Background 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. Results 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. Conclusion 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.

Context In the clinic it is important to differentiate primary hyperparathyroidism (PHPT) from the more benign, inherited disorder, familial hypocalciuric hypercalcemia (FHH). Since the conditions may sometimes overlap biochemically, identification of calcium-sensing receptor (CASR) gene variants causative of FHH (but not PHPT) is the most decisive diagnostic aid. When novel variants are identified, bioinformatics and functional assessment are required to establish pathogenicity. Objective We identified 3 novel CASR transmembrane domain missense variants, Thr699Asn, Arg701Gly, and Thr808Pro, in 3 probands provisionally diagnosed with FHH and examined the variants using bioinformatics and functional analysis. Methods Bioinformatics assessment utilized wANNOVAR software. For functional characterization, each variant was cloned into a mammalian expression vector; wild-type and variant receptors were transfected into HEK293 cells, and their expression and cellular localization were assessed by Western blotting and confocal immunofluorescence, respectively. Receptor activation in HEK293 cells was determined using an IP-One ELISA assay following stimulation with Ca++ ions. Results Bioinformatics analysis of the variants was unable to definitively assign pathogenicity. Compared with wild-type receptor, all variants demonstrated impaired expression of mature receptor reaching the cell surface and diminished activation at physiologically relevant Ca++ concentrations. Conclusion Three CASR missense variants identified in probands provisionally diagnosed with FHH result in receptor inactivation and are therefore likely causative of FHH. Inactivation may be due to inadequate processing/trafficking of mature receptor and/or conformational changes induced by the variants affecting receptor signaling. This study demonstrates the value of functional studies in assessing genetic variants identified in hypercalcemic patients.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of genes important to bone biology in general, and particularly for identifying components unique to intramembranous ossification, which cannot be captured at other skeletal sites. We assessed genetic determinants of SK-BMD in 43,800 individuals, identifying 59 genome-wide significant loci (4 novel), explaining 12.5% of its variance. Pathway and enrichment analyses of the association signals resulted in clustering within gene-sets involved in regulating the development of the skeleton; overexpressed in the musculoskeletal system; and enriched in enhancer and transcribed regions in osteoblasts. From the four novel loci (mapping to ZIC1, PRKAR1A, ATP6V1C1, GLRX3), two (ZIC1 and PRKAR1A) have previously been related to craniofacial developmental defects. Functional validation of skull development in zebrafish revealed abnormal cranial bone initiation that culminated in ectopic sutures and reduced BMD in mutated zic1 and atp6v1c1 fish and asymmetric bone growth and elevated BMD in mutated prkar1a fish. We confirmed a role of ZIC1 loss-of-function in suture patterning and discovered ATP6V1C1 gene associated with suture development. In light of the evidence presented suggesting that SK-BMD is genetically related to craniofacial abnormalities, our study opens new avenues to the understanding of the pathophysiology of craniofacial defects and towards the effective pharmacological treatment of bone diseases.

To examine differences in the genetic architecture of BMI between tall and short people, we conducted genome-wide and follow-up analyses using UK Biobank data. We identify 57 loci as height-specific, detect differences in SNP-based heritability between tall and short people and show how genetic correlations between the two rises during the lifespan. Using phenome-wide analyses (PHEWAS), a significant association between a short people-specific locus on MC4R and energy portion size was detected. We identify one locus (GPC5-GPC6) with different effect directions on BMI in short and tall people. PHEWAS indicates this locus is associated with bone mineral density. Transcriptome-wide analyses hint that genes differentially associated with BMI in short vs tall people are enriched in brain tissue. Our findings highlight the role of height in the genetic underpinnings of BMI, provide biological insight into mechanisms underlying height-dependent differences in BMI and show that in short and tall people obesity is a risk factor that differentially increases susceptibility for disease.

The bone marrow microenvironment (BMM) plays a key role in leukemia progression, but its molecular complexity in pre-B cell acute lymphoblastic leukemia (B-ALL), the most common cancer in children, remains poorly understood. To gain further insight, we used single-cell RNA sequencing to characterize the kinetics of the murine BMM during B-ALL progression. Normal pro- and pre-B cells were found to be the most affected at the earliest stages of disease and this was associated with changes in expression of genes regulated by the AP1-transcription factor complex and regulatory factors NELFE, MYC and BCL11A. Granulocyte–macrophage progenitors show reduced expression of the tumor suppressor long non-coding RNA Neat1 and disruptions in the rate of transcription. Intercellular communication networks revealed monocyte-dendritic precursors to be consistently active during B-ALL progression, with enriched processes including cytokine-mediated signaling pathway, neutrophil-mediated immunity and regulation of cell migration and proliferation. In addition, we confirmed that the hematopoietic stem and progenitor cell compartment was perturbed during leukemogenesis. These findings extend our understanding of the complexity of changes and molecular interactions among the normal cells of the BMM during B-ALL progression.

Osteosarcoma (OS) is a primary malignant bone tumour which usually occurs in children and adolescents. OS is primarily a result of chromosomal aberrations, a combination of acquired genetic changes and, hereditary, resulting in the dysregulation of cellular functions. The Hippo signalling pathway regulates cell and tissue growth by modulating cell proliferation, differentiation, and migration in developing organs. Mammalian STE20-like 1/2 (MST1/2) protein kinases are activated by neurofibromatosis type 2, Ras association domain family member 2, kidney and brain protein, or other factors. Interactions between MST1/2 and salvador family WW domain-containing protein 1 activate large tumour suppressor kinase 1/2 proteins, which in turn phosphorylate the downstream Yes-associated protein 1/transcriptional coactivator with PDZ-binding motif (YAP/TAZ). Moreover, dysregulation of this pathway can lead to aberrant cell growth, resulting in tumorigenesis. Interestingly, small molecules targeting the Hippo signalling pathways, through affecting YAP/TAZ cellular localisation and their interaction with members of the TEA/ATTS domain family of transcriptional enhancers are being developed and hold promise for the treatment of OS. This review discusses the existing knowledge about the involvement of the Hippo signalling cascade in OS and highlights several small molecule inhibitors as potential novel therapeutics.

© 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc. We aimed to report the first genomewide association study (GWAS) meta-analysis of dual-energy X-ray absorptiometry (DXA)-derived hip shape, which is thought to be related to the risk of both hip osteoarthritis and hip fracture. Ten hip shape modes (HSMs) were derived by statistical shape modeling using SHAPE software, from hip DXA scans in the Avon Longitudinal Study of Parents and Children (ALSPAC; adult females), TwinsUK (mixed sex), Framingham Osteoporosis Study (FOS; mixed), Osteoporotic Fractures in Men study (MrOS), and Study of Osteoporotic Fractures (SOF; females) (total N = 15,934). Associations were adjusted for age, sex, and ancestry. Five genomewide significant (p < 5 × 10−9, adjusted for 10 independent outcomes) single-nucleotide polymorphisms (SNPs) were associated with HSM1, and three SNPs with HSM2. One SNP, in high linkage disequilibrium with rs2158915 associated with HSM1, was associated with HSM5 at genomewide significance. In a look-up of previous GWASs, three of the identified SNPs were associated with hip osteoarthritis, one with hip fracture, and five with height. Seven SNPs were within 200 kb of genes involved in endochondral bone formation, namely SOX9, PTHrP, RUNX1, NKX3-2, FGFR4, DICER1, and HHIP. The SNP adjacent to DICER1 also showed osteoblast cis-regulatory activity of GSC, in which mutations have previously been reported to cause hip dysplasia. For three of the lead SNPs, SNPs in high LD (r2 > 0.5) were identified, which intersected with open chromatin sites as detected by ATAC-seq performed on embryonic mouse proximal femora. In conclusion, we identified eight SNPs independently associated with hip shape, most of which were associated with height and/or mapped close to endochondral bone formation genes, consistent with a contribution of processes involved in limb growth to hip shape and pathological sequelae. These findings raise the possibility that genetic studies of hip shape might help in understanding potential pathways involved in hip osteoarthritis and hip fracture. © 2018 The Authors. Journal of Bone and Mineral Research Published by Wiley Periodicals, Inc.

Publisher's version (útgefin grein), Polycystic ovary syndrome (PCOS) is a disorder characterized by hyperandrogenism, ovulatory dysfunction and polycystic ovarian morphology. Affected women frequently have metabolic disturbances including insulin resistance and dysregulation of glucose homeostasis. PCOS is diagnosed with two different sets of diagnostic criteria, resulting in a phenotypic spectrum of PCOS cases. The genetic similarities between cases diagnosed based on the two criteria have been largely unknown. Previous studies in Chinese and European subjects have identified 16 loci associated with risk of PCOS. We report a fixed-effect, inverse-weighted-variance meta-analysis from 10,074 PCOS cases and 103,164 controls of European ancestry and characterisation of PCOS related traits. We identified 3 novel loci (near PLGRKT, ZBTB16 and MAPRE1), and provide replication of 11 previously reported loci. Only one locus differed significantly in its association by diagnostic criteria; otherwise the genetic architecture was similar between PCOS diagnosed by self-report and PCOS diagnosed by NIH or non-NIH Rotterdam criteria across common variants at 13 loci. Identified variants were associated with hyperandrogenism, gonadotropin regulation and testosterone levels in affected women. Linkage disequilibrium score regression analysis revealed genetic correlations with obesity, fasting insulin, type 2 diabetes, lipid levels and coronary artery disease, indicating shared genetic architecture between metabolic traits and PCOS. Mendelian randomization analyses suggested variants associated with body mass index, fasting insulin, menopause timing, depression and male-pattern balding play a causal role in PCOS. The data thus demonstrate 3 novel loci associated with PCOS and similar genetic architecture for all diagnostic criteria. The data also provide the first genetic evidence for a male phenotype for PCOS and a causal link to depression, a previously hypothesized comorbid disease. Thus, the genetics provide a comprehensive view of PCOS that encompasses multiple diagnostic criteria, gender, reproductive potential and mental health., This work has been supported by MRC grant MC_U106179472 (FD, KO, JRBP), Samuel Oschin Comprehensive Cancer Institute Developmental Funds, Center for Bioinformatics and Functional Genomics and Department of Biomedical Sciences Developmental Funds (MRJ), NCI P30CA177558 (CH), NCI UM1CA186107 (PK), European Regional Development Fund (Project No. 2014-2020.4.01.15-0012) and the European Union’s Horizon 2020 research and innovation program under grant agreements No 692065 (TL, RM, AS) and 692145 (RM), NICHD R01HD065029 (RS), Estonian Ministry of Education and Research (grant IUT34-16 to TL), NICHD R01HD057450 (MU), NICHD P50HD044405 (AD), NICHD R01HD057223 (AD), R01HD085227 (MGH, AD), deCode Genetics (GT, UT, KS, US), Raine Medical Research Foundation Priming Grant (BHM), SCGOPHCG RAC 2015-16/034 (SGW, BGAS), 2016-17/018 (BGAS), NIHR BRC, Wellcome Trust, MRC (TDS), Eris M. Field Chair in Diabetes Research (MOG), NIDDK P30 DK063491 (MOG), NIDDK U01DK094431, U01DK048381 (DE), NICHD U10HD38992 (RL), Estonian Ministry of Education and Research (grant IUT34-16), Enterprise Estonia (grant EU48695); the EU-FP7 Marie Curie Industry-Academia Partnerships and Pathways (IAPP, grant SARM, EU324509 to AS), Wellcome (090532, 098381, 203141); European Commission (ENGAGE: HEALTH-F4-2007-201413 to MIM), MRC G0802782, MR/M012638/1 (SF), Li Ka Shing Foundation, WT-SSI/John Fell Funds, NIHR Biomedical Research Centre, Oxford, Widenlife and NICHD 5P50HD028138-27 (CML), NICHD R01HD065029, ADA 1-10-CT-57, Harvard Clinical and Translational Science Center, from the National Center for Research Resources 1UL1 RR025758 (CKW). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

Osteoporosis is a complex disease with a strong genetic component. Genomewide association studies (GWAS) have been very successful at identifying common genetic variants associated with bone parameters. A recently published study documented the results of the largest GWAS for bone mineral density (BMD) performed to date (n = 142,487), identifying 307 conditionally independent single-nucleotide polymorphisms (SNPs) as associated with estimated BMD (eBMD) at the genomewide significance level. The vast majority of these variants are non-coding SNPs. Expression quantitative trait locus (eQTL) studies using disease-specific cell types have increasingly been integrated with the results from GWAS to identify genes through which the observed GWAS associations are likely mediated. We generated a unique human osteoclast-specific eQTL data set using cells differentiated in vitro from 158 participants. We then used this resource to characterize the 307 recently identified BMD GWAS SNPs for association with nearby genes (±500 kb). After correction for multiple testing, 24 variants were found to be significantly associated with the expression of 32 genes in the osteoclast-like cells. Bioinformatics analysis suggested that these variants and those in strong linkage disequilibrium with them are enriched in regulatory regions. Several of the eQTL associations identified are relevant to genes that present strongly as having a role in bone, particularly IQGAP1, CYP19A1, CTNNB1, and COL6A3. Supporting evidence for many of the associations was obtained from publicly available eQTL data sets. We have also generated strong evidence for the presence of a regulatory region on chromosome 15q21.2 relevant to both the GLDN and CYP19A1 genes. In conclusion, we have generated a unique osteoclast-specific eQTL resource and have used this to identify 32 eQTL associations for recently identified BMD GWAS loci, which should inform functional studies of osteoclast biology. © 2018 American Society for Bone and Mineral Research.

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.

Osteoporosis is a complex condition with contributions from, and interactions between, multiple genetic loci and environmental factors. This review summarizes key advances in the application of genetic approaches for the identification of osteoporosis susceptibility genes. Genome-wide linkage analysis (GWLA) is the classical approach for identification of genes that cause monogenic diseases; however, it has shown limited success for complex diseases like osteoporosis. In contrast, genome-wide association studies (GWAS) have successfully identified over 200 osteoporosis susceptibility loci with genome-wide significance, and have provided most of the candidate genes identified to date. Phenome-wide association studies (PheWAS) apply a phenotype-to-genotype approach which can be used to complement GWAS. PheWAS is capable of characterizing the association between osteoporosis and uncommon and rare genetic variants. Another alternative approach, whole genome sequencing (WGS), will enable the discovery of uncommon and rare genetic variants in osteoporosis. Meta-analysis with increasing statistical power can offer greater confidence in gene searching through the analysis of combined results across genetic studies. Recently, new approaches to gene discovery include animal phenotype based models such as the Collaborative Cross and ENU mutagenesis. Site-directed mutagenesis and genome editing tools such as CRISPR/Cas9, TALENs and ZNFs have been used in functional analysis of candidate genes in vitro and in vivo. These resources are revolutionizing the identification of osteoporosis susceptibility genes through the use of genetically defined inbred mouse libraries, which are screened for bone phenotypes that are then correlated with known genetic variation. Identification of osteoporosis-related susceptibility genes by genetic approaches enables further characterization of gene function in animal models, with the ultimate aim being the identification of novel therapeutic targets for osteoporosis.

Paget’s disease of bone (PDB) is characterised by focal abnormalities of bone remodelling, with increased osteoclastic resorption the primary feature of the disease. Genetic factors have been shown to play an important role in PDB, and genome-wide association studies (GWAS) have identified 7 genetic loci as associated with PDB at the genome-wide level. Expression quantitative trait locus (eQTL) studies using cell types that are directly relevant to the disease of interest are increasingly being used to identify putative effector genes for GWAS loci. We have recently constructed a unique osteoclast-specific eQTL resource using cells differentiated in vitro from 158 subjects for study of the genetics of bone disease. Considering the major role osteoclasts have in PDB, we used this resource to investigate potential genetic regulatory effects for the 7 PDB genome-wide significant loci on genes located within 500 kb of each locus. After correction for multiple testing, we observed statistically significant associations for rs4294134 with expression of the gene STMP1, and rs2458413 with expression of the genes DPYS and DCSTAMP. The eQTL associations observed for rs4294134 with STMP1, and rs2458413 with DCSTAMP were further supported by eQTL data from other tissue types. The product of the STMP1 gene has not been extensively studied, however the DCSTAMP gene has an established role in osteoclast differentiation and the associations seen between rs2458413 and PDB are likely mediated through regulatory effects on this gene. This study highlights the value of eQTL data in determining which genes are relevant to GWAS loci.

Bone mineral density (BMD) assessed by DXA is used to evaluate bone health. In children, total body (TB) measurements are commonly used; in older individuals, BMD at the lumbar spine (LS) and femoral neck (FN) is used to diagnose osteoporosis. To date, genetic variants in more than 60 loci have been identified as associated with BMD. To investigate the genetic determinants of TB-BMD variation along the life course and test for age-specific effects, we performed a meta-analysis of 30 genome-wide association studies (GWASs) of TB-BMD including 66,628 individuals overall and divided across five age strata, each spanning 15 years. We identified variants associated with TB-BMD at 80 loci, of which 36 have not been previously identified; overall, they explain approximately 10% of the TB-BMD variance when combining all age groups and influence the risk of fracture. Pathway and enrichment analysis of the association signals showed clustering within gene sets implicated in the regulation of cell growth and SMAD proteins, overexpressed in the musculoskeletal system, and enriched in enhancer and promoter regions. These findings reveal TB-BMD as a relevant trait for genetic studies of osteoporosis, enabling the identification of variants and pathways influencing different bone compartments. Only variants in ESR1 and close proximity to RANKL showed a clear effect dependency on age. This most likely indicates that the majority of genetic variants identified influence BMD early in life and that their effect can be captured throughout the life course.

Osteoporosis is a complex disease with a strong genetic component. Genomewide association studies (GWAS) have been very successful at identifying common genetic variants associated with bone parameters. A recently published study documented the results of the largest GWAS for bone mineral density (BMD) performed to date (n = 142,487), identifying 307 conditionally independent single-nucleotide polymorphisms (SNPs) as associated with estimated BMD (eBMD) at the genomewide significance level. The vast majority of these variants are non-coding SNPs. Expression quantitative trait locus (eQTL) studies using disease-specific cell types have increasingly been integrated with the results from GWAS to identify genes through which the observed GWAS associations are likely mediated. We generated a unique human osteoclast-specific eQTL data set using cells differentiated in vitro from 158 participants. We then used this resource to characterize the 307 recently identified BMD GWAS SNPs for association with nearby genes (±500 kb). After correction for multiple testing, 24 variants were found to be significantly associated with the expression of 32 genes in the osteoclast-like cells. Bioinformatics analysis suggested that these variants and those in strong linkage disequilibrium with them are enriched in regulatory regions. Several of the eQTL associations identified are relevant to genes that present strongly as having a role in bone, particularly IQGAP1, CYP19A1, CTNNB1, and COL6A3. Supporting evidence for many of the associations was obtained from publicly available eQTL data sets. We have also generated strong evidence for the presence of a regulatory region on chromosome 15q21.2 relevant to both the GLDN and CYP19A1 genes. In conclusion, we have generated a unique osteoclast-specific eQTL resource and have used this to identify 32 eQTL associations for recently identified BMD GWAS loci, which should inform functional studies of osteoclast biology. © 2018 American Society for Bone and Mineral Research.

Osteoporosis is a common and debilitating bone disease that is characterised by a low bone mineral density (BMD), a highly heritable trait. Genome-wide association studies (GWAS) have proven to be very successful in identifying common genetic variants associated with BMD adjusted for age, gender and weight, however a large portion of the genetic variance for this trait remains unexplained. There is evidence to suggest significant genetic correlation between body size traits and BMD. It has also recently been suggested that unintended bias can be introduced as a result of adjusting a phenotype for a correlated trait. We performed a GWAS meta-analysis in two populations (total n = 6,696) using BMD data adjusted for only age and gender, in an attempt to identify genetic variants associated with BMD including those that may have potential pleiotropic effects on BMD and body size traits. We observed a single variant, rs2566752, associated with spine BMD at the genome-wide significance level in the meta-analysis (P = 3.36 × 10 −09 ). Logistic regression analysis also revealed an association between rs2566752 and fracture rate in one of our study cohorts (P = 0.017, n = 5,654). This is an intronic variant located in the wntless Wnt ligand secretion mediator (WLS) gene (1p31.3), a known BMD locus which encodes an integral component of the Wnt ligand secretion pathway. Bioinformatics analyses of variants in moderate LD with rs2566752 produced strong evidence for a regulatory role for the variants rs72670452, rs17130567 and rs1430738. Expression quantitative trait locus (eQTL) analysis suggested that the variants rs12568456 and rs17130567 are associated with expression of the WLS gene in whole blood, cerebellum and temporal cortex brain tissue (P = 0.034–1.19 × 10 −23 ). Gene-wide association testing using the VErsatile Gene-based Association Study 2 (VEGAS2) software revealed associations between the coiled-coil domain containing 170 (CCDC170) gene, located adjacent to the oestrogen receptor 1 (ESR1) gene, and BMD at the spine, femoral neck and total hip sites (P = 1.0 × 10 −06 , 2.0 × 10 −06 and 2.0 × 10 −06 respectively). Genetic variation at the WLS and CCDC170/ESR1 loci were found to be significantly associated with BMD adjusted for only age and gender at the genome-wide level in this meta-analysis.

17β-Estradiol is important in maintaining bone structure, and regulation of its synthesis plays an important role in the development of postmenopausal osteoporosis. We and others have demonstrated associations between variation in the CYP19A1 gene (encoding aromatase) and areal bone mineral density (aBMD) phenotypes in women. In the present study 33 tag polymorphisms were genotyped across the CYP19A1 gene in a population of 1,185 Caucasian postmenopausal women to test the association between sequence variations, total DXA hip aBMD, and circulating 17β-estradiol levels. An in silico bioinformatics analysis was performed for single nucleotide polymorphisms (SNPs) associated with aBMD to identify putative functional effects, while linkage disequilibrium analysis of these SNPs was undertaken with previously published sequence variants. Five SNPs located in the central third of the gene were strongly associated with total-hip aBMD after adjustment for age (P = 0.006-0.013). A haplotype analysis of these five SNPs revealed an association between the haplotype C-G-G-G-C and increased aBMD (P = 0.008) and the haplotype A-A-A-A-A and a decreased aBMD (P = 0.021). The haplotype frequency was 9.0% for C-G-G-G-C and 15.4% for A-A-A-A-A, with the variation in mean total-hip aBMD explained by the haplotype analyses being 5% and 7%, respectively. None of these polymorphisms was significantly associated with circulating 17β-estradiol levels. In conclusion, common genetic variations within the CYP19A1 gene are significantly associated with aBMD in postmenopausal Caucasian women.

Previously, homozygous deletion of the UGT2B17 gene has shown association with hip fracture. Using a high-throughput qRT-PCR assay, we genotyped UGT2B17 copy number variation (CNV) in 1,347 elderly Caucasian women and examined for effects on bone phenotypes. We found no evidence of association between UGT2B17 CNV and osteoporosis risk in this population.Genetic studies of osteoporosis commonly examine SNPs in candidate genes or whole genome analyses, but insertions and deletions of DNA, collectively called CNV, also comprise a large amount of the genetic variability between individuals. Previously, homozygous deletion of the UGT2B17 gene in CNV 4q13.2, which encodes an enzyme that mediates the glucuronidation of steroid hormones, has shown association with the risk of hip fracture.We used a quantitative real-time PCR assay for genotyping the UGT2B17 CNV in a well-characterized population study of 1,347 Caucasian women aged 75.2 ± 2.7 years (mean ± SD), to assess the effect of the CNV on bone mass density (BMD) at the total hip site and osteoporosis risk.The UGT2B17 CNV distribution was consistent with the expected Hardy-Weinberg distribution and not different from frequencies previously reported in a Caucasian population. Data from ANCOVA of age- and weight-adjusted BMD for UGT2B17 CNV genotype showed no significant difference between genotype groups. Individuals with homozygous or heterozygous deletion of the UGT2B17 gene showed no increased risk of incident fragility fracture.These data suggest that quantitative real-time PCR is a rapid and efficient technique for determination of candidate CNVs, including the UGT2B17 CNV; however, we found no evidence of an effect of UGT2B17 CNV on osteoporosis risk in elderly Caucasian women.

Neuritin 1 (NRN1), an activity-regulated gene with multiple roles in neurodevelopment and synaptic plasticity, is located within the 6p24-p25 interval on chromosome 6, previously identified as linked to a subtype of schizophrenia (SZ) characterized by pervasive cognitive deficit (CD). We have tested the effect of NRN1 sequence variation on susceptibility to SZ and on general cognitive ability in patients and non-psychiatric control subjects by re-sequencing the coding regions of NRN1 and its flanking sequences, and genotyping 19 single-nucleotide polymorphisms (SNPs) in 336 SZ patients and 172 healthy control individuals. All participants completed comprehensive neurocognitive assessment, including tests estimating premorbid/prior IQ and current IQ. Logistic regression analyses found no significant association for any of the 19 SNPs with SZ or its CD subtype. However, linear regression analysis gave significant association (P = 0.024 and P = 0.027 after correction for multiple testing) for polymorphisms rs1475157 and rs9405890 with current IQ in the patient group. In SZ, the rs1475157-rs9405890 haplotypes revealed a highly significant association with the abstraction component of current ("fluid") intelligence (P = 0.0014), and with percentage loss of IQ points between premorbid and current intelligence (P = 0.0041). Results in the control group were not significant after correction. This is the first analysis of association between variation in NRN1 and SZ. The findings suggest a role of NRN1 as a modifier of cognitive functioning in SZ, with implications for future research into the impact of the environment on the development and maintenance of "fluid" intelligence.

Summary Background Osteoporosis is diagnosed by the measurement of bone mineral density, which is a highly heritable and multifactorial trait. We aimed to identify genetic loci that are associated with bone mineral density. Methods In this genome-wide association study, we identified the most promising of 314 075 single nucleotide polymorphisms (SNPs) in 2094 women in a UK study. We then tested these SNPs for replication in 6463 people from three other cohorts in western Europe. We also investigated allelic expression in lymphoblast cell lines. We tested the association between the replicated SNPs and osteoporotic fractures with data from two studies. Findings We identified genome-wide evidence for an association between bone mineral density and two SNPs (p −8 ). The SNPs were rs4355801, on chromosome 8, near to the TNFRSF11B (osteoprotegerin) gene, and rs3736228, on chromosome 11 in the LRP5 (lipoprotein-receptor-related protein) gene. A non-synonymous SNP in the LRP5 gene was associated with decreased bone mineral density (rs3736228, p=6·3×10 −12 for lumbar spine and p=1·9×10 −4 for femoral neck) and an increased risk of both osteoporotic fractures (odds ratio [OR] 1·3, 95% CI 1·09–1·52, p=0·002) and osteoporosis (OR 1·3, 1·08–1·63, p=0·008). Three SNPs near the TNFRSF11B gene were associated with decreased bone mineral density (top SNP, rs4355801: p=7·6×10 −10 for lumbar spine and p=3·3×10 −8 for femoral neck) and increased risk of osteoporosis (OR 1·2, 95% CI 1·01–1·42, p=0·038). For carriers of the risk allele at rs4355801, expression of TNFRSF11B in lymphoblast cell lines was halved (p=3·0×10 −6 ). 1883 (22%) of 8557 people were at least heterozygous for these risk alleles, and these alleles had a cumulative association with bone mineral density (trend p=2·3×10 −17 ). The presence of both risk alleles increased the risk of osteoporotic fractures (OR 1·3, 1·08–1·63, p=0·006) and this effect was independent of bone mineral density. Interpretation Two gene variants of key biological proteins increase the risk of osteoporosis and osteoporotic fracture. The combined effect of these risk alleles on fractures is similar to that of most well-replicated environmental risk factors, and they are present in more than one in five white people, suggesting a potential role in screening. Funding Wellcome Trust, European Commission, NWO Investments, Arthritis Research Campaign, Chronic Disease Research Foundation, Canadian Institutes of Health Research, European Society for Clinical and Economic Aspects of Osteoporosis, Genome Canada, Genome Quebec, Canada Research Chairs, National Health and Medical Research Council of Australia, and European Union.

Normal thyroid function is essential for health, but its genetic architecture remains poorly understood. Here, for the heritable thyroid traits thyrotropin (TSH) and free thyroxine (FT4), we analyse whole-genome sequence data from the UK10K project (N=2,287). Using additional whole-genome sequence and deeply imputed data sets, we report meta-analysis results for common variants (MAF≥1%) associated with TSH and FT4 (N=16,335). For TSH, we identify a novel variant in SYN2 (MAF=23.5%, P=6.15 × 10(-9)) and a new independent variant in PDE8B (MAF=10.4%, P=5.94 × 10(-14)). For FT4, we report a low-frequency variant near B4GALT6/SLC25A52 (MAF=3.2%, P=1.27 × 10(-9)) tagging a rare TTR variant (MAF=0.4%, P=2.14 × 10(-11)). All common variants explain ≥20% of the variance in TSH and FT4. Analysis of rare variants (MAF1%) using sequence kernel association testing reveals a novel association with FT4 in NRG1. Our results demonstrate that increased coverage in whole-genome sequence association studies identifies novel variants associated with thyroid function.

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

BACKGROUND: Low bone mineral density (BMD) is a primary risk factor for osteoporosis and is a highly heritable trait, but appears to be influenced by many genes. Genome-wide linkage studies have highlighted the chromosomal region 3p14-p22 as a quantitative trait locus for BMD (LOD 1.1 - 3.5). The FLNB gene, which is thought to have a role in cytoskeletal actin dynamics, is located within this chromosomal region and presents as a strong candidate for BMD regulation. We have previously identified significant associations between four SNPs in the FLNB gene and BMD in women. We have also previously identified associations between five SNPs located 5' of the transcription start site (TSS) and in intron 1 of the FLNB gene and expression of FLNB mRNA in osteoblasts in vitro. The latter five SNPs were genotyped in this study to test for association with BMD parameters in a family-based population of 769 Caucasian women. RESULTS: Using FBAT, significant associations were seen for femoral neck BMD Z-score with the SNPs rs11720285, rs11130605 and rs9809315 (P = 0.004 - 0.043). These three SNPs were also found to be significantly associated with total hip BMD Z-score (P = 0.014 - 0.026). We then combined the genotype data for these three SNPs with the four SNPs we previously identified as associated with BMD and performed a conditional analysis to determine whether they represent multiple independent associations with BMD. The results from this analysis suggested that these variants represent a single association signal. CONCLUSIONS: The SNPs identified in our studies as associated with BMD appear to be part of a single association signal between the FLNB gene and BMD in our data. FLNB is one of several genes located in 3p14-p22 that has been identified as significantly associated with BMD in Caucasian women.

Background: Osteoporosis is a highly heritable trait. Many candidate genes have been proposed as being involved in regulating bone mineral density (BMD). Few of these findings have been replicated in independent studies.Objective: To assess the relationship between BMD and fracture and all common single-nucleotide polymorphisms (SNPs) in previously proposed osteoporosis candidate genes.Design: Large-scale meta-analysis of genome-wide association data.Setting: 5 international, multicenter, population-based studies.Participants: Data on BMD were obtained from 19 195 participants (14 277 women) from 5 populations of European origin. Data on fracture were obtained from a prospective cohort (n = 5974) from the Netherlands.Measurements: Systematic literature review using the Human Genome Epidemiology Navigator identified autosomal genes previously evaluated for association with osteoporosis. We explored the common SNPs arising from the haplotype map of the human genome (HapMap) across all these genes. BMD at the femoral neck and lumbar spine was measured by dual-energy x-ray absorptiometry. Fractures were defined as clinically apparent, site-specific, validated nonvertebral and vertebral low-energy fractures.Results: 150 candidate genes were identified and 36 016 SNPs in these loci were assessed. SNPs from 9 gene loci (ESR1, LRP4, ITGA1, LRP5, SOST, SPP1, TNFRSF11A, TNFRSF11B, and TNFSF11) were associated with BMD at either site. For most genes, no SNP was statistically significant. For statistically significant SNPs (n = 241), effect sizes ranged from 0.04 to 0.18 SD per allele. SNPs from the LRP5, SOST, SPP1, and TNFRSF11A loci were significantly associated with fracture risk; odds ratios ranged from 1.13 to 1.43 per allele. These effects on fracture were partially independent of BMD at SPP1 and SOST.Limitation: Only common polymorphisms in linkage disequilibrium with SNPs in HapMap could be assessed, and previously reported associations for SNPs in some candidate genes could not be excluded.Conclusion: In this large-scale collaborative genome-wide meta-analysis, 9 of 150 candidate genes were associated with regulation of BMD, 4 of which also significantly affected risk for fracture. However, most candidate genes had no consistent association with BMD.

The 1p36 region of the human genome has been identified as containing a QTL for BMD in multiple studies. We analysed the TNFRSF1B gene from this region, which encodes the TNF receptor 2, in two large population-based cohorts. Our results suggest that variation in TNFRSF1B is associated with BMD. The TNFRSF1B gene, encoding the TNF receptor 2, is a strong positional and functional candidate gene for impaired bone structure through the role that TNF has in bone cells. The aims of this study were to evaluate the role of variations in the TNFRSF1B gene on bone structure and osteoporotic fracture risk in postmenopausal women. Six SNPs in TNFRSF1B were analysed in a cohort of 1,190 postmenopausal Australian women, three of which were also genotyped in an independent cohort of 811 UK postmenopausal women. Differences in phenotypic means for genotype groups were examined using one-way ANOVA and ANCOVA. Significant associations were seen for IVS1+5580A>G with BMD and QUS parameters in the Australian population (P = 0.008 − 0.034) and with hip BMD parameters in the UK population (P = 0.005 − 0.029). Significant associations were also observed between IVS1+6528G>A and hip BMD parameters in the UK cohort (P = 0.0002 − 0.003). We then combined the data from the two cohorts and observed significant associations between both IVS1+5580A>G and IVS1+6528G>A and hip BMD parameters (P = 0.002 − 0.033). Genetic variation in TNFRSF1B plays a role in the determination of bone structure in Caucasian postmenopausal women, possibly through effects on osteoblast and osteoclast differentiation.

Follistatin has been reported as a candidate gene for polycystic ovary syndrome (PCOS) from linkage and association studies. Acting to regulate the development of ovarian follicles and as an antagonist to aromatase activity, alterations in follistatin function or expression may result in key features of PCOS such as reduced serum FSH, impaired ovarian follicle development and augmented ovarian androgen production. We investigated polymorphisms in the FST gene to determine if genetic variation is associated with susceptibility to PCOS or key phenotypic features of PCOS patients in a case –control association study. One hundred and seventy-three PCOS patients of Caucasian descent (mean age 30.0 + 4.8 years), conforming to the NIH diagnostic criteria, were recruited from a clinical practice database and 107 normal ovulating women (mean age 38.8+ 13.4 years) were recruited from the general community as control subjects. Morphometric data, biochemistry and genomic DNA were collected from study subjects and genotyping was performed on seven Single nucleotide polymorphisms (SNPs) in the FST gene region. Allele frequencies of the SNPs were rs1423560 G/C (0.99/0.01), rs3797297 C/A (0.80/0.20), rs11745088 C/G (0.98/0.02), rs3203788 A/T (0.98/0.02) and rs1062809 G/C (1.00/—), rs1127760 A/T (0.98/0.02) and rs1127761 A/T (0.98/0.02), and these were not significantly different between the PCOS and control groups ( P< 0.05). Statistical analysis revealed significant associations between the SNP rs3797297 and sex hormone-binding globulin (P 5 0.04) and free androgen index (FAI) ( P< 0.01). We conclude that FST is not a susceptibility locus for PCOS; however, the SNP rs3797297 from FST gene was associated with androgenic markers for PCOS and may be of importance in the hyperandrogenaemia of the disease.

The murine arachidonate 15-lipoxygenase gene (Alox15) has recently been identified as a negative regulator of peak bone mineral density (BMD). The human ALOX15 gene shares significant sequence homology with the murine Alox15 gene; however, the human arachidonate 12-lipoxygenase gene (ALOX12) is functionally more similar to the mouse gene. Multiple single-nucleotide polymorphisms (SNPs) in the human ALOX15 and ALOX12 genes have previously been reported to be significantly associated with BMD in humans. On the basis of these data, we carried out our own investigation of the human ALOX15 and ALOX12 genes and their relationship with hip and spine BMD parameters. The study population consisted of 779 postmenopausal women with a mean (+/- standard deviation) age of 62.5 +/- 5.9 years at BMD measurement and was recruited from a single large general practice in Chingford, northeast London. Three SNPs from ALOX15 and five from ALOX12 were analyzed. None of the SNPs that we analyzed in ALOX15 were significantly associated with any of the BMD parameters or fracture data. However, we found that three SNPs from ALOX12, all previously associated with spine BMD in women, were significantly associated with spine and various hip BMD parameters in our cohort (P = 0.029-0.049). In conclusion, we found no association between polymorphism in ALOX15 and BMD phenotypes but were able to replicate previous findings that genetic variation in ALOX12 seems to play a role in determining bone structure in Caucasian women.

Bone mass is the single most important risk factor for osteoporotic fractures in the elderly and is mainly influenced by genetic factors accounting for 40-75% of the inter-individual variation. Critical for the bone remodeling process is the balance between the newly discovered members of the tumor necrosis factor ligand and receptor superfamilies, osteoprotegerin (OPG) and receptor activator of nuclear factor-kappaB ligand, which mediate the effects of many upstream regulators of bone metabolism. In the present study, we evaluated the impact of sequence variations in the OPG gene on bone mass, bone-related biochemistry including serum OPG and fracture frequency in elderly Australian women. A total of 1101 women were genotyped for 3 different single nucleotide polymorphisms (SNP) within the OPG gene (G1181C, T950C and A163G). The effects of these SNPs and serum OPG on calcaneal quantitative ultrasound measurements, osteodensitometry of the hip and bone-related biochemistry were examined. We found no significant relationship between sequence variations in the OPG gene or serum OPG and bone mass, bone-related biochemistry or fracture frequency. Our findings confirm some recent publications investigating the same SNPs but diverge from others, indicating that generalization of the relationships found in this type of study must be done with caution and signify the importance of determining associations between polymorphisms and osteoporosis in different ethnic groups.

The CYP19A1 gene encodes the enzyme aromatase, which is responsible for the final step in the biosynthesis of estrogens. Work from our group has shown that a microsatellite in intron 4 is associated with bone density and estrogen levels. Moreover a threonine to methionine missense mutation at amino acid position 201 in exon 6 of the aromatase gene causes a 5 fold increase in aromatase activity in an in-vitro assay. In this study we have characterized the effect of this variant on a large population of elderly women. A cohort of 1257 post-menopausal Western Australian women aged 75 ± 2.5 years in 1998 were recruited from the Western Australian electoral roll. Bone mineral density (BMD) at the hip and spine sites was measured by DXA (Hologic QDR 4500). Serum estradiol was measured by RIA (Orion Diagnostica, Finland). Genotyping of the T201M site was by Taqman (Applied Biosystems, Foster City, CA). Spinal deformity was calculated from morphometric x-ray calculations. The frequencies of the genotypes in the population were CC = 0.83, CT = 0.16, TT = 0.01. Using a dominant model the variant (T) was associated with increased estradiol for the heterozygote and rarer homozygote combined (CC: 25.5 ± 14.9 pmol/L vs. CT, TT: 38.3 ± 27.4 pmol/L; P = 4.5 × 10− 21). In addition, the variant allele was associated with increased BMD at the total hip (CC: 806 ± 126 mg/cm2 vs. CT, TT: 831 ± 129 mg/cm2; P = 0.015), femoral neck (CC: 686 ± 104 mg/cm2 vs. CT, TT: 706 ± 109 mg/cm2; P = 0.021). This variant was also associated with a reduction in spinal deformity (L4–T12) (CC: 37.6% vs. CT, TT: 21.1%; P = 0.012). We found that carriers of the T allele (17% of the population) show increased circulating estrogen, increased DXA BMD and reduced spine deformity. The observed associations are in keeping with the molecular effects of the nucleotide change previously characterised in vitro. These findings serve as the basis for further examination of the regulatory mechanisms of aromatase expression in human bone tissue and suggest that this variant may be of clinical use in predicting propensity to disease.