Genome-wide pleiotropy of osteoporosis-related phenotypes: the Framingham Study

Role of NEL‑like molecule‑1 in osteogenesis/chondrogenesis (Review)

A dynamic balance exists between osteogenesis and osteoclastogenesis in bone tissue, which can lead to several bone diseases, such as osteoporosis, osteoarthritis, bone necrosis and bone defects, in cases of insufficient osteogenesis or excessive osteoclastogenesis. NEL‑like molecule‑1 (NELL‑1) was first discovered in 1999 as an osteogenic factor that can prevent or treat bone diseases by increasing osteogenic levels. To date, research has identified multiple signaling pathways involved in improving osteogenic levels. Furthermore, to apply NELL‑1 in clinical practice, researchers have optimized its osteogenic effect by combining it with other molecules, changing its molecular structure and performing bone tissue engineering. Currently, research on NELL‑1 is gaining increasing attention. In the near future, it will definitely be applied in clinical practice to eliminate diseases. Thus, the present study provides a comprehensive review of NELL‑1 in enhancing osteogenic levels from the perspectives of the molecular mechanism, interactions with other molecules/cells, molecular‑level changes, applications in bone tissue engineering and its expression in tumors, providing a solid theoretical basis for its clinical application.

Neural EGFL like 1 as a novel gene for Trabecular Bone Score in older adults: The Bushehr Elderly Health (BEH) program

Neural EGFL like 1 (NELL-1), is a secreted glycoprotein and stimulates osteogenic cell differentiation and bone mineralization. This study aimed to explore the relationship between NELL-1 and Trabecular Bone Score (TBS) as a novel tool for the evaluation of osteoporosis in an elderly population-based cohort study in Iran. A single-locus analysis was performed on TBS using data from 2,071 participants in the Bushehr Elderly Health (BEH) Program. The study investigated 376 independent single nucleotide polymorphisms (SNPs) within the NELL-1 on chromosome 11p15.1. The association between SNPs and the mean TBS L1 to L4 was analyzed through an additive model. Significant variants in the additive model (PFDR<0.05) were further examined within dominant, recessive, over-dominant, and co-dominant models. Multiple linear regression was employed to assess the relationship between the genetic risk score (GRS) derived from significant SNPs and TBS. Three SNPs within the NELL-1 showed a statistically significant association with TBS after adjusting for age and sex. The associations for rs1901945 (β = 0.013, PFDR = 0.0007), rs1584851 (β = -0.011, PFDR = 0.0003), and rs58028601 (β = 0.011, PFDR = 0.0003) were significant in the additive model. Additionally, significant results were observed for rs1901945 and rs58028601 in the dominant model (P<0.05). The GRS showed a statistically significant relationship with TBS, considering adjustments for age, sex, Body Mass Index, type 2 diabetes, and smoking (β = 0.077, P = 1.7×10-5). This study highlights the association of NELL-1 with TBS, underscoring its potential as a candidate for further research and personalized medicine concerning the impact of this gene on bone quality.

CD300e: Emerging role and mechanism as an immune-activating receptor

As a transmembrane protein, CD300e is primarily expressed in myeloid cells. It belongs to the CD300 glycoprotein family, functioning as an immune-activating receptor. Dysfunction of CD300e has been suggested in many diseases, such as infections, immune disorders, obesity, and diabetes, signifying its potential as a key biomarker for disease diagnosis and treatment. This review is aimed to explore the roles and potential mechanisms of CD300e in regulating oxidative stress, immune cell activation, tissue damage and repair, and lipid metabolism, shedding light on its role as a diagnostic marker or a therapeutic target, particularly for infections and autoimmune disorders.

NELL1 membranous nephropathy: clinical associations provide mechanistic clues

Neural epidermal growth factor-like 1 (NELL1) membranous nephropathy (MN) is notable for its segmental deposit distribution, IgG1 dominant deposits, and comparatively high rate of spontaneous remission. It has been associated with a variety of exposures and secondary conditions, specifically use of thiol-containing medications - including lipoic acid, bucillamine, and tiopronin - as well as traditional indigenous medications (TIM) particularly those with high mercury content, and non-steroid anti-inflammatory drugs (NSAIDs). Malignancies, graft vs. host disease (GVHD), infection, and autoimmune conditions have also been associated with NELL1 MN. Herein, we provide a detailed summary of the clinicopathologic features of NELL1 and associations with underlying conditions, with a focus on treatment and outcomes. Rare cases of dual NELL1 and phospholipase A2 receptor (PLA2R) positive MN are reviewed. Genome-wide association study of NELL1, role of NELL1 in other physiologic and pathologic processes, and connection between NELL1 MN and malignancy with relevance of NELL1 tumor staining are examined. Finally, relationships and potential disease mechanisms of thiol- and mercury- associated NELL1 MN are discussed.

Bisphosphonate conjugation enhances the bone-specificity of NELL-1-based systemic therapy for spaceflight-induced bone loss in mice

Microgravity-induced bone loss results in a 1% bone mineral density loss monthly and can be a mission critical factor in long-duration spaceflight. Biomolecular therapies with dual osteogenic and anti-resorptive functions are promising for treating extreme osteoporosis. We previously confirmed that NELL-like molecule-1 (NELL-1) is crucial for bone density maintenance. We further PEGylated NELL-1 (NELL-polyethylene glycol, or NELL-PEG) to increase systemic delivery half-life from 5.5 to 15.5 h. In this study, we used a bio-inert bisphosphonate (BP) moiety to chemically engineer NELL-PEG into BP-NELL-PEG and specifically target bone tissues. We found conjugation with BP improved hydroxyapatite (HA) binding and protein stability of NELL-PEG while preserving NELL-1's osteogenicity in vitro. Furthermore, BP-NELL-PEG showed superior in vivo bone specificity without observable pathology in liver, spleen, lungs, brain, heart, muscles, or ovaries of mice. Finally, we tested BP-NELL-PEG through spaceflight exposure onboard the International Space Station (ISS) at maximal animal capacity (n = 40) in a long-term (9 week) osteoporosis therapeutic study and found that BP-NELL-PEG significantly increased bone formation in flight and ground control mice without obvious adverse health effects. Our results highlight BP-NELL-PEG as a promising therapeutic to mitigate extreme bone loss from long-duration microgravity exposure and musculoskeletal degeneration on Earth, especially when resistance training is not possible due to incapacity (e.g., bone fracture, stroke).

Deciphering the chromatin spatial organization landscapes during BMMSC differentiation

The differentiation imbalance in bone marrow mesenchymal stem cells (BMMSCs) is critical for the development of bone density diseases as the population ages. BMMSCs are precursor cells for osteoblasts and adipocytes; however, the chromatin organization landscapes during BMMSC differentiation remain elusive. In this study, we systematically delineate the four-dimensional (4D) genome and dynamic epigenetic atlas of BMMSCs by RNA sequencing (RNA-seq), assay for transposase-accessible chromatin sequencing (ATAC-seq), and high-throughput chromosome conformation capture (Hi-C). The structure analyses reveal 17.5% common and 28.5%-30% specific loops among BMMSCs, osteoblasts, and adipocytes. The subsequent correlation of genome-wide association studies (GWAS) and expression quantitative trait locus (eQTL) data with multi-omics analysis reveal 274 genes and 3634 single nucleotide polymorphisms (SNPs) associated with bone degeneration and osteoporosis (OP). We hypothesize that SNP mutations affect transcription factor (TF) binding sites, thereby affecting changes in gene expression. Furthermore, 26 motifs, 260 TFs, and 291 SNPs are identified to affect the eQTL. Among these genes, DAAM2, TIMP2, and TMEM241 were found to be essential for diseases such as bone degeneration and OP and may serve as potential drug targets.

Private Spaceflight: A New Landscape for Dealing with Medical Risk

