Skeletal Site-Related Variation in Human Trabecular Bone Transcriptome and Signaling

BACKGROUND: The skeletal site-specific influence of multiple genes on bone morphology is recognised, but the question as to how these influences may be exerted at the molecular and cellular level has not been explored. METHODOLOGY: To address this question, we have compared global gene expression profiles of human trabecular bone from two different skeletal sites that experience vastly different degrees of mechanical loading, namely biopsies from iliac crest and lumbar spinal lamina. PRINCIPAL FINDINGS: In the lumbar spine, compared to the iliac crest, the majority of the differentially expressed genes showed significantly increased levels of expression; 3406 transcripts were up- whilst 838 were down-regulated. Interestingly, all gene transcripts that have been recently demonstrated to be markers of osteocyte, as well as osteoblast and osteoclast-related genes, were markedly up-regulated in the spine. The transcriptome data is consistent with osteocyte numbers being almost identical at the two anatomical sites, but suggesting a relatively low osteocyte functional activity in the iliac crest. Similarly, osteoblast and osteoclast expression data suggested similar numbers of the cells, but presented with higher activity in the spine than iliac crest. This analysis has also led to the identification of expression of a number of transcripts, previously known and novel, which to our knowledge have never earlier been associated with bone growth and remodelling. CONCLUSIONS AND SIGNIFICANCE: This study provides molecular evidence explaining anatomical and micro-architectural site-related changes in bone cell function, which is predominantly attributable to alteration in cell transcriptional activity. A number of novel signaling molecules in critical pathways, which have been hitherto not known to be expressed in bone cells of mature vertebrates, were identified

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.Funding: ORCADES was supported by the Chief Scientist Office of the Scottish Government (CZB/4/276, CZB/4/710), the Royal Society, the MRC Human Genetics Unit, Arthritis Research UK and the European Union framework programme 6 EUROSPAN project (contract no. LSHG-CT-2006-018947). DNA extractions were performed at the Wellcome Trust Clinical Research Facility in Edinburgh. We would like to acknowledge the invaluable contributions of Lorraine Anderson and the research nurses in Orkney, the administrative team in Edinburgh and the people of Orkney. CABRIO was supported by the Instituto de Salud Carlos III and Fondos FEDER from the EU (PI 11/1092 and PI12/615). The AOGC study was funded by the Australian National Health and Medical Research Council (Project grant 511132). Lothian Birth Cohort 1921 phenotype collection was supported by the UK’s Biotechnology and Biological Sciences Research Council (BBSRC), The Royal Society and The Chief Scientist Office of the Scottish Government. Phenotype collection in the Lothian Birth Cohort 1936 was supported by Age UK (The Disconnected Mind project). Genotyping of the cohorts was funded by the BBSRC. The work was undertaken by the University of Edinburgh Centre for Cognitive Ageing and Cognitive Epidemiology, part of the cross council Lifelong Health and Wellbeing Initiative (MR/K026992/1). Funding from the BBSRC and Medical Research Council (MRC) is gratefully acknowledged. Research work on Slovenian case and control samples was funded by Slovenian Research Agency (project no. P3-0298 and J3-2330). The Danish National Birth Cohort (DNBC) is a result of major grants from the Danish National Research Foundation, the Danish Pharmacists’Fund, the Egmont Foundation, the March of Dimes Birth Defects Foundation, the Augustinus Foundation and the Health Fund of the Danish Health Insurance Societies. The DNBC biobank is a part of the Danish National Biobank resource, which is supported by the Novo Nordisk Foundation. Dr Bjarke Feenstra is supported by an Oak Foundation Fellowship. The Framingham Study was funded by grants from the US National Institute for Arthritis, Musculoskeletal and Skin Diseases and National Institute on Aging (R01 AR 41398 and R01 AR061162; DPK and R01 AR 050066; DK). The Framingham Heart Study of the National Heart, Lung, and Blood Institute of the National Institutes of Health and Boston University School of Medicine were supported by the National Heart, Lung, and Blood Institute’s Framingham Heart Study (N01-HC-25195) and its contract with Affymetrix, Inc. for genotyping services (N02-HL-6-4278). Analyses reflect intellectual input and resource development from the Framingham Heart Study investigators participating in the SNP Health Association Resource (SHARe) project. A portion of this research was conducted using the Linux Cluster for Genetic Analysis (LinGA-II) funded by the Robert Dawson Evans Endowment of the Department of Medicine at Boston University School of Medicine and Boston Medical Center. This research was performed within the Genetic Factors for Osteoporosis (GEFOS) consortium, funded by the European Commission (HEALTH-F2-2008-201865-GEFOS).Acknowledgments: The authors are grateful to the patients and controls from the different centres who agreed to participate in this study. We would like to thank Ms Dilruba Kabir at the Rheumatology and Bone Disease Unit, CGEM-IGMM, Edinburgh, UK; Mr Matt Sims at the MRC Epidemiology Unit, University of Cambridge, UK; Ms Mila Jhamai and Ms Sarah Higgins at the Genetics Laboratory of Erasmus MC, Rotterdam, The Netherlands; Ms Johanna Hadler, Ms Kathryn A Addison and Ms Karena Pryce of the University of Queensland Centre for Clinical Genomics, Brisbane, Australia, for technical support on the genotyping stage; and Mr Marijn Verkerk and Dr Anis Abuseiris at the Genetics Laboratory of Erasmus MC, Rotterdam, for assistance on the data analysis. We would like to acknowledge the invaluable contributions of Lorraine Anderson and the research nurses in Orkney, the administrative team in Edinburgh and the people of Orkney. We would also like to thank Professor Nick Gilbert and Dr Giovanny Rodriguez-Blanco for their comments and advice on the manuscript preparation. This study makes use of data generated by the Wellcome Trust Case Control Consortium. A full list of the investigators who contributed to the generation of the data is available at www.wtccc.org.uk

