Body mass index stratified meta-analysis of genome-wide association studies of polycystic ovary syndrome in women of European ancestry.

BACKGROUND: 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. RESULTS: 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/m2). Two genome-wide significant loci were identified for lean PCOS, led by rs12000707 within DENND1A (P = 1.55 × 10-12) and rs2228260 within XBP1 (P = 3.68 × 10-8). One additional locus, LINC02905, was highlighted as significantly associated with lean PCOS through gene-based analyses (P = 1.76 × 10-6). 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-9) and a gene-based association was identified with ERBB4 (P = 1.59 × 10-6). 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. CONCLUSIONS: 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.

Bone Trans-omics: Integrating Omics to Unveil Mechanistic Molecular Networks Regulating Bone Biology and Disease.

PURPOSE OF REVIEW: Recent advancements in "omics" technologies and bioinformatics have afforded researchers new tools to study bone biology in an unbiased and holistic way. The purpose of this review is to highlight recent studies integrating multi-omics data gathered from multiple molecular layers (i.e.; trans-omics) to reveal new molecular mechanisms that regulate bone biology and underpin skeletal diseases. RECENT FINDINGS: Bone biologists have traditionally relied on single-omics technologies (genomics, transcriptomics, proteomics, and metabolomics) to profile measureable differences (both qualitative and quantitative) of individual molecular layers for biological discovery and to investigate mechanisms of disease. Recently, literature has grown on the implementation of integrative multi-omics to study bone biology, which combines computational and informatics support to connect multiple layers of data derived from individual "omic" platforms. This emerging discipline termed "trans-omics" has enabled bone biologists to identify and construct detailed molecular networks, unveiling new pathways and unexpected interactions that have advanced our mechanistic understanding of bone biology and disease. While the era of trans-omics is poised to revolutionize our capacity to answer more complex and diverse questions pertinent to bone pathobiology, it also brings new challenges that are inherent when trying to connect "Big Data" sets. A concerted effort between bone biologists and interdisciplinary scientists will undoubtedly be needed to extract physiologically and clinically meaningful data from bone trans-omics in order to advance its implementation in the field.

Molecular structure, expression, and the emerging role of Siglec-15 in skeletal biology and cancer.

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.

Relationships between attention to emotion and anxiety among a community sample of adolescents.

BACKGROUND: Attentional bias to threat has been implicated as a cognitive mechanism in anxiety disorders for youth. Yet, prior studies documenting this bias have largely relied on a method with questionable reliability (i.e. dot-probe task) and small samples, few of which included adolescents. The current study sought to address such limitations by examining relations between anxiety - both clinically diagnosed and dimensionally rated - and attentional bias to threat. METHODS: The study included a community sample of adolescents and employed eye-tracking methodology intended to capture possible biases across the full range of both automatic (i.e. vigilance bias) and controlled attentional processes (i.e. avoidance bias, maintenance bias). We examined both dimensional anxiety (across the full sample; n = 215) and categorical anxiety in a subset case-control analysis (n = 100) as predictors of biases. RESULTS: Findings indicated that participants with an anxiety disorder oriented more slowly to angry faces than matched controls. Results did not suggest a greater likelihood of initial orienting to angry faces among our participants with anxiety disorders or those with higher dimensional ratings of anxiety. Greater anxiety severity was associated with greater dwell time to neutral faces. CONCLUSIONS: This is the largest study to date examining eye-tracking metrics of attention to threat among healthy and anxious youth. Findings did not support the notion that anxiety is characterized by heightened vigilance or avoidance/maintenance of attention to threat. All effects detected were extremely small. Links between attention to threat and anxiety among adolescents may be subtle and highly dependent on experimental task dimensions.

Correction: Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

[This corrects the article DOI: 10.1371/journal.pgen.1007813.].

The Hippo in the room: Targeting the Hippo signalling pathway for osteosarcoma therapies.

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.