As NASA and other space agencies make plans to proceed with human exploration missions beyond low earth orbit (LEO), the private sector, including Space X, Virgin Galactic, Blue Origin, Space Adventures and others, echo these plans with initiatives of their own to send humans further into space. Development of more sub-orbital flight opportunities, orbital flight opportunities to LEO and even higher risk endeavors will certainly result in exposure to medical risks for an expanding and heterogeneous population of civilians. To date, a handful of "space tourists" have flown to the International Space Station (ISS), at their own expense, ushering in a new era in which anyone with reasonably good health and even those with physical disability may consider becoming space travelers. Indeed, medical and behavioral issues of healthy, professional astronauts, have not been problematic on short orbital flights. However, recent attempts to test the potential limitations in astronauts on extended duration orbital flights in preparation for future missions beyond LEO raise concern about individual differences in ability to tolerate the hazardous spaceflight environment. Given the rapid development of opportunities for non-professionals and the employees of private companies to travel into space, this is an appropriate time to consider the development of selection strategies for non-government space travelers, including the development of genomic and other modern tools to assess susceptibility to spaceflight risk.

Isolation of NELL 1 Aptamers for Rhabdomyosarcoma Targeting

NELL1 (Neural epidermal growth factor-like (EGFL)-like protein) is an important biomarker associated with tissue and bone development and regeneration. NELL1 upregulation has been linked with metastasis and negative prognosis in rhabdomyosarcoma (RMS). Furthermore, multiple recent studies have also shown the importance of NELL1 in inflammatory bowel disease and membranous nephropathy, amongst other diseases. In this study, several anti-NELL1 DNA aptamers were selected from a randomized ssDNA pool using a fluorescence-guided method and evaluated for their binding affinity and selectivity. Several other methods such as a metabolic assay and confocal microscopy were also applied for the evaluation of the selected aptamers. The top three candidates were evaluated further, and AptNCan3 was shown to have a binding affinity up to 959.2 nM. Selectivity was examined in the RH30 RMS cells that overexpressed NELL1. Both AptNCan2 and AptNCan3 could significantly suppress metabolic activity in RMS cells. AptNCan3 was found to locate on the cell membrane and also on intracellular vesicles, which matched the location of NELL1 shown by antibodies in previous research. These results indicate that the selected anti-NELL1 aptamer showed strong and highly specific binding to NELL1 and therefore has potential to be used for in vitro or in vivo studies and treatments.

NELL-1 in Genome-Wide Association Studies across Human Diseases

Neural epidermal growth factor-like (EGFL)-like protein (NELL)-1 is a potent and key osteogenic factor in the development and regeneration of skeletal tissues. Intriguingly, accumulative data from genome-wide association studies (GWASs) have started unveiling potential broader roles of NELL-1 beyond its functions in bone and cartilage. With exploration of the genetic variants of the entire genome in large-scale disease cohorts, GWASs have been used for establishing the connection between specific single-nucleotide polymorphisms of NELL1, in addition to osteoporosis, metabolic diseases, inflammatory conditions, neuropsychiatric diseases, neurodegenerative disorders, and malignant tumors. This review summarizes the findings from GWASs on the manifestation, significance level, implications on function, and correlation of specific NELL1 single-nucleotide polymorphisms in various disorders in humans. By offering a unique and comprehensive correlation between genetic variants and plausible functions of NELL1 in GWASs, this review illustrates the wide range of potential effects of a single gene on the pathogenesis of multiple disorders in humans.

Inactivation of Nell-1 in Chondrocytes Significantly Impedes Appendicular Skeletogenesis

NELL-1, an osteoinductive protein, has been shown to regulate skeletal ossification. Interestingly, an interstitial 11p14.1-p15.3 deletion involving the Nell-1 gene was recently reported in a patient with short stature and delayed fontanelle closure. Here we sought to define the role of Nell-1 in endochondral ossification by investigating Nell-1-specific inactivation in Col2α1-expressing cell lineages. Nell-1flox/flox ; Col2α1-Cre+ (Nell-1Col2α1 KO) mice were generated for comprehensive analysis. Nell-1Col2α1 KO mice were born alive but displayed subtle femoral length shortening. At 1 and 3 months postpartum, Nell-1 inactivation resulted in dwarfism and premature osteoporotic phenotypes. Specifically, Nell-1Col2α1 KO femurs and tibias exhibited significantly reduced length, bone mineral density (BMD), bone volume per tissue volume (BV/TV), trabecular number/thickness, cortical volume/thickness/density, and increased trabecular separation. The decreased bone formation rate revealed by dynamic histomorphometry was associated with altered numbers and/or function of osteoblasts and osteoclasts. Furthermore, longitudinal observations by in vivo micro-CT showed delayed and reduced mineralization at secondary ossification centers in mutants. Histologically, reduced staining intensities of Safranin O, Col-2, Col-10, and fewer BrdU-positive chondrocytes were observed in thinner Nell-1Col2α1 KO epiphyseal plates along with altered distribution and weaker expression level of Ihh, Patched-1, PTHrP, and PTHrP receptor. Primary Nell-1Col2α1 KO chondrocytes also exhibited decreased proliferation and differentiation, and its downregulated expression of the Ihh-PTHrP signaling molecules can be partially rescued by exogenous Nell-1 protein. Moreover, intranuclear Gli-1 protein and gene expression of the Gli-1 downstream target genes, Hip-1 and N-Myc, were also significantly decreased with Nell-1 inactivation. Notably, the rescue effects were diminished/reduced with application of Ihh signaling inhibitors, cyclopamine or GANT61. Taken together, these findings suggest that Nell-1 is a pivotal modulator of epiphyseal homeostasis and endochondral ossification. The cumulative chondrocyte-specific Nell-1 inactivation significantly impedes appendicular skeletogenesis resulting in dwarfism and premature osteoporosis through inhibiting Ihh signaling and predominantly altering the Ihh-PTHrP feedback loop. © 2018 American Society for Bone and Mineral Research.

Clinical advantages and disadvantages of anabolic bone therapies targeting the WNT pathway

The WNT signalling pathway is a key regulator of bone metabolism, particularly bone formation, which has helped to define the role of osteocytes - the most abundant bone cells - as orchestrators of bone remodelling. Several molecules involved in the control of the WNT signalling pathway have been identified as potential targets for the development of bone-building therapeutics for patients with osteoporosis. Several of these molecules have been investigated in animal models, but only inhibitors of sclerostin (which is produced by osteocytes) have been investigated in phase III clinical studies. Here, we review the rationale for these developments and the specificity and potential off-target actions of WNT-based therapeutics. We also describe the available preclinical and clinical studies and discuss the benefits and risks of using sclerostin inhibitors for the management of patients with osteoporosis.

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.

Bone Strength Estimated by Micro-Finite Element Analysis (µFEA) Is Heritable and Shares Genetic Predisposition With Areal BMD: The Framingham Study

