Identification of Transcripts with Shared Roles in the Pathogenesis of Postmenopausal Osteoporosis and Cardiovascular Disease.

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.

Mechanical-Stress-Related Epigenetic Regulation of ZIC1 Transcription Factor in the Etiology of Postmenopausal Osteoporosis.

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.

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.

Identifying Novel Gene Variants in Coronary Artery Disease and Shared Genes With Several Cardiovascular Risk Factors.

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.

Genetic Polymorphism of miR-196a-2 is Associated with Bone Mineral Density (BMD).

MicroRNAs (miRNAs) are small non-coding RNA molecules that post-transcriptionally regulate the translation of messenger RNAs. Given the crucial role of miRNAs in gene expression, genetic variants within miRNA-related sequences may affect miRNA function and contribute to disease risk. Osteoporosis is characterized by reduced bone mass, and bone mineral density (BMD) is a major diagnostic proxy to assess osteoporosis risk. Here, we aimed to identify miRNAs that are involved in BMD using data from recent genome-wide association studies (GWAS) on femoral neck, lumbar spine and forearm BMD. Of 242 miRNA-variants available in the GWAS data, we found rs11614913:C > T in the precursor miR-196a-2 to be significantly associated with femoral neck-BMD (p-value = 9.9 × 10-7, ß = -0.038) and lumbar spine-BMD (p-value = 3.2 × 10-11, ß = -0.061). Furthermore, our sensitivity analyses using the Rotterdam study data showed a sex-specific association of rs11614913 with BMD only in women. Subsequently, we highlighted a number of miR-196a-2 target genes, expressed in bone and associated with BMD, that may mediate the miRNA function in BMD. Collectively, our results suggest that miR-196a-2 may contribute to variations in BMD level. Further biological investigations will give more insights into the mechanisms by which miR-196a-2 control expression of BMD-related genes.

Increased ID2 levels in adult precursor B cells as compared with children is associated with impaired Ig locus contraction and decreased bone marrow output.

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.