Genetic sharing with cardiovascular disease risk factors and diabetes reveals novel bone mineral density loci

Bone Mineral Density (BMD) is a highly heritable trait, but genome-wide association studies have identified few genetic risk factors. Epidemiological studies suggest associations between BMD and several traits and diseases, but the nature of the suggestive comorbidity is still unknown. We used a novel genetic pleiotropy-informed conditional False Discovery Rate (FDR) method to identify single nucleotide polymorphisms (SNPs) associated with BMD by leveraging cardiovascular disease (CVD) associated disorders and metabolic traits. By conditioning on SNPs associated with the CVD-related phenotypes, type 1 diabetes, type 2 diabetes, systolic blood pressure, diastolic blood pressure, high density lipoprotein, low density lipoprotein, triglycerides and waist hip ratio, we identified 65 novel independent BMD loci (26 with femoral neck BMD and 47 with lumbar spine BMD) at conditional FDR &lt; 0.01. Many of the loci were confirmed in genetic expression studies. Genes validated at the mRNA levels were characteristic for the osteoblast/osteocyte lineage, Wnt signaling pathway and bone metabolism. The results provide new insight into genetic mechanisms of variability in BMD, and a better understanding of the genetic underpinnings of clinical comorbidity.</p

Assessment of the genetic and clinical determinants of fracture risk: genome wide association and mendelian randomisation study.