Life-Course Genome-wide Association Study Meta-analysis of Total Body BMD and Assessment of Age-Specific Effects.

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.

Large-scale genome-wide meta-analysis of polycystic ovary syndrome suggests shared genetic architecture for different diagnosis criteria.

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.

Genome-wide association study meta-analysis for quantitative ultrasound parameters of bone identifies five novel loci for broadband ultrasound attenuation.

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.

The relationship between depression and chronotype: A longitudinal assessment during childhood and adolescence.

BACKGROUND/OBJECTIVE: During adolescence, chronotype shifts toward "eveningness." "Eveningness" is related to negative physical and mental health outcomes. Little is known about what influences the shift in chronotype beyond pubertal status. The current study examined the influence of earlier depression predicting later individual differences in adolescent chronotype, accounting for pubertal status, and the prospective prediction of later increases in depression from earlier chronotype. METHODS: Youth (age M = 12.06, SD = 2.35; 56.5% girls) from the community completed repeated assessments of depression, including both self-reports (14 assessments) and diagnostic interviews (eight assessments), over a 48-month period. At the 36-month timepoint, participants completed chronotype and pubertal development measures. Regression and ANOVA analyses examined: (1) the influence of earlier depression levels (baseline to 36 months) upon chronotype, and (2) chronotype (at 36 months) upon later depression (48 months). RESULTS: Youth with higher earlier depression symptoms (ß = -0.347, P < .001) and history of depression diagnosis (ß = -0.13, P = .045) showed a greater eveningness preference controlling for pubertal status, age, and gender. Further, depression diagnosis history interacted with pubertal status to predict chronotype: (F(1,243) = 4.171, P = .045) such that the influence of depression on chronotype was greatest among postpubertal youth (t = 3.271, P = .002). Chronotype (greater eveningness preference) predicted prospective increases in depression symptoms (ß = -0.16, P = .03) and onset of depressive episode (b = -0.085, OR = 0.92, P = .03) 1 year later. CONCLUSION: Depression, experienced earlier in life, predicts greater preference for eveningness, especially among postpubertal youth. In turn, later depression is predicted by evening preference. These findings suggest the reciprocal interplay between mood and biological rhythms, especially depression and chronotype, during adolescence.

A Preliminary Multimethod Comparison of Sleep Among Adolescents With and Without Generalized Anxiety Disorder.

Previous studies suggest that youth with anxiety disorders experience their sleep as more disrupted and unsatisfying than their healthy peers. However, it is unclear whether these subjective complaints align with objective measures of sleep quantity and quality. The purpose of this preliminary study was to assess subjective and objective sleep parameters, and their relationships with anxiety symptomatology, among adolescents (62.8% female, 81.4% Caucasian), ages 12 to 18 (M = 15.29 years), with generalized anxiety disorder (n = 26) and controls without any psychopathology (n = 17). We measured sleep over 7 nights using sleep diaries and actigraphy and collected self- and parent-report questionnaires pertaining to sleep, anxiety, and depression. Repeated-measures mixed models were used to examine relationships between nightly sleep duration and morning anxiety. We found a number of differences in sleep between our anxious and healthy participants. Via sleep diary, our anxious participants had longer sleep onset latencies and lower satisfaction with sleep relative to controls, whereas via actigraphy we found longer sleep onset latencies but greater overall sleep duration among anxious versus control participants. Actigraphic measures of sleep disturbance were associated with parent-report of anxiety and depression. Our mixed-model analyses revealed that decreases in nightly sleep duration were associated with increased morning anxiety, but only among our participants with generalized anxiety disorder. Findings suggest that sleep disturbance among anxious adolescents can be detected using both subjective and objective measures and that, for these individuals, fluctuations in sleep duration may have real consequences for daytime anxiety.

Genome-wide association study using family-based cohorts identifies the WLS and CCDC170/ESR1 loci as associated with bone mineral density.

BACKGROUND: 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. RESULTS: 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). CONCLUSIONS: 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.

