Search

Genome-wide association studies (GWAS) have repeatedly identified genetic variants associated with bone mineral density (BMD) and osteoporotic fracture in non-coding regions of C7ORF76, a poorly studied gene of unknown function. The aim of the present study was to elucidate the causality and molecular mechanisms underlying the association. We re-sequenced the genomic region in two extreme BMD groups from the BARCOS cohort of postmenopausal women to search for functionally relevant variants. Eight selected variants were tested for association in the complete cohort and 2 of them (rs4342521 and rs10085588) were found significantly associated with lumbar spine BMD and nominally associated with osteoporotic fracture. cis-eQTL analyses of these 2 SNPs, together with SNP rs4727338 (GWAS lead SNP in Estrada et al., Nat Genet. 44:491-501, 2012), performed in human primary osteoblasts, disclosed a statistically significant influence on the expression of the proximal neighbouring gene SLC25A13 and a tendency on the distal SHFM1. We then studied the functionality of a putative upstream regulatory element (UPE), containing rs10085588. Luciferase reporter assays showed transactivation capability with a strong allele-dependent effect. Finally, 4C-seq experiments in osteoblastic cell lines showed that the UPE interacted with different tissue-specific enhancers and a lncRNA (LOC100506136) in the region. In summary, this study is the first one to analyse in depth the functionality of C7ORF76 genomic region. We provide functional regulatory evidence for the rs10085588, which may be a causal SNP within the 7q21.3 GWAS signal for osteoporosis., (Copyright © 2019 Elsevier Inc. All rights reserved.)

MicroRNAs (miRNAs) are important regulators of many cellular processes, including the differentiation and activity of osteoblasts, and therefore, of bone turnover. MiR-320a is overexpressed in osteoporotic bone tissue but its role in osteoblast function is unknown. In the present study, functional assays were performed with the aim to elucidate the mechanism of miR-320a action in osteoblastic cells. MiR-320a was either overexpressed or inhibited in human primary osteoblasts (hOB) and gene expression changes were evaluated through microarray analysis. In addition, the effect of miR-320a on cell proliferation, viability, and oxidative stress in hOB was evaluated. Finally, matrix mineralization and alkaline phosphatase activity were assessed in order to evaluate osteoblast functionality. Microarray results showed miR-320a regulation of a number of key osteoblast genes and of genes involved in oxidative stress. Regulation of osteoblast differentiation and ossification appeared as the best significant biological processes (PANTHER P value = 3.74E-05; and P value = 3.06E-04, respectively). The other enriched pathway was that of the cellular response to cadmium and zinc ions, mostly by the overexpression of metallothioneins. In hOBs, overexpression of miR-320a increased cell proliferation and oxidative stress levels whereas mineralization capacity was reduced. In conclusion, overexpression of miR-320a increased stress oxidation levels and was associated with reduced osteoblast differentiation and functionality, which could trigger an osteoporotic phenotype., Competing Interests: The authors have declared that no competing interests exist.

