Abstract

The human aging is marked by several body changes, including in bone mineral density (BMD). Research shows that microRNAs are important modulators of bone metabolism. The present research aims to analyze the whole blood concentration of 10 selected microRNAs (miRs) and their association with absolute and relative scores of BMD in specific osseous site of Brazilian very-old adults. Forty noninstitutionalized and apparently healthy, very old (≥80 years) outpatients were eligible for research. Anthropometry, biochemistry, and densitometry measurements were performed along with coronary artery calcification (CAC) scores and tested across total circulating levels of microRNAs. As expected, the relative BMD scores for the lumbosacral region (L1 to S5) and for the femoral head and neck observed in the sample denote weakened bone architecture, compatible with prevalent osteopenia and osteoporosis. In this context, one single significant association was found, and negatively implicated the miR-34a-5p with both absolute (β = −0.36, for BMD) and relative (β = −0.43, for T-score) densitometry indexes of the femoral head (adjusted to sex and physical activity practice), but not with the other sites. No difference in total blood concentrations of the miRs was found according to CAC scores. Our findings indicate greater circulating levels for miR-34a-5p among very-old adults who display the lowest scores of BMD, being a finding consistent with a modest contribution of the miR (along with co-variables) to the mineralization of that site. Attesting clinical relevance of our findings demands forthcoming studies.

1. Introduction

The human aging is marked by several body changes and greater risk for disorders as hormonal and metabolic imbalances as well as sarcopenia and osteoporosis. This last entity is characterized by progressive bone deterioration and a concomitant increase in risk for fractures among very-old men and women [1, 2]. Osteoporosis is mostly (but not solely) derived from loss of bone mineral density (BMD) and is increasingly potentiated by meno/andropausal-related hormonal changes, sedentary lifestyle, and other age-related conditions [3, 4].

BMD is maintained by the precise balance between formation and reabsorption of calcium in all bones, a balance accomplished by osteoclasts and osteoblasts [5]. Current research shows that microRNAs (miRs) are important modulators of bone metabolism [6, 7], with miRs as small non-protein-encoding RNA molecules containing 18 to 23 nucleotides in average, in one single-stranded chain [8, 9].

Being responsible for the posttranscriptional silencing of thousands of genes, mostly by the inhibition of protein translation, miRs are associated with virtually all biological process, including the regulation of the mineral content in osseous structures [610]. Since the onset of the bone mineralization process, several in vitro studies point out developmental roles for miRs [1114], whereas experimental studies show differential in situ expression of these elements in osteoporotic bones [1517]. By now, a body of evidence shows that a specific set of miRs (including miR-1, miR-92a, miR-126, miR-130a, miR-146a, miR-155, and miR-221, for instance) directly associate with the age-associated decline in BMD [1822]. In this context, the miR-34a is suggested as a powerful suppressor of osteoclastogenesis and of calcium reabsorption in bones [23].

In this context, there is report of overexpression of the miR-21 and miR-100 in blood samples of patients with osteoporosis, in a manner that is proportional to scores of bone mineral density [24]. This evidence brings about a plausible role of miRs as potential biomarkers for decline of BMD and, if implicated in the osteoporosis pathophysiology, as eligible targets for preventive and therapeutic strategies of the condition. Thus, the present research has main objective to analyze the association of these 10 selected miRs with scores of bone mineral density of Brazilian very-old adults.

2. Materials and Methods

2.1. Subjects

All patients resided in the Brazilian Federal District, located in the Midwest Brazil. This cross-sectional report analyzed data from consecutive noninstitutionalized apparently healthy ambulatory patients, aged 80 years or more, participants in the ongoing Brazilian Study on Healthy Aging, a prospective study designed to identify markers for primary prevention of vascular events and other prevalent health outcomes in geriatrics settings [25, 26]. For this report, a subset that has never manifested bone fractures or recurrent falls was analyzed. Additional selection criteria were the absence of autoimmune disease (including rheumatic disorders), chronic or recurrent infections, prior or current neoplastic disease, or use of steroidal or nonsteroidal anti-inflammatory drugs in the past 30 days. Use of BMD-related drugs was investigated for each patient. The present report presents analyses of admission data of 40 subjects who fulfilled inclusion criteria. None were under any particular dietary regimen or were regular practitioners of physical exercises.

