Abstract

Background

Osteoporosis individual susceptibility is determined by the interaction of multiple genetic variants and environmental factors. The aim of this study was to conduct SNP-SNP interaction analyses in candidate genes influencing heel quantitative ultrasound (QUS) parameter in early adulthood to identify novel insights into the mechanism of disease.

Methods

The study population included 575 healthy subjects (mean age 20.41; SD 2.36). To assess bone mass QUS was performed to determine Broadband ultrasound attenuation (BUA, dB/MHz). A total of 32 SNPs mapping to loci that have been characterized as genetic markers for QUS and/or BMD parameters were selected as genetic markers in this study. The association of all possible SNP pairs with QUS was assessed by linear regression and a SNP-SNP interaction was defined as a significant departure from additive effects.

Conclusion

Our analysis of SNP-SNP interaction in candidate genes of QUS in Caucasian young adults reveal several interactions, especially between ESR1 and LRP5 genes, that did not reach statistical significance. Although our results do not support a relevant genetic contribution of SNP-SNP epistatic interactions to QUS in young adults, further studies in larger independent populations would be necessary to support these preliminary findings.

Keywords

Gene interactionQuantitative ultrasoundCandidate gene

Background

Osteoporosis represents an important public health problem worldwide that triggers near 9 million fractures annually [1]. Its prevalence is increasing due to the ageing of the population [2]. Osteoporosis is characterized by a deterioration of bone tissue microarchitecture that leads to low bone mineral density (BMD) and an increased risk of bone fragility fractures [3]. As occurs in many complex diseases, it is widely accepted that osteoporosis individual susceptibility is determined by the interaction of multiple genetic variants and environmental factors [4, 5].

An important determinant of osteoporosis risk later in life is the peak bone mass (PBM) [6]. Enhancing bone mass accrual to maximize PBM, which is attained by early adulthood, could help in the elderly to reduce the risk of fracture [6]. Current evidence suggests that genetic factors are major contributors to regulation of PBM, accounting for 50% and 80% of the variance in BMD [7, 8]. Therefore, identifying genetic factors affecting bone mass in early adulthood would be worthwhile.

Recently, it has been demonstrated that QUS predicts fracture independently of BMD [9, 10]. QUS provides information about bone mass and other quality aspects which may be of relevance in determining fracture risk such as bone microstructure, elasticity or connectivity [11, 12]. Besides this, previous studies have shown the ability of heel QUS to predict independently fracture risk [11, 13]. Compared to dual-energy x-ray absorptiometry (DXA), QUS is more easily accessible, low-cost, non-invasive and non-ionizing [14].

To identify the genetic factors underlying osteoporosis susceptibility, extensive genetic studies have been conducted reporting association of several genetic variants with different bone phenotypes [15–18]. However, these explain little of the heritability in complex phenotypes [19]. The missing heritability problem refers to the observation that the number of significant associations discovered does not form a substantial proportion of heritability for most traits. Thus, SNP-SNP interactions could help explain the missing heritability of common complex traits and therefore, their identification is considered of relevance. To date, only a few studies have been performed to identify SNP-SNP interactions influencing osteoporosis-related traits [20–26]. Most of them have investigated BMD determined by DXA, and therefore the possible influence of genetic interactions on QUS parameters is still unknown [25–27].

The independent association between single nucleotide polymorphisms (SNPs) in WNT16 (rs2908007, rs2908004, and rs2707466) [18, 28, 29], RSPO3 (rs774121) [18] and LRP5 (rs3736228) [30] genes and QUS parameters have been well established. In previous replication studies from our group we investigated the implication of 32 genetic markers reported to influence QUS and/or BMD parameters, and confirmed that WTN16, RSPO3 and LRP5 are genetic factors that determine bone mineralization in young adults [28, 30].

Taking into consideration all this evidence, and bearing in mind that it is likely that interactions of common variants at different loci are influencing QUS traits, the aim of this study was to identify SNP-SNP interactions between SNPs independently associated with QUS, that could contribute to QUS traits variation in early adulthood.

Methods

Study subjects