Numerous GWAS and candidate gene studies have highlighted the role of the Wnt pathway in bone biology. Our objective has been to study in detail the allelic architecture of three Wnt pathway genes: WNT16, DKK1 and SOST, in the context of osteoporosis. We have resequenced the coding and some regulatory regions of these three genes in two groups with extreme bone mineral density (BMD) (n = ∼50, each) from the BARCOS cohort. No interesting novel variants were identified. Thirteen predicted functional variants have been genotyped in the full cohort (n = 1490), and for ten of them (with MAF > 0.01), the association with BMD has been studied. We have found six variants nominally associated with BMD, of which 2 WNT16 variants predicted to be eQTLs for FAM3C (rs55710688, in the Kozak sequence and rs142005327, within a putative enhancer) withstood multiple-testing correction. In addition, two rare variants in functional regions (rs190011371 in WNT16b 3'UTR and rs570754792 in the SOST TATA box) were found only present in three women each, all with BMD below the mean of the cohort. Our results reinforce the higher importance of regulatory versus coding variants in these Wnt pathway genes and open new ways for functional studies of the relevant variants.

Background: MicroRNAs (miRNAs) are important regulators of gene expression, with documented roles in bone metabolism and osteoporosis, suggesting potential therapeutic targets. Our aim was to identify miRNAs differentially expressed in fractured vs nonfractured bones. Additionally, we performed a miRNA profiling of primary osteoblasts to assess the origin of these differentially expressed miRNAs., Methods: Total RNA was extracted from (a) fresh femoral neck trabecular bone from women undergoing hip replacement due to either osteoporotic fracture (OP group, n = 6) or osteoarthritis in the absence of osteoporosis (Control group, n = 6), matching the two groups by age and body mass index, and (b) primary osteoblasts obtained from knee replacement due to osteoarthritis (n = 4). Samples were hybridized to a microRNA array containing more than 1900 miRNAs. Principal component analysis (PCA) plots and heat map hierarchical clustering were performed. For comparison of expression levels, the threshold was set at log fold change > 1.5 and a p-value < 0.05 (corrected for multiple testing)., Results: Both PCA and heat map analyses showed that the samples clustered according to the presence or absence of fracture. Overall, 790 and 315 different miRNAs were detected in fresh bone samples and in primary osteoblasts, respectively, 293 of which were common to both groups. A subset of 82 miRNAs was differentially expressed (p < 0.05) between osteoporotic and control osteoarthritic samples. The eight miRNAs with the lowest p-values (and for which a validated miRNA qPCR assay was available) were assayed, and two were confirmed: miR-320a and miR-483-5p. Both were over-expressed in the osteoporotic samples and expressed in primary osteoblasts. miR-320a is known to target CTNNB1 and predicted to regulate RUNX2 and LEPR, while miR-483-5p down-regulates IGF2. We observed a reduction trend for this target gene in the osteoporotic bone., Conclusions: We identified two osteoblast miRNAs over-expressed in osteoporotic fractures, which opens novel prospects for research and therapy.

Osteoporosis is a complex disease involving many putative genetic factors. Association analysis of functional SNPs in candidate genes is an important tool for their identification. However, this approach is affected by limited power, population stratification, and other drawbacks that lead to discordant results. Replication in independent cohorts is essential. We performed association analyses of three functional polymorphisms previously associated with bone phenotypes--namely, Ala222Val in MTHFR, Ile1062Val in LRP6, and -13910C>T in LCT--in a cohort of 944 postmenopausal Spanish women, all of them with lumbar spine (LS) bone mineral density (BMD) data and most with femoral neck (FN) BMD and fracture data. We found significant differences between genotypes only for the MTHFR polymorphism and vertebral factures, with an OR of 2.27 (95% CI 1.17-4.38) for the TT vs. CC/CT genotypes, P = 0.018. We present genotype and allele frequency data for LCT -13910C>T for a Spanish population, where the T allele (conferring lactase persistence) has a frequency of 38.6%. Genotype frequencies were consistent with observed clines in Europe and with the prevalence of lactase nonpersistence. The LCT -13910C>T polymorphism was significantly associated with height and weight, such that T allele carriers were 0.88 cm taller (95% CI 0.08-1.59 cm, P = 0.032, adjusted by age) than CC individuals and TT homozygotes were 1.91 kg heavier than CC/CT individuals (95% CI 0.11-3.71 kg, P = 0.038, adjusted by age). In conclusion, no significant association was observed between the studied polymorphisms and LS BMD or FN BMD in postmenopausal Spanish women, and only MTHFR Ala222Val was associated with vertebral fractures.

The RANKL/RANK/OPG pathway is essential for bone remodeling regulation. Many hormones and cytokines are involved in regulating gene expression in most of the pathway components. Moreover, any deregulation of this pathway can alter bone metabolism, resulting in loss or gain of bone mass. Whether osteoblasts from osteoporotic and nonosteoporotic patients respond differently to cytokines is unknown. The aim of this study was to compare the effect of interleukin (IL)-1beta, proftaglandin E(2) (PGE(2)), and transforming growth factor-beta1 (TGF-beta1) treatments on OPG and RANKL gene expression in normal (n = 11) and osteoporotic (n = 8) primary osteoblasts. OPG and RANKL mRNA levels of primary human osteoblastic (hOB) cell cultures were assessed by real-time PCR. In all cultures, OPG mRNA increased significantly in response to IL-1beta treatment and decreased in response to TGF-beta1 whereas PGE(2) treatment had no effect. RANKL mRNA levels were significantly increased by all treatments. Differences in OPG and RANKL responses were observed between osteoporotic and nonosteoporotic hOB: in osteoporotic hOB, the OPG response to IL-1beta treatment was up to three times lower (P = 0.009), whereas that of RANKL response to TGF-beta1 was five times higher (P = 0.002) after 8 h of treatment, as compared with those in nonosteoporotic hOBs. In conclusion, osteoporotic hOB cells showed an anomalous response under cytokine stimulation, consistent with an enhanced osteoclastogenesis resulting in high levels of bone resorption., ((c) 2010 Wiley-Liss, Inc.)

LRP5 encodes the low-density lipoprotein receptor-related protein 5, a transmembrane protein involved in Wnt signaling. LRP5 is an important regulator of osteoblast growth and differentiation, affecting bone mass in vertebrates. Whether common variations in LRP5 are associated with normal BMD variation or osteoporotic phenotypes is of great relevance. We used a haplotype-based approach to search for common disease-associated variants in LRP5 in a cohort of 964 Spanish postmenopausal women. Twenty-four SNPs were selected, covering the LRP5 region, including the missense changes p.V667M and p.A1330V. The SNPs were genotyped and evaluated for association with BMD at the lumbar spine (LS) or femoral neck (FN) and with osteoporotic fracture, at single SNP and haplotype levels, by regression methods. Association with LS BMD was found for SNP 1, rs312009, located in the 5'-flanking region (p = 0.011, recessive model). SNP 6, rs2508836, in intron 1, was also associated with BMD, both at LS (p = 0.025, additive model) and FN (p = 0.031, recessive model). Two polymorphisms were associated with fracture: SNP 11, rs729635, in intron 1, and SNP 15, rs643892, in intron 5 (p = 0.007 additive model and p = 0.019 recessive model, respectively). Haplotype analyses did not provide additional information, except for haplotype "GC" of the block located at the 3'end of the gene. This haplotype spans intron 22 and the 3' untranslated region and was associated with FN BMD (p = 0.029, one copy of the haplotype versus none). In silico analyses showed that SNP 1 (rs312009) lies in a putative RUNX2 binding site. Electro-mobility shift assays confirmed RUNX2 binding to this site.

Context: The COLIA1 gene is a strong candidate for susceptibility to osteoporosis. The causal genetic variants are currently unclear, but the most likely are functional polymorphisms in the promoter and intron 1 of COLIA1., Objective: The objective of the study was to determine whether promoter and intron 1 polymorphisms of COLIA1 or haplotypes defined by these polymorphisms regulate bone mineral density (BMD) in women., Design: This was a population-based association study involving 3270 women from the United Kingdom who took part in a regional osteoporosis screening program., Main Outcome Measures: BMD at the lumbar spine (LS-BMD) and femoral neck (FN-BMD) was measured on two occasions approximately 6 yr apart, in relation to polymorphisms and haplotypes defined by polymorphisms within the COLIA1 intron 1 (+1245G/T; rs1800012) and promoter (-1997G/T; rs1107946; -1663IndelT; rs2412298)., Results: The polymorphisms were in strong linkage disequilibrium, and three haplotypes accounted for more than 95% of alleles at the COLIA1 locus. The individual polymorphisms were associated with BMD, but the most consistent associations were with haplotypes defined by all three polymorphisms. Homozygote carriers of haplotype 2 (-1997G/-1663delT/+1245T) had reduced BMD at baseline (P = 0.007 for LS-BMD; P = 0.008 for FN-BMD), whereas homozygotes for haplotype 3 (-1997T/-1663insT/+1245G) had increased BMD (P = 0.007 for LS-BMD). Similar associations were observed at follow-up for haplotype 3, but the association with haplotype 2 was weaker due to increased uptake of hormone replacement therapy in homozygotes for this haplotype., Conclusions: Two haplotypes defined by polymorphisms in the 5' flank of the COLIA1 regulate BMD in a bidirectional manner in women.

Background: Polymorphisms of the vitamin D receptor (VDR) gene have been implicated in the genetic regulation of bone mineral density (BMD). However, the clinical impact of these variants remains unclear., Objective: To evaluate the relation between VDR polymorphisms, BMD, and fractures., Design: Prospective multicenter large-scale association study., Setting: The Genetic Markers for Osteoporosis consortium, involving 9 European research teams., Participants: 26,242 participants (18,405 women)., Measurements: Cdx2 promoter, FokI, BsmI, ApaI, and TaqI polymorphisms; BMD at the femoral neck and the lumbar spine by dual x-ray absorptiometry; and fractures., Results: Comparisons of BMD at the lumbar spine and femoral neck showed nonsignificant differences less than 0.011 g/cm2 for any genotype with or without adjustments. A total of 6067 participants reported a history of fracture, and 2088 had vertebral fractures. For all VDR alleles, odds ratios for fractures were very close to 1.00 (range, 0.98 to 1.02) and collectively the 95% CIs ranged from 0.94 (lowest) to 1.07 (highest). For vertebral fractures, we observed a 9% (95% CI, 0% to 18%; P = 0.039) risk reduction for the Cdx2 A-allele (13% risk reduction in a dominant model)., Limitations: The authors analyzed only selected VDR polymorphisms. Heterogeneity was detected in some analyses and may reflect some differences in collection of fracture data across cohorts. Not all fractures were related to osteoporosis., Conclusions: The FokI, BsmI, ApaI, and TaqI VDR polymorphisms are not associated with BMD or with fractures, but the Cdx2 polymorphism may be associated with risk for vertebral fractures.

Context: Both bone mineral density (BMD) and fracture risk have a strong genetic component. Estrogen receptor alpha (ESR1) is a candidate gene for osteoporosis, but previous studies of ESR1 polymorphisms in this field were hampered by small sample size, lack of standardization, and inconclusive results., Objective: To generate large-scale evidence on whether 3 common ESR1 polymorphisms (intron 1 polymorphisms XbaI [dbSNP: rs9340799] and PvuII [dbSNP: rs2234693] and promoter TA repeats microsatellite) and haplotypes thereof are associated with BMD and fractures., Design and Setting: Meta-analysis of individual-level data involving standardized genotyping of 18 917 individuals in 8 European centers., Main Outcome Measures: BMD of femoral neck and lumbar spine; all fractures and vertebral fractures by genotype., Results: No between-center heterogeneity was observed for any outcome in any genetic contrast. None of the 3 polymorphisms or haplotypes had any statistically significant effect on BMD in adjusted or unadjusted analyses, and estimated differences between genetic contrasts were 0.01 g/cm2 or less. Conversely, we found significant reductions in fracture risk. In women homozygous for the absence of an XbaI recognition site, the adjusted odds of all fractures were reduced by 19% (odds ratio, 0.81 [95% CI, 0.71-0.93]; P = .002) and vertebral fractures by 35% (odds ratio, 0.65 [95% CI, 0.49-0.87]; P = .003). Effects on fractures were independent of BMD and unaltered in adjusted analyses. No significant effects on fracture risk were seen for PvuII and TA repeats., Conclusions: ESR1 is a susceptibility gene for fractures, and XbaI determines fracture risk by mechanisms independent of BMD. Our study demonstrates the value of adequately powered studies with standardized genotyping and clinical outcomes in defining effects of common genetic variants on complex diseases.

Introducción: WNT16 es un gen importante en la homeostasis ósea, que se encuentra en un locus muy complejo que también incluye los genes vecinos: ING3, FAM3C y CPED1. Además del claro papel de WNT16 sobre la determinación de la densidad mineral ósea (DMO), también se han encontrado pruebas sobre la importancia de estos tres genes vecinos en el metabolismo óseo. Queda por tanto esclarecer si las variantes en WNT16 asociadas a la DMO realizan su efecto propiamente sobre WNT16 o si lo hacen modificando la expresión de estos genes vecinos. Material y métodos: Hemos determinado los niveles de expresión de CPED1 y FAM3C en osteoblastos primarios y hemos comprobado si variantes de WNT16 se comportan como loci de rasgos cuantitativos de expresión (expresion quantitative trait loci; eQTL) de estos genes. Resultados: La variante de cambio de aminoácido rs2908004 en WNT16 actúa como eQTL de FAM3C en osteoblastos primarios bajo la hipótesis de modelo dominante. Discusión: Es posible que el efecto de esta variante sobre la DMO sea debido a la modificación de los niveles de expresión de FAM3C además o en vez de un efecto directo de la proteína WNT16 mutante resultante del cambio de aminoácido.

Objetivos: Identificar microRNAs (miRNAs) diferencialmente expresados en muestras óseas con fractura osteoporótica respecto a huesos sanos.Material y métodos: Se extrajo RNA total a partir de hueso trabecular fresco del cuello femoral de mujeres sometidas a reemplazo de cadera, ya sea debido a fractura osteoporótica (n=6) o por artrosis en ausencia de osteoporosis (según la DMO) (n=6). Las muestras se hibridaron en un array de miRNAs y se realizaron diagramas de PCA y de mapa de calor. Para la comparación de los niveles de expresión, se fijó como significativo un umbral de cambio de >1,5 veces y un valor p

Objetivos: El sistema RANK/RANKL/OPG está implicado en la determinación de la densidad mineral ósea (DMO) y de la microarquitectura ósea. Nuestro estudio pretende evaluar si existen SNPs en la región 3’UTR del gen RANK asociados con fenotipos osteoporóticos. Material y métodos: Se genotiparon siete variantes genéticas en 1.098 mujeres de la cohorte BARCOS, y se evaluó su asociación con la DMO y las fracturas osteoporóticas. Se testó una interacción con el SNP rs9594738 en el gen RANKL el cual fue previamente asociado con la DMO. Resultados: Ninguno de los SNPs se asoció significativamente con la DMO. El SNP rs78326403 se asoció con fracturas de muñeca/antebrazo (Modelo Log-aditivo odds ratio (OR)=3,12 [IC 95%: 1,69 ; 5,75]; p=7,16×10-4), mientras que el SNP rs884205 se asoció con fracturas de columna vertebral (OR=4,05 Recesivo; [IC 95%: 1,59 ; 10,35]; p=8,24×10-3). Por último, se detectó una interacción entre el SNP rs9594738 del RANKL y el rs78326403 del RANK sobre la prevalencia de fractura (p=0,039). El análisis del efecto de los genotipos compuestos rs9594738 y rs78326403 dio un aumento de la prevalencia de fracturas en sujetos con un mayor número de alelos desfavorables, siendo las OR 2,76 [IC 95%: 1,30 ; 5,81]; p=0,007) y 5,14 [IC 95%: 1,37 ; 15,67]; p=0,007) para 2 y ≥3 alelos desfavorables, respectivamente, en comparación con ninguno/1. Conclusiones: Dos SNPs en el 3’UTR del gen RANK predisponen a la fractura osteoporótica sitio-dependiente. Una interacción con el SNP rs9594738 del RANKL sugiere un efecto aditivo de la DMO y la resistencia ósea.