This study was approved by the institutional ethics committee and conducted following the Helsinki Declaration. Participation was voluntary, and written informed consent was obtained from each participant.

2.2. Anthropometric Characteristics

Anthropometric assessments were body mass index (BMI); weight (kg)/height (m2); and waist circumference (WC), according to the standards recommended by the World Health Organization [27].

2.3. Biochemical Analysis

At baseline, and after 12 h of overnight fasting, each participant underwent blood drawing for biochemical analysis and freezing at −80°C of serum and whole blood samples. For metabolic analysis, the following measures were performed using reagents from Roche Diagnostics (Mannheim, USA) and carried out by the same certified clinical laboratory according to routine clinical biochemistry, being expressed in standard units: blood glucose, glycated hemoglobin (HbA1c), total cholesterol and fractions, triglycerides, total serum protein, high sensitivity C-reactive protein, aspartate aminotransferase, alanine aminotransferase, alkaline phosphatase, gamma-glutamyl transferase, uric acid, and serum creatinine. Creatinine clearance was estimated using serum creatinine as described by Cockcroft and Gault [28].

Also, each individual underwent a detailed clinical examination. Type-2 diabetes was characterized by fasting blood glucose ≥126 mg/dl or use of insulin or oral antidiabetic drugs [29]. Hypertensive cases were defined following to the VII Brazilian Guidelines of Hypertension [30]. Metabolic syndrome was defined by the criteria of the National Cholesterol Education Program Adult Treatment Panel III [31]. Consumption of drugs to treat BMD-related conditions was investigated and recorded. Subjects were considered physically fit according to standards set by the World Health Organization for prescription of exercises for older adults [32].

2.4. Dual Energy X-Ray Absorptiometry Measurements