A total of 575 Caucasian healthy individuals (400 females and 175 males, median age 20,41 ± 2,69 recruited from different centres of Granada (Spain) composed the study population after giving written informed consent. Local ethics committees approved the study that was conducted following the Declaration of Helsinki. Were excluded from the study those individuals who reported therapy with hormonal contraceptives, diseases of metabolic or endocrine systems and history of bone disease.

Covariates

A body composition analyzer (TANITA BC-418MA) was used to estimate weight, fat mass and lean mass to the nearest 0,11 Kg. Height measurement was performed with a Harpenden stadiometer (Holtain 602VR®) to the nearest 0,1 cm. Body mass index (BMI) was calculated as weight divided by height squared (kg/m2). Physical activity (PA) was quantified using the short form and self-administered version of the International Physical Activity Questionnaire (IPAQ) [31]. Dietary calcium intake (DCI) was obtained from a 72 h recall method including intakes on Thursday, Friday, and Saturday. During the interviews and for a better precision of the food records, pictorial food models and standard household measures were used. All the 72 h recall were computerized with a nutrient analysis program (Nutriber 1.1.5) to convert food records into amounts of nutrients intake.

QUS of the heel

Ultrasound measurements (BUA, dB/MHz) were performed using the CUBA clinical ultrasound bone densitometer (McCue Ultrasonic Limited, Compton, Winchester, UK) at the right calcaneus to assess bone mass. This localization was selected on the basis of its accessibility and its high content of trabecular bone [32]. To guarantee the long-term stability of the ultrasound bone densitometer, daily calibrations were performed using a with physical phantom.

SNPs selection and genotyping process

A total of 32 SNPs mapping to loci that have been characterized as genetic markers for QUS and/or BMD parameters were selected.

Genotyping was carried out using a custom design TaqMan OpenArray (Life Technologies, Carlsbad, CA, USA). Genotyping plate containing a predesigned TaqMan genotyping assay for every polymorphism was conducted in the Genomic and Genotyping unit of GENYO centre (Pfizer-University of Granada-Junta de Andalucía Centre for Genomics and Oncological Research). The amplification reaction was performed following the conditions recommended by the manufacturer. For detection of the fluorescence and interpretation of the genotyping results the QuantStudio 12 K Flex Real-Time PCR System (Applied Biosystems) was used.

All the Taqman assays included in the array showed a genotyping call rate above 95%. The accuracy of the genotyping process was verified including in all plates negative controls and duplicate samples. We observed 100% reproducibility.

Statistical analysis

For all analyses the genotype of each SNP was encoded as 0, 1, 2, where 0 and 2 denotes homozygotes with major allele and minor allele respectively and 1 codes heterozygotes.

The Pearson’s goodness-of-fit χ2 test was used to assess the Hardy-Weinberg equilibrium (HWE). First, to reduce redundancy between SNPs, only one SNP out of all SNPs in strong pairwise linkage disequilibrium (r2 > 0.90) was kept, leading to the final inclusion of 30 SNPs (rs2707466 and rs7988338 were excluded) (Table 1). All possible pairwise SNP-SNP interactions between the selected 30 SNPs (435 possible interactions) were assessed by bivariate linear regression modelling of QUS measurements. Linear regression analyses were adjusted for the following covariates: sex, age, BMI, physical activity and calcium intake. The interaction between SNP1 and SNP2 was defined as the additional effect of their concomitant carriage on QUS over the addition of their independent effects (departure from additivity): βinteraction = βobserved – βexpected with βexpected = βSNP1 + βSNP2, where βobserved is the observed effect of the concomitant carriage of SNP1 and SNP2 (versus no carriage) in a bivariate linear regression, βSNP1 is the effect (β coefficient) of SNP1 in the absence of SNP2, βSNP2 is the effect of SNP2 in the absence of SNP1 in the same bivariate regression framework (bivariate refers only to the number of genetic markers (SNPs), as non-genetic covariates have systematically been included for all analyses). The results were reported as βinteraction (change of the outcome variable QUS with 95% confidence intervals (CIs). To all associations, the highly conservative Bonferroni correction considering the number of tested SNPs was applied. The cut-off value for significance was set at p<0.0001 (0.05/435). The STATA software was used to perform statistical analyses (version 11.0; STATA Corporation, College Station, TX).

Table 1

General information for the studied single nucleotide polymorphisms (SNPs)

Results

The characteristics of the 575 study subjects have been published previously [28]. Table 1 shows marker information, including rsID and minor allele frequency (MAF) for 32 SNPs genotyped. All the SNPs were in HWE and none failed the frequency (MAF < 0.01) or missingess (genotyping >0.05) tests.

The pairwise SNP-SNP interactions tested by linear regression showed multiple interactions (Table 2). The interaction comprising SNPs rs9340799 and rs3736228 that map in the ESR1 and LRP5 genes respectively, revealed the lowest p value (0.001) after adjusting for sex, age, BMI, physical activity and calcium intake. This interaction amounted to a β (95% CI) = 14.289 (5.548, 23.029). Interestingly, we found several interactions between different polymorphisms in these genes, further supporting the presence of epistasis between these two loci (rs556442 in LRP5 and rs9340799 in ESR1; rs2234693 in ESR1 and rs3736228 in LRP5; rs2306862 in LRP5 and rs9340799 in ESR1; rs556442 in LRP5 and rs2234693 in ESR1). However, none of the pairwise SNP-SNP tests reached the significant cut off p value considering the Bonferroni correction for multiple testing (p = 0.0001).

On the other hand, our analyses revealed SNP-SNP interactions within the same genes (rs2908004 and rs3801387 in WNT16; rs3736228 and rs2306862 in LRP5; rs2982552 and rs3020331 in ESR1; rs2982552 and rs9340799 in ESR1).

Discussion

In this study, we aimed to provide a more comprehensive view of the genetic architecture underlying QUS by analysing two-way interactions between common SNPs in candidate genes. A number of SNPs in candidate genes were analysed considering potentially important covariates for bone mineralization process such as, age, sex, calcium intake, physical activity and BMI.

This study is the first to investigate SNP-SNP interactions with heel QUS trait. Our findings identified several SNP-SNP interactions, highlighting that observed between rs9340799 in ESR1 and rs3736228 in LRP5. However, the results from these analyses should be interpreted with caution since SNP-SNP interactions did not reach the highly conservative Bonferroni significance threshold for multiple testing (p<0.0001). Even though these data suggest that epistatic SNP-SNP interactions does not influence heel ultrasound measurements in young Caucasian adults, the possibility of a false negative due to a limited statistical power cannot be excluded.

Oestrogen action in bone are mediated mainly through ERα codified by the ESR1 gene, which has been associated with different osteoporosis related traits in previous studies [33]. The positive effects of the oestrogens on the skeleton have been well established; oestrogens play a major role in the aetiology of osteoporosis by the regulation of bone turnover and inhibition of bone loss [34]. These biological effects are mediated by binding and activation of specific oestrogen receptors (ERα and ERβ) [35]. The rs9340799 variant is located in the first intron of the ESR1 gene and although this intron may contain regulatory elements, the functional implication of this genetic variant remains unknown. In addition, low-density lipoprotein receptor-related protein 5 (LRP5) codifies a protein belonging to the Wnt canonical signal, which regulates osteoblast and osteocyte function [36]. Similarly, the functional relevance of the rs3736228 SNP located in LRP5 is still unknown [37, 38].

The results of this study could be suggestive of a well-regulated cross-talk between ESR1 and LRP5 genes in bone physiology. Although ESR1 and LRP5 develop different roles in skeletal homeostasis, they could be interacting in bone cell biology through synergistic or antagonistic actions. It could be hypothesised that individuals carrying risk alleles at rs9340799 and rs3736228 SNPs have a lower expression of both ESR1 and LRP5, which could lead to an impaired bone mass. As we have analysed a cohort including only young adults, our findings also may suggest that the interaction between rs9340799 and rs3736228 SNPs might be implicated in bone mass accrual. Although previous studies have reported the individual contribution of these polymorphisms to osteoporosis-related phenotypes in early adulthood [30, 38–41], the molecular mechanisms through which these polymorphisms might interact are still under investigation. Future studies are required to elucidate the molecular mechanisms by this interaction is implicated in the bone mass acquisition process during early adulthood.

Besides the interaction between LRP5 and ESR1, our model suggested others such as TMEM135-WNT16, ESR1-DKK1 or OPG-LRP5. Interestingly, all of these genes belong to well characterize pathways implicated in the complex mechanisms that regulate bone formation. Thus, a novel line of research would be very interesting to elucidate how they may be interacting at the molecular level.

On the other hand, SNP-SNP interactions between variants that were not in LD within the same gene were observed. That was the case of genes WTN16, LPR5 and ESR1. Thus, beyond the individual effect of SNPs in a gene, epistatic SNP-SNP interactions in the same gene occur and could be a potential factor contributing to the unexplained heritability of bone mass acquisition.

On the other hand, it could be hypothesised that the effects of SNP-SNP interactions would be weak and therefore, a very large population size would be needed to detect associations surpassing the p value threshold for multiple testing (p<0.0001). Thus, a limitation of the present study could be the population size included that would not reach enough statistical power to detect epistatic SNP-SNP interactions with weak effect. Therefore, in order to support our preliminary results, further studies including larger population size and meta-analysis are needed. On the other hand, the current approach considered only two-locus interactions and therefore more complex interactions between three or more genetic markers are unknown. Again, larger studies are required to confirm our findings and to elucidate mechanisms by which the genetic interaction between these genes influences quantitative bone phenotypes in early adulthood.

Conclusions

Our analysis of SNP-SNP interaction in candidate genes of QUS in Caucasian young adults reveal several interactions, especially between ESR1 and LRP5 genes, that did not reach statistical significance. Although our results do not support a relevant genetic contribution of SNP-SNP epistatic interactions to QUS in young adults, further studies in larger independent populations would be necessary to support these preliminary findings.

Abbreviations

BMD:

Bone mineral density

BMI:

Body mass index

BUA:

Broadband ultrasound attenuation

CIs:

Confidence intervals

DCI:

Dietary calcium intake

DXA:

Dual-energy x-ray absorptiometry

GWAS:

Genome wide association studies

HWE:

Hardy-Weinberg equilibrium

IPAQ:

International Physical Activity Questionnaire

MAF:

Minor allele frequency

PBM:

Peak bone mass

QUS:

quantitative ultrasound

Declarations

Acknowledgements

Not applicable.

Funding

This study was supported by a grant PI-0414-2014 from Consejería de Salud (Junta de Andalucía, Spain). Correa-Rodríguez M is a predoctoral fellow (FPU13/00143) from the Ministerio de Educación, Cultura y Deporte (Programa de Formación del Profesorado Universitario).

Availability data and materials

The datasets during and/or analysed during the current study available from the corresponding author on reasonable request.

Authors’ contributions

M C-R and J S-RV monitored data collection and drafted and revised the paper. S V and J M wrote the statistical analysis plan, cleaned and analyzed the data. B R-M and G O analyzed the data, and drafted and revised the paper. All authors have approved the manuscript.

Ethics approval and consent to participate

This study met the institutional and/or national research committee as well as 1964 Helsinki declaration and its later amendments ethical standards for procedures involving human participants. All study participants gave their written informed consent and the local ethics committee of the University of Granada (Spain) approved the study (78/CEIH/2015).

Consent for publication

Not applicable.

Competing interests

Gisela Orozco is a member of the Editorial Board. The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Authors’ Affiliations

(1)

Faculty of Health Sciences, University of Granada, Granada, Spain

(2)

Arthritis Research UK Centre for Genetics and Genomics, Centre for Musculoskeletal Research, Faculty of Biology, Medicine and Health, Manchester Academic Health Science Centre, The University of Manchester, Manchester, UK

Correa-Rodríguez M, Schmidt-RioValle J, Rueda-Medina B. The rs3736228 polymorphism in the LRP5 gene is associated with calcaneal ultrasound parameter but not with body composition in a cohort of young Caucasian adults. J Bone Miner Metab. 2016;27Google Scholar