Objectives: To identify microRNAs (miRNAs) differentially expressed in bone samples with osteoporotic fracture compared with healthy bones. Methods: Total RNA was extracted from fresh trabecular bone of the femoral neck of women undergoing hip replacement surgery, either because to osteoporotic fracture (n=6) or in the absence of osteoarthritis osteoporosis (based on BMD) (n=6). The samples were hybridized on an array of miRNAs and PCA diagrams and heat map were made. To compare expression levels, >1.5 times and a value p<0.05 Student's T test (corrected for multiple testing) was set as a threshold of significant change. Results: Both PCA analysis and the heat map showed a samples grouping whether there was fracture or not. 790 were detected miRNAs in bone samples, 82 of which were altered in the osteoporotic samples. After validation in another panel of 6 samples 6 osteoporotic and non-osteoporotic by PCR real time of the most significant miRNAs, and for which there was a test available, the miRNAs, miR-320a and miR-22-3p were confirmed. These two miRNAs were detected in cultures of primary osteoblasts, although they did not maintain the same pattern of expression in total bone samples. Conclusions: We have shown that there are differences in the expression of miRNAs in samples with osteoporotic fracture. This opens prospects for research and design of new therapies. [ABSTRACT FROM AUTHOR]

Objectives: The RANK/RANKL/OPG system is involved in the determination of bone mineral density (BMD) and bone microarchitecture. Our study seeks to evaluate if there are SNPs in the 3'UTR region of the RANK gene associated with osteoporotic phenotypes. Material and methods: Seven genetic variants in 1,098 women from the BARCOS cohort were genotyped, and their association with BMD and osteoporotic fractures evaluated. An interaction with SNP rs9594738 in the RANKL gene which was previously associated with BMD was tested. Results: None of the SNPs were associated significantly with BMD. SNP rs78326403 was associated with wrist/forearm fractures (Log-additive model odds ratio (OR)=3.12 [IC 95%: 1.69 ; 5.75]; p=7.16x10-4), while SNP rs884205 was associated with fractures of the spinal column (OR=4.05 Recessive; [95% CI: 1.59 ; 10.35]; p=8.24x10-3). Lastly, an interaction was detected between SNP rs9594738 from RANKL and rs78326403 from RANK on the presence of fracture (p=0.039). The analysis of the effects of combined genotypes rs9594738 and rs78326403 pointed to an increase in the prevalence of fractures in subjects with a greater number of unfavourable alleles, the ORs being 2.76 [95% CI: 1.30 ; 5.81]; p=0.007) and 5.14 [95% CI: 1.37 ; 15.67]; p=0.007) for 2 and ≥3 unfavourable alleles respectively, in comparison with none/1. Conclusions: Two SNPs in 3'UTR from the RANK gene predispose to site-dependent osteoporotic fracture. An interaction with SNP rs9594738 from RANKL suggests an additive effect of BMD and bone strength. [ABSTRACT FROM AUTHOR]