Genetic factors contribute to the risk of bone fractures, partly because of effects on bone strength. High-resolution peripheral quantitative computed tomography (HR-pQCT) estimates bone strength using micro-finite element analysis (µFEA). The goal of this study was to investigate if the bone failure load estimated by HR-pQCT-based µFEA is heritable and to what extent it shares genetic regulation with areal bone mineral density (aBMD). Bone microarchitecture was measured by HR-pQCT at the ultradistal tibia and ultradistal radius in adults from the Framingham Heart Study (n = 1087, mean age 72 years; 57% women). Radial and tibial failure load in compression were estimated by µFEA. Femoral neck (FN) and ultradistal forearm (UD) aBMD were measured by dual-energy X-ray absorptiometry (DXA). Heritability (h2 ) of failure load and aBMD and genetic correlations between them was estimated adjusting for covariates (age and sex). Failure load values at the non-weight-bearing ultradistal radius and at the weight-bearing ultradistal tibia were highly correlated (r = 0.906; p < 0.001). Estimates of h2 adjusted for covariates were 0.522 for the radius and 0.497 for the tibia. Additional adjustment for height did not impact on the h2 results, but adjustment for aBMD at the UD and FN somewhat decreased h2 point estimates: 0.222 and 0.380 for radius and tibia, respectively. In bivariate analysis, there was a high phenotypic and genetic correlation between covariate-adjusted failure load at the radius and UD aBMD (ρP  = 0.826, ρG  = 0.954, respectively), whereas environmental correlations were lower (ρE  = 0.696), all highly significant (p < 0.001). Similar correlations were observed between tibial failure load and femoral neck aBMD (ρP  = 0.577, ρG  = 0.703, both p < 0.001; ρE  = 0.432, p < 0.05). These data from adult members of families from a population-based cohort suggest that bone strength of distal extremities estimated by micro-finite element analysis is heritable and shares some genetic composition with areal BMD, regardless of the skeletal site. © 2017 American Society for Bone and Mineral Research.

The rs3736228 polymorphism in the LRP5 gene is associated with calcaneal ultrasound parameter but not with body composition in a cohort of young Caucasian adults

