Novel Genetic Variants Associated With Increased Vertebral Volumetric BMD, Reduced Vertebral Fracture Risk, and Increased Expression of SLC1A3 and EPHB2

Genome-wide association study reveals the candidate genes of humerus quality in laying duck

Bone plays a crucial role in poultry's health and production. However, during the selection and cage farming, there has been a decline in bone quality. As the development of breeding theory, researchers find that it's possible to enhance bone quality through selective breeding.This study measure 8 humerus quality in 260 samples of the 350-day-old female duck. By descripting the basic characteristic traits, mechanical property traits we found that all the bone quality traits had a large variable coefficient, especially mechanical properties trait (20-70%), indicating that there was a large difference in bone health status among laying ducks. The phenotypic correlations showed a high correlation between weight and density, diameter and perimeter, breaking and toughness (r = 0.52-0.68). And then, we performed the Genome-wide association study (GWAS) to reveal the candidate genes of humerus quality in ducks. Seven candidate protein-coding genes were identified with perimeter trait, and 52 protein-coding genes were associated with toughness trait. We also analysed the candidate region and performed KEGG and GO analyse for 75 candidate genes. Furthermore, the expression analyse of the above candidate genes in different stage of humerus and different tissues were performed. Finally, AP2A2, SMAD3, SMNDC1, NFIA, EPHB2, PMEPA1, UNC5C, ESR1, VAV3, NFATC2 deserve further focus. The obtained results can contribute to new insight into bone quality and provide new genetic biomarkers for application in duck breeding programs.

GREM2 inactivation increases trabecular bone mass in mice

Osteoporosis is a common skeletal disease affecting millions of individuals world-wide, with an increased risk of fracture, and a decreased quality of life. Despite its well-known consequences, the etiology of osteoporosis and optimal treatment methods are not fully understood. Human genetic studies have identified genetic variants within the FMN2/GREM2 locus to be associated with trabecular volumetric bone mineral density (vBMD) and vertebral and forearm fractures, but not with cortical bone parameters. GREM2 is a bone morphogenetic protein (BMP) antagonist. In this study, we employed Grem2-deficient mice to investigate whether GREM2 serves as the plausible causal gene for the fracture signal at the FMN2/GREM2 locus. We observed that Grem2 is moderately expressed in bone tissue and particularly in osteoblasts. Complete Grem2 gene deletion impacted mouse survival and body growth. Partial Grem2 inactivation in Grem2+/- female mice led to increased trabecular BMD of femur and increased trabecular bone mass in tibia due to increased trabecular thickness, with an unchanged cortical thickness, as compared with wildtype littermates. Furthermore, Grem2 inactivation stimulated osteoblast differentiation, as evidenced by higher alkaline phosphatase (Alp), osteocalcin (Bglap), and osterix (Sp7) mRNA expression after BMP-2 stimulation in calvarial osteoblasts and osteoblasts from the long bones of Grem2-/- mice compared to wildtype littermates. These findings suggest that GREM2 is a possible target for novel osteoporotic treatments, to increase trabecular bone mass and prevent osteoporotic fractures.

Large-scale circulating proteome association study (CPAS) meta-analysis identifies circulating proteins and pathways predicting incident hip fractures

Hip fractures are associated with significant disability, high cost, and mortality. However, the exact biological mechanisms underlying susceptibility to hip fractures remain incompletely understood. In an exploratory search of the underlying biology as reflected through the circulating proteome, we performed a comprehensive Circulating Proteome Association Study (CPAS) meta-analysis for incident hip fractures. Analyses included 6430 subjects from two prospective cohort studies (Cardiovascular Health Study and Trøndelag Health Study) with circulating proteomics data (aptamer-based 5 K SomaScan version 4.0 assay; 4979 aptamers). Associations between circulating protein levels and incident hip fractures were estimated for each cohort using age and sex-adjusted Cox regression models. Participants experienced 643 incident hip fractures. Compared with the individual studies, inverse-variance weighted meta-analyses yielded more statistically significant associations, identifying 23 aptamers associated with incident hip fractures (conservative Bonferroni correction 0.05/4979, P < 1.0 × 10-5). The aptamers most strongly associated with hip fracture risk corresponded to two proteins of the growth hormone/insulin growth factor system (GHR and IGFBP2), as well as GDF15 and EGFR. High levels of several inflammation-related proteins (CD14, CXCL12, MMP12, ITIH3) were also associated with increased hip fracture risk. Ingenuity pathway analysis identified reduced LXR/RXR activation and increased acute phase response signaling to be overrepresented among those proteins associated with increased hip fracture risk. These analyses identified several circulating proteins and pathways consistently associated with incident hip fractures. These findings underscore the usefulness of the meta-analytic approach for comprehensive CPAS in a similar manner as has previously been observed for large-scale human genetic studies. Future studies should investigate the underlying biology of these potential novel drug targets.

Association between AXIN1 gene polymorphism (rs9921222) of WNT signaling pathway and susceptibility to osteoporosis in Egyptian patients: a case-control study

BACKGROUND

Osteoporosis (OP) is the most prevalent metabolic bone disease. Numerous genetic loci are strongly related to OP. AXIN1 is a significant gene that serves an important role in the WNT signaling pathway. The aim of this study was to explore the association between the AXIN1 genetic polymorphism (rs9921222) and OP susceptibility.

METHODS

A total of 101 subjects were enrolled in the study (50 patients with OP and 51 healthy individuals). Genomic DNA was extracted from whole blood using the QIAamp DNA Blood Mini Kit, and the AXIN1 gene polymorphism (rs9921222) was genotyped by TaqMan allelic discrimination assays. A logistic regression analysis was used to assess the association between genotypes and OP risk.

RESULTS

We found that AXIN1 rs9921222 had a significant association with the susceptibility of OP under the homozygote model (TT vs. CC: OR = 16.6, CI = 2.03-136.4, p = 0.009), (CT vs. CC: OR = 6.3, CI = 1.23-31.8, p = 0.027), recessive genetic model (TT vs.TC-CC: OR = 13.6, CI = 1.7-110.4, p = 0.015), and the dominant model (TT-TC vs. CC: OR = 9.7, CI = 2.6-36.3, p < 0.001). Allele T was significantly associated with OP risk (T vs. C: OR = 10.5, CI = 3.5-31.15, p = 0.001). There was a statistically significant difference between genotypes in mean platelet volume (p = 0.004), and platelet distribution width (p = 0.025). In addition, lumbar spine bone density, and femur neck bone density were significantly different between genotypes (p < 0.001).

CONCLUSION

AXIN1 rs9921222 was associated with OP susceptibility in the Egyptian population and should be considered a potential determinant risk for OP.

Risk Variants in or Near ZBTB40 AND NFATC1 Increase the Risk of Both IBD and Adverse Bone Health Outcomes Highlighting Common Genetic Underpinnings Across Both Diseases

BACKGROUND

Inflammatory bowel disease (IBD) is associated with an increased risk of osteoporosis and bone fracture. The aims of this study were to (1) confirm the association between IBD and low bone density and (2) test for shared risk variants across diseases.

