- Research article
- Open Access
Age-related microRNAs in older breast cancer patients: biomarker potential and evolution during adjuvant chemotherapy
BMC Cancer volume 18, Article number: 1014 (2018)
MicroRNAs (miRNAs) are important regulators of cellular function and have been associated with both aging and cancer, but the impact of chemotherapy on age-related miRNAs has barely been studied.
Our aim was to examine whether chemotherapy accelerates the aging process in elderly breast cancer patients using miRNA expression profiling.
We monitored age-related miRNAs in blood of women, aged 70 or older, receiving adjuvant chemotherapy (docetaxel and cyclophosphamide, TC) for invasive breast cancer (chemo group, CTG, n = 46). A control group of older breast cancer patients without chemotherapy was included for comparison (control group, CG, n = 43). All patients underwent geriatric assessment at inclusion (T0), after 3 months (T1) and 1 year (T2). Moreover, we analysed the serum expression of nine age-related miRNAs (miR-20a, miR-30b, miR-34a, miR-106b, miR-191, miR-301a, miR-320b, miR-374a, miR-378a) at each timepoint.
Except for miR-106b, which behaved slightly different in CTG compared to CG, all miRNAs showed moderate fluctuations during the study course with no significant differences between groups. Several age-related miRNAs correlated with clinical frailty (miR-106b, miR-191, miR-301a, miR-320b, miR-374a), as well as with other biomarkers of aging, particularly Interleukin-6 (IL-6) and Monocyte Chemoattractant Protein-1 (MCP-1) (miR-106b, miR-301a, miR-374a-5p, miR-378a-3p). Moreover, based on their ‘aging miRNA’ profiles, patients clustered into two distinct groups exhibiting significantly different results for several biological/clinical aging parameters.
These results further corroborate our earlier report, stating that adjuvant TC chemotherapy does not significantly boost aging progression in elderly breast cancer patients. Our findings also endorsed specific age-related miRNAs as promising aging/frailty biomarkers in oncogeriatric populations.
ClinicalTrials.gov, NCT00849758. Registered on 20 February 2009. This clinical trial was registered prospectively.
(Breast) cancer treatment in the elderly represents a major challenge in clinical oncology. Despite their important representation within the overall breast cancer (BC) population, older women are often excluded from standard BC medical treatment regimens, or are offered less aggressive (and possibly less effective) therapies. This is attributable to concerns about increased risk of side effects, as well as decreased medical fitness of older cancer patients, either perceived or assumed by subjective clinical evaluation. According to international guidelines for patient management in geriatric oncology, it is not justified to base treatment choices on chronological age only . In addition to the already available geriatric assessment (GA) tools, new clinical and biological tools, including circulating biomarkers of aging, are currently being developed to better assess global health and individualize therapeutic strategies [2,3,4,5,6,7,8].
Besides the risks that chemotherapy imposes to frail patients, the impact that it may have on the aging process of more fit elderly patients is still not well understood. Based on the observation of cellular senescence induced by cytotoxic agents (including anticancer molecules) , researchers have been investigating possible effects of chemotherapy on aging progression. Data obtained from follow-up studies have shown that adult survivors of childhood cancers develop degenerative diseases (typical of old age) earlier in life and with a higher incidence compared to the general population . Chemotherapy-induced aging and frailty are assumed to be caused, among other factors, by production of free radical intermediates, persistent DNA damage not counterbalanced by adequate DNA repair mechanisms, and/or senescence triggered by telomere instability . In a small cohort of patients with head and neck cancer, telomeres of peripheral blood mononuclear cells were reported to be severely shortened after combined chemo-radiotherapy, especially in older patients . Similar results were also obtained in patients who underwent chemotherapy for BC  and non-Hodgkin lymphoma .
To investigate the hypothesized acceleration of the aging process by cancer treatment we have recently conducted a prospective clinical study in older (70+) BC patients, that monitored the evolution during adjuvant chemotherapy of both clinical aging parameters (GA) and various biological markers described in the literature as potential biomarkers of aging. This biomarker panel included leukocyte telomere length, the inflammation-related plasma cytokines/chemokines IL-6, Tumor Necrosis Factor-alpha (TNFα), Interleukin-10 (IL-10), Regulated on Activation, Normal T-Cell Expressed and Secreted (RANTES)/ C-C motif chemokine-5 (CCL5) and MCP-1/C-C motif chemokine-2 (CCL2) as markers of ‘inflammaging’, and the aging-related protein Insulin-like Growth Factor-1 (IGF-1) [5, 15,16,17,18,19,20,21,22,23].
Clinical and biomarker assessments were accomplished at three different time points, i.e. prior to initiation of chemotherapy, after 3 months (last cycle of chemotherapy) and after 1 year. The study comprised a cohort of BC patients receiving docetaxel/cyclophosphamide treatment and a control cohort receiving only hormone treatment. The results of the above-mentioned work have been published recently .
In the present study, we attempted to further expand our insight into the potential connections between adjuvant chemotherapy and biological aging, by analysing the behaviour of circulating aging-related miRNAs previously identified by our group .
Indeed, in addition to the above-mentioned aging biomarkers, several miRNAs have also been implicated in the aging process [25, 26]. MiRNAs are small molecules (18–22 nucleotides) belonging to the wider class of non-coding RNAs, which exert non-permanent epigenetic functions via the post-transcriptional regulation of gene expression . Mutations and/or altered expression of miRNAs have been associated with several pathological conditions, including cancer [28,29,30,31,32,33,34,35,36,37,38,39,40]. Moreover, miRNAs may play a role in the aging process, as some of them are part of molecular pathways that regulate cellular senescence. The expression of miR-34a, for instance, is induced by p53, and high levels of this miRNA result in growth arrest [41, 42]. This miRNA has also been positively associated with myocardial aging. Its primary target is the longevity-associated deacetylase SIRT-1  that regulates chromatin remodelling, stress responses, DNA repair and insulin regulation.
Several of the miRNAs reported in the literature to be associated with aging have also been found to be over- or underexpressed in cancer patients [29, 34, 35, 39, 40]. To identify plasma miRNAs with aging biomarker potential within a cancer population, we have first carried out a plasma miRNA screening study involving young and older BC patients . From this pilot study, we have established a panel of miRNAs that were found to be differentially expressed according to age among BC patients: miR-20a, miR-30b, miR-34a, miR-106b, miR-191, miR-301a, miR-320b, miR374aand miR-378a-3p. Here, we have examined the quantitative fluctuations of these ‘aging miRNAs’ during and after adjuvant chemotherapy.
From 2009 to 2012, women at least 70 years old who were affected by locally-advanced, non-metastatic BC and eligible for adjuvant systemic chemotherapy were enrolled at 5 hospitals in Belgium, henceforth referred to as “chemotherapy group” (CTG). This group consisted originally of 57 patients, and miRNA data were available at the 3 time points in 46 patients.
In parallel, a comparable series of patients, not eligible for systemic chemotherapy, but only for endocrine treatment were included as a “control group” (CG) . This group consisted originally of 52 patients, and miRNA data were available at the 3 time points in 43 patients.
The antineoplastic therapy administered to the CTG consisted of docetaxel at a dose of 75 mg/m2 plus cyclophosphamide at 600 mg/m2 every 3 weeks for a total of 4 cycles (TC scheme)  . G-CSF (granulocyte-colony stimulating factor) was administered at each cycle, according to the National Comprehensive Cancer Network (NCCN) guidelines . The planned adjuvant treatment also included an aromatase inhibitor (to be administered after chemotherapy completion) in case of hormone-sensitive tumors, and trastuzumab administration in case of Her2 positive tumors.
Conversely, patients in the CG received an aromatase inhibitor as sole medical treatment. In both groups, radiation therapy was either or not administered according to institutional practice. Enrolment took place after breast surgery. Blood samples were collected at three time points: T0: between 3 and 6 weeks after surgery, always before the first cycle of chemotherapy; T1: 3 months after inclusion (in principle the day of the fourth and last cycle of chemotherapy for patients in the CTG); T2: 1 year after inclusion. At each time point, patients also underwent clinical geriatric evaluation.
This study was approved by the local ethics committees of Clinique Sainte Elisabeth (Namur), Imelda Hospital (Bonheiden), Jessa Hospital (Hasselt), Zol Ziekenhuis Oost-Limburg (Genk), Jules Bordet Institute (Brussels), and by the University Hospitals Leuven central ethics committee. All patients signed a written informed consent in accordance to the Helsinki Declaration. This clinical trial was registered prospectively.
This article adheres to CONSORT guidelines, where applicable.
Clinical geriatric evaluation
Detailed information on GA tools and results of clinical evaluation accomplished at the different time points have been extensively documented in our primary study on aging biomarkers and chemotherapy . Briefly, patients were screened at baseline with the G8 screening tool and the Flemish version of the Triage Risk Screening Tool (fTRST), and social data were collected (age, living situation, marital status and educational level).
At each time point also a geriatric assessment (GA) was performed, as well as a frailty assessment with Balducci Frailty score and Leuven Oncogeriatric Frailty Score (LOFS) . Moreover, we assessed performance status according to the Eastern Cooperative Oncology Group - Performance Status (ECOG-PS), and quality of life (QoL) using the EORTC QLQ-C30 questionnaire.
At both T1 and T2, adverse events (using the CTCAE v4.0 classification) and unexpected hospitalizations were also monitored.
Our primary endpoint was to evaluate whether aging-related miRNAs changed during the study period and, if so, whether relevant differences could be detected over time between CTG and CG.
As secondary endpoints, we also assessed (i) potential correlations of miRNAs measured at inclusion (T0) with chronological age, clinical geriatric assessment parameters and aging biomarkers reported in our primary paper ; (ii) whether miRNAs at inclusion had a predictive value towards acute and/or irreversible decline in functionality and in QoL; (iii) whether miRNAs at inclusion predicted toxicity and unexpected hospitalizations during and after chemotherapy; (iv) correlation patterns between the 9 miRNAs and possible formation of patient clusters based on miRNA expression profiles.
Blood sample collection and processing
At each time point, 4-mL whole blood specimens were collected from each patient in BD Vacutainer SST II Advance serum tubes. After incubation at room temperature for 20 to 60 min, the blood samples were centrifuged at 1300×g for 10 min at 4 °C and supernatants (serum) were aliquoted and stored at − 80 °C.
Aging biomarker analysis
Methodology and results of biomarker analyses (i.e. leukocyte telomere length, circulating IL-6, IL-10, TNFα, RANTES/CCL5, MCP-1/CCL2 and IGF-1) performed at the 3 time points were described in detail in our previous publication .
Isolation of miRNAs from serum
For each sample, 250 μL of serum were thawed on ice and then centrifuged at 3000 x g for 5 min to remove debris. To 200 μL of supernatant, 1 μg of carrier MS2 RNA (Roche) was added in order to stabilize RNA during extraction and cDNA synthesis. Moreover, 1 μL of the synthetic RNA spike-in UniSp6 was added to allow evaluation of the efficiency and uniformity of the entire RNA extraction/cDNA synthesis procedure. Then, miRNAs were isolated with the miRCURY™ RNA Isolation Kit–Biofluids (Exiqon), following the manufacturer’s instructions. Spin columns were finally eluted twice with 25 μL DNase/RNase-free water each time. Both eluates were pooled and stored at − 80 °C.
cDNA synthesis and qPCR
On each purified miRNA sample, cDNA synthesis was performed in duplicate. 2 μL of miRNA extract were processed using the Universal cDNA synthesis kit II (Exiqon), according to the manufacturer’s instructions. cDNA samples were stored at − 80 °C until PCR analysis.
Measurement of relative amounts of transcripts was carried out by real-time qPCR analysis using Pick-&-Mix microRNA PCR panels (96 well Ready-to-Use custom plates) with Exilent SYBR® Green Master Mix (Exiqon). For each RNA sample, both duplicate cDNAs were assessed in a single plate. Every plate included primers for: 5 reference miRNAs (hsa-miR-23a-3p, hsa-miR-29a-3p, hsa-miR-29c-3p, hsa-miR-140-3p, hsa-miR-484, further referred to as miR-23a, miR-29a, miR-29c, miR-140 and miR-484) used for data normalization, the 9 aging-related miRNAs of interest selected for the study (see below) and the synthetic spike-in UniSp6 to allow evaluation of miRNA extraction/reverse transcription efficiency. The 9 aging miRNAs included 3 that were previously shown to increase with aging (hsa-miR-34a-5p, hsa-miR-320b and hsa-miR-378a-3p, further referred to as miR-34a, miR-320b and miR-378a), and 6 that previously showed negative association with aging (hsa-miR-20a-3p, hsa-miR-30b-5p, hsa-miR-106b-5p, hsa-miR-191-5p, hsa-miR-301a-3p and hsa-miR-374a-5p, further referred to as miR-20a, miR-30b, miR-106b, miR-191, miR-301a and miR-374a) . All plates additionally included the interplate calibrator UniSp3, in order to allow detection of global amplification differences due to inter-run variation.
Assays were carried out following the manufacturer’s specifications. Briefly, cDNA was diluted 50x in nuclease-free water and mixed with an equal volume of 2x Exilent SYBR Green master mix (Exiqon). Final reaction volume was 10 μL. Plates were run on a LightCycler 480 (LC480, Roche) instrument applying the following thermal cycling protocol: activation step (10 min at 95 °C); 45 amplification cycles (10 s at 95 °C, 1 min at 60 °C, ramp rate 1,6 °C/s); melting curve analysis.
Quality control and processing of PCR data
As haemolysis can alter the relative amounts of different serum miRNAs through the release of intracellular miRNAs from erythrocytes, a quality control was performed using the miR-451/miR-23a-3p hemolysis test . Following qPCR analysis of the expression of miR-451 (highly expressed in erythrocytes) and miR-23a-3p (stably expessed in biofluids), all samples with a Delta Cp value (Cp miR-451 minus Cp miR-23a-3p) higher than 5 were excluded from further analysis. Samples with borderline results (Delta Cp between 4 and 5) were double-checked for haemolysis using a second method. Thawed serum samples were briefly spun down to remove debris, and then the absorbance spectrum was assessed on a Nanodrop ND-1000. Samples showing an apparent absorption peak at 415 nm (the hemoglobin absorption maximum) were excluded from the study. The qPCR data were processed using the MultiD GenEx software. We visually inspected expression profiles of all miRNAs and the UniSp6 across all samples on a bidimensional line plot. Samples with a clearly deviating expression for the entire miRNA panel were excluded from further analysis. Normalization was performed using the 5 reference transcripts miR-23a-3p, miR-29a-3p, miR-29c-3p, miR-140-3p and miR-484. These were validated in our previous paper  as suitable reference miRNAs for serum/plasma samples by the algorithm tools GeNorm and NormFinder and were now again confirmed to be stably expressed across all serum samples. Technical repeats (duplicate cDNAs per sample) were averaged and finally, all values were converted to relative quantities and then log-transformed (Log2 scale).
For the primary endpoint, miRNAs were modelled as response variables in linear models for repeated measures with time, group and their interaction as explanatory variables. An unstructured residual covariance matrix was modelled to account for clustering.
The evolution over time in the CTG was assessed by estimating the change in miRNA level between inclusion (T0) and 3 months (T1) and between inclusion (T0) and 12 months (T2). Results were presented by the mean change between time points with 95% confidence interval (CI). The difference in evolution between chemotherapy and control patients was assessed by a test for group by time interaction.
For the secondary endpoints, Spearman correlations were used for studying univariable association of miRNAs with continuous or ordinal variables. Kruskal-Wallis tests were used to compare miRNA levels between more than 2 groups, and Mann-Whitney U tests for comparisons between two groups. Multivariable models: a backward selection procedure was applied for selecting a set of miRNA as independent predictors of response variables (age, clinical aging parameters and aging biomarkers). Linear regression was used for continuous variables, logistic regression for binary variables, and proportional odds models for ordinal variables. Mann-Whitney U tests were used for comparing miRNA levels between patients with and without decline in functionality, unexpected hospitalization or grade II-III-IV toxicity.
The association between the miRNAs were studied by means of Pearson correlations. To identify groups (clusters) of patients with similar miRNA profiles, a disjoint cluster analysis was performed based on minimizing the sum of squared (euclidian) distances from the cluster means; miRNA values were standardized for this analysis. To decide upon the number of clusters, we took into account the pseudo F statistic (larger means better fit) and the number of patients per cluster. The SAS procedure FASTCLUS was used for this analysis. Analyses were performed for data measured at inclusion. Mann-Whitney U tests were used for comparing patients within two clusters on ordinal or continuous variables. Fisher exact tests were used for comparing clusters on categorical or binary outcomes (decline, hospitalization, toxicity).
All tests were two sided, and a 5% significance level was considered for all tests.
All analyses have been performed using SAS software, version 9.4 of the SAS System for Windows. Copyright © 2002 SAS Institute Inc. SAS and all other SAS Institute Inc. product or service names are registered trademarks or trademarks of SAS Institute Inc., Cary, NC, USA.
Figures were performed using using GraphPad Prism version 6.00 for Windows, GraphPad Software, La Jolla, CA, USA.
Evolution of aging miRNAs over time during BC treatment
For each miRNA, time evolution in both CTG and CG is shown in Fig. 1 and the corresponding statistics are summarized in Table 1. Several miRNAs showed significant changes in patients of CTG during the course of the study: miR-34a was increased at T1 (p = 0.0039) while miR-30b and miR-374a were decreased at T1 (p = 0.0003 and 0.0394, respectively). For miR-374a, these changes appeared to be transient: the initial miRNA levels measured at inclusion (T0) were restored after 1 year (T2). In contrast, the observed changes of miR-30b and miR34a still persisted after 1 year, albeit somewhat less pronounced (mean differences were − 0.22 at T2 versus − 0.39 at T1 for miR-30b and 0.60 at T2 versus 0.76 at T1 for miR-34a). However, for none of these three miRNAs, a significant difference in evolution over time could be demonstrated when comparing CTG with CG, as indicated by the lack of a statistically significant group by time interaction. Plasma levels of miR-106b and miR-191 were found to be slightly increased (p = 0.0497), respectively decreased (p = 0.0272) in CTG at T2 but not at T1. Moreover, a significant group by time interaction (p = 0.024) for miR-106b may point to a different evolution of this miRNA in CTG versus CG. No significant modifications were observed for the other miRNAs (miR-20a, miR-301a, miR-320b, miR-374a, miR-378a) during the time course of the study (Fig. 1, Table 1). Overall, the observed differences over time did not appear to depend on the type of administered treatment, hinting toward the lack of an effect of chemotherapy on aging in the analysed population. Of note, 3 miRNAs (miR-20a, miR-301a, miR-320b) were significantly different at baseline (in the direction of increased aging) in CG compared to CTG, corresponding to the fact that clinical aging was also slightly more pronounced in CG (see Table 3 from our recent publication ).
Association of aging miRNAs with patient’s chronological age
In initial univariable analyses, patient age at inclusion, strongly tended to be associated with several previously identified aging related miRNAs  measured at T0 (i.e. miR-30b, miR-374a, miR-106b, miR-301a, miR-320b), even within this ‘purely elderly population’ (Table 2).
In a next step, a backward multivariable model selection procedure was applied, resulting in a model with miR-301a as the only independent explanatory variable for age. The model revealed a negative association between miRNA-301a and age: higher age is associated with lower miR-301a levels (p = 0.0006) (Table 3), which is in accordance with our previous findings .
Association of miRNAs with clinical aging
Next, we investigated whether the ‘aging miRNAs’, measured at inclusion, correlated with the patient’s clinical aging status, also assessed at T0.
The output of the univariable analysis (Table 2) reported a positive correlation of miR-191, miR-301a and miR-374a with LOFS score. In the subsequent multivariable model (Table 3), however, miR-320b and miR-374a emerged as independent predictors for LOFS, with higher LOFS scores being associated with higher miR374a and lower miR-320b levels. Note that different statistical techniques used for univariable and multivariable analysis may account for the apparent discrepancy in miRNAs arising as significant predictors.
Conversely, none of the miRNAs showed an association with the Balducci frailty score as ordinal outcome: the preliminary observed difference of miR-301a levels between the 3 categories fit, vulnerable and frail (p = 0.037) (Table 2) was not confirmed in the multivariable analysis (Table 3).
Although none of the miRNAs showed a significant individual association with G8 total score (Table 2), 4 miRNAs (miR-106b, miR-191, miR-320b, and miR-374a) resulted as independent predictors of total G8 from the backward multivariable model selection procedure (Table 3). This may point out toward a coordinated set of molecular changes reflecting the parameters assessed by this score. Concerning fTRST, which returns a higher score (scale 0–6) with increasing frailty, a negative correlation was observed for miR-301a, which was confirmed by the multivariable model selection (Tables 2 and 3).
Association of miRNAs at inclusion with other aging biomarkers
We also examined possible correlations between the ‘aging miRNAs’ and other aging biomarkers measured at T0.
Mean leukocyte telomere length (T/S ratio) was negatively correlated with miR-320b (Table 4), which was confirmed to be the only explanatory variable in the multivariable model selection (Table 3). Two other miRNAs, miR-34a and miR-106b, were borderline significant in the univariable analysis (Table 4), but were not retained in the multivariable model selection procedure (Table 3).
In univariable analysis, 5 miRNAs were found to be negatively correlated with IL-6, a cytokine well known to be increased during aging, particularly in frail individuals. Those were miR-30b, miR-106b, miR-191, miR301a and miR-374a. Conversely, miR-378a-3p showed a positive correlation with IL-6 (Table 4). These correlation trends are in line with our previous findings that miR-30b, miR-106b, miR-191 and miR-374 are all decreased, while miR-378a is increased, in elderly versus young patients, and with the widely documented age-related increase in plasma IL-6 levels [15, 16]. Of these 6 miRNAs that were significantly associated with IL-6, 3 resulted as significant independent predictors of IL-6 in the multivariable model: miR-106b, miR-374a and miR-378a (Table 3). Furthermore, TNF-α and MCP-1, which are also known to gradually increase in plasma during aging, showed pronounced associations with several ‘aging miRNAs’. As expected, negative associations were found for miRNAs showing decreased expression with higher age (miR-106b, miR-191, miR-301a, miR-374a) whereas positive associations were found for miRNAs showing increased expression with higher age (miR-378a, miR-320b) (Table 3). Higher miR-320b values were confirmed to be associated with higher TNF-α levels (p = 0.0140) in the subsequent multivariable selection model (Table 3). For MCP-1, the multivariable selection model further corroborated miR-301a (p = 0.0121) and miR-378a (p = 0.0025) as independent explanatory variables (Table 2). For the other plasma biomarkers (i.e. IL-10, RANTES, IGF-1), no consistent correlations were established in univariable and/or multivariable analyses (Tables 4 and 3).
Association of miRNAs with adverse effects of chemotherapy: decline of functionality and QoL, unexpected hospitalization and toxicity
In CTG, none of the individual miRNAs measured at inclusion (T0) was predictive of decline in functionality or decline in QoL at 3 months (T1) or at 1 year (T2): initial miRNA levels at T0 did not significantly differ between patients who experienced a decline in functionality and/or QoL during the course of the study and patients who did not (all p ≥ 0.05). Also, in patients exhibiting a decline in functionality and/or QoL at 3 months or at 1 year, miRNA plasma levels at the corresponding time point were not significantly altered (all p ≥ 0.05). Moreover, miRNAs plasma level at inclusion neither predicted grade II-III-IV toxicity at 3 months, nor unexpected hospitalization during the whole time course (all p > 0.1) (results not shown).
Correlations and cluster analysis of the 9 miRNAs
We have also examined the interrelationship between the 9 miRNAs included in the study; Table 5 summarizes Spearman’s correlation coefficients and associated p-values, based on miRNA measurements at inclusion. As expected, strong correlations exist between several miRNAs of the ‘aging miRNA’ panel, most particularly miR-30b, miR-191, miR-301a and miR-374a. Accordingly, a disjoint cluster analysis based on T0 miRNA measurements, revealed two main patient clusters of which one (cluster A) consistently scores lower on miR-20a, miR-30b, miR-191, miR-301a and miR-374a and higher on miR-378a compared to the other (cluster B), which is consistent with an ‘older’ aging profile for cluster A (Fig. 2). For the remaining miRNAs (miR-34a, miR106b and miR-320b), differences between patient clusters were either small or inconsistent, as shown in Fig. 2. One patient did not fit either of both clusters and was excluded from the cluster analysis.
In a next step, we compared both patient groups to determine whether they also showed differences with respect to aging biomarkers and/or clinical variables at inclusion. Interestingly, patients from cluster A indeed exhibited significantly higher fTRST, IL-6, TNFα and MCP-1 (Table 6), along with a higher mean age. LOFS was also apparently decreased in these patients (mean LOFS were 6.9 and 7.7 for clusters A and B, respectively), but this difference was not statistically significant (Table 6). Moreover, cluster A patients showed a markedly higher tendency to experience a decline in QoL during chemotherapy: 31.9% of cluster A patients, versus only 7.4% of cluster B patients, scored lower on QoL at 3 months (i.e. at the end of chemotherapy treatment) as compared to inclusion (p = 0.051).
We have recently published a scientific article reporting on the evolution of clinical and biological aging markers in older BC patients receiving chemotherapy . This study demonstrated that adjuvant TC chemotherapy had basically no impact on aging and frailty during a one-year period; we only detected a modest and temporary alteration of clinical aging indicators, while established aging biomarkers such as IL-6 did not show significant fluctuations during the a one-year period.
Aiming to corroborate these findings with additional aging biomarkers, we decided to analyse, on the same patient cohort, a panel of 9 miRNAs that were clearly differentially expressed in elderly BC patients compared to their younger counterparts in a previous exploratory miRNA pilot study also carried out by our group . This panel included miR-34a, miR-320b and miR378a, which were previously found to be positively correlated with age, and miR-20a, miR-30b, miR-106b, miR-191, miR-301a and miR-374a, which previously showed negative associations with age.
We have found that several of these ‘aging miRNAs’ changed during the course of the study, but these changes seemed rather transient and, most importantly, no striking differences in miRNA evolution between CTG and CG were revealed. Although an apparent discrepancy between CTG and CG was noted for miR-106b, the relevance of this observation can be questioned since the statistical significance was only modestly persuasive and the direction of the change is not in line with logic expectancies: given the previously established negative correlation of miR-106b with age, one would rather expect a decrease of miR-106b in chemotherapy-treated patients. Thus, the present study basically supports our previous statement that chemotherapy treatment is not associated with marked acceleration of the aging process in elderly BC patients.
Besides this primary research question, we have also further explored the relationships existing between the 9 selected ‘aging miRNAs’ and other parameters reflecting the patient’s clinical and biological age. In line with our previous findings, the present results show a decrease of miR-301a with calendar age at inclusion. Four other miRNAs of the panel, namely miR-30b, miR-106b, miR-320b and miR-374a-5p, showed trends consistent with our pilot data, but were not retained as independent explanatory variables for age in a multivariable model, probably because of the small age range of this study cohort (all ≥70).
Of relevance, miR-320b and miR-374a, which we previously found to be associated only with chronological age, now showed a clear correlation with frailty (LOFS) as well. Frail patients, scoring lower on the LOFS scale, exhibited higher miR-320b and lower miR-374a levels, as could be anticipated from the respectively positive and negative correlations of these miRNAs with age. Thus, these two miRNAs that evolve oppositely with aging are actually also oppositely correlated with LOFS. It should be noted that the basic Balducci frailty score, unlike the more refined LOFS, did not show significant associations with any of the 9 miRNAs in multivariable analysis. For the G8 screening tool, scored at inclusion and yielding a lower score with increasing frailty, 4 miRNAs of the aging panel (i.e. miR-106b, miR-191, miR-320b and miR-374a) were identified as independent predictors in the multivariable model selection. This finding should, however, be interpreted with caution, since 2 of these miRNAs (miR-106a and miR-191) were correlated in the opposite direction as expected and none of the four miRNAs showed a significant correlation with G8 score in univariable analysis. For the fTRST, yielding a higher score with increasing frailty, we found a consistent negative association with miR-301a, which is in full agreement with the age-related decrease of this miRNA . Thus, nearly all associations found are consistent with the known age-related trends. Overall, these data are suggestive of a link between the aging miRNAs and clinical frailty, whereas our previous miRNA study only uncovered strong associations of the 9 miRNAs with chronological age but not with frailty . The larger number of patients included in the present study probably afforded a greater statistical power, allowing detection of biological/clinical correlations that remained unrevealed in our initial miRNA pilot study.
With regard to the correlations between the aging miRNA panel and other biomarkers of aging, 6 of the selected miRNAs (though 3 of them failed to reach statistical significance in the multivariable analysis), correlated with circulating levels of IL-6, a well-known aging biomarker  that we previously also found to be predictive of frailty . Moreover, lower mean leukocyte telomere length, a well-described hallmark of biological aging, and increased plasma TNFα were both associated with higher miR-320b levels. This is once more in accordance with the age-related increase in miR-320b, as previously reported . Finally, the plasma inflammatory chemokine MCP-1, previously reported to increase with aging , was positively correlated with miR-378a, which also increases with age, and negatively correlated with miR-301a, which decreases with age.
When studying interconnections between the 9 miRNAs, we found that patients actually clustered into two major distinct subgroups with different ‘aging miRNA’ expression profiles: patients from one cluster showed significantly lower expression of miR-30b, miR-191, miR-301a and miR-374a and higher expression of miR-378a, as compared to the other cluster, which is consistent with an ‘older’ miRNA profile. In perfect agreement, these patients exhibited increased IL-6, TNFα, MCP-1 and also significantly higher fTRST and a markedly higher propensity towards decline in functionality during chemotherapy treatment. Apparently, aging miRNA-based patient clustering actually goes along with signs of different biological and clinical age between clusters.
Taken together, our findings point to marked correlations of (1°) chronological age with miR-301a, (2°) biomarkers of aging with miR-106b, miR-301a, miR-320b, miR-374a and miR-378a, and (3°) clinical frailty with miR-320b and miR-374a. Experiments on mice have previously suggested that miR-320b may be involved in cardiac degeneration due to hyperglycemia, as high levels of this miRNA are found in the diabetic heart . For miR-374a, however, this is the first study, to the best of our knowledge, that reveals an association with aging/frailty. We understand that our finding is limited by its correlative nature, so that further research is needed to assess whether a causative relationship links these particular miRNAs and physical involution.
Some miRNAs previously found over- or under-expressed with age, like miR-301a, or also with aging-related health conditions (such as miR-34a and miR-20a) [25, 26, 48] confirmed a trend in the same direction in our present study, albeit without reaching a formal statistical significance. One possible explanation could be that there is no actual association between the examined variables. However, since we are studying aging miRNAs in the context of elderly BC patients, we also have to bear in mind the multifaceted nature of the neoplastic process: indeed several miRNAs overexpressed during aging are underexpressed in patients affected by neoplasms or vice versa. In fact, an association with cancer opposite to that observed during aging has previously been demonstrated for miR-34a, miR-106b, miR-191-5p, miR-301a, although some studies generated conflicting results [29, 32, 34, 37,38,39,40,41, 49]. In particular, experiments on murine cell cultures revealed that miR-34a, showing age-related increase, acts as an oncosuppressor molecule within the p-53 regulatory pathway . Evidence of reduced levels of this miRNA in hepatocellular carcinoma samples suggests that these functions may be exerted in humans as well [35, 50]. On the other hand, miR-106b and miR-20b, which both decrease at higher age, belong to the so-called onco-miR cluster 17~ 92 and are known to promote cancer growth by different mechanisms, including the escape of senescence [28, 33, 36, 38, 40]. In particular, miR-106b, which promotes cell cycle progression by targeting the regulatory protein p21, is upregulated both in cancer tissue and plasma of BC patients, and may be responsible of chemoresistance in BC cell cultures .
Moreover, several mechanisms and molecular pathways link cancer and aging in a not straightforwardly understandable way. Aging is a multifactorial process, in which a major role is played by senescent cells. Cellular senescence is a protective anti-cancer mechanism at young age, arresting the cell cycle in case of acquired irreversible DNA damage. However, the drawbacks of an accumulation of senescent cells, most particularly the altered secretion of inflammatory mediators, are mirrored at higher age by nervous, cardiovascular, immune-endocrine degeneration and a higher susceptibility to cancer itself. Based on all these data, we cannot exclude a confounding role of the oncological pathology itself. It is thus important to study aging biomarkers specifically in cancer populations, and not only in non-cancer older populations.
The lack of convincing statistical significance for certain ‘aging miRNAs’ towards particular aging/frailty-related response variables may also be explained by the notion that several miRNAs may act in a weak but subtly coordinated fashion scarcely appreciable by conventional statistics. As a consequence, we were not surprised to observe general trends in the directions we had foreseen, yet without individual miRNAs standing out as ‘the ultimate’ aging/frailty-related miRNA(s). What we found is rather an intricate network of weakly acting miRNAs; in fact each of the 9 miRNAs from our ‘aging miRNA panel’ somehow showed significant, or at least borderline, association(s) with one or more of the biological and/or clinical aging/frailty markers we previously evaluated , indicating that our ‘aging miRNA’ panel indeed consists of miRNAs that are truly relevant in the context of the aging process.
Of course, as stated previously, this study has its limitations due to the intrinsic peculiarities and heterogeneity of this particular oncogeriatric population, and should be considered as a starting point towards a thorough understanding of the complex relationships between cancer, aging and chemotherapy.
In conclusion, the uppermost result of this study is the absence of any difference in the serum aging miRNA profile between chemotherapy-treated and control patients during and after adjuvant chemotherapy. This confirms our earlier findings with regard to the lack of a major impact of chemotherapy on a variety of biological and clinical aging markers.
Moreover, we created a miRNA-based clustering model resulting into two distinct groups exhibiting significantly different results for several biological/clinical aging parameters.
Although this is only a beginning, gaining knowledge about the interaction between the aging process and cancer treatment is a crucial step in the process of tailoring our treatments to older patients.
Monocyte Chemoattractant Protein-1
Docetaxel and cyclophosphamide
Tumor Necrosis Factor-alpha
Regulated on Activation, Normal T-Cell Expressed and Secreted
C-C motif chemokine-5
C-C motif chemokine-2
Insulin-like Growth Factor-1
Granulocyte-colony stimulating factor
National Comprehensive Cancer Network
Human epidermal growth factor receptor 2
Flemish version of the Triage Risk Screening Tool
Leuven Oncogeriatric Frailty Score
Eastern Cooperative Oncology Group – Performance Status
Quality of Life
Wildiers H, Kunkler I, Biganzoli L, et al. Management of breast cancer in elderly individuals: recommendations of the International Society of Geriatric Oncology. Lancet Oncol. 2007;8:1101–15.
Vandenberk B, Brouwers B, Hatse S, Wildiers H. p16INK4a: a central player in cellular senescence and a promising aging biomarker in elderly cancer patients. J Geriatr Oncol. 2011;2:259–69.
Pallis AG, Hatse S, Brouwers B, et al. Evaluating the physiological reserves of older patients with cancer: the value of potential biomarkers of aging? J Geriatr Oncol. 2014;5:204–18.
Decoster L, Van Puyvelde K, Mohile S, et al. Screening tools for multidimensional health problems warranting a geriatric assessment in older cancer patients: an update on SIOG recommendationsdagger. Ann Oncol. 2015;26:288–300.
Brouwers B, Dalmasso B, Hatse S, et al. Biological ageing and frailty markers in breast cancer patients. Aging. 2015;7:319–33.
Wildiers H, Heeren P. Puts M, et al. International Society of Geriatric Oncology consensus on geriatric assessment in older patients with cancer. J Clin Oncol. 2014;32:2595–603.
Balducci L, Extermann M. Management of the frail person with advanced cancer. Crit Rev Oncol Hematol. 2000;33:143–8.
Balducci L, Extermann M. Management of cancer in the older person: a practical approach. Oncologist. 2000;5:224–37.
Buttiglieri S, Ruella M, Risso A, et al. The aging effect of chemotherapy on cultured human mesenchymal stem cells. Exp Hematol. 2011;39:1171–81.
Franco S, Ozkaynak MF, Sandoval C, et al. Telomere dynamics in childhood leukemia and solid tumors: a follow-up study. Leukemia. 2003;17:401–10.
Maccormick RE. Possible acceleration of aging by adjuvant chemotherapy: a cause of early onset frailty? Med Hypotheses. 2006;67:212–5.
Unryn BM, Hao D, Gluck S, Riabowol KT. Acceleration of telomere loss by chemotherapy is greater in older patients with locally advanced head and neck cancer. Clin Cancer Res. 2006;12:6345–50.
Schroder CP, Wisman GB, de Jong S, et al. Telomere length in breast cancer patients before and after chemotherapy with or without stem cell transplantation. Br J Cancer. 2001;84:1348–53.
Lee JJ, Nam CE, Cho SH, Park KS, Chung IJ, Kim HJ. Telomere length shortening in non-Hodgkin's lymphoma patients undergoing chemotherapy. Ann Hematol. 2003;82:492–5.
Cohen HJ, Pieper CF, Harris T, Rao KM, Currie MS. The association of plasma IL-6 levels with functional disability in community-dwelling elderly. J Gerontol A Biol Sci Med Sci. 1997;52:M201–8.
Ferrucci L, Harris TB, Guralnik JM, et al. Serum IL-6 level and the development of disability in older persons. J Am Geriatr Soc. 1999;47:639–46.
Leng SX, Cappola AR, Andersen RE, et al. Serum levels of insulin-like growth factor-I (IGF-I) and dehydroepiandrosterone sulfate (DHEA-S), and their relationships with serum interleukin-6, in the geriatric syndrome of frailty. Aging Clin Exp Res. 2004;16:153–7.
Saurwein-Teissl M, Blasko I, Zisterer K, Neuman B, Lang B, Grubeck-Loebenstein B. An imbalance between pro- and anti-inflammatory cytokines, a characteristic feature of old age. Cytokine. 2000;12:1160–1.
Stumpf C, Lehner C, Yilmaz A, Daniel WG, Garlichs CD. Decrease of serum levels of the anti-inflammatory cytokine interleukin-10 in patients with advanced chronic heart failure. Clin Sci. 2003;105:45–50.
Antonelli A, Rotondi M, Fallahi P, et al. Increase of CXC chemokine CXCL10 and CC chemokine CCL2 serum levels in normal ageing. Cytokine. 2006;34:32–8.
Gerli R, Monti D, Bistoni O, et al. Chemokines, sTNF-Rs and sCD30 serum levels in healthy aged people and centenarians. Mech Ageing Dev. 2000;121:37–46.
Mansfield AS, Nevala WK, Dronca RS, Leontovich AA, Shuster L, Markovic SN. Normal ageing is associated with an increase in Th2 cells, MCP-1 (CCL1) and RANTES (CCL5), with differences in sCD40L and PDGF-AA between sexes. Clin Exp Immunol. 2012;170:186–93.
Cawthon RM, Smith KR, O'Brien E, Sivatchenko A, Kerber RA. Association between telomere length in blood and mortality in people aged 60 years or older. Lancet. 2003;361:393–5.
Brouwers B, Hatse S, Dal Lago L, Neven P, Vuylsteke P, Dalmasso B, et al. The impact of adjuvant chemotherapy in older breast cancer patients on clinical and biological aging parameters. Oncotarget. 2016;7:29977–88. https://doi.org/10.18632/oncotarget.8796.
Hatse S, Brouwers B, Dalmasso B, et al. Circulating MicroRNAs as easy-to-measure aging biomarkers in older breast cancer patients: correlation with chronological age but not with fitness/frailty status. PLoS One. 2014;9:e110644.
Khee SG, Yusof YA, Makpol S. Expression of senescence-associated microRNAs and target genes in cellular aging and modulation by tocotrienol-rich fraction. Oxidative Med Cell Longev. 2014;2014:725929.
Bartel DP. MicroRNAs: genomics, biogenesis, mechanism, and function. Cell. 2004;116:281–97.
Borgdorff V, Lleonart ME, Bishop CL, et al. Multiple microRNAs rescue from Ras-induced senescence by inhibiting p21(Waf1/Cip1). Oncogene 2010;29:2262–2271.
Cole KA, Attiyeh EF, Mosse YP, et al. A functional screen identifies miR-34a as a candidate neuroblastoma tumor suppressor gene. Mol Cancer Res. 2008;6:735–42.
Grassmann F, Schoenberger PG, Brandl C, et al. A circulating microrna profile is associated with late-stage neovascular age-related macular degeneration. PLoS One. 2014;9:e107461.
Hafez MM, Hassan ZK, Zekri AR, et al. MicroRNAs and metastasis-related gene expression in Egyptian breast cancer patients. Asian Pac J Cancer Prev. 2012;13:591–8.
Ichikawa T, Sato F, Terasawa K, et al. Trastuzumab produces therapeutic actions by upregulating miR-26a and miR-30b in breast cancer cells. PLoS One. 2012;7:e31422.
Ivanovska I, Ball AS, Diaz RL, et al. MicroRNAs in the miR-106b family regulate p21/CDKN1A and promote cell cycle progression. Mol Cell Biol. 2008;28:2167–74.
Li JY, Zhang Y, Zhang WH, Jia S, Kang Y, Tian R. Effects of differential distribution of microvessel density, possibly regulated by miR-374a, on breast cancer prognosis. Asian Pac J Cancer Prev. 2013;14:1715–20.
Li N, Fu H, Tie Y, et al. miR-34a inhibits migration and invasion by down-regulation of c-met expression in human hepatocellular carcinoma cells. Cancer Lett. 2009;275:44–53.
Liu C, Kelnar K, Vlassov AV, Brown D, Wang J, Tang DG. Distinct microRNA expression profiles in prostate cancer stem/progenitor cells and tumor-suppressive functions of let-7. Cancer Res. 2012;72:3393–404.
McDermott AM, Miller N, Wall D, et al. Identification and validation of oncologic miRNA biomarkers for luminal A-like breast cancer. PLoS One. 2014;9:e87032.
Waters PS, Dwyer RM, Brougham C, et al. Impact of tumour epithelial subtype on circulating microRNAs in breast cancer patients. PLoS One. 2014;9:e90605.
Yu H, Li H, Qian H, et al. Upregulation of miR-301a correlates with poor prognosis in triple-negative breast cancer. Med Oncol. 2014;31:283.
Zheng R, Pan L, Gao J, et al. Prognostic value of miR-106b expression in breast cancer patients. J Surg Res. 2015;195:158–65.
He L, He X, Lim LP, et al. A microRNA component of the p53 tumour suppressor network. Nature. 2007;447:1130–4.
Kumamoto K, Spillare EA, Fujita K, et al. Nutlin-3a activates p53 to both down-regulate inhibitor of growth 2 and up-regulate mir-34a, mir-34b, and mir-34c expression, and induce senescence. Cancer Res. 2008;68:3193–203.
Tabuchi T, Satoh M, Itoh T, Nakamura M. MicroRNA-34a regulates the longevity-associated protein SIRT1 in coronary artery disease: effect of statins on SIRT1 and microRNA-34a expression. Clin Sci. 2012;123:161–71.
Jones SE, Collea R, Paul D, et al. Adjuvant docetaxel and cyclophosphamide plus trastuzumab in patients with HER2-amplified early stage breast cancer: a single-group, open-label, phase 2 study. Lancet Oncol. 2013;14:1121–8.
Crawford J, Armitage J, Balducci L, et al. Myeloid growth factors. J Natl Compr Cancer Netw. 2009;7:64–83.
Blondal T, Jensby Nielsen S, Baker A, et al. Assessing sample and miRNA profile quality in serum and plasma or other biofluids. Methods. 2013;59:S1–6.
Costantino S, Paneni F, Luscher TF, Cosentino F. MicroRNA profiling unveils hyperglycaemic memory in the diabetic heart. Eur Heart J. 2016;37:572–6.
ElSharawy A, Keller A, Flachsbart F, et al. Genome-wide miRNA signatures of human longevity. Aging Cell. 2012;11:607–16.
Freres P, Josse C, Bovy N, et al. Neoadjuvant chemotherapy in breast Cancer patients induces miR-34a and miR-122 expression. J Cell Physiol. 2015;230:473–81.
Xu X, Chen W, Miao R, et al. miR-34a induces cellular senescence via modulation of telomerase activity in human hepatocellular carcinoma by targeting FoxM1/c-Myc pathway. Oncotarget. 2015;6:3988–4004.
Hans Wildiers is a recipient of the ‘Fonds voor Wetenschappelijk Onderzoek – Vlaanderen (FWO).’ The study was supported by a grant from the Belgian ‘Stichting tegen Kanker’.
FWO and Stichting tegen Kanker had no role in the design of the study, as well as in the collection, analysis, interpretation of data and in writing the manuscript.
Availability of data and materials
The datasets generated and analysed during the current study are not publicly available because there is no formal, broadly recognized platform for these data, but can be provided the corresponding author on reasonable request.
Ethics approval and consent to participate
The study was approved by the local ethics committees of Clinique Sainte Elisabeth (Namur), Imelda Hospital (Bonheiden), Jessa Hospital (Hasselt), Zol Ziekenhuis Oost-Limburg (Genk), Jules Bordet Institute (Brussels). Final approval was given by the University Hospitals Leuven central ethics committee on 20 April 2009. All patients signed a written informed consent in accordance to the Helsinki Declaration.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Dalmasso, B., Hatse, S., Brouwers, B. et al. Age-related microRNAs in older breast cancer patients: biomarker potential and evolution during adjuvant chemotherapy. BMC Cancer 18, 1014 (2018). https://doi.org/10.1186/s12885-018-4920-6
- Breast cancer
- Adjuvant chemotherapy