OBJECTIVES: To identify the genetic determinants of fracture risk and assess the role of 15 clinical risk factors on osteoporotic fracture risk. DESIGN: Meta-analysis of genome wide association studies (GWAS) and a two-sample mendelian randomisation approach. SETTING: 25 cohorts from Europe, United States, east Asia, and Australia with genome wide genotyping and fracture data. PARTICIPANTS: A discovery set of 37 857 fracture cases and 227 116 controls; with replication in up to 147 200 fracture cases and 150 085 controls. Fracture cases were defined as individuals (>18 years old) who had fractures at any skeletal site confirmed by medical, radiological, or questionnaire reports. Instrumental variable analyses were performed to estimate effects of 15 selected clinical risk factors for fracture in a two-sample mendelian randomisation framework, using the largest previously published GWAS meta-analysis of each risk factor. RESULTS: Of 15 fracture associated loci identified, all were also associated with bone mineral density and mapped to genes clustering in pathways known to be critical to bone biology (eg, SOST, WNT16, and ESR1) or novel pathways (FAM210A, GRB10, and ETS2). Mendelian randomisation analyses showed a clear effect of bone mineral density on fracture risk. One standard deviation decrease in genetically determined bone mineral density of the femoral neck was associated with a 55% increase in fracture risk (odds ratio 1.55 (95% confidence interval 1.48 to 1.63; P=1.5×10-68). Hand grip strength was inversely associated with fracture risk, but this result was not significant after multiple testing correction. The remaining clinical risk factors (including vitamin D levels) showed no evidence for an effect on fracture. CONCLUSIONS: This large scale GWAS meta-analysis for fracture identified 15 genetic determinants of fracture, all of which also influenced bone mineral density. Among the clinical risk factors for fracture assessed, only bone mineral density showed a major causal effect on fracture. Genetic predisposition to lower levels of vitamin D and estimated calcium intake from dairy sources were not associated with fracture risk

Phenotypic Dissection of Bone Mineral Density Reveals Skeletal Site Specificity and Facilitates the Identification of Novel Loci in the Genetic Regulation of Bone Mass Attainment

Heritability of bone mineral density (BMD) varies across skeletal sites, reflecting different relative contributions of genetic and environmental influences. To quantify the degree to which common genetic variants tag and environmental factors influence BMD, at different sites, we estimated the genetic (rg) and residual (re) correlations between BMD measured at the upper limbs (UL-BMD), lower limbs (LL-BMD) and skull (SK-BMD), using total-body DXA scans of ~4,890 participants recruited by the Avon Longitudinal Study of Parents and their Children (ALSPAC). Point estimates of rg indicated that appendicular sites have a greater proportion of shared genetic architecture (LL-/UL-BMD rg = 0.78) between them, than with the skull (UL-/SK-BMD rg = 0.58 and LL-/SK-BMD rg = 0.43). Likewise, the residual correlation between BMD at appendicular sites (re = 0.55) was higher than the residual correlation between SK-BMD and BMD at appendicular sites (re = 0.20-0.24). To explore the basis fo

Orexin-A and Orexin-B During the Postnatal Development of the Rat Brain

Orexin-A and orexin-B are hypothalamic neuropeptides isolated from a small group of neurons in the hypothalamus, which project their axons to all major parts of the central nervous system. Despite the extensive information about orexin expression and function at different parts of the nervous system in adults, data about the development and maturation of the orexin system in the brain are a bit contradictory and insufficient. A previous study has found expression of orexins in the hypothalamus after postnatal day 15 only, while others report orexins detection at embryonic stages of brain formation. In the present study, we investigated the distribution of orexin-A and orexin-B neuronal cell bodies and fibers in the brain at three different postnatal stages: 1-week-, 2-week-old and adult rats. By means of immunohistochemical techniques, we demonstrated that a small subset of cells in the lateral hypothalamus, and the perifornical and periventricular areas were orexin-A and orexin-B positive not only in 2-week-old and adult rats but also in 1-week-old animals. In addition, orexin-A and orexin-B expressing neuronal varicosities were found in many other brain regions. These results suggest that orexin-A and orexin-B play an important role in the early postnatal brain development. The widespread distribution of orexinergic projections through all these stages may imply an involvement of the two neurotransmitters in a large variety of physiological and behavioral processes also including higher brain functions like learning and memory

Characterisation of CART-containing neurons and cells in the porcine pancreas, gastro-intestinal tract, adrenal and thyroid glands

<p>Abstract</p> <p>Background</p> <p>The peptide CART is widely expressed in central and peripheral neurons, as well as in endocrine cells. Known peripheral sites of expression include the gastrointestinal (GI) tract, the pancreas, and the adrenal glands. In rodent pancreas CART is expressed both in islet endocrine cells and in nerve fibers, some of which innervate the islets. Recent data show that CART is a regulator of islet hormone secretion, and that CART null mutant mice have islet dysfunction. CART also effects GI motility, mainly via central routes. In addition, CART participates in the regulation of the hypothalamus-pituitary-adrenal-axis. We investigated CART expression in porcine pancreas, GI-tract, adrenal glands, and thyroid gland using immunocytochemistry.</p> <p>Results</p> <p>CART immunoreactive (IR) nerve cell bodies and fibers were numerous in pancreatic and enteric ganglia. The majority of these were also VIP IR. The finding of intrinsic CART containing neurons indicates that pancreatic and GI CART IR nerve fibers have an intrinsic origin. No CART IR endocrine cells were detected in the pancreas or in the GI tract. The adrenal medulla harboured numerous CART IR endocrine cells, most of which were adrenaline producing. In addition CART IR fibers were frequently seen in the adrenal cortex and capsule. The capsule also contained CART IR nerve cell bodies. The majority of the adrenal CART IR neuronal elements were also VIP IR. CART IR was also seen in a substantial proportion of the C-cells in the thyroid gland. The majority of these cells were also somatostatin IR, and/or 5-HT IR, and/or VIP IR.</p> <p>Conclusion</p> <p>CART is a major neuropeptide in intrinsic neurons of the porcine GI-tract and pancreas, a major constituent of adrenaline producing adrenomedullary cells, and a novel peptide of the thyroid C-cells. CART is suggested to be a regulatory peptide in the porcine pancreas, GI-tract, adrenal gland and thyroid.</p

Bivariate genome-wide association meta-analysis of pediatric musculoskeletal traits reveals pleiotropic effects at the SREBF1/TOM1L2 locus

Bone mineral density is known to be a heritable, polygenic trait whereas genetic variants contributing to lean mass variation remain largely unknown. We estimated the shared SNP heritability and performed a bivariate GWAS meta-analysis of total-body lean mass (TB-LM) and total-body less head bone mineral density (TBLH-BMD) regions in 10,414 children. The estimated SNP heritability is 43% for TBLH-BMD, and 39% for TB-LM, with a shared genetic component of 43%. We identify variants with pleiotropic effects in eight loci, including seven established bone mineral density loci: _WNT4, GALNT3, MEPE, CPED1/WNT16, TNFSF11, RIN3, and PPP6R3/LRP5_. Variants in the _TOM1L2/SREBF1_ locus exert opposing effects TB-LM and TBLH-BMD, and have a stronger association with the former trait. We show that _SREBF1_ is expressed in murine and human osteoblasts, as well as in human muscle tissue. This is the first bivariate GWAS meta-analysis to demonstrate genetic factors with pleiotropic effects on bone mineral density and lean mass