METHODS

The study cohort included patients from the Michigan Genomics Initiative. Student's t tests (continuous) and chi-square tests (categorical) were used for univariate analyses. Multivariable logistic regression was performed to test the effect of IBD on osteoporosis or osteopenia. Publicly available genome-wide association summary statistics were used to identify variants that alter the risk of IBD and bone density, and Mendelian randomization (MR) was used to identify causal effects of genetically predicted IBD on bone density.

RESULTS

There were 51 405 individuals in the Michigan Genomics Initiative cohort including 10 378 (20.2%) cases of osteoporosis or osteopenia and 1404 (2.7%) cases of IBD. Patients with osteoporosis or osteopenia were more likely to be older (64 years of age vs 56 years of age; P < .001), female (67% vs 49%; P < .001), and have a lower body mass index (29 kg/m2 vs 30 kg/m2; P < .001). IBD patients with (odds ratio, 4.60; 95% confidence interval, 3.93-5.37) and without (odds ratio, 1.77; 95% confidence interval, 1.42-2.21) steroid use had a significantly higher risk of osteoporosis or osteopenia. Twenty-one IBD variants associated with reduced bone mineral density at P ≤ .05 and 3 IBD risk variants associated with reduced bone mineral density at P ≤ 5 × 10-8. Of the 3 genome-wide significant variants, 2 increased risk of IBD (rs12568930-T: MIR4418;ZBTB40; rs7236492-C: NFATC1). MR did not reveal a causal effect of genetically predicted IBD on bone density (MR Egger, P = .30; inverse variance weighted, P = .63).

CONCLUSIONS

Patients with IBD are at increased risk for low bone density, independent of steroid use. Variants in or near ZBTB40 and NFATC1 are associated with an increased risk of IBD and low bone density.

Aberrant paracrine signalling for bone remodelling underlies the mutant histone-driven giant cell tumour of bone

Oncohistones represent compelling evidence for a causative role of epigenetic perturbations in cancer. Giant cell tumours of bone (GCTs) are characterised by a mutated histone H3.3 as the sole genetic driver present in bone-forming osteoprogenitor cells but absent from abnormally large bone-resorbing osteoclasts which represent the hallmark of these neoplasms. While these striking features imply a pathogenic interaction between mesenchymal and myelomonocytic lineages during GCT development, the underlying mechanisms remain unknown. We show that the changes in the transcriptome and epigenome in the mesenchymal cells caused by the H3.3-G34W mutation contribute to increase osteoclast recruitment in part via reduced expression of the TGFβ-like soluble factor, SCUBE3. Transcriptional changes in SCUBE3 are associated with altered histone marks and H3.3^G34W enrichment at its enhancer regions. In turn, osteoclasts secrete unregulated amounts of SEMA4D which enhances proliferation of mutated osteoprogenitors arresting their maturation. These findings provide a mechanism by which GCTs undergo differentiation in response to denosumab, a drug that depletes the tumour of osteoclasts. In contrast, hTERT alterations, commonly found in malignant GCT, result in the histone-mutated neoplastic cells being independent of osteoclasts for their proliferation, predicting unresponsiveness to denosumab. We provide a mechanism for the initiation of GCT, the basis of which is dysfunctional cross-talk between bone-forming and bone-resorbing cells. The findings highlight the role of tumour/microenvironment bidirectional interactions in tumorigenesis and how this is exploited in the treatment of GCT.

Inbreeding depression and the probability of racing in the Thoroughbred horse

Small effective population sizes and active inbreeding can lead to inbreeding depression due to deleterious recessive mutations exposed in the homozygous state. The Thoroughbred racehorse has low levels of population genetic diversity, but the effects of genomic inbreeding in the population are unknown. Here, we quantified inbreeding based on runs of homozygosity (ROH) using 297 K SNP genotypes from 6128 horses born in Europe and Australia, of which 13.2% were unraced. We show that a 10% increase in inbreeding (F_ROH) is associated with a 7% lower probability of ever racing. Moreover, a ROH-based genome-wide association study identified a haplotype on ECA14 which, in its homozygous state, is linked to a 32.1% lower predicted probability of ever racing, independent of F_ROH. The haplotype overlaps a candidate gene, EFNA5, that is highly expressed in cartilage tissue, which when damaged is one of the most common causes of catastrophic musculoskeletal injury in racehorses. Genomics-informed breeding aiming to reduce inbreeding depression and avoid damaging haplotype carrier matings will improve population health and racehorse welfare.

GATA4 and estrogen receptor alpha bind at SNPs rs9921222 and rs10794639 to regulate AXIN1 expression in osteoblasts

Osteoporosis is a serious public health problem that affects 200 million people worldwide. Genome-wide association studies have revealed the association between several single nucleotide polymorphisms (SNPs) near WNT/β-catenin signaling genes and bone mineral density (BMD). The activation of β-catenin by WNT ligands is required for osteoblast differentiation. SNP rs9921222 is an intronic variant of AXIN1 (a scaffold protein in the destruction complex that regulates β-catenin signaling) that correlates with BMD. However, the biological mechanism of SNP rs9921222 has never been reported. Here, we show that the genotype of SNP rs9921222 correlates with the expression of AXIN1 in human osteoblasts. RNA and genomic DNA were analyzed from primary osteoblasts from 111 different individuals. Homozygous TT at rs9921222 correlates with a higher expression of AXIN1 than homozygous CC. Regional association analysis showed that rs9921222 is in high linkage disequilibrium (LD) with SNP rs10794639. In silico transcription factor analysis predicted that rs9921222 is within a GATA4 motif and rs10794639 is adjacent to an estrogen receptor alpha (ERα) motif. Mechanistically, GATA4 and ERα bind at SNPs rs9921222 and rs10794639 as detected by ChIP-qPCR. Luciferase assays demonstrate that rs9921222 is the causal SNP to alter ERα and GATA4 binding. GATA4 promoted the expression, and in contrast, ERα suppressed the expression of AXIN1 via the histone deacetylase complex member SIN3A. Functionally, the level of AXIN1 negatively correlates with the level of transcriptionally active β-catenin. In summary, we have discovered a molecular mechanism of the SNP rs9921222 to regulate AXIN1 through GATA4 and ERα binding in human osteoblasts.

Associations of osteoprotegerin (OPG) TNFRSF11B gene polymorphisms with risk of fractures in older adult populations: meta-analysis of genetic and genome-wide association studies

The meta-analysis of osteoprotegerin (OPG) (TNFRSF11B) polymorphisms from genetic association studies and genome-wide association studies was performed in order to test the hypothesis of association between OPG polymorphisms and fracture. The findings showed a significant 13% to 37% protective effect of OPG on fractures in postmenopausal women (PSM) (rs2073618), overall, ≥ 60y and Western subjects (rs3134069 and rs3134070).

PURPOSE

Fractures in older people usually result from compromised bone integrity. The multifactorial aetiology of fractures includes both genetic and environmental factors. Inconsistency of reported associations of osteoprotegerin (OPG) (TNFRSF11B) polymorphisms with fracture in the older adult population warranted a meta-analysis to determine more precise estimates.