The aim of the present study was to investigate the possible influence of low-density lipoprotein receptor-related protein 5 (LRP5) and sclerostin (SOST) genes as genetic factors contributing to calcaneal quantitative ultrasound (QUS) and body composition variables in a population of young Caucasian adults. The study population comprised a total of 575 individuals (mean age 20.41years; SD 2.36) whose bone mass was assessed through QUS to determine broadband ultrasound attenuation (BUA, dB/MHz). Body composition measurements were performed using a body composition analyser. Seven single-nucleotide polymorphisms (SNPs) of LRP5 (rs2306862, rs599083, rs556442 and rs3736228) and SOST (rs4792909, rs851054 and rs2023794) were selected as genetic markers and genotyped using TaqMan OpenArray^® technology. Linear regression analysis was used to test the possible association of the tested SNPs with QUS and body composition parameters. Linear regression analysis revealed that the rs3736228 SNP of LPR5 was significantly associated with BUA after adjustment for age, sex, weight, height, physical activity and calcium intake (P = 0.028, β (95% CI) = 0.089 (0.099-1.691). For the remaining SNPs, no significant association with the QUS measurement was observed. Regarding body composition, no significant association was found between LRP5 and SOST polymorphisms and body mass index, total fat mass and total lean mass after adjustment for age and sex as covariates. We concluded that the rs3736228 LRP5 genetic polymorphism influences calcaneal QUS parameter in a population of young Caucasian adults. This finding suggests that LRP5 might be an important genetic marker contributing to bone mass accrual early in life.

Neural EGF-like protein 1 (NELL-1): Signaling crosstalk in mesenchymal stem cells and applications in regenerative medicine

Bone tissue regeneration holds the potential to solve both osteoporosis and large skeletal defects, two problems associated with significant morbidity. The differentiation of mesenchymal stem cells into the osteogenic lineage requires a specific microenvironment and certain osteogenic growth factors. Neural EGF Like-Like molecule 1 (NELL-1) is a secreted glycoprotein that has proven, both in vitro and in vivo, to be a potent osteo-inductive factor. Furthermore, it has been shown to repress adipogenic differentiation and inflammation. NELL-1 can work synergistically with other osteogenic factors such as Bone Morphogenic Protein (BMP) -2 and -9, and has shown promise for use in tissue engineering and as a systemically administered drug for the treatment of osteoporosis. Here we provide a comprehensive up-to-date review on the molecular signaling cascade of NELL-1 in mesenchymal stem cells and potential applications in bone regenerative engineering.

RANKL/RANK/OPG Polymorphisms and Heel Quantitative Ultrasound in Young Adults

BACKGROUND

The receptor activator of the nuclear factor-kappa B ligand (RANKL), the receptor activator of nuclear factor-kappa B (RANK), and the osteoprotegerin (OPG) signaling pathway play an important role in the regulation of bone remodeling and osteoclast differentiation. Quantitative ultrasound (QUS) is a relatively recent and noninvasive method providing structural information on microstructure, bone elasticity, and connectivity. However, in contrast to bone mineral density measurements, the possible association of the RANKL/RANK/OPG pathway with heel QUS has not been analyzed.

OBJECTIVES

The aim of this study was to assess, for the first time, the contribution of the RANKL/RANK/OPG pathway genes in the genetic background of heel QUS parameters.

METHODS

Ten single-nucleotide polymorphisms (SNPs) of RANKL (rs9594759, rs12585014, rs7988338, rs2148073), RANK (rs1805034, rs12458117, rs3018362), and OPG (rs4355801, rs3102735, rs2073618) were selected as genetic markers and genotyped using Open Array technology in 575 self-reported Caucasian individuals aged 18-25. Bone mass in the right calcaneus was estimated with QUS to obtain the broadband ultrasound attenuation (BUA) measurement (dB/MHz). Linear regression analyses were performed to test the possible association between the SNPs and BUA.

RESULTS

Linear regression analysis of all the tested SNPs revealed no significant association with the BUA parameter after adjusting for age, gender, weight, height, physical activity, and calcium intake. The lowest p-value was observed for the rs9594759 RANKL polymorphism and heel QUS (p = .06; b* = -.075, 95% CI [-0.960, 0.028]).

CONCLUSION

Our results suggest that the polymorphism of the RANKL, RANK, and OPG genes does not make a significant genetic contribution to heel ultrasound measurements in a population of young Caucasian adults. Further studies replicating the results in independent populations are needed to support these initial findings.

Heritability and Genetic Correlations for Bone Microarchitecture: The Framingham Study Families

High-resolution peripheral quantitative computed tomography (HR-pQCT) measures bone microarchitecture and volumetric bone mineral density (vBMD), important risk factors for osteoporotic fractures. We estimated the heritability (h² ) of bone microstructure indices and vBMD, measured by HR-pQCT, and genetic correlations (ρ_G ) among them and between them and regional aBMD measured by dual-energy X-ray absorptiometry (DXA), in adult relatives from the Framingham Heart Study. Cortical (Ct) and trabecular (Tb) traits were measured at the distal radius and tibia in up to 1047 participants, and ultradistal radius (UD) aBMD was obtained by DXA. Heritability estimates, adjusted for age, sex, and estrogenic status (in women), ranged from 19.3% (trabecular number) to 82.8% (p < 0.01, Ct.vBMD) in the radius and from 51.9% (trabecular thickness) to 98.3% (cortical cross-sectional area fraction) in the tibia. Additional adjustments for height, weight, and radial aBMD had no major effect on h² estimates. In bivariate analyses, moderate to high genetic correlations were found between radial total vBMD and microarchitecture traits (ρ_G from 0.227 to 0.913), except for cortical porosity. At the tibia, a similar pattern of genetic correlations was observed (ρ_G from 0.274 to 0.948), except for cortical porosity. Environmental correlations between the microarchitecture traits were also substantial. There were high genetic correlations between UD aBMD and multivariable-adjusted total and trabecular vBMD at the radius (ρ_G  = 0.811 and 0.917, respectively). In summary, in related men and women from a population-based cohort, cortical and trabecular microarchitecture and vBMD at the radius and tibia were heritable and shared some h² with regional aBMD measured by DXA. These findings of high heritability of HR-pQCT traits, with a slight attenuation when adjusting for aBMD, supports further work to identify the specific variants underlying volumetric bone density and fine structure of long bones. Knowledge that some of these traits are genetically correlated can serve to reduce the number of traits for genetic association studies. © 2016 American Society for Bone and Mineral Research.

Efficacy of Intraperitoneal Administration of PEGylated NELL-1 for Bone Formation

Systemically delivered NEL-like molecule-1 (NELL-1), a potent pro-osteogenic protein, promotes bone formation in healthy and osteoporotic mouse models. PEGylation of NELL-1 (NELL-PEG) increases the half-life of the protein in a mouse model without compromising its osteogenic potential, thereby improving its pharmacokinetics upon systemic delivery. This study consists of a twofold approach: a biodistribution test and an in vivo osteogenic potential test. The biodistribution test compared two commonly used administration methods for drug delivery other than intravenous-intraperitoneal (IP) and subcutaneous (SC)-to examine NELL-PEG biodistribution in mice. Compared to a single-dose SC injection (1.25 mg/kg), a single-dose IP administration yielded a higher protein uptake in the targeted bone sites. When the IP injection dose was doubled to 2.5 mg/kg, the protein remained in the femurs, tibias, and vertebrae for up to 72 h. Next, based on the results of the biodistribution study, IP administration was selected to further investigate the in vivo osteogenic effects of weekly NELL-PEG injection (q7d). In vivo, the IP administered NELL-PEG group showed significantly greater bone mineral density, bone volume fraction, and trabecular bone formation in the targeted bone sites compared to the phosphate-buffered saline control. In summary, weekly NELL-PEG injection via IP administration successfully enhanced the overall bone quality. These findings demonstrate that systemic delivery of NELL-PEG via IP administration may serve as an effective osteogenic therapy for preventing and treating osteoporosis.

Vertebral Implantation of NELL-1 Enhances Bone Formation in an Osteoporotic Sheep Model

BACKGROUND

Vertebral compression fractures related to osteoporosis greatly afflict the aging population. One of the most commonly used therapy today is balloon kyphoplasty. However, this treatment is far from ideal and is associated with significant side effects. NELL-1, an osteoinductive factor that possesses both pro-osteogenic and anti-osteoclastic properties, is a promising candidate for an alternative to current treatment modalities. This study utilizes the pro-osteogenic properties of recombinant human NELL-1 (rhNELL-1) in lumbar spine vertebral defect model in osteoporotic sheep.

METHODS

Osteoporosis was induced through ovariectomy, dietary depletion of calcium and vitamin D, and steroid administration. After osteoporotic induction, lumbar vertebral body defect creation was performed. Sheep were randomly implanted with the control vehicle, comprised of hyaluronic acid (HA) with hydroxyapatite-coated β-tricalcium phosphate (β-TCP), or the treatment material of rhNELL-1 protein lyophilized onto β-TCP mixed with HA. Analysis of lumbar spine defect healing was performed by radiographic, histologic, and computer-simulated biomechanical testing.

RESULTS

rhNELL-1 treatment significantly increased lumbar spine bone formation, as determined by bone mineral density, % bone volume, and mean cortical width as assessed by micro-computed tomography. Histological analysis revealed a significant increase in bone area and osteoblast number and decrease in osteoclast number around the implant site. Computer-simulated biomechanical analysis of trabecular bone demonstrated that rhNELL-1-treatment resulted in a significantly more stress-resistant composition.

CONCLUSION

Our findings suggest rhNELL-1-based vertebral implantation successfully improved cortical and cancellous bone regeneration in the lumbar spine of osteoporotic sheep. rhNELL-1-based bone graft substitutes represent a potential new local therapy.

The genetics of bone mass and susceptibility to bone diseases

Osteoporosis is characterized by low bone mass and an increased risk of fracture. Genetic factors, environmental factors and gene-environment interactions all contribute to a person's lifetime risk of developing an osteoporotic fracture. This Review summarizes key advances in understanding of the genetics of bone traits and their role in osteoporosis. Candidate-gene approaches dominated this field 20 years ago, but clinical and preclinical genetic studies published in the past 5 years generally utilize more-sophisticated and better-powered genome-wide association studies (GWAS). High-throughput DNA sequencing, large genomic databases and improved methods of data analysis have greatly accelerated the gene-discovery process. Linkage analyses of single-gene traits that segregate in families with extreme phenotypes have led to the elucidation of critical pathways controlling bone mass. For example, components of the Wnt-β-catenin signalling pathway have been validated (in both GWAS and functional studies) as contributing to various bone phenotypes. These notable advances in gene discovery suggest that the next decade will witness cataloguing of the hundreds of genes that influence bone mass and osteoporosis, which in turn will provide a roadmap for the development of new drugs that target diseases of low bone mass, including osteoporosis.

NELL-1 in the treatment of osteoporotic bone loss

NELL-1 is a secreted, osteoinductive protein whose expression rheostatically controls skeletal ossification. Overexpression of NELL-1 results in craniosynostosis in humans and mice, whereas lack of Nell-1 expression is associated with skeletal undermineralization. Here we show that Nell-1-haploinsufficient mice have normal skeletal development but undergo age-related osteoporosis, characterized by a reduction in osteoblast:osteoclast (OB:OC) ratio and increased bone fragility. Recombinant NELL-1 binds to integrin β1 and consequently induces Wnt/β-catenin signalling, associated with increased OB differentiation and inhibition of OC-directed bone resorption. Systemic delivery of NELL-1 to mice with gonadectomy-induced osteoporosis results in improved bone mineral density. When extended to a large animal model, local delivery of NELL-1 to osteoporotic sheep spine leads to significant increase in bone formation. Altogether, these findings suggest that NELL-1 deficiency plays a role in osteoporosis and demonstrate the potential utility of NELL-1 as a combination anabolic/antiosteoclastic therapeutic for bone loss.

Genetic regulation of bone metabolism in the chicken: similarities and differences to Mammalian systems

Birds have a unique bone physiology, due to the demands placed on them through egg production. In particular their medullary bone serves as a source of calcium for eggshell production during lay and undergoes continuous and rapid remodelling. We take advantage of the fact that bone traits have diverged massively during chicken domestication to map the genetic basis of bone metabolism in the chicken. We performed a quantitative trait locus (QTL) and expression QTL (eQTL) mapping study in an advanced intercross based on Red Junglefowl (the wild progenitor of the modern domestic chicken) and White Leghorn chickens. We measured femoral bone traits in 456 chickens by peripheral computerised tomography and femoral gene expression in a subset of 125 females from the cross with microarrays. This resulted in 25 loci for female bone traits, 26 loci for male bone traits and 6318 local eQTL loci. We then overlapped bone and gene expression loci, before checking for an association between gene expression and trait values to identify candidate quantitative trait genes for bone traits. A handful of our candidates have been previously associated with bone traits in mice, but our results also implicate unexpected and largely unknown genes in bone metabolism. In summary, by utilising the unique bone metabolism of an avian species, we have identified a number of candidate genes affecting bone allocation and metabolism. These findings can have ramifications not only for the understanding of bone metabolism genetics in general, but could also be used as a potential model for osteoporosis as well as revealing new aspects of vertebrate bone regulation or features that distinguish avian and mammalian bone.

In this work we seek to further the understanding of bone genetics by mapping bone traits and gene expression in the chicken. Bone in female birds is special due to egg production. In this study, we combine the genetic mapping of bone traits with bone gene expression to find candidate quantitative trait genes that explain the differences between wild and domestic chickens in terms of bone production. The concept of combining genetic mapping and gene expression mapping is not new, and has already been successful in isolating bone-related genes in mammals, however this is the first time it has been applied to an avian system with such unique bone modelling processes. We aim to reveal new molecular mechanisms of bone regulation, and many of the candidates we find are new, highlighting the potential this technique has to identify the potential differences between avian and mammalian bone biology.

Pharmacokinetics and osteogenic potential of PEGylated NELL-1 in vivo after systemic administration

Osteoporosis is a skeletal disorder attributable to an imbalance in osteoblast and osteoclast activity. NELL-1, a secretory protein that promotes osteogenesis while suppressing osteoclastic activity, holds potential as an osteoporosis therapy. Recently, we demonstrated that PEGylation of NELL-1 significantly improves its thermostability while preserving its bioactivity in vitro. However, the effect of PEGylation on the pharmacokinetics and osteogenic potential of NELL-1 in vivo have yet to be investigated. The present study demonstrated that PEGylation of NELL-1 significantly increases the elimination half-life time of the protein from 5.5 h to 15.5 h while distributing more than 2-3 times the amount of protein to bone tissues (femur, tibia, vertebrae, calvaria) in vivo when compared to naked NELL-1. In addition, microCT and DXA analyses demonstrated that systemic NELL-PEG therapy administered every 4 or 7 days significantly increases not only femoral and lumbar BMD and percent bone volume, but also new bone formation throughout the overall skeleton after four weeks of treatment. Furthermore, immunohistochemistry revealed increased osteocalcin expression, while TRAP staining showed reduced osteoclast numbers in NELL-PEG groups. Our findings suggest that the PEGylation technique presents a viable and promising approach to further develop NELL-1 into an effective systemic therapeutic for the treatment of osteoporosis.

Bioactivity and circulation time of PEGylated NELL-1 in mice and the potential for osteoporosis therapy

Osteoporosis is a progressive bone disease due to low osteoblast activity and/or high osteoclast activity. NELL-1 is a potential therapy for osteoporosis because it specifically increases osteoblast differentiation. However, similar to other protein drugs, the bioavailability of NELL-1 may be limited by its in vivo half-life and rapid clearance from body. The purpose of the present study is to prolong NELL-1 circulation time in vivo by PEGylation with three monomeric PEG sizes (5, 20, 40 kDa). While linear PEG 5k yielded the most efficient PEGylation and the most thermally stable conjugate, linear PEG 20k resulted in the conjugate with the highest Mw and longest in vivo circulation. Compared to non-modified NELL-1, all three PEGylated conjugates showed enhanced thermal stability and each prolonged the in vivo circulation time significantly. Furthermore, PEGylated NELL-1 retained its osteoblastic activity without any appreciable cytotoxicity. These findings motivate further studies to evaluate the efficacy of PEGylated NELL-1 on the prevention and treatment of osteoporosis.

Heritability of bone mineral density in a multivariate family-based study

There is evidence for a genetic contribution to bone mineral density (BMD×). Different loci affecting BMD have been identified by diverse linkage and genome-wide association studies. We studied the heritability of and the correlations among six densitometric phenotypes and four bone mass/fracture phenotypes. For this purpose, we used a family-based study of the genetics of osteoporosis, the Genetic Analysis of Osteoporosis Project. The primary aim of our study was to examine the roles of genetic and environmental factors in determining osteoporosis-related phenotypes. The project consisted of 11 extended families from Spain. All of them were selected through a proband with osteoporosis. BMD was measured using dual-energy X-ray absorptiometry. The proportion of variance of BMD attributable to significant covariates ranged from 25% (for femoral neck BMD) to 48% (for whole-body total BMD). The vast majority of the densitometric phenotypes had highly significant heritability, ranging from 0.252 (whole-body total BMD) to 0.537 (trochanteric BMD) after correcting for covariate effects. All of the densitometric phenotypes showed high and significant genetic correlations (from -0.772 to -1.000) with a low bone mass/osteopenia condition (Affected 3). Our findings provide additional evidence on the heritability of BMD and a strong genetic correlation between BMD and bone mass/fracture phenotypes in a Spanish population. Our results emphasize the importance of detecting genetic risk factors and the benefit of early diagnosis and especially therapeutic and preventive strategies.

Endocrine crosstalk between muscle and bone

The musculoskeletal system is a complex organ comprised of the skeletal bones, skeletal muscles, tendons, ligaments, cartilage, joints, and other connective tissue that physically and mechanically interact to provide animals and humans with the essential ability of locomotion. This mechanical interaction is undoubtedly essential for much of the diverse shape and forms observed in vertebrates and even in invertebrates with rudimentary musculoskeletal systems such as fish. It makes sense from a historical point of view that the mechanical theories of musculoskeletal development have had tremendous influence of our understanding of biology, because these relationships are clear and palpable. Less visible to the naked eye or even to the microscope is the biochemical interaction among the individual players of the musculoskeletal system. It was only in recent years that we have begun to appreciate that beyond this mechanical coupling of muscle and bones, these 2 tissues function at a higher level through crosstalk signaling mechanisms that are important for the function of the concomitant tissue. Our brief review attempts to present some of the key concepts of these new concepts and is outline to present muscles and bones as secretory/endocrine organs, the evidence for mutual genetic and tissue interactions, pathophysiological examples of crosstalk, and the exciting new directions for this promising field of research aimed at understanding the biochemical/molecular coupling of these 2 intimately associated tissues.

A multi-trait, meta-analysis for detecting pleiotropic polymorphisms for stature, fatness and reproduction in beef cattle

Polymorphisms that affect complex traits or quantitative trait loci (QTL) often affect multiple traits. We describe two novel methods (1) for finding single nucleotide polymorphisms (SNPs) significantly associated with one or more traits using a multi-trait, meta-analysis, and (2) for distinguishing between a single pleiotropic QTL and multiple linked QTL. The meta-analysis uses the effect of each SNP on each of n traits, estimated in single trait genome wide association studies (GWAS). These effects are expressed as a vector of signed t-values (t) and the error covariance matrix of these t values is approximated by the correlation matrix of t-values among the traits calculated across the SNP (V). Consequently, t'V⁻¹t is approximately distributed as a chi-squared with n degrees of freedom. An attractive feature of the meta-analysis is that it uses estimated effects of SNPs from single trait GWAS, so it can be applied to published data where individual records are not available. We demonstrate that the multi-trait method can be used to increase the power (numbers of SNPs validated in an independent population) of GWAS in a beef cattle data set including 10,191 animals genotyped for 729,068 SNPs with 32 traits recorded, including growth and reproduction traits. We can distinguish between a single pleiotropic QTL and multiple linked QTL because multiple SNPs tagging the same QTL show the same pattern of effects across traits. We confirm this finding by demonstrating that when one SNP is included in the statistical model the other SNPs have a non-significant effect. In the beef cattle data set, cluster analysis yielded four groups of QTL with similar patterns of effects across traits within a group. A linear index was used to validate SNPs having effects on multiple traits and to identify additional SNPs belonging to these four groups.

We describe novel methods for finding significant associations between a genome wide panel of SNPs and multiple complex traits, and further for distinguishing between genes with effects on multiple traits and multiple linked genes affecting different traits. The method uses a meta-analysis based on estimates of SNP effects from independent single trait genome wide association studies (GWAS). The method could therefore be widely used to combine already published GWAS results. The method was applied to 32 traits that describe growth, body composition, feed intake and reproduction in 10,191 beef cattle genotyped for approximately 700,000 SNP. The genes found to be associated with these traits can be arranged into 4 groups that differ in their pattern of effects and hence presumably in their physiological mechanism of action. For instance, one group of genes affects weight and fatness in the opposite direction and can be described as a group of genes affecting mature size, while another group affects weight and fatness in the same direction.

Genetic determinants of heel bone properties: genome-wide association meta-analysis and replication in the GEFOS/GENOMOS consortium

Quantitative ultrasound of the heel captures heel bone properties that independently predict fracture risk and, with bone mineral density (BMD) assessed by X-ray (DXA), may be convenient alternatives for evaluating osteoporosis and fracture risk. We performed a meta-analysis of genome-wide association (GWA) studies to assess the genetic determinants of heel broadband ultrasound attenuation (BUA; n = 14 260), velocity of sound (VOS; n = 15 514) and BMD (n = 4566) in 13 discovery cohorts. Independent replication involved seven cohorts with GWA data (in silico n = 11 452) and new genotyping in 15 cohorts (de novo n = 24 902). In combined random effects, meta-analysis of the discovery and replication cohorts, nine single nucleotide polymorphisms (SNPs) had genome-wide significant (P < 5 × 10(-8)) associations with heel bone properties. Alongside SNPs within or near previously identified osteoporosis susceptibility genes including ESR1 (6q25.1: rs4869739, rs3020331, rs2982552), SPTBN1 (2p16.2: rs11898505), RSPO3 (6q22.33: rs7741021), WNT16 (7q31.31: rs2908007), DKK1 (10q21.1: rs7902708) and GPATCH1 (19q13.11: rs10416265), we identified a new locus on chromosome 11q14.2 (rs597319 close to TMEM135, a gene recently linked to osteoblastogenesis and longevity) significantly associated with both BUA and VOS (P < 8.23 × 10(-14)). In meta-analyses involving 25 cohorts with up to 14 985 fracture cases, six of 10 SNPs associated with heel bone properties at P < 5 × 10(-6) also had the expected direction of association with any fracture (P < 0.05), including three SNPs with P < 0.005: 6q22.33 (rs7741021), 7q31.31 (rs2908007) and 10q21.1 (rs7902708). In conclusion, this GWA study reveals the effect of several genes common to central DXA-derived BMD and heel ultrasound/DXA measures and points to a new genetic locus with potential implications for better understanding of osteoporosis pathophysiology.

Genome-wide approaches for identifying genetic risk factors for osteoporosis

Osteoporosis, the most common type of bone disease worldwide, is clinically characterized by low bone mineral density (BMD) and increased susceptibility to fracture. Multiple genetic and environmental factors and gene-environment interactions have been implicated in its pathogenesis. Osteoporosis has strong genetic determination, with the heritability of BMD estimated to be as high as 60%. More than 80 genes or genetic variants have been implicated in risk of osteoporosis by hypothesis-free genome-wide studies. However, these genes or genetic variants can only explain a small portion of BMD variation, suggesting that many other genes or genetic variants underlying osteoporosis risk await discovery. Here, we review recent progress in genome-wide studies of osteoporosis and discuss their implications for medicine and the major challenges in the field.

NELL-1 injection maintains long-bone quantity and quality in an ovariectomy-induced osteoporotic senile rat model

Over 10 million Americans have osteoporosis, and is the predominant cause of fractures in the elderly. Treatment of fractures in the setting of osteoporosis is complicated by a suboptimal bone regenerative response due to a decline in the number of osteoblasts, their function, and survival. Consequently, an osteogenic therapeutic to prevent and treat fractures in patients with osteoporosis is needed. Nel-like molecule-1 (NELL-1), a novel osteoinductive growth factor, has been shown to promote bone regeneration. In this study, we aim to demonstrate the capacity of recombinant NELL-1 to prevent ovariectomy (OVX)-induced osteoporosis in a senile rat model. Ten-month-old female Sprague-Dawley rats underwent either sham surgery or OVX. Subsequently, 50 μL of 600 μg/mL NELL-1 lyophilized onto a 0-50-μm tricalcium phosphate (TCP) carrier was injected into the femoral bone marrow cavity while phosphate-buffered saline (PBS) control was injected into the contralateral femur. Our microcomputed tomography results showed that OVX+PBS/TCP control femurs showed a continuous decrease in the bone volume (BV) and bone mineral density (BMD) from 2 to 8 weeks post-OVX. In contrast, OVX+NELL-1/TCP femurs showed resistance to OVX-induced bone resorption showing BV and BMD levels similar to that of SHAM femurs at 8 weeks post-OVX. Histology showed increased endosteal-woven bone, as well as decreased adipocytes in the bone marrow of NELL-1-treated femurs compared to control. NELL-1-treated femurs also showed increased immunostaining for bone differentiation markers osteopontin and osteocalcin. These findings were validated in vitro, in which addition of NELL-1 in OVX bone marrow stem cells resulted in increased osteogenic differentiation. Thus, NELL-1 effectively enhances in situ osteogenesis in the bone marrow, making it potentially useful in the prevention and treatment of osteoporotic fractures.

CLIP Test: a new fast, simple and powerful method to distinguish between linked or pleiotropic quantitative trait loci in linkage disequilibria analysis

An important question arises when mapping quantitative trait loci (QTLs) for genetically correlated traits: is the correlation due to pleiotropy (a single QTL affecting more than one trait) and/or close linkage (different QTLs that are physically close to each other and influence the traits)? In this article, we propose the Close Linkage versus Pleiotropism (CLIP) test, a fast, simple and powerful method to distinguish between these two situations. The CLIP test is based on the comparison of the square of the observed correlation between a combination of apparent effects at the marker level to the minimal value it can take under the pleiotropic assumption. A simulation study was performed to estimate the power and alpha risk of the CLIP test and compare it to a test that evaluated whether the confidence intervals of the two QTLs overlapped or not (CI test). On average, the CLIP test showed a higher power (68%) to detect close-linked QTLs than the CI test (43%) and a same alpha risk (4%).

Assessing bone health in children and adolescents

During normal childhood and adolescence, the skeleton undergoes tremendous change. Utilizing the processes of modeling and remodeling, the skeleton acquires its adult configuration and ultimately achieves peak bone mass. Optimization of peak bone mass requires the proper interaction of environmental, dietary, hormonal, and genetic influences. A variety of acute and chronic conditions, as well as genetic polymorphisms, are associated with reduced bone density, which can lead to an increased risk of fracture both in childhood and later during adulthood. Bone densitometry has an established role in the evaluation of adults with bone disorders, and the development of suitable reference ranges for children now permits the application of this technology to younger individuals. We present a brief overview of the factors that determine bone density and the emerging role of bone densitometry in the assessment of bone mass in growing children and adolescents.

Body size and pubertal development explain ethnic differences in structural geometry at the femur in Asian, Hispanic, and white early adolescent girls living in the U.S

Variation in structural geometry is present in adulthood, but when this variation arises and what influences this variation prior to adulthood remains poorly understood. Ethnicity is commonly the focus of research of skeletal integrity and appears to explain some of the variation in quantification of bone tissue. However, why ethnicity explains variation in skeletal integrity is unclear.

METHODS

Here we examine predictors of bone cross sectional area (CSA) and section modulus (Z), measured using dual-energy X-ray absorptiometry (DXA) and the Advanced Hip Analysis (AHA) program at the narrow neck of the femur in adolescent (9-14 years) girls (n=479) living in the United States who were classified as Asian, Hispanic, or white if the subject was 75% of a given group based on parental reported ethnicity. Protocols for measuring height and weight follow standardized procedures. Total body lean mass (LM) and total body fat mass (FM) were quantified in kilograms using DXA. Total dietary and total dairy calcium intakes from the previous month were estimated by the use of an electronic semi-quantitative food frequency questionnaire (eFFQ). Physical activity was estimated for the previous year by a validated self-administered modifiable activity questionnaire for adolescents with energy expenditure calculated from the metabolic equivalent (MET) values from the Compendium of Physical Activities. Multiple regression models were developed to predict CSA and Z.

RESULTS

Age, time from menarche, total body lean mass (LM), total body fat mass (FM), height, total calcium, and total dairy calcium all shared a significant (p<0.05), positive relationship with CSA. Age, time from menarche, LM, FM, and height shared significant (p<0.05), positive relationships with Z. For both CSA and Z, LM was the most important covariate. Physical activity was not a significant predictor of geometry at the femoral neck (p≥0.339), even after removing LM as a covariate. After adjusting for covariates, ethnicity was not a significant predictor in regression models for CSA and Z.

CONCLUSION

Variability in bone geometry at the narrow neck of the femur is best explained by body size and pubertal maturation. After controlling for these covariates there were no differences in bone geometry between ethnic groups.

Polymorphisms in the inflammatory genes CIITA, CLEC16A and IFNG influence BMD, bone loss and fracture in elderly women

Osteoclast activity and the fine balance between bone formation and resorption is affected by inflammatory factors such as cytokines and T lymphocyte activity, mediated by major histocompatibility complex (MHC) molecules, in turn regulated by the MHC class II transactivator (MHC2TA). We investigated the effect of functional polymorphisms in the MHC2TA gene (CIITA), and two additional genes; C-type lectin domain 16A (CLEC16A), in linkage disequilibrium with CIITA and Interferon-γ (IFNG), an inducer of CIITA; on bone density, bone resorption markers, bone loss and fracture risk in 75 year-old women followed for up to 10 years (OPRA n = 1003) and in young adult women (PEAK-25 n = 999). CIITA was associated with BMD at age 75 (lumbar spine p = 0.011; femoral neck (FN) p = 0.049) and age 80 (total body p = 0.015; total hip p = 0.042; FN p = 0.028). Carriers of the CIITA rs3087456(G) allele had 1.8–3.4% higher BMD and displayed increased rate of bone loss between age 75 and 80 (FN p = 0.013; total hip p = 0.030; total body p = 3.8E⁻⁵). Despite increasing bone loss, the rs3087456(G) allele was protective against incident fracture overall (p = 0.002), osteoporotic fracture and hip fracture. Carriers of CLEC16A and IFNG variant alleles had lower BMD (p<0.05) and ultrasound parameters and a lower risk of incident fracture (CLEC16A, p = 0.011). In 25-year old women, none of the genes were associated with BMD. In conclusion, variation in inflammatory genes CIITA, CLEC-16A and INFG appear to contribute to bone phenotypes in elderly women and suggest a role for low-grade inflammation and MHC class II expression for osteoporosis pathogenesis.

Clinical review: Genome-wide association studies of skeletal phenotypes: what we have learned and where we are headed

CONTEXT

The primary goals of genome-wide association studies (GWAS) are to discover new molecular and biological pathways involved in the regulation of bone metabolism that can be leveraged for drug development. In addition, the identified genetic determinants may be used to enhance current risk factor profiles.

EVIDENCE ACQUISITION

There have been more than 40 published GWAS on skeletal phenotypes, predominantly focused on dual-energy x-ray absorptiometry-derived bone mineral density (BMD) of the hip and spine.

EVIDENCE SYNTHESIS

Sixty-six BMD loci have been replicated across all the published GWAS, confirming the highly polygenic nature of BMD variation. Only seven of the 66 previously reported genes (LRP5, SOST, ESR1, TNFRSF11B, TNFRSF11A, TNFSF11, PTH) from candidate gene association studies have been confirmed by GWAS. Among 59 novel BMD GWAS loci that have not been reported by previous candidate gene association studies, some have been shown to be involved in key biological pathways involving the skeleton, particularly Wnt signaling (AXIN1, LRP5, CTNNB1, DKK1, FOXC2, HOXC6, LRP4, MEF2C, PTHLH, RSPO3, SFRP4, TGFBR3, WLS, WNT3, WNT4, WNT5B, WNT16), bone development: ossification (CLCN7, CSF1, MEF2C, MEPE, PKDCC, PTHLH, RUNX2, SOX6, SOX9, SPP1, SP7), mesenchymal-stem-cell differentiation (FAM3C, MEF2C, RUNX2, SOX4, SOX9, SP7), osteoclast differentiation (JAG1, RUNX2), and TGF-signaling (FOXL1, SPTBN1, TGFBR3). There are still 30 BMD GWAS loci without prior molecular or biological evidence of their involvement in skeletal phenotypes. Other skeletal phenotypes that either have been or are being studied include hip geometry, bone ultrasound, quantitative computed tomography, high-resolution peripheral quantitative computed tomography, biochemical markers, and fractures such as vertebral, nonvertebral, hip, and forearm.

CONCLUSIONS

Although several challenges lie ahead as GWAS moves into the next generation, there are prospects of new discoveries in skeletal biology. This review integrates findings from previous GWAS and provides a roadmap for future directions building on current GWAS successes.

The transcriptional profile of mesenchymal stem cell populations in primary osteoporosis is distinct and shows overexpression of osteogenic inhibitors

Primary osteoporosis is an age-related disease characterized by an imbalance in bone homeostasis. While the resorptive aspect of the disease has been studied intensely, less is known about the anabolic part of the syndrome or presumptive deficiencies in bone regeneration. Multipotent mesenchymal stem cells (MSC) are the primary source of osteogenic regeneration. In the present study we aimed to unravel whether MSC biology is directly involved in the pathophysiology of the disease and therefore performed microarray analyses of hMSC of elderly patients (79–94 years old) suffering from osteoporosis (hMSC-OP). In comparison to age-matched controls we detected profound changes in the transcriptome in hMSC-OP, e.g. enhanced mRNA expression of known osteoporosis-associated genes (LRP5, RUNX2, COL1A1) and of genes involved in osteoclastogenesis (CSF1, PTH1R), but most notably of genes coding for inhibitors of WNT and BMP signaling, such as Sclerostin and MAB21L2. These candidate genes indicate intrinsic deficiencies in self-renewal and differentiation potential in osteoporotic stem cells. We also compared both hMSC-OP and non-osteoporotic hMSC-old of elderly donors to hMSC of ∼30 years younger donors and found that the transcriptional changes acquired between the sixth and the ninth decade of life differed widely between osteoporotic and non-osteoporotic stem cells. In addition, we compared the osteoporotic transcriptome to long term-cultivated, senescent hMSC and detected some signs for pre-senescence in hMSC-OP.

Our results suggest that in primary osteoporosis the transcriptomes of hMSC populations show distinct signatures and little overlap with non-osteoporotic aging, although we detected some hints for senescence-associated changes. While there are remarkable inter-individual variations as expected for polygenetic diseases, we could identify many susceptibility genes for osteoporosis known from genetic studies. We also found new candidates, e.g. MAB21L2, a novel repressor of BMP-induced transcription. Such transcriptional changes may reflect epigenetic changes, which are part of a specific osteoporosis-associated aging process.

Genetic influence on bone mineral density in Korean twins and families: the healthy twin study

Bone mineral density (BMD), a representative marker of osteoporosis risk, is found to be highly heritable in this Korean study, which is very consistent with the findings in Western populations. This finding strongly supports that genetic factors are significant determinants of osteoporosis risk along with individual biological and behavioral factors.

INTRODUCTION

Although genetic factors are known to contribute significantly to variations in BMD in Western populations, such an association has not been fully evaluated in an Asian population. This study was conducted to determine the role of genetic factors on BMD in Korean population.

METHODS

The study participants were 2,728 men and women consisting of 497 monozygotic (MZ) twin pairs, 119 dizygotic (DZ) twin pairs, and 1,496 first-degree relatives from the Healthy Twin Study. BMD was measured using dual-energy X-ray absorptiometry. Quantitative genetic analysis based on a variance decomposition model was performed.

RESULTS

Age and the measured covariates accounted for 17~61% of the variation in BMD, depending on the sites of measurement. After accounting for the covariate effects, the heritability of BMD at the whole body, thoracic and lumbar spine, whole ribs, whole pelvis, whole arms, and whole legs were 0.76, 0.72, 0.73, 0.71, 0.51, and 0.75, respectively. The pair-wise correlation of BMD was the highest within MZ twin pairs, followed by DZ twin pairs, sibling pairs, and parents-child pairs. Cross-trait correlation analysis revealed a positive genetic correlation between BMDs at different sites, ranging from 0.80 (arm and leg BMD) to 0.50 (pelvis and arm BMD).

CONCLUSIONS

The high heritability of BMD in this Korean population similar to those found in Western populations and the significant common genetic basis between BMDs at different sites strongly supports a significant role of genetic determinants on the risk of osteoporosis.

Moving toward System Genetics through Multiple Trait Analysis in Genome-Wide Association Studies

Association studies are a staple of genotype–phenotype mapping studies, whether they are based on single markers, haplotypes, candidate genes, genome-wide genotypes, or whole genome sequences. Although genetic epidemiological studies typically contain data collected on multiple traits which themselves are often correlated, most analyses have been performed on single traits. Here, I review several methods that have been developed to perform multiple trait analysis. These methods range from traditional multivariate models for systems of equations to recently developed graphical approaches based on network theory. The application of network theory to genetics is termed systems genetics and has the potential to address long-standing questions in genetics about complex processes such as coordinate regulation, homeostasis, and pleiotropy.

Insights into the genetics of osteoporosis from recent genome-wide association studies

Osteoporosis, which is characterised by reduced bone mineral density (BMD) and an increased risk of fragility fractures, is the result of a complex interaction between environmental factors and genetic variants that confer susceptibility. Heritability studies have shown that BMD and other osteoporosis-related traits such as ultrasound properties of bone, skeletal geometry and bone turnover have significant inheritable components. Although previous linkage and candidate gene studies have provided few replicated loci for osteoporosis, genome-wide association approaches have produced clear and reproducible findings. To date, 20 genome-wide association studies (GWASs) for osteoporosis and related traits have been conducted, identifying dozens of genes. Further meta-analyses of GWAS data and deep resequencing of rare variants will uncover more novel susceptibility loci and ultimately provide possible therapeutic targets for fracture prevention.

No association between polymorphisms and haplotypes of COL1A1 and COL1A2 genes and osteoporotic fracture in postmenopausal Chinese women

AIM

To study whether genetic polymorphisms of COL1A1 and COL1A2 genes affected the onset of fracture in postmenopausal Chinese women.

METHODS

SNPs in COL1A1 and COL1A2 genes were identified via direct sequencing in 32 unrelated postmenopausal Chinese women. Ten SNPs were genotyped in 1252 postmenopausal Chinese women. The associations were examined using both single-SNP and haplotype tests using logistic regression.

RESULTS

Twenty four (4 novel) and 28 (7 novel) SNPs were identified in COL1A1 and COL1A2 gene, respectively. The distribution frequencies of 2 SNPs in COL1A1 (rs2075554 and rs2586494) and 3 SNPs in COL1A2 (rs42517, rs1801182, and rs42524) were significantly different from those documented for the European Caucasian population. No significant difference was observed between fracture and control groups with respect to allele frequency or genotype distribution in 9 selected SNPs and haplotype. No significant association was found between fragility fracture and each SNP or haplotype. The results remained the same after additional corrections for other risk factors such as weight, height, and bone mineral density.

CONCLUSION

Our results show no association between common genetic variations of COL1A1 and COL1A2 genes and fracture, suggesting the complex genetic background of osteoporotic fractures.

Identification of homogeneous genetic architecture of multiple genetically correlated traits by block clustering of genome-wide associations

Genome-wide association studies (GWAS) using high-density genotyping platforms offer an unbiased strategy to identify new candidate genes for osteoporosis. It is imperative to be able to clearly distinguish signal from noise by focusing on the best phenotype in a genetic study. We performed GWAS of multiple phenotypes associated with fractures [bone mineral density (BMD), bone quantitative ultrasound (QUS), bone geometry, and muscle mass] with approximately 433,000 single-nucleotide polymorphisms (SNPs) and created a database of resulting associations. We performed analysis of GWAS data from 23 phenotypes by a novel modification of a block clustering algorithm followed by gene-set enrichment analysis. A data matrix of standardized regression coefficients was partitioned along both axes--SNPs and phenotypes. Each partition represents a distinct cluster of SNPs that have similar effects over a particular set of phenotypes. Application of this method to our data shows several SNP-phenotype connections. We found a strong cluster of association coefficients of high magnitude for 10 traits (BMD at several skeletal sites, ultrasound measures, cross-sectional bone area, and section modulus of femoral neck and shaft). These clustered traits were highly genetically correlated. Gene-set enrichment analyses indicated the augmentation of genes that cluster with the 10 osteoporosis-related traits in pathways such as aldosterone signaling in epithelial cells, role of osteoblasts, osteoclasts, and chondrocytes in rheumatoid arthritis, and Parkinson signaling. In addition to several known candidate genes, we also identified PRKCH and SCNN1B as potential candidate genes for multiple bone traits. In conclusion, our mining of GWAS results revealed the similarity of association results between bone strength phenotypes that may be attributed to pleiotropic effects of genes. This knowledge may prove helpful in identifying novel genes and pathways that underlie several correlated phenotypes, as well as in deciphering genetic and phenotypic modularity underlying osteoporosis risk.

How pleiotropic genetics of the musculoskeletal system can inform genomics and phenomics of aging

Genetic study can provide insight into the biologic mechanisms underlying inter-individual differences in susceptibility to (or resistance to) organisms' aging. Recent advances in molecular genetics and genetic epidemiology provide the necessary tools to perform a study of the genetic sources of biological aging. However, to be successful, the genetic study of a complex condition requires a heritable phenotype to be developed and validated. Genome-wide association studies offer an unbiased approach to identify new candidate genes for human diseases. It is hypothesized that convergent results from multiple aging-related traits will point out the genes responsible for the general aging of the organism. This perspective focuses on the musculoskeletal aging as an example of an approach to identify a downstream common pathway that summarizes aging processes. Since the musculoskeletal traits are linked to the state of many vital functions, disability, and ultimately survival rates, we postulate that there is significance in studying musculoskeletal aging. Construction of an integrated phenotype of aging can be achieved based on shared genetics among multiple musculoskeletal biomarkers. Valid biomarkers from other systems of the organism should be similarly explored. The new composite aging score needs to be validated by determining whether it predicts all-cause mortality, incidences of major chronic diseases, and disability late in life. Comprehensive databases on biomarkers of musculoskeletal aging in multiple large cohort studies, along with information on various health outcomes, are needed to validate the proposed measure of biological aging.

Refined QTLs of osteoporosis-related traits by linkage analysis with genome-wide SNPs: Framingham SHARe

Genome-wide association studies (GWAS) using high-density array of single-nucleotide polymorphisms (SNPs) offer an unbiased strategy to identify new candidate genes for osteoporosis. We used a subset of autosomal SNPs from the Affymetrix 500K+50K SNP GeneChip marker set to examine genetic linkage with multiple highly heritable osteoporosis-related traits, including BMD of the hip and spine, heel ultrasound (attenuation and speed of sound), and geometric indices of the hip, in two generations from the Framingham Osteoporosis Study. Variance component linkage analysis was performed using normalized residuals (adjusted for age, height, BMI, and estrogen status in women). Multipoint linkage analyses produced LOD scores > or =3.0 for BMD on chromosomes (chr.) 9 and 11 and for ultrasound speed of sound on chr. 5. Hip geometric traits were linked with higher LOD scores, such as with shaft width on chr. 4 (LOD=3.9) and chr. 16 (LOD=3.8) and with shaft section modulus on chr. 22 (LOD=4.0). LOD score > or =5.0 was obtained for femoral neck width on chr. 7. In conclusion, with an SNP-based linkage approach, we identified several novel potential QTLs and confirmed previously identified chromosomal regions linked to bone mass and geometry. Subsequent focus on the spectrum of genetic polymorphisms in these refined regions may contribute to finding variants predisposing to osteoporosis.