Participants underwent bone mineral density (BMD) assessments of the femoral neck, femoral head, and lumbarsacral (L1 to S5) regions with dual energy X-ray absorptiometry (DXA; (Lunar Prodigy Advance, GE Healthcare, USA) according to standard protocol to measure BMD (g/cm2). The difference between an individual’s BMD and the mean BMD for a reference population was expressed in standard deviation term (T-score). Also, fat mass (kg) was incidentally measured during the procedure.

2.5. Cardiac Computed Tomography

Coronary artery calcification (CAC) was assessed by computed tomography, performed in a 64-slice scanner (Aquillion 64, Toshiba, Ottawara, Japan). Axial slices of 3 mm thickness with 3 mm table-feed were acquired at 70% of R-R interval with prospective electrocardiography triggering. CAC was defined as a minimum of 3 contiguous pixels with a peak Hounsfield unit density >130, being scored by a certified radiologist using the Agatston score to express its extent.

2.6. MicroRNA Selection

To select the set of microRNAs to be investigated, we have searched the DIANA Tools (TarBase v7.0) [33] for miRs validated as interacting with target genes encoding key noncollagenous, mineralization-associated proteins of the bone matrix, as described by Javed et al. [34]. A manual search was performed by entering Ensembl code individually, with each human(hsa)-miR being selected for analyses when the database reported an association (positive result) and the source tissue specified as bone. From this strategy, 10 miRs were rendered eligible for analyses, as follows: miR-1-3p, miR-21-5p, miR-34a-5p, miR-92a-3p, miR-100-5p, miR-126-3p, miR-130a-3p, miR-146a-5p, miR-155-5p, and miR-221-3p.

2.7. MicroRNA Analysis

All blood samples were obtained by puncture of the antecubital vein with 1 ml aliquots immediately frozen and stored at −80°C until further processing. Before starting the extraction procedure, samples were thawed on ice. Total RNA (including small RNAs) was purified from 700 µl blood samples using the miRVana™ PARIS™ kit (Thermo Fisher Scientific, Wilmington, DE, USA) adapted by adding multiple acid-phenol : chloroform extraction steps to produce clearer aqueous phases. Concentration and optical quality of samples were determined using the NanoDrop Lite. Reactions of cDNA synthesis for each miR were performed with the TaqMan® microRNA Reverse Transcription Kit (Thermo Fisher Scientific) and the appropriate TaqMan® RT primer, with 30 min incubations at 16°C and then at 42°C, terminated at 85°C for 5 min. The quantitative real-time PCR (qPCR) step was performed in 10 μl reaction mixtures using a TaqMan Universal PCR Master Mix (Thermo Fisher Scientific) with specific TaqMan® assays (primers and probes) and the ABI 7500 equipment (Applied Biosystems, Foster City, CA, USA), following the manufacturer’s instructions. RNU43 was used as endogenous control to normalize Ct values obtained for each target miR, with expression compared by means of the 2−ΔCt method [35].

2.8. Statistical Analysis

To address the aim of evaluating the occurrence and strength of the association between whole blood, circulating levels of miRs, and scores (continuous and categorical) of anthropometric, clinical, and biochemical traits of potential confounding effect in the main model, our statistical analyses started at obtaining correlation coefficients between these variables. For that, normal distribution of all variables was assessed using the Kolmogorov–Smirnov test.

The association between continuous, normally distributed variables was evaluated using Pearson’s correlation test, whereas the involvement of a least one categorical or nonnormally distributed continuous variable was dealt using a linear regression model. For these latter analyses, absence or presence of a given feature was represented by 0 or 1, respectively. When a variable was nonnormally distributed, a logarithmic transformation was applied. Whenever an interaction was noticed, regression analyses were run using adjustment for the confounding variable(s) or condition(s). For these analyses, a value (two-tailed) was rendered significant following Bonferroni’s principle of correction for multiple comparisons whenever a given trait is tested across k independent variables (e.g., if k = 10 tests, α = 0.005).

Also, Student’s t-test was used to verify if miRs concentrations varied across the subjects grouped according to the median scores of CAC and BMD in the sample, with adjustment to covariates when appropriate. Finally, linear multivariate regression analysis, the stepwise method, was performed to assess at which extent the whole blood miRs concentrations explained the variability in the calcification indexes. All analyses were performed with the Statistical Package for the Social Sciences (SPSS) for Windows (version 17.0).

3. Results

After applying the exclusion criteria, 40 very-old subjects (men and women) were eligible for research. The common clinical, anthropometric, and metabolic features of the subjects investigated are summarized in Table 1. Most notably, average relative bone mineral density scores for the lumbosacral and the femoral head and neck were expressed in negative values, denoting weakened bone architecture among the subjects, compatible with the expected prevalence of osteopenia and osteoporosis.

Table 2 presents correlation and regression analyses of circulating levels of miRs with common clinical, anthropometric, and metabolic features of the very-old adults investigated. It became evident that the overall serum levels of miRs investigated vary neither according to the common characteristics as age, sex, or body fat nor across conditions as type-2 diabetes and hypertension (among other traits) presented by the volunteers. Complementarily, average miRs levels were also compared using the t-test across groups of dichotomic variables (sex, physical activity, T2DM, hypertension, and metabolic syndrome), and no difference was found (not shown). On the other hand, similar analyses found that average scores of BMD at all osseous sites varied according to sex () and physical activity practice () (not shown in Tables), with these traits being included as covariables in the forthcoming inferential analyses.

The linear regression analyses of circulating levels of miRs with the mineral content in the coronary and the osseous territories are summarized in Table 3. According to the significance threshold adopted herein, one single relevant relationship was found and negatively implicated the miR-34a with both absolute (β = −0.36, for BMD) and relative (β = −0.43, for T-score) mineral density scores in the femoral head. Similar associations were observed neither for the lumbosacral region nor the femoral neck.

Figure 1 presents comparison of log-transformed levels of all three microRNAs with greater trends for an association (from the prior analyses) with the femoral head BMD (miR-34a, miR-126 and miR-146) across individuals grouped according to the median value (0.82 g/cm2) of bone mineral density in the femoral head. Results reassure that a significant difference could only be devised for miR-34a, with a greater concentration of circulating level for miR-34a among carriers of lower BMD. No difference was found in total blood concentrations of any of the miR tested across the sample split according to median CAC scores, with or without adjustments.

4. Discussion

A baseline finding in our setting consisted in mean negative values for relative scores of BMD on all osseous sites investigated (lumbosacral and femoral head and neck) in the sample of very-old adults. Regardless of this expected scenario of an important age-related decrease of BMD, a remarkable and specific negative association of miR-34a with the bone density scores of the femoral head could be devised. Also, as expected, mineral density scores were significantly higher for men than for women at all osseous sites investigated (). But, our results seem not biased by admitting both sexes since the male/female ratio did not vary across the sample slip according to median scores of BMD (χ2 = 2.157; ; for femoral head as example) or CAC (χ2 = 2.837; ). Thus, we suggest that the miR-34a (among all miRs studied) exhibited a consistent relationship with the mineralization status of the femoral bone, so to be detected out of a sampling of the peripheral bloodstream. Stepwise multivariate regression showed that the miR-34a levels along with sex and physical activity accounted for roughly 58% of the variance in the mineral density of the femoral head.

It has been demonstrated that an in vitro inhibition of miR-34a by means of antisense-miR was responsible for an increased differentiation of osteoblasts and, consequently, an increase in cultured bone matrix formation [11, 23]. On the other hand, overexpression of miR-34a caused an inverse effect. Thereby, research elsewhere already suggests that in situ expression of mirR-34a it is negatively correlated to BMD, rendering levels particularly high of this element as a plausible biomarker for osteoporosis [15]. Nonetheless, the actual physiological mechanism(s) on how miR-34a acts to suppress matrix mineralization is yet to be fully described.

In this respect, there is evidence that the expression of transforming growth factor-β-induced factor 2 (Tgif2) during in vitro formation and differentiation of both animal and human osteoclasts is targeted by miR34a [11]. Tgif2 is soundly acknowledged as a proosteoclastogenic factor and miR34a inhibits its expression, consequently leading to decreased calcium reabsorption and bone matrix renewal [23]. These results are consistent with our findings on elderly with lower BMD showing more pronounced total circulating levels of miR-34a.

Despite the homogeneity of our subjects on what concerns age strata, clinical conditions, and biochemical profile, our investigation has limitations greatly related to factors inherent to a developing country as Brazil and that were not investigated herein. The remarkable interethnic variation owing to genetic admixture is certainly one of these factors [36]. Also, diversity in functional capacity and in food intake cannot be ruled out as confounding factors. Besides all, the most remarkable limitation may be represented by the difficulty of relating findings in a peripheral blood sample with pathophysiological processes in specific tissues/organs. In this respect, one should note that bone-derived exosomes are released into the extracellular environment [37], and it is plausible that at least part of the microRNA content in the whole blood may reflect concentrations (and imbalances) existing elsewhere, including the bone millieu.

So far, the authors are unaware of other studies that have correlated total circulating levels of microRNAs in the blood with bone mineral density of very-old women and men. Our aged human model suggests that total concentration of miR-34a explains at least in part the complex phenotype of bone structure. However, attesting clinical relevance of the finding presented herein still depends on forthcoming studies, so as to assess its diagnostic value (if any) for future clinical practice.

Data Availability

Data can be made available upon reasonable request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Authors’ Contributions

OT Nóbrega performed all microRNAs assessments and analyzed and interpreted the results. NI Viana and ST Reis assisted and supervised, respectively, the laboratory component. WM Freitas and AC Sposito executed and supervised, respectively, the medical and radiological components. Diego IV Perez assisted with statistical analyses. KRM Leite and M Srougi designed and coordinated the study. OT Nóbrega and GS Morais Junior participated in the preparation of the original manuscript.

Acknowledgments

This research was supported by CNPq (grant #445692/2014-6), FAPDF (grant #193.000.967-2015), and CAPES (Finance Code #001). A.C. Sposito received a fellowship for productivity in research and O. T. Nóbrega received a senior post-doctoral fellowship, both from CNPq, Brazil. G.S. Morais-Junior had a stipend from CAPES, Brazil.