METHODS

We searched for all available literature on OPG (TNFRSF11B) and fracture. Four polymorphisms were examined, one exonic (rs2073618) and three intronic (rs3134069, rs3134070 and rs3102735). The first two intron polymorphisms were combined (OPGI: osteoprotegerin intron) on account of complete linkage disequilibrium. Risks were estimated with odds ratios (ORs) and 95% confidence intervals (CIs) using the allele-genotype model that included variant (var), wild-type (wt) and heterozygote (het). Multiple comparisons were Bonferroni-corrected. We used meta-regression to examine sources of heterogeneity. Zero heterogeneity (homogeneity: I² = 0%) and high significance (P^a < 0.00001) were the criteria for strength of evidence. Significant outcomes were subjected to sensitivity analysis and publication bias assessment.

RESULTS

From 13 articles (11 genetic association and two genome-wide), this meta-analysis generated five significant pooled ORs, all indicating reduced risks (ORs 0.44-0.87). Of the five, four highly significant comparisons (P^a ≤ 0.00001-0.002) survived the Bonferroni correction, one in rs2073618 het model of the postmenopausal women (OR 0.87, 95% CI 0.81-0.92, I² = 0%) and the other three in OPGI wt model of the overall analysis, ≥ 60 y and Western subjects (ORs 0.63-0.71, 95% CI 0.47-0.86, I² = 97-99%). These findings were consistent, had high significance and high statistical power and were robust and without evidence of publication bias. Four covariates (year of publication, study quality, fracture type/site and sample size) were the sources of heterogeneity in the OPGI overall outcomes (P^a = 0.0001-0.03).

CONCLUSION

Evidence showed that the OPG (TNFRSF11B) polymorphisms reduced the risk for fracture in older adults, particularly protective among postmenopausal women, ≥ 60 y and Western subjects.

A transcriptome-wide association study to detect novel genes for volumetric bone mineral density

Transcriptome-wide association studies (TWAS) systematically investigate the association of genetically predicted gene expression with disease risk, providing an effective approach to identify novel susceptibility genes. Osteoporosis is the most common metabolic bone disease, associated with reduced bone mineral density (BMD) and increased risk of osteoporotic fractures, whereas genetic factors explain approximately 70% of the variance in phenotypes associated with bone. BMD is commonly assessed using dual-energy X-ray absorptiometry (DXA) to obtain measurements (g/cm2) of areal BMD. However, quantitative computed tomography (QCT) measured 3D volumetric BMD (vBMD) (g/cm3) has important advantages compared with DXA since it can evaluate cortical and trabecular microstructural features of bone quality, which can be used to directly predict fracture risk. Here, we performed the first TWAS for volumetric BMD (vBMD) by integrating genome-wide association studies (GWAS) data from two independent cohorts, namely the Framingham Heart Study (FHS, n = 3298) and the Osteoporotic Fractures in Men (MrOS, n = 4641), with tissue-specific gene expression data from the Genotype-Tissue Expression (GTEx) project. We first used stratified linkage disequilibrium (LD) score regression approach to identify 12 vBMD-relevant tissues, for which vBMD heritability is enriched in tissue-specific genes of the given tissue. Focusing on these tissues, we subsequently leveraged GTEx expression reference panels to predict tissue-specific gene expression levels based on the genotype data from FHS and MrOS. The associations between predicted gene expression levels and vBMD variation were then tested by MultiXcan, an innovative TWAS method that integrates information available across multiple tissues. We identified 70 significant genes associated with vBMD, including some previously identified osteoporosis-related genes such as LYRM2 and NME8, as well as some novel loci such as DNAAF2 and SPAG16. Our findings provide novel insights into the pathophysiological mechanisms of osteoporosis and highlight several novel vBMD-associated genes that warrant further investigation.

WNT4 Balances Development vs Disease in Gynecologic Tissues and Women's Health

The WNT family of proteins is crucial in numerous developmental pathways and tissue homeostasis. WNT4, in particular, is uniquely implicated in the development of the female phenotype in the fetus, and in the maintenance of müllerian and reproductive tissues. WNT4 dysfunction or dysregulation can drive sex-reversal syndromes, highlighting the key role of WNT4 in sex determination. WNT4 is also critical in gynecologic pathologies later in life, including several cancers, uterine fibroids, endometriosis, and infertility. The role of WNT4 in normal decidualization, implantation, and gestation is being increasingly appreciated, while aberrant activation of WNT4 signaling is being linked both to gynecologic and breast cancers. Notably, single-nucleotide polymorphisms (SNPs) at the WNT4 gene locus are strongly associated with these pathologies and may functionally link estrogen and estrogen receptor signaling to upregulation and activation of WNT4 signaling. Importantly, in each of these developmental and disease states, WNT4 gene expression and downstream WNT4 signaling are regulated and executed by myriad tissue-specific pathways. Here, we review the roles of WNT4 in women's health with a focus on sex development, and gynecologic and breast pathologies, and our understanding of how WNT4 signaling is controlled in these contexts. Defining WNT4 functions provides a unique opportunity to link sex-specific signaling pathways to women's health and disease.

Level and change in bone microarchitectural parameters and their relationship with previous fracture and established bone mineral density loci

BACKGROUND

Osteoporosis is characterised by a reduction of bone mineral density (BMD) and predisposition to fracture. Bone microarchitecture, measured by high resolution peripheral quantitative computed tomography (HR-pQCT), has been related to fragility fractures and BMD and has been the subject of large-scale genome-wide analysis. We investigated whether fracture was related to baseline values and longitudinal changes in bone microarchitecture and whether bone microarchitecture was associated with established BMD loci.

METHODS

115 males and 99 females (aged 72-81 at baseline) from the Hertfordshire Cohort Study (HCS) were analysed. Fracture history was determined in 2011-2012 by self-report and vertebral fracture assessment. Participants underwent HR-pQCT scans of the distal radius and tibia in 2011-2012 and 2017. Previous fracture in relation to baseline values and changes in tibial HR-pQCT parameters was examined using sex-adjusted logistic regression with and without adjustment for age, sociodemographic, lifestyle and clinical characteristics; baseline values and changes in parameters associated with previous fracture were then examined in relation to four established BMD loci after adjustment for sex and age.

RESULTS

Previous fracture was related to: higher trabecular area (fully-adjusted odds ratio [95% CI] per SD greater baseline value: 2.18 [1.27,3.73], p = 0.005); lower total volumetric BMD (0.53 [0.34,0.84], p = 0.007), cortical area (0.53 [0.30,0.95], p = 0.032), cortical BMD (0.56 [0.36,0.88], p = 0.011) and cortical thickness (0.45 [0.27,0.77], p = 0.004); and greater declines in trabecular BMD (p = 0.001). Associations were robust in sex- and fully-adjusted analysis. Relationships between BMD loci and these HR-pQCT parameters were weak: rs3801387 (WNT16) was related to decline in trabecular BMD (p = 0.011) but no other associations were significant (p > 0.05).