The Influence of DNA Methylation on Bone Cells.

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.

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.

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 (r(2) = 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.

Global miRNA expression analysis of serous and clear cell ovarian carcinomas identifies differentially expressed miRNAs including miR-200c-3p as a prognostic marker.

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.

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.

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.

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

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

Social stratification reflected in bone mineral density and stature: Spectral imaging and osteoarchaeological findings from medieval Norway.

This study presents skeletal material from five medieval burial sites in Eastern Norway, confined to one royal burial church, one Dominican monastery, and three burial sites representing parish populations. We combine osteological analysis and Dual Energy X-Ray Absorptiometry, studying the remains of 227 individuals (102 females and 125 males) employing young, middle, and old adult age categories. The aim is to assess bone mineral density as a skeletal indicator of socioeconomic status including stature as a variable. We detected that socioeconomic status significantly affected bone mineral density and stature. Individuals of high status had higher bone mineral density (0.07 g/cm2, p = 0.003) and taller stature (1.85 cm, p = 0.017) than individuals from the parish population. We detected no significant relationship between young adult bone mineral density and socioeconomic status (p = 0.127 and 0.059 for females and males, respectively). For males, high young adult bone mineral density and stature varied concordantly in both status groups. In contrast, females of high status were significantly taller than females in the parish population (p = 0.011). Our findings indicate quite different conditions during growth and puberty for the two groups of females. The age-related pattern of bone variation also portrayed quite different trajectories for the two socioeconomic status groups of both sexes. We discuss sociocultural practices (living conditions during childhood and puberty, as well as nutritional and lifestyle factors in adult life), possibly explaining the differences in bone mineral density between the high-status and parish population groups. The observation of greater differences in bone mineral density and stature for females than males in the medieval society of Norway is also further discussed.

Heavy-load exercise in older adults activates vasculogenesis and has a stronger impact on muscle gene expression than in young adults.

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.

In silico discovery of blood cell macromolecular associations.

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.

Zic1 transcription factor in bone: neural developmental protein regulates mechanotransduction in osteocytes.

A transcriptome analysis compared gene expression in human bone biopsy samples taken from lumbar spine and iliac crest, sites that experience high and low levels of mechanical stress, respectively. The analysis revealed that the zinc finger protein of cerebellum (Zic) family member transcription factor Zic1 was the most up-regulated gene in the lumbar spine (202-fold; P<10(-7)) in comparison with the iliac crest. Software analysis of differential gene expression in the biopsy samples identified the ciliary-related proteins PATCH1 and GLI-Kruppel family members Gli1 and Gli3 as part of a potential molecular network associated with Zic1. RT-PCR confirmed the expression of Zic1, Gli1, and Gli3 and other related key signaling mediators in osteoblastic cells and osteocytes in vitro. Zic1 was immunolocalized in the cytosol and nucleus of the murine osteocyte cell line MLO-Y4 and osteoblast-like cells MC3T3-E1 and in primary rat osteoblasts. MLO-Y4 cells subjected to prolonged oscillatory fluid flow showed increased localization of Zic1 in the nucleus with diminished levels in the cytosol, but no such changes were seen in MC3T3-E1 cells. A shear stress-induced increase in T-cell factor/lymphoid enhancer factor transcriptional activity was abolished by Zic1 gene silencing. These results suggest that Zic1, perhaps together with Gli1 and Gli3, may act as a link between mechanosensing and Wnt signaling. We conclude that Zic1, a neural developmental transcription factor, plays an important role in shear flow mechanotransduction in osteocytes.

Lipocalin 2 serum levels correlate with age and bone turnover biomarkers in healthy subjects but not in postmenopausal osteoporotic women.

PURPOSE: Lipocalin 2 (LCN2) is an adipokine involved in many physiological functions, including bone metabolism. We previously demonstrated its implication in mouse models of mechanical unloading-induced osteoporosis and in a cohort of bed rest volunteers. We therefore aimed at studying its involvement in postmenopausal osteoporosis. METHODS: We measured serum LCN2 and correlated its levels to Dickkopf WNT Signaling Pathway Inhibitor 1 (DKK1), Tartrate Resistant Acid Phosphatase 5B (TRAcP5B), sclerostin, urinary N-terminal telopeptide of type I collagen (NTX), serum C-terminal telopeptide of type I collagen (CTX), parathyroid hormone and vitamin K by ELISA performed in a cohort of younger (50-65 years) and older (66-90 years) osteoporotic women in comparison to healthy subjects. A cohort of male healthy and osteoarthritic patients was also included. Sobel mediation analysis was used to test indirect associations among age, LCN2 and DKK1 or NTX. RESULTS: LCN2 levels were unchanged in osteoporotic and in osteoarthritis patients when compared to healthy subjects and did not correlate with BMD. However, serum LCN2 correlated with age in healthy women (R = 0.44; P = 0.003) and men (R = 0.5; P = 0.001) and serum concentrations of DKK1 (R = 0.47; P = 0.003) and urinary NTX (R = 0.34; P = 0.04). Sobel mediation analysis showed that LCN2 mediates an indirect relationship between age and DKK1 (P = 0.02), but not with NTX, in healthy subjects. CONCLUSIONS: Taken together, the results suggest a hitherto unknown association between LCN2, DKK1 and age in healthy individuals, but not in postmenopausal osteoporotic women.

Striking decrease in the total precursor B-cell compartment during early childhood as evidenced by flow cytometry and gene expression changes.

The number of circulating B-cells in peripheral blood plateaus between 2 and 24 months of age, and thereafter declines gradually. How this reflects the kinetics of the precursor B-cell pool in the bone marrow is of clinical interest, but has not been studied thoroughly in humans. The authors analyzed bone marrow (n = 37) from healthy children and adults (flow cytometry) searching for age-related changes in the total precursor B-cell compartment. In an age-matched cohort (n = 25) they examined age-related global gene expression changes (Affymetrix) in unsorted bone marrow with special reference to the recombination activating gene 1, RAG1. Subsequently, they searched the entire gene set for transcripts correlating to the RAG1 profile to discover other known and possibly new precursor B-cell related transcripts. Both methods disclosed a marked, transient increase of total precursor B-cells at 6-20 months, followed by a rapid decrease confined to the first 2 years. The decline thereafter was considerably slower, but continued until adulthood. The relative composition of total precursor B-cells, however, did not change significantly with age. The authors identified 54 genes that were highly correlated to the RAG1 profile (r >or= .9, p < 1 x 10(-8)). Of these 54 genes, 15 were characteristically B-lineage associated like CD19, CD79, VPREB, EBF1, and PAX5; the remaining 39 previously not described as distinctively B-lineage related. The marked, transient increase in precursor B-cells and RAG1 transcriptional activity is not reflected by a similar peak in B-cells in peripheral blood, whereas the sustained plateau concurs in time.

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

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

Mendelian Randomization Analysis Reveals a Causal Influence of Circulating Sclerostin Levels on Bone Mineral Density and Fractures.

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.

Calmodulin-dependent kinase 1beta is expressed in the epiphyseal growth plate and regulates proliferation of mouse calvarial osteoblasts in vitro.

The Ca(2+)/Calmodulin-dependent protein kinase (CaMK) family is activated in response to elevation of intracellular Ca(2+), and includes CaMK1 (as well as CaMK2 and CaMK4), which exists as different isoforms (alpha, beta, gamma and delta). CaMK1 is present in several cell types and may be involved in various cellular processes, but its role in bone is unknown. In situ hybridization was used to determine the spatial and temporal expression of CaMK1beta during endochondral bone development in mouse embryos and newborn pups. The cellular and subcellular distribution of CaMK1 was assessed by quantitative immunogold electron microscopy (EM). The role of CaMK1beta in mouse calvarial osteoblasts was investigated by using small interfering RNA (siRNA) to silence its expression, while in parallel monitoring cell proliferation and levels of skeletogenic transcripts. cRNA in situ hybridization and EM studies show that CaMK1beta is mainly located in developing long bones and vertebrae (from ED14.5 until day 10 after birth), with highest expression in epiphyseal growth plate hypertrophic chondrocytes. By RT-PCR, we show that CaMK1beta2 (but not beta1) is expressed in mouse hind limbs (in vivo) and mouse calvarial osteoblasts (in vitro), and also in primary human articular chondrocyte cultures. Silencing of CaMK1beta in mouse calvarial osteoblasts by siRNA significantly decreases osteoblast proliferation and c-Fos gene expression (approx. 50%), without affecting skeletogenic markers for more differentiated osteoblasts (i.e. Cbfa1/Runx2, Osterix (Osx), Osteocalcin (Oc), Alkaline phosphatase (Alp) and Osteopontin (Opn)). These results identify CaMK1beta as a novel regulator of osteoblast proliferation, via mechanisms that may at least in part involve c-Fos, thus implicating CaMK1beta in the regulation of bone and cartilage development.