Osteoprotegerin plays a key role in bone remodelling. We studied the association between 24 polymorphisms and haplotypes on the OPG gene and bone mineral density and fractures. After multiple-testing correction, one SNP and two block-haplotypes were significantly associated with FN BMD. Two other block-haplotypes were associated with fracture. Osteoprotegerin (OPG) plays a key role in bone remodelling. Here we studied the association between polymorphisms and haplotypes on the OPG gene and bone mineral density (BMD) and fractures. Twenty-four single nucleotide polymorphisms (SNPs) were selected to cover six haplotypic blocks and were genotyped in 964 postmenopausal Spanish women. Haplotypes were established with HaploStats. Association was analysed by GLM (for BMD) and logistic regression (for fractures) both at single SNP and haplotype levels. Upon adjustment for multiple testing ( p < 0.0073), one of the SNPs (SNP #17, rs1032129) remained significantly associated with FN BMD ( p = 0.001). Four block-haplotypes stood multiple-testing correction. Two remained associated with FN BMD and two with fracture. The association of block-4 haplotype “AC” (of SNPs #18 and #17) with FN BMD ( p = 0.0002) was stronger than that of SNP#17 alone and was the best result overall. A global assessment of the results indicated that all the alleles and haplotypes with a protective effect, at p < 0.05, belonged to a frequent long-range haplotype. In conclusion, these results provide a detailed picture of the involvement of common variants and haplotypes of the OPG gene in bone phenotypes. [ABSTRACT FROM AUTHOR]