CONCLUSION

Baseline values of HR-pQCT parameters and greater decline in trabecular BMD were associated with fracture. Change in trabecular BMD was associated with WNT16 which has been demonstrated to influence bone health in murine models and human genome-wide association studies (GWAS).

Genetic variants modify the associations of concentrations of methylmalonic acid, vitamin B-12, vitamin B-6, and folate with bone mineral density

BACKGROUND

Elevated plasma homocysteine has been found to be associated with an increased risk of osteoporosis, especially hip and vertebral fractures. The plasma concentration of homocysteine is dependent on the activities of several B vitamin-dependent enzymes, such as methylenetetrahydrofolate reductase (MTHFR), methionine synthase (MTR), methionine synthase reductase (MTRR), and cystathionine β-synthase (CBS).

OBJECTIVES

We investigated whether genetic variants in some of the genes involved in 1 carbon metabolism modify the association of B vitamin-related measures with bone mineral density (BMD) and strength.

METHODS

We measured several B vitamins and biomarkers in participants of the Framingham Offspring Study, and performed analyses of methylmalonic acid (MMA) continuously and <210 nmol/L; pyridoxal-5'-phosphate; vitamin B-12 continuously and ≥258 pmol/L; and folate. The outcomes of interest included areal and volumetric BMD, measured by DXA and quantitative computed tomography (QCT), respectively. We evaluated associations between the bone measures and interactions of single nucleotide polymorphism with a B vitamin or biomarker in Framingham participants (n = 4310 for DXA and n = 3127 for QCT). For analysis of DXA, we validated the association results in the B-PROOF cohort (n = 1072). Bonferroni-corrected locus-wide significant thresholds were defined to account for multiple testing.

RESULTS

The interactions between rs2274976 and vitamin B-12 and rs34671784 and MMA <210 nmol/L were associated with lumbar spine BMD, and the interaction between rs6586281 and vitamin B-12 ≥258 pmol/L was associated with femoral neck BMD. For QCT-derived traits, 62 interactions between genetic variants and B vitamins and biomarkers were identified.

CONCLUSIONS

Some genetic variants in the 1-carbon methylation pathway modify the association of B vitamin and biomarker concentrations with bone density and strength.  These interactions require further replication and functional validation for a mechanistic understanding of the role of the 1-carbon metabolism pathway on BMD and risks of fracture.

Epidemiological transition to mortality and refracture following an initial fracture

This study sought to redefine the concept of fracture risk that includes refracture and mortality, and to transform the risk into "skeletal age". We analysed data obtained from 3521 women and men aged 60 years and older, whose fracture incidence, mortality, and bone mineral density (BMD) have been monitored since 1989. During the 20-year follow-up period, among 632 women and 184 men with a first incident fracture, the risk of sustaining a second fracture was higher in women (36%) than in men (22%), but mortality risk was higher in men (41%) than in women (25%). The increased risk of mortality was not only present with an initial fracture, but was accelerated with refractures. Key predictors of post-fracture mortality were male gender (hazard ratio [HR] 2.4; 95% CI, 1.79–3.21), advancing age (HR 1.67; 1.53–1.83), and lower femoral neck BMD (HR 1.16; 1.01–1.33). A 70-year-old man with a fracture is predicted to have a skeletal age of 75. These results were incorporated into a prediction model to aid patient-doctor discussion about fracture vulnerability and treatment decisions.

The Musculoskeletal Knowledge Portal: Making Omics Data Useful to the Broader Scientific Community