Differences in sleep disturbances among offspring of parents with and without bipolar disorder: association with conversion to bipolar disorder.

OBJECTIVES: Disruptions in sleep and dysregulation in circadian functioning may represent core abnormalities in the pathophysiology of bipolar disorder (BP). However, it is not clear whether these dysfunctions are state or trait markers of BP. This report compared sleep and circadian phenotypes among three groups: offspring of parents with BP diagnosed with BP at intake (BP/OB; n = 47), offspring of parents with BP without BP at intake (non-BP/OB; n = 386), and offspring of matched control parents who did not have BP (controls; n = 301). We also examined the association of baseline sleep parameters with subsequent development of BP among the non-BP/OB group. METHODS: Pittsburgh Bipolar Offspring Study youth (ages 6-18 years) and their parents completed assessments every two years pertaining to the child's sleep and circadian phenotypes and current psychopathology. Mixed-effects models examined differences in baseline sleep and circadian variables among the three groups. RESULTS: BP/OB offspring who were in a mood episode differed significantly on sleep parameters from the non-BP/OB and the offspring of controls, such as having inadequate sleep. Mixed logistic regression procedures showed that baseline sleep and circadian variables, such as frequent waking during the night, significantly predicted the development of BP among non-BP/OB over longitudinal follow-up. CONCLUSIONS: While lifetime diagnostic status accounted for differences among the groups in sleep and circadian disturbances, psychopathology explained the differences even further. Additionally, sleep disturbance may be a prognostic indicator of the development of BP in high-risk youth. Future studies are required to further disentangle whether sleep and circadian disruption are state or trait features of BP.

Epigenome-wide Association Study Shows Differential DNA Methylation of MDC1, KLF9, and CUTA in Autoimmune Thyroid Disease.

CONTEXT: Autoimmune thyroid disease (AITD) includes Graves disease (GD) and Hashimoto disease (HD), which often run in the same family. AITD etiology is incompletely understood: Genetic factors may account for up to 75% of phenotypic variance, whereas epigenetic effects (including DNA methylation [DNAm]) may contribute to the remaining variance (eg, why some individuals develop GD and others HD). OBJECTIVE: This work aimed to identify differentially methylated positions (DMPs) and differentially methylated regions (DMRs) comparing GD to HD. METHODS: Whole-blood DNAm was measured across the genome using the Infinium MethylationEPIC array in 32 Australian patients with GD and 30 with HD (discovery cohort) and 32 Danish patients with GD and 32 with HD (replication cohort). Linear mixed models were used to test for differences in quantile-normalized ß values of DNAm between GD and HD and data were later meta-analyzed. Comb-p software was used to identify DMRs. RESULTS: We identified epigenome-wide significant differences (P < 9E-8) and replicated (P < .05) 2 DMPs between GD and HD (cg06315208 within MDC1 and cg00049440 within KLF9). We identified and replicated a DMR within CUTA (5 CpGs at 6p21.32). We also identified 64 DMPs and 137 DMRs in the meta-analysis. CONCLUSION: Our study reveals differences in DNAm between GD and HD, which may help explain why some people develop GD and others HD and provide a link to environmental risk factors. Additional research is needed to advance understanding of the role of DNAm in AITD and investigate its prognostic and therapeutic potential.

Conditional testing of multiple variants associated with bone mineral density in the FLNB gene region suggests that they represent a single association signal.

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.

Identification of Differentially Expressed Genes and Molecular Pathways Involved in Osteoclastogenesis Using RNA-seq.

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.

Leveraging osteoclast genetic regulatory data to identify genes with a role in osteoarthritis.

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.

Age-dependent genetic regulation of osteoarthritis: independent effects of immune system genes.

OBJECTIVES: 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. METHODS: 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. RESULTS: 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. CONCLUSIONS: 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.

Sugar transporter Slc37a2 regulates bone metabolism in mice via a tubular lysosomal network in osteoclasts.

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.

Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis.

Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.