1
|
Reppe S, Gundersen S, Sandve GK, Wang Y, Andreassen OA, Medina-Gomez C, Rivadeneira F, Utheim TP, Hovig E, Gautvik KM. Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease. Int J Mol Sci 2024; 25:5554. [PMID: 38791593 PMCID: PMC11121938 DOI: 10.3390/ijms25105554] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 04/11/2024] [Revised: 05/11/2024] [Accepted: 05/13/2024] [Indexed: 05/26/2024] Open
Abstract
Epidemiological evidence suggests existing comorbidity between postmenopausal osteoporosis (OP) and cardiovascular disease (CVD), but identification of possible shared genes is lacking. The skeletal global transcriptomes were analyzed in trans-iliac bone biopsies (n = 84) from clinically well-characterized postmenopausal women (50 to 86 years) without clinical CVD using microchips and RNA sequencing. One thousand transcripts highly correlated with areal bone mineral density (aBMD) were further analyzed using bioinformatics, and common genes overlapping with CVD and associated biological mechanisms, pathways and functions were identified. Fifty genes (45 mRNAs, 5 miRNAs) were discovered with established roles in oxidative stress, inflammatory response, endothelial function, fibrosis, dyslipidemia and osteoblastogenesis/calcification. These pleiotropic genes with possible CVD comorbidity functions were also present in transcriptomes of microvascular endothelial cells and cardiomyocytes and were differentially expressed between healthy and osteoporotic women with fragility fractures. The results were supported by a genetic pleiotropy-informed conditional False Discovery Rate approach identifying any overlap in single nucleotide polymorphisms (SNPs) within several genes encoding aBMD- and CVD-associated transcripts. The study provides transcriptional and genomic evidence for genes of importance for both BMD regulation and CVD risk in a large collection of postmenopausal bone biopsies. Most of the transcripts identified in the CVD risk categories have no previously recognized roles in OP pathogenesis and provide novel avenues for exploring the mechanistic basis for the biological association between CVD and OP.
Collapse
|
2
|
Zheng J, Wheeler E, Pietzner M, Andlauer TFM, Yau MS, Hartley AE, Brumpton BM, Rasheed H, Kemp JP, Frysz M, Robinson J, Reppe S, Prijatelj V, Gautvik KM, Falk L, Maerz W, Gergei I, Peyser PA, Kavousi M, de Vries PS, Miller CL, Bos M, van der Laan SW, Malhotra R, Herrmann M, Scharnagl H, Kleber M, Dedoussis G, Zeggini E, Nethander M, Ohlsson C, Lorentzon M, Wareham N, Langenberg C, Holmes MV, Davey Smith G, Tobias JH. Lowering of Circulating Sclerostin May Increase Risk of Atherosclerosis and Its Risk Factors: Evidence From a Genome-Wide Association Meta-Analysis Followed by Mendelian Randomization. Arthritis Rheumatol 2023; 75:1781-1792. [PMID: 37096546 PMCID: PMC10586470 DOI: 10.1002/art.42538] [Citation(s) in RCA: 6] [Impact Index Per Article: 6.0] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/27/2022] [Revised: 03/22/2023] [Accepted: 04/18/2023] [Indexed: 04/26/2023]
Abstract
OBJECTIVE In this study, we aimed to establish the causal effects of lowering sclerostin, target of the antiosteoporosis drug romosozumab, on atherosclerosis and its risk factors. METHODS A genome-wide association study meta-analysis was performed of circulating sclerostin levels in 33,961 European individuals. Mendelian randomization (MR) was used to predict the causal effects of sclerostin lowering on 15 atherosclerosis-related diseases and risk factors. RESULTS We found that 18 conditionally independent variants were associated with circulating sclerostin. Of these, 1 cis signal in SOST and 3 trans signals in B4GALNT3, RIN3, and SERPINA1 regions showed directionally opposite signals for sclerostin levels and estimated bone mineral density. Variants with these 4 regions were selected as genetic instruments. MR using 5 correlated cis-SNPs suggested that lower sclerostin increased the risk of type 2 diabetes mellitus (DM) (odds ratio [OR] 1.32 [95% confidence interval (95% CI) 1.03-1.69]) and myocardial infarction (MI) (OR 1.35 [95% CI 1.01-1.79]); sclerostin lowering was also suggested to increase the extent of coronary artery calcification (CAC) (β = 0.24 [95% CI 0.02-0.45]). MR using both cis and trans instruments suggested that lower sclerostin increased hypertension risk (OR 1.09 [95% CI 1.04-1.15]), but otherwise had attenuated effects. CONCLUSION This study provides genetic evidence to suggest that lower levels of sclerostin may increase the risk of hypertension, type 2 DM, MI, and the extent of CAC. Taken together, these findings underscore the requirement for strategies to mitigate potential adverse effects of romosozumab treatment on atherosclerosis and its related risk factors.
Collapse
|
3
|
Medina-Gomez C, Mullin BH, Chesi A, Prijatelj V, Kemp JP, Shochat-Carvalho C, Trajanoska K, Wang C, Joro R, Evans TE, Schraut KE, Li-Gao R, Ahluwalia TS, Zillikens MC, Zhu K, Mook-Kanamori DO, Evans DS, Nethander M, Knol MJ, Thorleifsson G, Prokic I, Zemel B, Broer L, McGuigan FE, van Schoor NM, Reppe S, Pawlak MA, Ralston SH, van der Velde N, Lorentzon M, Stefansson K, Adams HHH, Wilson SG, Ikram MA, Walsh JP, Lakka TA, Gautvik KM, Wilson JF, Orwoll ES, van Duijn CM, Bønnelykke K, Uitterlinden AG, Styrkársdóttir U, Akesson KE, Spector TD, Tobias JH, Ohlsson C, Felix JF, Bisgaard H, Grant SFA, Richards JB, Evans DM, van der Eerden B, van de Peppel J, Ackert-Bicknell C, Karasik D, Kague E, Rivadeneira F. Bone mineral density loci specific to the skull portray potential pleiotropic effects on craniosynostosis. Commun Biol 2023; 6:691. [PMID: 37402774 PMCID: PMC10319806 DOI: 10.1038/s42003-023-04869-0] [Citation(s) in RCA: 4] [Impact Index Per Article: 4.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 07/20/2022] [Accepted: 04/25/2023] [Indexed: 07/06/2023] Open
Abstract
Skull bone mineral density (SK-BMD) provides a suitable trait for the discovery of key genes in bone biology, particularly to intramembranous ossification, not captured at other skeletal sites. We perform a genome-wide association meta-analysis (n ~ 43,800) of SK-BMD, identifying 59 loci, collectively explaining 12.5% of the trait variance. Association signals cluster within gene-sets involved in skeletal development and osteoporosis. Among the four novel loci (ZIC1, PRKAR1A, AZIN1/ATP6V1C1, GLRX3), there are factors implicated in intramembranous ossification and as we show, inherent to craniosynostosis processes. Functional follow-up in zebrafish confirms the importance of ZIC1 on cranial suture patterning. Likewise, we observe abnormal cranial bone initiation that culminates in ectopic sutures and reduced BMD in mosaic atp6v1c1 knockouts. Mosaic prkar1a knockouts present asymmetric bone growth and, conversely, elevated BMD. In light of this evidence linking SK-BMD loci to craniofacial abnormalities, our study provides new insight into the pathophysiology, diagnosis and treatment of skeletal diseases.
Collapse
Grants
- UL1 TR000128 NCATS NIH HHS
- U01 AG042124 NIA NIH HHS
- U01 AG042145 NIA NIH HHS
- U01 AG042168 NIA NIH HHS
- U01 AG042140 NIA NIH HHS
- U24 AG051129 NIA NIH HHS
- R01 AR051124 NIAMS NIH HHS
- U01 AG027810 NIA NIH HHS
- U01 AR066160 NIAMS NIH HHS
- MC_UU_00007/10 Medical Research Council
- R01 HD058886 NICHD NIH HHS
- RC2 AR058973 NIAMS NIH HHS
- Wellcome Trust
- M01 RR000240 NCRR NIH HHS
- U01 AG042143 NIA NIH HHS
- UL1 RR026314 NCRR NIH HHS
- U01 AG042139 NIA NIH HHS
- EC | EU Framework Programme for Research and Innovation H2020 | H2020 Priority Excellent Science | H2020 European Research Council (H2020 Excellent Science - European Research Council)
- European Cooperation in Science and Technology (COST)
- Wellcome Trust (Wellcome)
- Department of Health | National Health and Medical Research Council (NHMRC)
- U.S. Department of Health & Human Services | NIH | Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD)
- ZonMw (Netherlands Organisation for Health Research and Development)
- EC | EC Seventh Framework Programm | FP7 Ideas: European Research Council (FP7-IDEAS-ERC - Specific Programme: "Ideas" Implementing the Seventh Framework Programme of the European Community for Research, Technological Development and Demonstration Activities (2007 to 2013))
- Vetenskapsrådet (Swedish Research Council)
- U.S. Department of Health & Human Services | NIH | National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS)
- Gouvernement du Canada | Canadian Institutes of Health Research (Instituts de Recherche en Santé du Canada)
- Nederlandse Organisatie voor Wetenschappelijk Onderzoek (Netherlands Organisation for Scientific Research)
- NCHA (Netherlands Consortium Healthy Ageing) Leiden/ Rotterdam; Dutch Ministry of Economic Affairs, Agriculture and Innovation (project KB-15-004-003); the Research Institute for Diseases in the Elderly [Netherlands] (014-93-015; RIDE2)
- Clinical and Translational Research Center (5-MO1-RR-000240 and UL1 RR-026314); U.S. Department of Health & Human Services | NIH | Eunice Kennedy Shriver National Institute of Child Health and Human Development (NICHD) GrantRecipient="Au50"
- European Commission FP6 STRP grant number 018947 (LSHG-CT-2006-01947); Netherlands Organization for Scientific Research and the Russian Foundation for Basic Research (NWO-RFBR 047.017.043); Netherlands Brain Foundation (project number F2013(1)-28) GrantRecipient="Au40"
- Chief Scientist Office of the Scottish Government (CZB/4/276, CZB/4/710) GrantRecipient="Au28"
- Chief Scientist Office of the Scottish Government (CZB/4/276, CZB/4/710) GrantRecipient="Au38"
- The Pawsey Supercomputing Centre (with Funding from the Australian Government and the Government of Western Australia; PG 16/0162, PG 17/director2025) GrantRecipient="Au45”
- European Commission (EC)
- U.S. Department of Health & Human Services | NIH | National Institute of Arthritis and Musculoskeletal and Skin Diseases (NIAMS);NIH Roadmap for Medical Research [USA]: U01 AG027810, U01 AG042124, U01 AG042139, U01 AG042140, U01 AG042143, U01 AG042145, U01 AG042168, U01 AR066160, and UL1 TR000128 GrantRecipient="Au39”
- Versus Arthritis [USA] 21937 GrantRecipient="Au57”
Collapse
|
4
|
Gautvik KM, Sachse D, Hinton AC, Olstad OK, Kiel DP, Hsu YH, Utheim TP, Lary CW, Reppe S. In silico discovery of blood cell macromolecular associations. BMC Genom Data 2022; 23:57. [PMID: 35879676 PMCID: PMC9317115 DOI: 10.1186/s12863-022-01077-3] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 03/02/2022] [Accepted: 07/13/2022] [Indexed: 11/23/2022] Open
Abstract
Background Physical molecular interactions are the basis of intracellular signalling and gene regulatory networks, and comprehensive, accessible databases are needed for their discovery. Highly correlated transcripts may reflect important functional associations, but identification of such associations from primary data are cumbersome. We have constructed and adapted a user-friendly web application to discover and identify putative macromolecular associations in human peripheral blood based on significant correlations at the transcriptional level. Methods The blood transcriptome was characterized by quantification of 17,328 RNA species, including 341 mature microRNAs in 105 clinically well-characterized postmenopausal women. Intercorrelation of detected transcripts signal levels generated a matrix with > 150 million correlations recognizing the human blood RNA interactome. The correlations with calculated adjusted p-values were made easily accessible by a novel web application. Results We found that significant transcript correlations within the giant matrix reflect experimentally documented interactions involving select ubiquitous blood relevant transcription factors (CREB1, GATA1, and the glucocorticoid receptor (GR, NR3C1)). Their responsive genes recapitulated up to 91% of these as significant correlations, and were replicated in an independent cohort of 1204 individual blood samples from the Framingham Heart Study. Furthermore, experimentally documented mRNAs/miRNA associations were also reproduced in the matrix, and their predicted functional co-expression described. The blood transcript web application is available at http://app.uio.no/med/klinmed/correlation-browser/blood/index.php and works on all commonly used internet browsers. Conclusions Using in silico analyses and a novel web application, we found that correlated blood transcripts across 105 postmenopausal women reflected experimentally proven molecular associations. Furthermore, the associations were reproduced in a much larger and more heterogeneous cohort and should therefore be generally representative. The web application lends itself to be a useful hypothesis generating tool for identification of regulatory mechanisms in complex biological data sets. Supplementary Information The online version contains supplementary material available at 10.1186/s12863-022-01077-3.
Collapse
|
5
|
Gautvik KM, Olstad OK, Raue U, Gautvik VT, Kvernevik KJ, Utheim TP, Ravnum S, Kirkegaard C, Wiig H, Jones G, Pilling LC, Trappe S, Raastad T, Reppe S. Heavy-load exercise in older adults activates vasculogenesis and has a stronger impact on muscle gene expression than in young adults. Eur Rev Aging Phys Act 2022; 19:23. [PMID: 36182918 PMCID: PMC9526277 DOI: 10.1186/s11556-022-00304-1] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Abstract] [Key Words] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/14/2022] [Accepted: 09/19/2022] [Indexed: 11/18/2022] Open
Abstract
BACKGROUND A striking effect of old age is the involuntary loss of muscle mass and strength leading to sarcopenia and reduced physiological functions. However, effects of heavy-load exercise in older adults on diseases and functions as predicted by changes in muscle gene expression have been inadequately studied. METHODS Thigh muscle global transcriptional activity (transcriptome) was analyzed in cohorts of older and younger adults before and after 12-13 weeks heavy-load strength exercise using Affymetrix microarrays. Three age groups, similarly trained, were compared: younger adults (age 24 ± 4 years), older adults of average age 70 years (Oslo cohort) and above 80 years (old BSU cohort). To increase statistical strength, one of the older cohorts was used for validation. Ingenuity Pathway analysis (IPA) was used to identify predicted biological effects of a gene set that changed expression after exercise, and Principal Component Analysis (PCA) was used to visualize differences in muscle gene expressen between cohorts and individual participants as well as overall changes upon exercise. RESULTS Younger adults, showed few transcriptome changes, but a marked, significant impact was observed in persons of average age 70 years and even more so in persons above 80 years. The 249 transcripts positively or negatively altered in both cohorts of older adults (q-value < 0.1) were submitted to gene set enrichment analysis using IPA. The transcripts predicted increase in several aspects of "vascularization and muscle contractions", whereas functions associated with negative health effects were reduced, e.g., "Glucose metabolism disorder" and "Disorder of blood pressure". Several genes that changed expression after intervention were confirmed at the genome level by containing single nucleotide variants associated with handgrip strength and muscle expression levels, e.g., CYP4B1 (p = 9.2E-20), NOTCH4 (p = 9.7E-8), and FZD4 (p = 5.3E-7). PCA of the 249 genes indicated a differential pattern of muscle gene expression in young and elderly. However, after exercise the expression patterns in both young and old BSU cohorts were changed in the same direction for the vast majority of participants. CONCLUSIONS The positive impact of heavy-load strength training on the transcriptome increased markedly with age. The identified molecular changes translate to improved vascularization and muscular strength, suggesting highly beneficial health effects for older adults.
Collapse
|
6
|
Datta HK, Kringen MK, Tuck SP, Salpingidou G, Olstad OK, Gautvik KM, Cockell SJ, Gautvik VT, Prediger M, Wu JJ, Birch MA, Reppe S. Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis. Int J Mol Sci 2022; 23:ijms23062957. [PMID: 35328378 PMCID: PMC8955993 DOI: 10.3390/ijms23062957] [Citation(s) in RCA: 2] [Impact Index Per Article: 1.0] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 01/04/2022] [Revised: 03/02/2022] [Accepted: 03/05/2022] [Indexed: 12/21/2022] Open
Abstract
Mechanical loading exerts a profound influence on bone density and architecture, but the exact mechanism is unknown. Our study shows that expression of the neurological transcriptional factor zinc finger of the cerebellum 1 (ZIC1) is markedly increased in trabecular bone biopsies in the lumbar spine compared with the iliac crest, skeletal sites of high and low mechanical stress, respectively. Human trabecular bone transcriptome analyses revealed a strong association between ZIC1 mRNA levels and gene transcripts characteristically associated with osteoblasts, osteocytes and osteoclasts. This supposition is supported by higher ZIC1 expression in iliac bone biopsies from postmenopausal women with osteoporosis compared with age-matched control subjects, as well as strongly significant inverse correlation between ZIC1 mRNA levels and BMI-adjusted bone mineral density (BMD) (Z-score). ZIC1 promoter methylation was decreased in mechanically loaded vertebral bone compared to unloaded normal iliac bone, and its mRNA levels correlated inversely with ZIC1 promoter methylation, thus linking mechanical stress to epigenetic control of gene expression. The findings were corroborated in cultures of rat osteoblast progenitors and osteoblast-like cells. This study demonstrates for the first time how skeletal epigenetic changes that are affected by mechanical forces give rise to marked alteration in bone cell transcriptional activity and translate to human bone pathophysiology.
Collapse
|
7
|
Gautvik KM, Günther CC, Prijatelj V, Medina-Gomez C, Shevroja E, Rad LH, Yazdani M, Lindalen E, Valland H, Gautvik VT, Olstad OK, Holden M, Rivadeneira F, Utheim TP, Reppe S. Distinct Subsets of Noncoding RNAs Are Strongly Associated With BMD and Fracture, Studied in Weight-Bearing and Non-Weight-Bearing Human Bone. J Bone Miner Res 2020; 35:1065-1076. [PMID: 32017184 DOI: 10.1002/jbmr.3974] [Citation(s) in RCA: 7] [Impact Index Per Article: 1.8] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/26/2019] [Revised: 01/22/2020] [Accepted: 01/26/2020] [Indexed: 12/14/2022]
Abstract
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.
Collapse
|
8
|
Zheng J, Maerz W, Gergei I, Kleber M, Drechsler C, Wanner C, Brandenburg V, Reppe S, Gautvik KM, Medina-Gomez C, Shevroja E, Gilly A, Park YC, Dedoussis G, Zeggini E, Lorentzon M, Henning P, Lerner UH, Nilsson KH, Movérare-Skrtic S, Baird D, Elsworth B, Falk L, Groom A, Capellini TD, Grundberg E, Nethander M, Ohlsson C, Davey Smith G, Tobias JH. Mendelian Randomization Analysis Reveals a Causal Influence of Circulating Sclerostin Levels on Bone Mineral Density and Fractures. J Bone Miner Res 2019; 34:1824-1836. [PMID: 31170332 PMCID: PMC6899787 DOI: 10.1002/jbmr.3803] [Citation(s) in RCA: 19] [Impact Index Per Article: 3.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/06/2018] [Revised: 05/14/2019] [Accepted: 05/23/2019] [Indexed: 12/21/2022]
Abstract
In bone, sclerostin is mainly osteocyte-derived and plays an important local role in adaptive responses to mechanical loading. Whether circulating levels of sclerostin also play a functional role is currently unclear, which we aimed to examine by two-sample Mendelian randomization (MR). A genetic instrument for circulating sclerostin, derived from a genomewide association study (GWAS) meta-analysis of serum sclerostin in 10,584 European-descent individuals, was examined in relation to femoral neck bone mineral density (BMD; n = 32,744) in GEFOS and estimated bone mineral density (eBMD) by heel ultrasound (n = 426,824) and fracture risk (n = 426,795) in UK Biobank. Our GWAS identified two novel serum sclerostin loci, B4GALNT3 (standard deviation [SD]) change in sclerostin per A allele (β = 0.20, p = 4.6 × 10-49 ) and GALNT1 (β = 0.11 per G allele, p = 4.4 × 10-11 ). B4GALNT3 is an N-acetyl-galactosaminyltransferase, adding a terminal LacdiNAc disaccharide to target glycocoproteins, found to be predominantly expressed in kidney, whereas GALNT1 is an enzyme causing mucin-type O-linked glycosylation. Using these two single-nucleotide polymorphisms (SNPs) as genetic instruments, MR revealed an inverse causal relationship between serum sclerostin and femoral neck BMD (β = -0.12, 95% confidence interval [CI] -0.20 to -0.05) and eBMD (β = -0.12, 95% CI -0.14 to -0.10), and a positive relationship with fracture risk (β = 0.11, 95% CI 0.01 to 0.21). Colocalization analysis demonstrated common genetic signals within the B4GALNT3 locus for higher sclerostin, lower eBMD, and greater B4GALNT3 expression in arterial tissue (probability >99%). Our findings suggest that higher sclerostin levels are causally related to lower BMD and greater fracture risk. Hence, strategies for reducing circulating sclerostin, for example by targeting glycosylation enzymes as suggested by our GWAS results, may prove valuable in treating osteoporosis. © 2019 The Authors. Journal of Bone and Mineral Research published by Wiley Periodicals, Inc.
Collapse
|
9
|
Alonso N, Estrada K, Albagha OME, Herrera L, Reppe S, Olstad OK, Gautvik KM, Ryan NM, Evans KL, Nielson CM, Hsu YH, Kiel DP, Markozannes G, Ntzani EE, Evangelou E, Feenstra B, Liu X, Melbye M, Masi L, Brandi ML, Riches P, Daroszewska A, Olmos JM, Valero C, Castillo J, Riancho JA, Husted LB, Langdahl BL, Brown MA, Duncan EL, Kaptoge S, Khaw KT, Usategui-Martín R, Del Pino-Montes J, González-Sarmiento R, Lewis JR, Prince RL, D’Amelio P, García-Giralt N, NoguéS X, Mencej-Bedrac S, Marc J, Wolstein O, Eisman JA, Oei L, Medina-Gómez C, Schraut KE, Navarro P, Wilson JF, Davies G, Starr J, Deary I, Tanaka T, Ferrucci L, Gianfrancesco F, Gennari L, Lucas G, Elosua R, Uitterlinden AG, Rivadeneira F, Ralston SH. Identification of a novel locus on chromosome 2q13, which predisposes to clinical vertebral fractures independently of bone density. Ann Rheum Dis 2018; 77:378-385. [PMID: 29170203 PMCID: PMC5912156 DOI: 10.1136/annrheumdis-2017-212469] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/29/2017] [Accepted: 11/01/2017] [Indexed: 12/20/2022]
Abstract
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.
Collapse
|
10
|
Reppe S, Datta HK, Gautvik KM. Omics analysis of human bone to identify genes and molecular networks regulating skeletal remodeling in health and disease. Bone 2017; 101:88-95. [PMID: 28450214 DOI: 10.1016/j.bone.2017.04.012] [Citation(s) in RCA: 15] [Impact Index Per Article: 2.1] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 12/10/2016] [Revised: 04/13/2017] [Accepted: 04/22/2017] [Indexed: 12/11/2022]
Abstract
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.
Collapse
|
11
|
Reppe S, Lien TG, Hsu YH, Gautvik VT, Olstad OK, Yu R, Bakke HG, Lyle R, Kringen MK, Glad IK, Gautvik KM. Distinct DNA methylation profiles in bone and blood of osteoporotic and healthy postmenopausal women. Epigenetics 2017. [PMID: 28650214 DOI: 10.1080/15592294.2017.1345832] [Citation(s) in RCA: 45] [Impact Index Per Article: 6.4] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/22/2022] Open
Abstract
DNA methylation affects expression of associated genes and may contribute to the missing genetic effects from genome-wide association studies of osteoporosis. To improve insight into the mechanisms of postmenopausal osteoporosis, we combined transcript profiling with DNA methylation analyses in bone. RNA and DNA were isolated from 84 bone biopsies of postmenopausal donors varying markedly in bone mineral density (BMD). In all, 2529 CpGs in the top 100 genes most significantly associated with BMD were analyzed. The methylation levels at 63 CpGs differed significantly between healthy and osteoporotic women at 10% false discovery rate (FDR). Five of these CpGs at 5% FDR could explain 14% of BMD variation. To test whether blood DNA methylation reflect the situation in bone (as shown for other tissues), an independent cohort was selected and BMD association was demonstrated in blood for 13 of the 63 CpGs. Four transcripts representing inhibitors of bone metabolism-MEPE, SOST, WIF1, and DKK1-showed correlation to a high number of methylated CpGs, at 5% FDR. Our results link DNA methylation to the genetic influence modifying the skeleton, and the data suggest a complex interaction between CpG methylation and gene regulation. This is the first study in the hitherto largest number of postmenopausal women to demonstrate a strong association among bone CpG methylation, transcript levels, and BMD/fracture. This new insight may have implications for evaluation of osteoporosis stage and susceptibility.
Collapse
|
12
|
Nielson CM, Liu CT, Smith AV, Ackert-Bicknell CL, Reppe S, Jakobsdottir J, Wassel C, Register TC, Oei L, Alonso N, Oei EH, Parimi N, Samelson EJ, Nalls MA, Zmuda J, Lang T, Bouxsein M, Latourelle J, Claussnitzer M, Siggeirsdottir K, Srikanth P, Lorentzen E, Vandenput L, Langefeld C, Raffield L, Terry G, Cox AJ, Allison MA, Criqui MH, Bowden D, Ikram MA, Mellstrom D, Karlsson MK, Carr J, Budoff M, Phillips C, Cupples LA, Chou WC, Myers RH, Ralston SH, Gautvik KM, Cawthon PM, Cummings S, Karasik D, Rivadeneira F, Gudnason V, Orwoll ES, Harris TB, Ohlsson C, Kiel DP, Hsu YH. Novel Genetic Variants Associated With Increased Vertebral Volumetric BMD, Reduced Vertebral Fracture Risk, and Increased Expression of SLC1A3 and EPHB2. J Bone Miner Res 2016; 31:2085-2097. [PMID: 27476799 PMCID: PMC5477772 DOI: 10.1002/jbmr.2913] [Citation(s) in RCA: 38] [Impact Index Per Article: 4.8] [Reference Citation Analysis] [Abstract] [Key Words] [MESH Headings] [Grants] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/13/2016] [Revised: 06/22/2016] [Accepted: 07/08/2016] [Indexed: 12/26/2022]
Abstract
Genome-wide association studies (GWASs) have revealed numerous loci for areal bone mineral density (aBMD). We completed the first GWAS meta-analysis (n = 15,275) of lumbar spine volumetric BMD (vBMD) measured by quantitative computed tomography (QCT), allowing for examination of the trabecular bone compartment. SNPs that were significantly associated with vBMD were also examined in two GWAS meta-analyses to determine associations with morphometric vertebral fracture (n = 21,701) and clinical vertebral fracture (n = 5893). Expression quantitative trait locus (eQTL) analyses of iliac crest biopsies were performed in 84 postmenopausal women, and murine osteoblast expression of genes implicated by eQTL or by proximity to vBMD-associated SNPs was examined. We identified significant vBMD associations with five loci, including: 1p36.12, containing WNT4 and ZBTB40; 8q24, containing TNFRSF11B; and 13q14, containing AKAP11 and TNFSF11. Two loci (5p13 and 1p36.12) also contained associations with radiographic and clinical vertebral fracture, respectively. In 5p13, rs2468531 (minor allele frequency [MAF] = 3%) was associated with higher vBMD (β = 0.22, p = 1.9 × 10-8 ) and decreased risk of radiographic vertebral fracture (odds ratio [OR] = 0.75; false discovery rate [FDR] p = 0.01). In 1p36.12, rs12742784 (MAF = 21%) was associated with higher vBMD (β = 0.09, p = 1.2 × 10-10 ) and decreased risk of clinical vertebral fracture (OR = 0.82; FDR p = 7.4 × 10-4 ). Both SNPs are noncoding and were associated with increased mRNA expression levels in human bone biopsies: rs2468531 with SLC1A3 (β = 0.28, FDR p = 0.01, involved in glutamate signaling and osteogenic response to mechanical loading) and rs12742784 with EPHB2 (β = 0.12, FDR p = 1.7 × 10-3 , functions in bone-related ephrin signaling). Both genes are expressed in murine osteoblasts. This is the first study to link SLC1A3 and EPHB2 to clinically relevant vertebral osteoporosis phenotypes. These results may help elucidate vertebral bone biology and novel approaches to reducing vertebral fracture incidence. © 2016 American Society for Bone and Mineral Research.
Collapse
|
13
|
Reppe S, Datta H, Gautvik KM. The Influence of DNA Methylation on Bone Cells. Curr Genomics 2016; 16:384-92. [PMID: 27019613 PMCID: PMC4765525 DOI: 10.2174/1389202916666150817202913] [Citation(s) in RCA: 15] [Impact Index Per Article: 1.9] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/03/2014] [Revised: 04/19/2015] [Accepted: 06/26/2015] [Indexed: 01/14/2023] Open
Abstract
DNA methylation in eukaryotes invokes heritable alterations of the of the cytosine base in DNA without changing the underlying genomic DNA sequence. DNA methylation may be modified by environmental exposures as well as gene polymorphisms and may be a mechanistic link between environmental risk factors and the development of disease. In this review, we consider the role of DNA methylation in bone cells (osteoclasts/osteoblasts/osteocytes) and their progenitors with special focus on in vitro and ex vivo analyses. The number of studies on DNA methylation in bone cells is still somewhat limited, nevertheless it is getting increasingly clear that this type of the epigenetic changes is a critical regulator of gene expression. DNA methylation is necessary for proper development and function of bone cells and is accompanied by disease characteristic functional alterations as presently reviewed including postmenopausal osteoporosis and mechanical strain.
Collapse
|
14
|
LeBlanc M, Zuber V, Andreassen BK, Witoelar A, Zeng L, Bettella F, Wang Y, McEvoy LK, Thompson WK, Schork AJ, Reppe S, Barrett-Connor E, Ligthart S, Dehghan A, Gautvik KM, Nelson CP, Schunkert H, Samani NJ, Ridker PM, Chasman DI, Aukrust P, Djurovic S, Frigessi A, Desikan RS, Dale AM, Andreassen OA. Identifying Novel Gene Variants in Coronary Artery Disease and Shared Genes With Several Cardiovascular Risk Factors. Circ Res 2015; 118:83-94. [PMID: 26487741 DOI: 10.1161/circresaha.115.306629] [Citation(s) in RCA: 40] [Impact Index Per Article: 4.4] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 04/07/2014] [Accepted: 10/20/2015] [Indexed: 01/02/2023]
Abstract
RATIONALE Coronary artery disease (CAD) is a critical determinant of morbidity and mortality. Previous studies have identified several cardiovascular disease risk factors, which may partly arise from a shared genetic basis with CAD, and thus be useful for discovery of CAD genes. OBJECTIVE We aimed to improve discovery of CAD genes and inform the pathogenic relationship between CAD and several cardiovascular disease risk factors using a shared polygenic signal-informed statistical framework. METHODS AND RESULTS Using genome-wide association studies summary statistics and shared polygenic pleiotropy-informed conditional and conjunctional false discovery rate methodology, we systematically investigated genetic overlap between CAD and 8 traits related to cardiovascular disease risk factors: low-density lipoprotein cholesterol, high-density lipoprotein cholesterol, triglycerides, type 2 diabetes mellitus, C-reactive protein, body mass index, systolic blood pressure, and type 1 diabetes mellitus. We found significant enrichment of single-nucleotide polymorphisms associated with CAD as a function of their association with low-density lipoprotein, high-density lipoprotein, triglycerides, type 2 diabetes mellitus, C-reactive protein, body mass index, systolic blood pressure, and type 1 diabetes mellitus. Applying the conditional false discovery rate method to the enriched phenotypes, we identified 67 novel loci associated with CAD (overall conditional false discovery rate <0.01). Furthermore, we identified 53 loci with significant effects in both CAD and at least 1 of low-density lipoprotein, high-density lipoprotein, triglycerides, type 2 diabetes mellitus, C-reactive protein, systolic blood pressure, and type 1 diabetes mellitus. CONCLUSIONS The observed polygenic overlap between CAD and cardiometabolic risk factors indicates a pathogenic relation that warrants further investigation. The new gene loci identified implicate novel genetic mechanisms related to CAD.
Collapse
|
15
|
Reppe S, Noer A, Grimholt RM, Halldórsson BV, Medina-Gomez C, Gautvik VT, Olstad OK, Berg JP, Datta H, Estrada K, Hofman A, Uitterlinden AG, Rivadeneira F, Lyle R, Collas P, Gautvik KM. Methylation of bone SOST, its mRNA, and serum sclerostin levels correlate strongly with fracture risk in postmenopausal women. J Bone Miner Res 2015; 30:249-56. [PMID: 25155887 DOI: 10.1002/jbmr.2342] [Citation(s) in RCA: 70] [Impact Index Per Article: 7.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 03/24/2014] [Revised: 07/27/2014] [Accepted: 08/12/2014] [Indexed: 12/31/2022]
Abstract
Inhibition of sclerostin, a glycoprotein secreted by osteocytes, offers a new therapeutic paradigm for treatment of osteoporosis (OP) through its critical role as Wnt/catenin signaling regulator. This study describes the epigenetic regulation of SOST expression in bone biopsies of postmenopausal women. We correlated serum sclerostin to bone mineral density (BMD), fractures, and bone remodeling parameters, and related these findings to epigenetic and genetic disease mechanisms. Serum sclerostin and bone remodeling biomarkers were measured in two postmenopausal groups: healthy (BMD T-score > -1) and established OP (BMD T-score < -2.5, with at least one low-energy fracture). Bone specimens were used to analyze SOST mRNAs, single nucleotide polymorphisms (SNPs), and DNA methylation changes. The SOST gene promoter region showed increased CpG methylation in OP patients (n = 4) compared to age and body mass index (BMI) balanced controls (n = 4) (80.5% versus 63.2%, p = 0.0001) with replication in independent cohorts (n = 27 and n = 36, respectively). Serum sclerostin and bone SOST mRNA expression correlated positively with age-adjusted and BMI-adjusted total hip BMD (r = 0.47 and r = 0.43, respectively; both p < 0.0005), and inversely to serum bone turnover markers. Five SNPs, one of which replicates in an independent population-based genomewide association study (GWAS), showed association with serum sclerostin or SOST mRNA levels under an additive model (p = 0.0016 to 0.0079). Genetic and epigenetic changes in SOST influence its bone mRNA expression and serum sclerostin levels in postmenopausal women. The observations suggest that increased SOST promoter methylation seen in OP is a compensatory counteracting mechanism, which lowers serum sclerostin concentrations and reduces inhibition of Wnt signaling in an attempt to promote bone formation.
Collapse
|
16
|
Benestad HB, Fossum S, Gautvik KM, Storm-Mathisen J, Storm-Mathisen I, Vaage JT. MINNEORD. TIDSSKRIFT FOR DEN NORSKE LEGEFORENING 2015. [DOI: 10.4045/tidsskr.15.0846] [Citation(s) in RCA: 0] [Impact Index Per Article: 0] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/02/2022] Open
|
17
|
Capulli M, Olstad OK, Onnerfjord P, Tillgren V, Muraca M, Gautvik KM, Heinegård D, Rucci N, Teti A. The C-terminal domain of chondroadherin: a new regulator of osteoclast motility counteracting bone loss. J Bone Miner Res 2014; 29:1833-46. [PMID: 24616121 DOI: 10.1002/jbmr.2206] [Citation(s) in RCA: 16] [Impact Index Per Article: 1.6] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 07/24/2013] [Revised: 01/23/2014] [Accepted: 02/06/2014] [Indexed: 11/12/2022]
Abstract
Chondroadherin (CHAD) is a leucine-rich protein promoting cell attachment through binding to integrin α2 β1 and syndecans. We observed that CHAD mRNA and protein were lower in bone biopsies of 50-year-old to 65-year-old osteoporotic women and in bone samples of ovariectomized mice versus gender/age-matched controls, suggesting a role in bone metabolism. By the means of an internal cyclic peptide (cyclicCHAD), we observed that its integrin binding sequence impaired preosteoclast migration through a nitric oxide synthase 2-dependent mechanism, decreasing osteoclastogenesis and bone resorption in a concentration-dependent fashion, whereas it had no effect on osteoblasts. Consistently, cyclicCHAD reduced transcription of two nitric oxide downstream genes, migfilin and vasp, involved in cell motility. Furthermore, the nitric oxide donor, S-nitroso-N-acetyl-D,L-penicillamine, stimulated preosteoclast migration and prevented the inhibitory effect of cyclicCHAD. Conversely, the nitric oxide synthase 2 (NOS2) inhibitor, N5-(1-iminoethyl)-l-ornithine, decreased both preosteoclast migration and differentiation, confirming a role of the nitric oxide pathway in the mechanism of action triggered by cyclicCHAD. In vivo, administration of cyclicCHAD was well tolerated and increased bone volume in healthy mice, with no adverse effect. In ovariectomized mice cyclicCHAD improved bone mass by both a preventive and a curative treatment protocol, with an effect in line with that of the bisphosphonate alendronate, that was mimicked by the NOS2 inhibitor [L-N6-(1-Iminoethyl)-lysine.2 dihydrochloride]. In both mouse models, cyclicCHAD reduced osteoclast and bone resorption without affecting osteoblast parameters and bone formation. In conclusion, CHAD is a novel regulator of bone metabolism that, through its integrin binding domain, inhibits preosteoclast motility and bone resorption, with a potential translational impact for the treatment of osteoporosis.
Collapse
|
18
|
Heier MS, Skinningsrud A, Paus E, Gautvik KM. Increased cerebrospinal fluid levels of nerve cell biomarkers in narcolepsy with cataplexy. Sleep Med 2014; 15:614-8. [PMID: 24784789 DOI: 10.1016/j.sleep.2014.02.005] [Citation(s) in RCA: 12] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Received: 11/01/2013] [Revised: 01/22/2014] [Accepted: 02/05/2014] [Indexed: 01/15/2023]
Abstract
BACKGROUND The association between narcolepsy with cataplexy and the hypocretinergic system in the central nervous system is strong since up to 75-90% of all patients have cerebrospinal fluid (CSF) hypocretin-1 deficiency. The predominant occurrence of HLADQB1*0602 tissue type in narcolepsy patients and recent results from genome-wide association studies suggest an underlying immunological mechanism. The present study was initiated to clarify whether measurement of nerve cell biomarkers in CSF could give additional knowledge of the pathophysiological mechanisms causing narcolepsy with cataplexy. METHODS Two patient groups with narcolepsy, comprising 18 patients with low CSF hypocretin-1 concentrations and typical cataplexy, and 18 patients with normal CSF hypocretin-1 levels and mild cataplexy-like symptoms, were compared to 17 controls. We measured the nerve cell biomarkers beta-amyloid (Aβ42), total tau protein (T-tau), phosphorylated tau (P-tau) and neuron-specific enolase (NSE) in CSF. RESULTS The concentrations of all biomarkers were significantly elevated in both patient groups compared to the controls. The concentration of beta-amyloid was significantly higher in the patient group with normal CSF hypocretin-1 concentration than in those with low concentrations, whereas the other biomarkers showed no difference between the patient groups. CONCLUSION The findings of elevated levels of CSF biomarkers independent of CSF hypocretin-1 reduction may reflect alterations in cell metabolism. The results suggest a more extensive affection of the sleep regulating cellular network, affecting other neuronal sites important in the regulation of sleep, in addition to the hypocretin-producing neurons.
Collapse
|
19
|
Vilming Elgaaen B, Olstad OK, Haug KBF, Brusletto B, Sandvik L, Staff AC, Gautvik KM, Davidson B. Global miRNA expression analysis of serous and clear cell ovarian carcinomas identifies differentially expressed miRNAs including miR-200c-3p as a prognostic marker. BMC Cancer 2014; 14:80. [PMID: 24512620 PMCID: PMC3928323 DOI: 10.1186/1471-2407-14-80] [Citation(s) in RCA: 87] [Impact Index Per Article: 8.7] [Reference Citation Analysis] [Abstract] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/09/2013] [Accepted: 02/07/2014] [Indexed: 12/18/2022] Open
Abstract
BACKGROUND Improved insight into the molecular characteristics of the different ovarian cancer subgroups is needed for developing a more individualized and optimized treatment regimen. The aim of this study was to a) identify differentially expressed miRNAs in high-grade serous ovarian carcinoma (HGSC), clear cell ovarian carcinoma (CCC) and ovarian surface epithelium (OSE), b) evaluate selected miRNAs for association with clinical parameters including survival and c) map miRNA-mRNA interactions. METHODS Differences in miRNA expression between HGSC, CCC and OSE were analyzed by global miRNA expression profiling (Affymetrix GeneChip miRNA 2.0 Arrays, n = 12, 9 and 9, respectively), validated by RT-qPCR (n = 35, 19 and 9, respectively), and evaluated for associations with clinical parameters. For HGSC, differentially expressed miRNAs were linked to differentially expressed mRNAs identified previously. RESULTS Differentially expressed miRNAs (n = 78) between HGSC, CCC and OSE were identified (FDR < 0.01%), of which 18 were validated (p < 0.01) using RT-qPCR in an extended cohort. Compared with OSE, miR-205-5p was the most overexpressed miRNA in HGSC. miR-200 family members and miR-182-5p were the most overexpressed in HGSC and CCC compared with OSE, whereas miR-383 was the most underexpressed. miR-205-5p and miR-200 members target epithelial-mesenchymal transition (EMT) regulators, apparently being important in tumor progression. miR-509-3-5p, miR-509-5p, miR-509-3p and miR-510 were among the strongest differentiators between HGSC and CCC, all being significantly overexpressed in CCC compared with HGSC. High miR-200c-3p expression was associated with poor progression-free (p = 0.031) and overall (p = 0.026) survival in HGSC patients. Interacting miRNA and mRNA targets, including those of a TP53-related pathway presented previously, were identified in HGSC. CONCLUSIONS Several miRNAs differentially expressed between HGSC, CCC and OSE have been identified, suggesting a carcinogenetic role for these miRNAs. miR-200 family members, targeting EMT drivers, were mostly overexpressed in both subgroups, among which miR-200c-3p was associated with survival in HGSC patients. A set of miRNAs differentiates CCC from HGSC, of which miR-509-3-5p and miR-509-5p are the strongest classifiers. Several interactions between miRNAs and mRNAs in HGSC were mapped.
Collapse
|
20
|
Jensen K, Rother MB, Brusletto BS, Olstad OK, Dalsbotten Aass HC, van Zelm MC, Kierulf P, Gautvik KM. Increased ID2 levels in adult precursor B cells as compared with children is associated with impaired Ig locus contraction and decreased bone marrow output. THE JOURNAL OF IMMUNOLOGY 2013; 191:1210-9. [PMID: 23825313 DOI: 10.4049/jimmunol.1203462] [Citation(s) in RCA: 13] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Track Full Text] [Subscribe] [Scholar Register] [Indexed: 02/06/2023]
Abstract
Precursor B cell production from bone marrow in mice and humans declines with age. Because the mechanisms behind are still unknown, we studied five precursor B cell subsets (ProB, PreBI, PreBII large, PreBII small, immature B) and their differentiation-stage characteristic gene expression profiles in healthy individual toddlers and middle-aged adults. Notably, the composition of the precursor B cell compartment did not change with age. The expression levels of several transcripts encoding V(D)J recombination factors were decreased in adults as compared with children: RAG1 expression was significantly reduced in ProB cells, and DNA-PKcs, Ku80, and XRCC4 were decreased in PreBI cells. In contrast, TdT was 3-fold upregulated in immature B cells of adults. Still, N-nucleotides, P-nucleotides, and deletions were similar for IGH and IGK junctions between children and adults. PreBII large cells in adults, but not in children, showed highly upregulated expression of the differentiation inhibitor, inhibitor of DNA binding 2 (ID2), in absence of changes in expression of the ID2-binding partner E2A. Further, we identified impaired Ig locus contraction in adult precursor B cells as a likely mechanism by which ID2-mediated blocking of E2A function results in reduced bone marrow B cell output in adults. The reduced B cell production was not compensated by increased proliferation in adult immature B cells, despite increased Ki67 expression. These findings demonstrate distinct regulatory mechanisms in B cell differentiation between adults and children with a central role for transcriptional regulation of ID2.
Collapse
|
21
|
Reppe S, Sachse D, Olstad OK, Gautvik VT, Sanderson P, Datta HK, Berg JP, Gautvik KM. Identification of transcriptional macromolecular associations in human bone using browser based in silico analysis in a giant correlation matrix. Bone 2013. [PMID: 23195995 DOI: 10.1016/j.bone.2012.11.015] [Citation(s) in RCA: 8] [Impact Index Per Article: 0.7] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Journal Information] [Submit a Manuscript] [Subscribe] [Scholar Register] [Indexed: 11/25/2022]
Abstract
Intracellular signaling is critically dependent on gene regulatory networks comprising physical molecular interactions. Presently, there is a lack of comprehensive databases for most human tissue types to verify such macromolecular interactions. We present a user friendly browser which helps to identify functional macromolecular interactions in human bone as significant correlations at the transcriptional level. The molecular skeletal phenotype has been characterized by transcriptome analysis of iliac crest bone biopsies from 84 postmenopausal women through quantifications of ~23,000 mRNA species. When the signal levels were inter-correlated, an array containing >260 million correlations was generated, thus recognizing the human bone interactome at the RNA level. The matrix correlation and p values were made easily accessible by a freely available online browser. We show that significant correlations within the giant matrix are reproduced in a replica set of 13 male vertebral biopsies. The identified correlations differ somewhat from transcriptional interactions identified in cell culture experiments and transgenic mice, thus demonstrating that care should be taken in extrapolating such results to the in vivo situation in human bone. The current giant matrix and web browser are a valuable tool for easy access to the human bone transcriptome and molecular interactions represented as significant correlations at the RNA-level. The browser and matrix should be a valuable hypothesis generating tool for identification of regulatory mechanisms and serve as a library of transcript relationships in human bone, a relatively inaccessible tissue.
Collapse
|
22
|
Elgaaen BV, Olstad OK, Sandvik L, Odegaard E, Sauer T, Staff AC, Gautvik KM. ZNF385B and VEGFA are strongly differentially expressed in serous ovarian carcinomas and correlate with survival. PLoS One 2012; 7:e46317. [PMID: 23029477 PMCID: PMC3460818 DOI: 10.1371/journal.pone.0046317] [Citation(s) in RCA: 31] [Impact Index Per Article: 2.6] [Reference Citation Analysis] [Abstract] [MESH Headings] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 06/26/2012] [Accepted: 08/29/2012] [Indexed: 01/08/2023] Open
Abstract
BACKGROUND The oncogenesis of ovarian cancer is poorly understood. The aim of this study was to identify mRNAs differentially expressed between moderately and poorly differentiated (MD/PD) serous ovarian carcinomas (SC), serous ovarian borderline tumours (SBOT) and superficial scrapings from normal ovaries (SNO), and to correlate these mRNAs with clinical parameters including survival. METHODS Differences in mRNA expression between MD/PD SC, SBOT and SNO were analyzed by global gene expression profiling (n = 23), validated by RT-qPCR (n = 41) and correlated with clinical parameters. RESULTS Thirty mRNAs differentially expressed between MD/PD SC, SBOT and SNO were selected from the global gene expression analyses, and 21 were verified (p<0.01) by RT-qPCR. Of these, 13 mRNAs were differentially expressed in MD/PD SC compared with SNO (p<0.01) and were correlated with clinical parameters. ZNF385B was downregulated (FC = -130.5, p = 1.2×10(-7)) and correlated with overall survival (p = 0.03). VEGFA was upregulated (FC = 6.1, p = 6.0×10(-6)) and correlated with progression-free survival (p = 0.037). Increased levels of TPX2 and FOXM1 mRNAs (FC = 28.5, p = 2.7×10(-10) and FC = 46.2, p = 5.6×10(-4), respectively) correlated with normalization of CA125 (p = 0.03 and p = 0.044, respectively). Furthermore, we present a molecular pathway for MD/PD SC, including VEGFA, FOXM1, TPX2, BIRC5 and TOP2A, all significantly upregulated and directly interacting with TP53. CONCLUSIONS We have identified 21 mRNAs differentially expressed (p<0.01) between MD/PD SC, SBOT and SNO. Thirteen were differentially expressed in MD/PD SC, including ZNF385B and VEGFA correlating with survival, and FOXM1 and TPX2 with normalization of CA125. We also present a molecular pathway for MD/PD SC.
Collapse
|
23
|
Medina-Gomez C, Kemp JP, Estrada K, Eriksson J, Liu J, Reppe S, Evans DM, Heppe DHM, Vandenput L, Herrera L, Ring SM, Kruithof CJ, Timpson NJ, Zillikens MC, Olstad OK, Zheng HF, Richards JB, St. Pourcain B, Hofman A, Jaddoe VWV, Smith GD, Lorentzon M, Gautvik KM, Uitterlinden AG, Brommage R, Ohlsson C, Tobias JH, Rivadeneira F. Meta-analysis of genome-wide scans for total body BMD in children and adults reveals allelic heterogeneity and age-specific effects at the WNT16 locus. PLoS Genet 2012; 8:e1002718. [PMID: 22792070 PMCID: PMC3390371 DOI: 10.1371/journal.pgen.1002718] [Citation(s) in RCA: 129] [Impact Index Per Article: 10.8] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 12/10/2011] [Accepted: 04/04/2012] [Indexed: 12/31/2022] Open
Abstract
To identify genetic loci influencing bone accrual, we performed a genome-wide association scan for total-body bone mineral density (TB-BMD) variation in 2,660 children of different ethnicities. We discovered variants in 7q31.31 associated with BMD measurements, with the lowest P = 4.1×10−11 observed for rs917727 with minor allele frequency of 0.37. We sought replication for all SNPs located ±500 kb from rs917727 in 11,052 additional individuals from five independent studies including children and adults, together with de novo genotyping of rs3801387 (in perfect linkage disequilibrium (LD) with rs917727) in 1,014 mothers of children from the discovery cohort. The top signal mapping in the surroundings of WNT16 was replicated across studies with a meta-analysis P = 2.6×10−31 and an effect size explaining between 0.6%–1.8% of TB-BMD variance. Conditional analyses on this signal revealed a secondary signal for total body BMD (P = 1.42×10−10) for rs4609139 and mapping to C7orf58. We also examined the genomic region for association with skull BMD to test if the associations were independent of skeletal loading. We identified two signals influencing skull BMD variation, including rs917727 (P = 1.9×10−16) and rs7801723 (P = 8.9×10−28), also mapping to C7orf58 (r2 = 0.50 with rs4609139). Wnt16 knockout (KO) mice with reduced total body BMD and gene expression profiles in human bone biopsies support a role of C7orf58 and WNT16 on the BMD phenotypes observed at the human population level. In summary, we detected two independent signals influencing total body and skull BMD variation in children and adults, thus demonstrating the presence of allelic heterogeneity at the WNT16 locus. One of the skull BMD signals mapping to C7orf58 is mostly driven by children, suggesting temporal determination on peak bone mass acquisition. Our life-course approach postulates that these genetic effects influencing peak bone mass accrual may impact the risk of osteoporosis later in life. Genetic investigations on bone mineral density (BMD) variation in children allow the identification of factors determining peak bone mass and their influence on developing osteoporosis later in life. We ran a genome-wide association study (GWAS) for total body BMD based on 2,660 children of different ethnic backgrounds, followed by replication in an additional 12,066 individuals comprising children, young adults, and elderly populations. Our GWAS meta-analysis identified two independent signals in the 7q31.31 locus, arising from SNPs in the vicinity of WNT16, FAM3C, and C7orf58. These variants were also associated with skull BMD, a skeletal trait with much less environmental influence for which one of the signals displayed age-specific effects. Integration of functional studies in a Wnt16 knockout mouse model and gene expression profiles in human bone tissue provided additional evidence that WNT16 and C7orf58 underlie the described associations. All together our findings demonstrate the relevance of these factors for bone biology, the attainment of peak bone mass, and their likely impact on bone fragility later in life.
Collapse
|
24
|
Heier MS, Jansson TS, Gautvik KM. Cerebrospinal fluid hypocretin 1 deficiency, overweight, and metabolic dysregulation in patients with narcolepsy. J Clin Sleep Med 2012; 7:653-8. [PMID: 22171205 DOI: 10.5664/jcsm.1474] [Citation(s) in RCA: 27] [Impact Index Per Article: 2.3] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/21/2022]
Abstract
STUDY OBJECTIVES The possible relationship between cerebrospinal fluid (CSF) hypocretin and leptin levels, overweight, and association to risk factors for diabetes 2 in narcolepsy with cataplexy were compared to patients with idiopathic hypersomnia and controls. PATIENTS 26 patients with narcolepsy, cataplexy, and hypocretin deficiency; 23 patients with narcolepsy, cataplexy, and normal hypocretin values; 11 patients with idiopathic hypersomnia; and 43 controls. MEASUREMENTS AND RESULTS Body mass index (BMI), serum leptin, and HbA1C were measured in patients and controls; and CSF hypocretin 1 and leptin measured in all patients. Female and male patients with narcolepsy and hypocretin deficiency had the highest mean BMI (27.8 and 26.2, respectively), not statistically different from patients with narcolepsy and normal hypocretin or controls, but statistically higher than the patients with idiopathic hypersomnia (p < 0.001 and 0.011, respectively). The number of obese patients (BMI > 30) was increased in both narcolepsy groups. Serum and CSF leptin levels correlated positively to BMI in patients and controls, but not to CSF hypocretin concentrations. HbA1C was within normal levels and similar in all groups. CONCLUSIONS The study confirms a moderate tendency to obesity (BMI > 30) and overweight in patients with narcolepsy and cataplexy. Obesity was not correlated to hypocretin deficiency or reduced serum or CSF leptin concentrations. We suggest that overweight and possible metabolic changes previously reported in narcolepsy, may be caused by other mechanisms.
Collapse
|
25
|
Estrada K, Styrkarsdottir U, Evangelou E, Hsu YH, Duncan EL, Ntzani EE, Oei L, Albagha OME, Amin N, Kemp JP, Koller DL, Li G, Liu CT, Minster RL, Moayyeri A, Vandenput L, Willner D, Xiao SM, Yerges-Armstrong LM, Zheng HF, Alonso N, Eriksson J, Kammerer CM, Kaptoge SK, Leo PJ, Thorleifsson G, Wilson SG, Wilson JF, Aalto V, Alen M, Aragaki AK, Aspelund T, Center JR, Dailiana Z, Duggan DJ, Garcia M, Garcia-Giralt N, Giroux S, Hallmans G, Hocking LJ, Husted LB, Jameson KA, Khusainova R, Kim GS, Kooperberg C, Koromila T, Kruk M, Laaksonen M, Lacroix AZ, Lee SH, Leung PC, Lewis JR, Masi L, Mencej-Bedrac S, Nguyen TV, Nogues X, Patel MS, Prezelj J, Rose LM, Scollen S, Siggeirsdottir K, Smith AV, Svensson O, Trompet S, Trummer O, van Schoor NM, Woo J, Zhu K, Balcells S, Brandi ML, Buckley BM, Cheng S, Christiansen C, Cooper C, Dedoussis G, Ford I, Frost M, Goltzman D, González-Macías J, Kähönen M, Karlsson M, Khusnutdinova E, Koh JM, Kollia P, Langdahl BL, Leslie WD, Lips P, Ljunggren Ö, Lorenc RS, Marc J, Mellström D, Obermayer-Pietsch B, Olmos JM, Pettersson-Kymmer U, Reid DM, Riancho JA, Ridker PM, Rousseau F, Slagboom PE, Tang NLS, Urreizti R, Van Hul W, Viikari J, Zarrabeitia MT, Aulchenko YS, Castano-Betancourt M, Grundberg E, Herrera L, Ingvarsson T, Johannsdottir H, Kwan T, Li R, Luben R, Medina-Gómez C, Palsson ST, Reppe S, Rotter JI, Sigurdsson G, van Meurs JBJ, Verlaan D, Williams FMK, Wood AR, Zhou Y, Gautvik KM, Pastinen T, Raychaudhuri S, Cauley JA, Chasman DI, Clark GR, Cummings SR, Danoy P, Dennison EM, Eastell R, Eisman JA, Gudnason V, Hofman A, Jackson RD, Jones G, Jukema JW, Khaw KT, Lehtimäki T, Liu Y, Lorentzon M, McCloskey E, Mitchell BD, Nandakumar K, Nicholson GC, Oostra BA, Peacock M, Pols HAP, Prince RL, Raitakari O, Reid IR, Robbins J, Sambrook PN, Sham PC, Shuldiner AR, Tylavsky FA, van Duijn CM, Wareham NJ, Cupples LA, Econs MJ, Evans DM, Harris TB, Kung AWC, Psaty BM, Reeve J, Spector TD, Streeten EA, Zillikens MC, Thorsteinsdottir U, Ohlsson C, Karasik D, Richards JB, Brown MA, Stefansson K, Uitterlinden AG, Ralston SH, Ioannidis JPA, Kiel DP, Rivadeneira F. Genome-wide meta-analysis identifies 56 bone mineral density loci and reveals 14 loci associated with risk of fracture. Nat Genet 2012; 44:491-501. [PMID: 22504420 PMCID: PMC3338864 DOI: 10.1038/ng.2249] [Citation(s) in RCA: 887] [Impact Index Per Article: 73.9] [Reference Citation Analysis] [Abstract] [MESH Headings] [Grants] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/07/2011] [Accepted: 03/16/2012] [Indexed: 12/15/2022]
Abstract
Bone mineral density (BMD) is the most widely used predictor of fracture risk. We performed the largest meta-analysis to date on lumbar spine and femoral neck BMD, including 17 genome-wide association studies and 32,961 individuals of European and east Asian ancestry. We tested the top BMD-associated markers for replication in 50,933 independent subjects and for association with risk of low-trauma fracture in 31,016 individuals with a history of fracture (cases) and 102,444 controls. We identified 56 loci (32 new) associated with BMD at genome-wide significance (P < 5 × 10(-8)). Several of these factors cluster within the RANK-RANKL-OPG, mesenchymal stem cell differentiation, endochondral ossification and Wnt signaling pathways. However, we also discovered loci that were localized to genes not known to have a role in bone biology. Fourteen BMD-associated loci were also associated with fracture risk (P < 5 × 10(-4), Bonferroni corrected), of which six reached P < 5 × 10(-8), including at 18p11.21 (FAM210A), 7q21.3 (SLC25A13), 11q13.2 (LRP5), 4q22.1 (MEPE), 2p16.2 (SPTBN1) and 10q21.1 (DKK1). These findings shed light on the genetic architecture and pathophysiological mechanisms underlying BMD variation and fracture susceptibility.
Collapse
|