The development of high-throughput genotyping technologies and large biobank collections, complemented with rapid methodological advances in statistical genetics, has enabled hypothesis-free genome-wide association studies (GWAS), which have identified hundreds of genetic variants across many loci associated with musculoskeletal conditions. Similarly, basic scientists have valuable molecular cellular and animal data based on musculoskeletal disease that would be enhanced by being able to determine the human translation of their findings. By integrating these large-scale human genomic musculoskeletal datasets with complementary evidence from model organisms, new and existing genetic loci can be statistically fine-mapped to plausibly causal variants, candidate genes, and biological pathways. Genes and pathways identified using this approach can be further prioritized as drug targets, including side-effect profiling and the potential for new indications. To bring together these big data, and to realize the vision of creating a knowledge portal, the International Federation of Musculoskeletal Research Societies (IFMRS) established a working group to collaborate with scientists from the Broad Institute to create the Musculoskeletal Knowledge Portal (MSK-KP)(http://mskkp.org/). The MSK consolidates omics datasets from humans, cellular experiments, and model organisms into a central repository that can be accessed by researchers. The vision of the MSK-KP is to enable better understanding of the biological mechanisms underlying musculoskeletal disease and apply this knowledge to identify and develop new disease interventions. © 2020 American Society for Bone and Mineral Research (ASBMR).

Distinguishing patients with idiopathic epilepsy from solitary cysticercus granuloma epilepsy and biochemical phenotype assessment using a serum biomolecule profiling platform

A major source of epilepsy is Neurocysticercosis (NCC), caused by Taenia solium infection. Solitary cysticercus granuloma (SCG), a sub-group of NCC induced epilepsy, is the most common form of NCC in India. Current diagnostic criteria for SCG epilepsy require brain imaging which may not be available in communities where the disease is endemic. Identification of serum changes and potential biomolecules that could distinguish SCG epilepsy from idiopathic generalized epilepsy (IE), without the initial need for imaging, could assist in disease identification, understanding, and treatment. The objective here was to investigate, using mass spectrometry (MS), sera biomolecule differences between patients with SCG epilepsy or IE to help distinguish these disorders based on physiological differences, to understand underlying phenotypes and mechanisms, and to lay ground work for future therapeutic and biomarker analyses. Sera were obtained from patients with SCG or IE (N = 29 each group). Serum mass peak profiling was performed with electrospray ionization (ESI) MS, and mass peak area means in the two groups were compared using leave one [serum sample] out cross validation (LOOCV). Serum LOOCV analysis identified significant differences between SCG and IE patient groups (p = 10-20), which became non-significant (p = 0.074) when the samples were randomly allocated to the groups and reanalyzed. Tandem MS/MS peptide analysis of serum mass peaks from SCG or IE patients was performed to help identify potential peptide/protein biochemical and phenotypic changes involving these two forms of epilepsy. Bioinformatic analysis of these peptide/protein changes suggested neurological, inflammatory, seizure, blood brain barrier, cognition, ion channel, cell death, and behavior related biochemical systems were being altered in these disease states. This study provides groundwork for aiding in distinguishing SCG and IE patients in minimally invasive, lower-cost manners, for improving understanding of underlying epilepsy mechanisms, and for further identifying discriminatory biomarkers and potential therapeutic targets.

Distinct Subsets of Noncoding RNAs Are Strongly Associated With BMD and Fracture, Studied in Weight-Bearing and Non-Weight-Bearing Human Bone

We investigated mechanisms resulting in low bone mineral density (BMD) and susceptibility to fracture by comparing noncoding RNAs (ncRNAs) in biopsies of non-weight-bearing (NWB) iliac (n = 84) and weight bearing (WB) femoral (n = 18) postmenopausal bone across BMDs varying from normal (T-score > -1.0) to osteoporotic (T-score ≤ -2.5). Global bone ncRNA concentrations were determined by PCR and microchip analyses. Association with BMD or fracture, adjusted by age and body mass index, were calculated using linear and logistic regression and least absolute shrinkage and selection operator (Lasso) analysis. At 10% false discovery rate (FDR), 75 iliac bone ncRNAs and 94 femoral bone ncRNAs were associated with total hip BMD. Eight of the ncRNAs were common for the two sites, but five of them (miR-484, miR-328-3p, miR-27a-5p, miR-28-3p, and miR-409-3p) correlated positively to BMD in femoral bone, but negatively in iliac bone. Of predicted pathways recognized in bone metabolism, ECM-receptor interaction and proteoglycans in cancer emerged at both sites, whereas fatty acid metabolism and focal adhesion were only identified in iliac bone. Lasso analysis and cross-validations identified sets of nine bone ncRNAs correlating strongly with adjusted total hip BMD in both femoral and iliac bone. Twenty-eight iliac ncRNAs were associated with risk of fracture (FDR < 0.1). The small nucleolar RNAs, RNU44 and RNU48, have a function in stabilization of ribosomal RNAs (rRNAs), and their association with fracture and BMD suggest that aberrant processing of rRNAs may be involved in development of osteoporosis. Cis-eQTL (expressed quantitative trait loci) analysis of the iliac bone biopsies identified two loci associated with microRNAs (miRNAs), one previously identified in a heel-BMD genomewide association study (GWAS). In this comprehensive investigation of the skeletal genetic background in postmenopausal women, we identified functional bone ncRNAs associated to fracture and BMD, representing distinct subsets in WB and NWB skeletal sites. © 2020 The Authors. Journal of Bone and Mineral Research published by American Society for Bone and Mineral Research.

Pregnancy and lactation, a challenge for the skeleton

In this review we discuss skeletal adaptations to the demanding situation of pregnancy and lactation. Calcium demands are increased during pregnancy and lactation, and this is effectuated by a complex series of hormonal changes. The changes in bone structure at the tissue and whole bone level observed during pregnancy and lactation appear to largely recover over time. The magnitude of the changes observed during lactation may relate to the volume and duration of breastfeeding and return to regular menses. Studies examining long-term consequences of pregnancy and lactation suggest that there are small, site-specific benefits to bone density and that bone geometry may also be affected. Pregnancy- and lactation-induced osteoporosis (PLO) is a rare disease for which the pathophysiological mechanism is as yet incompletely known; here, we discuss and speculate on the possible roles of genetics, oxytocin, sympathetic tone and bone marrow fat. Finally, we discuss fracture healing during pregnancy and lactation and the effects of estrogen on this process.

Gene-gene and gene-lifestyle interactions of AKAP11, KCNMA1, PUM1, SPTBN1, and EPDR1 on osteoporosis risk in middle-aged adults

OBJECTIVES

Osteoporosis is associated with genetic and environmental factors. The aim of this article was to determine how the polygenic risk scores (PRS) of genetic variants that affect osteoporosis and its related signaling interact with the lifestyle of middle-aged adults.

METHODS

The study examined 8845 participants from Ansan/Ansung cohorts. Osteoporosis was defined as a T-score of bone mineral density ≤-2.5 in either the wrist or tibia; 1136 participants had osteoporosis. Genome-wide association studies of individuals 40 to 65 y of age were conducted and the best gene-gene interactions from the genetic variants related to osteoporosis were selected and explored using the generalized multifactor dimensionality reduction method. PRS for the best model (PRSBM) was calculated by weighted PRS that was divided into low, medium, and high groups.

RESULTS

The model that contributed the most influence on osteoporosis risk with gene-gene interactions included AKAP11_rs238340, KCNMA1_ rs628948, PUM1_rs7529390, SPTBN1_ rs6752877, and EPDR1_rs2722298. The risk for osteoporosis in the tibia was elevated by 1.71-fold in the high PRSBM group compared with the low PRSBM group. Energy and nutrient intake did not have any interaction with PRSBM and thus did not influence risk for osteoporosis. However, interestingly, only coffee and caffeine intake did interact with PRSBM and affected risk for osteoporosis. In patients with low coffee (<3 cup/wk) and caffeine(<60 mg/d) consumption, osteoporosis risk was higher in the high PRSBM group than the low PRSBM group by 2.27- and 2.29-fold, respectively. In the low coffee intake group, bone mineral density in the high PRSBM group was significantly higher than in the low PRSBM arm.

CONCLUSIONS

Carriers with high PRSBM increased susceptibility to osteoporosis, especially in low coffee and caffeine intake. The results can be applied to personalized nutrition for lowering the risk for osteoporosis.

BMD-Related Genetic Risk Scores Predict Site-Specific Fractures as Well as Trabecular and Cortical Bone Microstructure

CONTEXT

It is important to identify patients at highest risk of fractures.

OBJECTIVE

To compare the separate and combined performances of bone-related genetic risk scores (GRSs) for prediction of forearm, hip and vertebral fractures separately, as well as of trabecular and cortical bone microstructure parameters separately.

DESIGN, SETTING, AND PARTICIPANTS

Using 1103 single nucleotide polymorphisms (SNPs) independently associated with estimated bone mineral density of the heel (eBMD), we developed a weighted GRS for eBMD and determined its contribution to fracture prediction beyond 2 previously developed GRSs for femur neck BMD (49 SNPs) and lumbar spine BMD (48 SNPs). Associations between these GRSs and forearm (ncases = 1020; ncontrols = 2838), hip (ncases = 1123; ncontrols = 2630) and vertebral (ncases = 288; ncontrols = 1187) fractures were evaluated in 3 Swedish cohorts. Associations between the GRSs and trabecular and cortical bone microstructure parameters (n = 426) were evaluated in the MrOS Sweden cohort.

RESULTS

We found that eBMDGRS was the only significant independent predictor of forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS were significant independent predictors of hip fractures. The eBMDGRS was the major GRS contributing to prediction of trabecular bone microstructure parameters while both FN-BMDGRS and eBMDGRS contributed information for prediction of cortical bone microstructure parameters.

CONCLUSIONS

The eBMDGRS independently predicts forearm and vertebral fractures while both FN-BMDGRS and eBMDGRS contribute independent information for prediction of hip fractures. We propose that eBMDGRS captures unique information about trabecular bone microstructure useful for prediction of forearm and vertebral fractures. These findings may facilitate personalized medicine to predict site-specific fractures as well as cortical and trabecular bone microstructure separately.

Haplotype-based genome-wide association studies for carcass and growth traits in chicken

There have been several genome-wide association study (GWAS) reported for carcass, growth, and meat traits in chickens. Most of these studies have been based on single SNPs GWAS. In contrast, haplotype-based GWAS reports have been limited. In the present study, 2 Northeast Agricultural University broiler lines divergently selected for abdominal fat content (NEAUHLF) and genotyped with the chicken 60K SNP chip were used to perform a haplotype-based GWAS. The lean and fat chicken lines were selected for abdominal fat content for 11 yr. Abdominal fat weight was significantly different between the 2 lines; however, there was no difference for body weight between the lean and fat lines. A total of 132 haplotype windows were significantly associated with abdominal fat weight. These significantly associated haplotype windows were primarily located on chromosomes 2, 4, 8, 10, and 26. Seven candidate genes, including SHH, LMBR1, FGF7, IL16, PLIN1, IGF1R, and SLC16A1, were located within these associated regions. These genes may play important roles in the control of abdominal fat content. Two regions on chromosomes 3 and 10 were significantly associated with testis weight. These 2 regions were previously detected by the single SNP GWAS using this same resource population. TCF21 on chromosome 3 was identified as a potentially important candidate gene for testis growth and development based on gene expression analysis and the reported function of this gene. TCF12, which was previously detected in our SNP by SNP interaction analysis, was located in a region on chromosome 10 that was significantly associated with testis weight. Six candidate genes, including TNFRSF1B, PLOD1, NPPC, MTHFR, EPHB2, and SLC35A3, on chromosome 21 may play important roles in bone development based on the known function of these genes. In addition, several regions were significantly associated with other carcass and growth traits, but no candidate genes were identified. The results of the present study may be helpful in understanding the genetic mechanisms of carcass and growth traits in chickens.

A road map for understanding molecular and genetic determinants of osteoporosis

Osteoporosis is a highly prevalent disorder characterized by low bone mineral density and an increased risk of fracture, termed osteoporotic fracture. Notably, bone mineral density, osteoporosis and osteoporotic fracture are highly heritable; however, determining the genetic architecture, and especially the underlying genomic and molecular mechanisms, of osteoporosis in vivo in humans is still challenging. In addition to susceptibility loci identified in genome-wide association studies, advances in various omics technologies, including genomics, transcriptomics, epigenomics, proteomics and metabolomics, have all been applied to dissect the pathogenesis of osteoporosis. However, each technology individually cannot capture the entire view of the disease pathology and thus fails to comprehensively identify the underlying pathological molecular mechanisms, especially the regulatory and signalling mechanisms. A change to the status quo calls for integrative multi-omics and inter-omics analyses with approaches in 'systems genetics and genomics'. In this Review, we highlight findings from genome-wide association studies and studies using various omics technologies individually to identify mechanisms of osteoporosis. Furthermore, we summarize current studies of data integration to understand, diagnose and inform the treatment of osteoporosis. The integration of multiple technologies will provide a road map to illuminate the complex pathogenesis of osteoporosis, especially from molecular functional aspects, in vivo in humans.

Genetic Dissection of Femoral and Tibial Microarchitecture

Our understanding of the genetic control of bone strength has relied mainly on estimates of bone mineral density. Here we have mapped genetic factors that influence femoral and tibial microarchitecture using high‐resolution x‐ray computed tomography (8‐μm isotropic voxels) across a family of 61 BXD strains of mice, roughly 10 isogenic cases per strain and balanced by sex. We computed heritabilities for 25 cortical and trabecular traits. Males and females have well‐matched heritabilities, ranging from 0.25 to 0.75. We mapped 16 genetic loci most of which were detected only in females. There is also a bias in favor of loci that control cortical rather than trabecular bone. To evaluate candidate genes, we combined well‐established gene ontologies with bone transcriptome data to compute bone‐enrichment scores for all protein‐coding genes. We aligned candidates with those of human genome‐wide association studies. A subset of 50 strong candidates fell into three categories: (1) experimentally validated genes already known to modulate bone function (Adamts4, Ddr2, Darc, Adam12, Fkbp10, E2f6, Adam17, Grem2, Ifi204); (2) candidates without any experimentally validated function in bone (eg, Greb1, Ifi202b), but linked to skeletal phenotypes in human cohorts; and (3) candidates that have high bone‐enrichment scores, but for which there is not yet any functional link to bone biology or skeletal system disease (including Ifi202b, Ly9, Ifi205, Mgmt, F2rl1, Iqgap2). Our results highlight contrasting genetic architecture between sexes and among major bone compartments. The alignment of murine and human data facilitates function analysis and should prove of value for preclinical testing of molecular control of bone structure. © 2019 The Authors. JBMR Plus published by Wiley Periodicals, Inc. on behalf of American Society for Bone and Mineral Research.

The genetic architecture of osteoporosis and fracture risk

Osteoporosis and fracture risk are common complex diseases, caused by an interaction of numerous disease susceptibility genes and environmental factors. With the advances in genomic technologies, large-scale genome-wide association studies (GWAS) have been performed which have broadened our understanding of the genetic architecture and biological mechanisms of complex disease. Currently, more than ~90 loci have been found associated with DXA derived bone mineral density (BMD), over ~500 loci with heel estimated BMD and several others with other less widely available bone parameters such as bone geometry, shape, and microarchitecture. Notably, several of the pathways identified by the GWAS efforts correspond to pathways that are currently targeted for the treatment of osteoporosis. Overall, tremendous progress in the field of the genetics of osteoporosis has been achieved with the discovery of WNT16, EN1, DAAM2, and GPC6 among others. Assessment of the function and biological mechanisms of the remaining genes may further untangle the complex genetic landscape of osteoporosis and fracture risk. With this review we aimed to provide a general overview of the existing GWAS studies on osteoporosis traits and fracture risk.

Association between rs12742784 polymorphism and hip fracture, bone mineral density, and EPHB2 mRNA expression levels in elderly Chinese women

Objective: This study aimed to determine the association between rs12742784 polymorphism in the non-coding area and hip fracture, bone mineral density (BMD), and EPHB2 mRNA expression levels in elderly Chinese women.Methods: We investigated 250 Chinese women (mean age: 63.5 ± 8.3 years) including 123 hip fracture patients and 127 non-fracture controls. All participants underwent clinical examination to meet the inclusion criteria. Lumbar and hip BMD were detected by dual-energy X-ray absorptiometry. rs12742784 polymorphism was determined by restriction fragment length polymorphism and EPHB2 mRNA expression levels were measured by real-time polymerase chain reaction.Results: Distribution of rs12742784 genotypes agreed with Hardy-Weinberg equilibrium. The frequency of the CT + TT genotype was significantly associated with decreased risk of hip fracture (adjusted odds ratio = 0.57, p < 0.01) after adjusting for age and body mass index, and with increased BMD and EPHB2 mRNA expression levels. The T allele of the rs12742784 single nucleotide polymorphism (SNP) was a protective factor for hip fracture (adjusted odds ratio = 0.56, p < 0.01).Conclusion: rs12742784 polymorphism was associated with EPHB2 mRNA expression levels, BMD, and hip fracture in Chinese women. The T allele of the rs12742784 SNP was a protective factor for osteoporosis and hip fracture.

Are Bone Mineral Density and Fractures Related to the Incidence and Progression of Radiographic Osteoarthritis of the Knee, Hip, and Hand in Elderly Men and Women? The Rotterdam Study

OBJECTIVE

To examine the longitudinal relationship between bone mineral density (BMD) and the incidence and progression of knee, hip, and hand osteoarthritis (OA), and to examine the relationship between prevalent vertebral and nonvertebral fractures and the incidence and progression of OA in elderly men and women in the Rotterdam Study.

METHODS

Age- and sex-specific quartiles of baseline femoral neck BMD data were constructed for 4,154 subjects. Radiographs were scored for incidence and progression of knee and hip OA, and for incidence of hand OA. Prevalent vertebral fractures were scored using the McCloskey/Kanis method, and prevalent nonvertebral fractures were reported by baseline interview.

RESULTS

Subjects in the highest quartile of femoral neck BMD had an increased risk of incident radiographic knee OA (ROA) (odds ratio [OR] 1.58 [95% confidence interval (95% CI) 1.14-2.18]), and an increased risk of incident hip ROA (OR 1.57 [95% CI 1.06-2.32]), compared to the lowest quartile. No significant relationship was found between high femoral neck BMD and progression of knee or hip ROA or the incidence of hand ROA. Prevalent vertebral and nonvertebral fractures were not related to an increase in the incidence or progression of knee or hip ROA. However, vertebral fractures were associated with incident hand ROA (OR 1.74 [95% CI 1.02-2.98]).

CONCLUSION

Results from the present study confirm earlier findings and thus provide strong evidence that high femoral neck BMD is a prognostic risk factor for the development of knee and hip ROA. Vertebral fractures were found to be a risk factor for incident hand ROA.

Mouse Models and Online Resources for Functional Analysis of Osteoporosis Genome-Wide Association Studies

Osteoporosis is a complex genetic disease in which the number of loci associated with the bone mineral density, a clinical risk factor for fracture, has increased at an exponential rate in the last decade. The identification of the causative variants and candidate genes underlying these loci has not been able to keep pace with the rate of locus discovery. A large number of tools and data resources have been built around the use of the mouse as model of human genetic disease. Herein, we describe resources available for functional validation of human Genome Wide Association Study (GWAS) loci using mouse models. We specifically focus on large-scale phenotyping efforts focused on bone relevant phenotypes and repositories of genotype-phenotype data that exist for transgenic and mutant mice, which can be readily mined as a first step toward more targeted efforts designed to deeply characterize the role of a gene in bone biology.

Dissecting the Genetics of Osteoporosis using Systems Approaches

Osteoporosis is a condition characterized by low bone mineral density (BMD) and an increased risk of fracture. Traits contributing to osteoporotic fracture are highly heritable, indicating that a comprehensive understanding of bone requires a thorough understanding of the genetic basis of bone traits. Towards this goal, genome-wide association studies (GWASs) have identified over 500 loci associated with bone traits. However, few of the responsible genes have been identified, and little is known of how these genes work together to influence systems-level bone function. In this review, we describe how systems genetics approaches can be used to fill these knowledge gaps.

Genetic effects on bone health

PURPOSE OF REVIEW

In recent years, the lower costs of arrays and sequencing technologies, and the better availability of data from genome-wide association studies (GWASs) have led to more reports on genetic factors that are associated with bone health. However, there remains the need for a summary of the newly identified genetic targets that are associated with bone metabolism, and the status of their functional characterization.

RECENT FINDINGS

GWASs revealed dozens of novel genetic loci that are associated with bone mineral density (BMD). Some of these targets have been functionally characterized, although the vast majority have not. Glypican 6, a membrane surface proteoglycan involved in cellular growth control and differentiation, was identified as a novel determinant of BMD and represents a possible drug target for treatment of osteoporosis. Pathway analysis also showed that cell-growth pathways and the SMAD proteins associated with low BMD.

SUMMARY

Hits that were significantly associated with BMD in different studies represent likely candidates (e.g. SOST, WNT16, ESR1 and RANKL) for functional characterization and development of osteoporosis treatments. Indeed, currently available treatment for osteoporosis (antibody against RANKL) appeared a significant target in four recent GWAS studies indicating their applicability and importance for future treatment development.

AKAP11 gene polymorphism is associated with bone mass measured by quantitative ultrasound in young adults

Background: Due to the increased prevalence of osteoporosis and direct health care cost of osteoporosis-related fractures, there is a growing interest in identifying genetic markers associated with osteoporosis phenotypes in order to develop genetic screening strategies. We aimed to analyze the possible associations between calcaneal Quantitative ultrasound (QUS), a valuable screening tool for assessing bone status in clinical practice, and ZBTB40 (rs7524102, rs6426749), SP7 (rs2016266) and AKAP11 (rs9533090) genes. Methods: A cross-sectional study was conducted on 550 healthy individuals of Caucasian ancestry (381 females and 169 males, median age 20.46±2,69). Bone mass was assessed through QUS to determine broadband ultrasound attenuation (BUA, dB/MHz). Single-nucleotide polymorphisms (SNPs) in ZBTB40 (rs7524102, rs6426749), SP7 (rs2016266) and AKAP11 (rs9533090) were selected as genetic markers and genotyped using TaqMan OpenArray® technology. Results: Linear regression analysis revealed that rs7524102 and rs6426749 in ZBTB40, and rs9533090 in AKAP11 were significantly associated with the calcaneal QUS parameter after adjustments for age, sex, weight, height, physical activity, and calcium intake (p=0.038, p=0.012 and p=0.008, respectively). After applying the Bonferroni correction for multiple testing (p=0.012), only the association of rs9533090 in AKAP11 remained significant. Conclusion:AKAP11 gene (rs9533090) influences QUS trait in a population of Caucasian young adults. The rs9533090 SNP may be considered a factor affecting peak bone mass acquisition.

Identification of a novel locus on chromosome 2q13, which predisposes to clinical vertebral fractures independently of bone density

OBJECTIVES

To identify genetic determinants of susceptibility to clinical vertebral fractures, which is an important complication of osteoporosis.

METHODS

Here we conduct a genome-wide association study in 1553 postmenopausal women with clinical vertebral fractures and 4340 controls, with a two-stage replication involving 1028 cases and 3762 controls. Potentially causal variants were identified using expression quantitative trait loci (eQTL) data from transiliac bone biopsies and bioinformatic studies.

RESULTS

A locus tagged by rs10190845 was identified on chromosome 2q13, which was significantly associated with clinical vertebral fracture (P=1.04×10-9) with a large effect size (OR 1.74, 95% CI 1.06 to 2.6). Bioinformatic analysis of this locus identified several potentially functional SNPs that are associated with expression of the positional candidate genes TTL (tubulin tyrosine ligase) and SLC20A1 (solute carrier family 20 member 1). Three other suggestive loci were identified on chromosomes 1p31, 11q12 and 15q11. All these loci were novel and had not previously been associated with bone mineral density or clinical fractures.

CONCLUSION

We have identified a novel genetic variant that is associated with clinical vertebral fractures by mechanisms that are independent of BMD. Further studies are now in progress to validate this association and evaluate the underlying mechanism.

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.

A genome-wide association study identifies nucleotide variants at SIGLEC5 and DEFA1A3 as risk loci for periodontitis

Periodontitis is one of the most common inflammatory diseases, with a prevalence of 11% worldwide for the severe forms and an estimated heritability of 50%. The disease is characterized by destruction of the alveolar bone due to an aberrant host inflammatory response to a dysbiotic oral microbiome. Previous genome-wide association studies (GWAS) have reported several suggestive susceptibility loci. Here, we conducted a GWAS using a German and Dutch case-control sample of aggressive periodontitis (AgP, 896 cases, 7,104 controls), a rare but highly severe and early-onset form of periodontitis, validated the associations in a German sample of severe forms of the more moderate phenotype chronic periodontitis (CP) (993 cases, 1,419 controls). Positive findings were replicated in a Turkish sample of AgP (223 cases, 564 controls). A locus at SIGLEC5 (sialic acid binding Ig-like lectin 5) and a chromosomal region downstream of the DEFA1A3 locus (defensin alpha 1-3) showed association with both disease phenotypes and were associated with periodontitis at a genome-wide significance level in the pooled samples, with P = 1.09E-08 (rs4284742,-G; OR = 1.34, 95% CI = 1.21-1.48) and P = 5.48E-10 (rs2738058,-T; OR = 1.28, 95% CI = 1.18-1.38), respectively. SIGLEC5 is expressed in various myeloid immune cells and classified as an inhibitory receptor with the potential to mediate tyrosine phosphatases SHP-1/-2 dependent signaling. Alpha defensins are antimicrobial peptides with expression in neutrophils and mucosal surfaces and a role in phagocyte-mediated host defense. This study identifies the first shared genetic risk loci of AgP and CP with genome-wide significance and highlights the role of innate and adaptive immunity in the etiology of periodontitis.

Utility of quantitative micro-computed tomographic analysis in zebrafish to define gene function during skeletogenesis

The zebrafish is a powerful experimental model to investigate the genetic and morphologic basis of vertebrate development. Analysis of skeletogenesis in this fish is challenging as a result of the small size of the developing and adult zebrafish. Many of the bones of small fishes such as the zebrafish and medaka are quite thin, precluding many standard assays of bone quality and morphometrics commonly used on bones of larger animals. Microcomputed tomography (microCT) is a common imaging technique used for detailed analysis of the skeleton of the zebrafish and determination of mutant phenotypes. However, the utility of this modality for analysis of the zebrafish skeleton, and the effect of inherent variation among individual zebrafish, including variables such as sex, age and strain, is not well understood. Given the increased use and accessibility of microCT, we set out to define the sensitivity of microCT methods in developing and adult zebrafish. We assessed skeletal shape and density measures in the developing vertebrae and parasphenoid of the skull base. We found most skeletal variables are tightly correlated to standard length, but that at later growth stages (>3months) there are age dependent effects on some skeletal measures. Further we find modest strain but not sex differences in skeletal measures. These data suggest that the appropriate control for assessing mutant phenotypes should be age and strain matched, ideally a wild-type sibling. By analyzing two mutants exhibiting skeletal dysplasia, we show that microCT imaging can be a sensitive method to quantify distinct skeletal parameters of adults. Finally, as developing zebrafish skeletons remain difficult to resolve by radiographic means, we define a contrast agent specific for bone that enhances resolution at early stages, permitting detailed morphometric analysis of the forming skeleton. This increased capability for detection extends the use of this imaging modality to leverage the zebrafish model to understand the development causes of skeletal dysplasias.

Omics analysis of human bone to identify genes and molecular networks regulating skeletal remodeling in health and disease

The skeleton is a metabolically active organ throughout life where specific bone cell activity and paracrine/endocrine factors regulate its morphogenesis and remodeling. In recent years, an increasing number of reports have used multi-omics technologies to characterize subsets of bone biological molecular networks. The skeleton is affected by primary and secondary disease, lifestyle and many drugs. Therefore, to obtain relevant and reliable data from well characterized patient and control cohorts are vital. Here we provide a brief overview of omics studies performed on human bone, of which our own studies performed on trans-iliacal bone biopsies from postmenopausal women with osteoporosis (OP) and healthy controls are among the first and largest. Most other studies have been performed on smaller groups of patients, undergoing hip replacement for osteoarthritis (OA) or fracture, and without healthy controls. The major findings emerging from the combined studies are: 1. Unstressed and stressed bone show profoundly different gene expression reflecting differences in bone turnover and remodeling and 2. Omics analyses comparing healthy/OP and control/OA cohorts reveal characteristic changes in transcriptomics, epigenomics (DNA methylation), proteomics and metabolomics. These studies, together with genome-wide association studies, in vitro observations and transgenic animal models have identified a number of genes and gene products that act via Wnt and other signaling systems and are highly associated to bone density and fracture. Future challenge is to understand the functional interactions between bone-related molecular networks and their significance in OP and OA pathogenesis, and also how the genomic architecture is affected in health and disease.

Using GWAS to identify novel therapeutic targets for osteoporosis

Osteoporosis is a common, increasingly prevalent, global health burden characterized by low bone mineral density (BMD) and increased risk of fracture. Despite its significant impact on human health, there is currently a lack of highly effective treatments free of side effects for osteoporosis. Therefore, a major goal in the field is to identify new drug targets. Genetic discovery has been shown to be effective in the unbiased identification of novel drug targets and genome-wide association studies (GWASs) have begun to provide insight into genetic basis of osteoporosis. Over the last decade, GWASs have led to the identification of ∼100 loci associated with BMD and other bone traits related to risk of fracture. However, there have been limited efforts to identify the causal genes underlying the GWAS loci or the mechanisms by which GWAS loci alter bone physiology. In this review, we summarize the current state of the field and discuss strategies for causal gene discovery and the evidence that the novel genes underlying GWAS loci are likely to be a new source of drug targets.

Quantitative imaging methods in osteoporosis

Osteoporosis is characterized by a decreased bone mass and quality resulting in an increased fracture risk. Quantitative imaging methods are critical in the diagnosis and follow-up of treatment effects in osteoporosis. Prior radiographic vertebral fractures and bone mineral density (BMD) as a quantitative parameter derived from dual-energy X-ray absorptiometry (DXA) are among the strongest known predictors of future osteoporotic fractures. Therefore, current clinical decision making relies heavily on accurate assessment of these imaging features. Further, novel quantitative techniques are being developed to appraise additional characteristics of osteoporosis including three-dimensional bone architecture with quantitative computed tomography (QCT). Dedicated high-resolution (HR) CT equipment is available to enhance image quality. At the other end of the spectrum, by utilizing post-processing techniques such as the trabecular bone score (TBS) information on three-dimensional architecture can be derived from DXA images. Further developments in magnetic resonance imaging (MRI) seem promising to not only capture bone micro-architecture but also characterize processes at the molecular level. This review provides an overview of various quantitative imaging techniques based on different radiological modalities utilized in clinical osteoporosis care and research.