Abstract
Background/Aim: Circulating microRNAs (miRs) in blood have been highlighted as diagnostic, prognostic and predictive biomarkers for “Precision medicine”. This study aimed to explore the possibility of using circulating precursor miRs (pre-miRs) as clinical biomarkers of recurrence in breast cancer. Materials and Methods: We performed miR microarray analyses of circulating miRs in blood from patients with or without recurrence of breast cancer and identified miR-488-5p as a recurrence-related miR. Then, the expression levels of pre-miR-488 or miR-488-5p were measured by RT-qPCR and the clinicopathological and prognostic significance of circulating pre-miR-488 was assessed in the blood of breast cancer patients. Results: A positive correlation was noted between pre-miR-488 and miR-488-5p expression in blood. In 330 cases of surgically-treated breast cancer, high expression of circulating pre-miR-488 was an independent poor prognostic factor for recurrence-free survival. Conclusion: Circulating pre-miR-488 expression could be a novel prognostic biomarker for predicting recurrence in breast cancer patients.
“Precision medicine” is defined as personalized treatment of patients with similar clinical symptoms on the basis of differences in genetics, biomarkers or phenotypes. It is expected to improve the prognosis of various malignancies, including breast cancer (1). There is growing interest in the use of “liquid biopsies” to identify biomarkers for “precision medicine” (2). “Liquid biopsy” is defined as the identification of biomarkers in body fluids, mainly in blood, and can provide non-invasive “real-time” information regarding tumor characteristics (2).
Recently, circulating microRNAs (miRs) in blood have been highlighted as diagnostic, prognostic and predictive biomarkers because of their high tissue specificity, stability and altered expression in tumor development (3). In particular, miRs in serum exosomes (actively secreted extracellular vesicles) are more stable than other forms, and they are expected to have potential applications as biomarkers (4, 5). Thus, miR analysis in liquid biopsies could offer a less invasive and more cost-effective biomarker for precision medicine (6).
miRs, short (17-22 bases) non-coding RNAs, regulate various cell functions, including proliferation, differentiation, cell fate determination and apoptosis. This is achieved by binding to the 3’-untranslated region of one or multiple target mRNAs followed by an inhibited expression (7). miRs play key roles in various aspects of tumor progression (8) as we reported elsewhere (5, 9, 10). The biogenesis of the majority of miRs involves transcription by RNA polymerase II into primary miRs and subsequent processing by the nuclear enzyme Drosha into shorter hairpin structures, called pre-miRs. The pre-miRs are then transported into the cytoplasm where they are further processed into mature miRs by the enzyme Dicer (11).
Many reports have shown the clinical significance of mature miR expression in various cancers (8). However, the significance of pre-miR expression remains unclear. Pre-miRs, like mature miRs, are enveloped by exosomes and secreted into body fluids where they function as intercellular messengers (12, 13). We focused on the possible advantages of assessing pre-miRs rather than mature miRs, as detection is less costly and easier to conduct.
In this study, miR-488-5p was first identified as a recurrence-related exosomal miR and the clinicopathological and prognostic significance of pre-miR-488 expression was assessed in the blood of breast cancer patients. Furthermore, the correlation of pre-miR-488 and miR-488-5p expression in the blood of such patients was examined.
Materials and Methods
Clinical samples. In an effort to identify recurrence-related exosomal miRs in serum, 10 serum samples were obtained from breast cancer patients who had undergone surgery for a primary tumor. Surgery was conducted at the Department of Surgery, Kyushu University Beppu Hospital (Beppu, Japan) in 2013. Among the 10 samples, 5 were from patients with recurrence, and 5 were from those without recurrence. Serum samples were stored at −80°C until miR extraction.
To determine the clinical significance of pre-miR-488 expression in the blood, we used updated clinical data files from breast cancer patients described elsewhere (14). Briefly, a total of 594 patients with breast cancer underwent resection of the primary tumor at the Department of Breast Oncology, National Kyushu Cancer Center (Fukuoka, Japan) from 2000 to 2008. Of these, 356 female patients with breast cancer without distant metastases, preoperative therapy or previous treatment for various cancers were included in this study. Among them, 330 were invasive ductal carcinomas (IDC), and 26 were ductal carcinomas in situ (DCIS). Survival analysis was performed for the 330 patients with IDC. The control samples were obtained from 11 healthy volunteers at Kyushu University Beppu Hospital. Aspiration of blood was conducted immediately before operation under general anesthesia. Blood was obtained through a venous catheter. The first ml of blood was discarded because of possible contamination by epidermal cells. Each 1 ml sample of blood was immediately mixed with 4 ml of ISOGEN-LS (Nippon Gene, Toyama, Japan) and stored at −80°C until RNA extraction. Polymorphprep™ (Alere Technologies AS, Oslo, Norway) was used to separate 3 blood cell fractions (polymorphonuclear leukocytes (PMN), mononuclear cells (MC) and erythrocytes (RBC)). Samples were obtained from 9 patients with IDC who underwent primary tumor resection at Kyushu University Beppu Hospital in 2017. Each fraction was mixed with ISOGEN II (Nippon Gene, Toyama, Japan) and stored at −80°C until RNA extraction.
Tumor tissues and paired normal tissues were obtained from 82 breast cancer patients who had undergone primary tumor resection at Kyushu University Beppu Hospital and affiliated hospitals between 1992 and 1996. The frozen tissue specimens were homogenized, mixed with ISOGEN (Nippon Gene) and stored at −80°C until RNA extraction. From this cohort, we used 30 high-quality RNA samples for RT-qPCR analysis.
All patients gave written informed consent to participate in this study, which was approved by each institutional review board and the Ethics and Indications Committee of each institution. The observation periods ranged from 0.3 to 6.9 years (median 3.8 years). Postoperative adjuvant therapy was performed according to the St.Gallen Consensus Conference (15). The patients underwent clinical examinations at least every 3 months and annual mammographies and were further tested only if they had symptoms. The stages and grades of the tumors were classified according to the AJCC/UICC TNM classification and stage groupings. All data for the samples, including age, pathological tumor size, nuclear grade, venous involvement, lymphatic involvement, lymph node metastasis, and the status of ER, PgR and HER2 were obtained from the medical records. The status of ER, PgR and HER2 was determined by immunohistochemical analysis. 0/1+ was defined as negative and 2+/3+ as positive in the HER2 status of this study. Overall survival (OS) was defined as the time between the date of diagnosis and the date of death. Recurrence-free survival (RFS) was defined as the length of time after surgical treatment for cancer during which the patient survived without a sign of recurrence.
Isolation of exosomes in serum. Exosomes were isolated from serum by ultracentrifugation as described previously (5).
Extraction of total RNA. Extraction of total RNA from serum exosomes was performed using the miReasy serum/plasma kit (Qiagen, Venlo, Netherlands) according to the manufacturer's protocol. Frozen tissue specimens were homogenized, and total RNA was extracted with ISOGEN according to the manufacturer's instructions. Total RNA of blood was extracted according to the ISOGEN-LS manufacturer's instructions.
PCR cloning. Human pre-miR-488 and RNA18S5 were amplified by PCR with the primer pairs as follows: pre-miR-488, forward 5’-CTCTCCCAGATAATGGCACTCTCA-3’ and reverse 5’-AGAGTCATCTGACCAAGAAATAGCC-3’; RNA18S5, forward 5’-AGTCCCTGCCCTTTGTACACA-3’ and reverse 5’-CGATCCGAGGGCCTCACTA-3’. The PCR products were cloned into a pCR2.1 cloning vector (Invitrogen, Carlsbad, CA) according to the manufacturer's instructions. We used each plasmid as a reference for calibration of the expression level of each target gene (pre-miR-488 or RNA18S5 at quantitative PCR (qPCR) as described below.
RT-qPCR Assay. RT-qPCR assessments of pre-miR-488, RNA18S5 and GAPDH were performed as previously described (5). In brief, RT was performed with random hexamers using M-MLV reverse transcriptase (Invitrogen, Carlsbad, CA, USA). qPCR was performed with LightCycler® FastStart DNA Master SYBR Green I (Roche, Basel, Switzerland). The raw data are presented as the relative quantity of target genes, normalized with respect to RNA18S5 or GAPDH. The primer sequences for RT-PCR were as follows: pre-miR-488, forward 5’-CTCTCCCAGATAATGGCACTCTCA-3’ and reverse 5’-AGAGTCATCTGACCAAGAAATAGCC-3’; RNA18S5, forward 5’-AGTCCCTGCCCTTTGTACACA-3’ and reverse 5’-CGATCCGAGGGCCTCACTA-3’ and GAPDH, forward, 5’-TTGGTATCGTGGAAGGACTC-3’ and reverse 5’-AGTAGAGGCAGGGATGATGT-3’.
RT-qPCR assessments of miR-488-5p and RNU6B were performed as described previously (5). In brief, miR-488-5p and RNU6B-specific cDNAs were synthesized using gene-specific primer sets according to the TaqMan® Micro-RNA Reverse Transcription Kit protocol (Applied Biosystems, Foster City, CA, USA). qPCR was performed using PCR LightCycler® 480 Probes Master (Applied Biosystems). miR-488-5p and RNU6B were purchased from Applied Biosystems.
Relative quantification of miR expression was calculated using the 2−ΔΔCt method. The raw data are presented as the relative quantity of each cloned plasmid or cDNA from Human Universal Reference Total RNA (Clontech Laboratories, Palo Alto, CA, USA) and normalized by internal control to RNA18S5, GAPDH or RNU6B. The PCR product was electrophoresed with 3% agarose gel and stained with ethidium bromide.
miR microarray analysis. We performed miR microarray as described previously (10). In brief, extracted total RNA was labeled with Hy5 using the miRCURY LNA microRNA Array labeling kit (Exiqon, Vedbaek, Denmark). Labeled RNAs were hybridized onto 3D-Gene Human miRNA Oligo chips 17v 1.00 (Toray, Kamakura, Japan). Fluorescent signals were scanned and analyzed using the 3D-Gene Scanner (Toray). The raw data from each spot were normalized by subtraction of the background signal mean intensity, determined by the 95% confidence intervals of the signal intensities of all blank spots. Valid measurements were considered as those in which the signal intensity of both duplicate spots was >2 s.d. of the background signal intensity. Cluster analysis was performed to identify exosomal miRs in serum that were statistically differentially expressed between breast cancer patients with and without recurrence (p<0.05, Student's t-test). The heatmap of the miRs profile is illustrated in Figure 1. The expression levels of these miRNAs range from low (green) to high (red). Patients 75, 83, 60, 95 and 94 had recurrence and patients 14, 7, 40, 35 and 25 had no recurrence.
Statistical analysis. For clinical analysis, cases were divided into two groups using the minimum p-value approach based on the miR-488 expression level, which is a comprehensive method to identify the optimal risk separation cutoff point in continuous gene expression measurements for survival analysis in multiple datasets (16). Associations between the variables were tested with the Mann–Whitney U-test, Student's t-test or Fisher's exact test. The degree of linearity was estimated by Spearman's rank correlation co-efficient. Survival curves were drawn according to the Kaplan–Meier method, and survival analysis was carried out by the log-rank test when two curves were being compared. Cox proportional hazards regression was used to determine univariate and multivariate hazard ratios for survival. A two-sided p<0.05 was deemed statistically significant. Statistical analyses were performed using JMP Pro 13 software (SAS Institute, Cary, NC, USA).
Results
miR profiling and identification of recurrence-related exosomal miRs in serum. It is well recognized that circulating miRs in blood are mainly carried by exosomes (8, 9). To identify recurrence-related exosomal miRs in the serum of breast cancer patients, miR microarray analyses of exosomal miRs was performed from 5 patients with recurrence and 5 patients without recurrence. The heatmap of the miR profile is illustrated in Figure 1. The expression levels of these miRNAs ranged from low (green) to high (red). Various miRs were significantly down-regulated or up-regulated in the cases with recurrence compared to those without recurrence. Among them, we focused on miR-488-5p as an up-regulated miR in the cases with recurrence because its clinical significance is unknown in breast cancer.
Pre-miR-488 expression in blood or breast cancer tissues of breast cancer patients. The findings shown in Figure 1 indicated that it would be useful to characterize the clinical significance of pre-miR-488 expression as well as its utility as a clinical biomarker in the blood of breast cancer patients. The sequence of pre-miR-488 is shown in Figure 2A. The underlined regions are paired primers for the PCR. The expression levels of circulating pre-miR-488 from 330 patients with IDC, 26 patients with DCIS and 11 healthy controls are plotted in Figure 2B. The pre-miR-488 levels of the control cases, DCIS cases, and IDC cases ranged from 0.061 to 0.496 (median, 0.191), from 0.012 to 0.328 (median, 0.078) and from 0.001 to 1.072 (median, 0.075), respectively. Interestingly, pre-miR-488 expression in the blood was significantly lower in patients with IDC as well as DCIS compared with healthy controls. There was no significant difference in pre-miR-488 expression levels in blood between of DCIS and IDC. Next, pre-miR-488 expression was assessed in IDC tissues (Figure 2C). Pre-miR-488 expression in tumor tissues was lower than normal tissues in 83.3% of 30 breast cancer patients although in some cases, pre-miR-488 expression in tumor tissues was much higher than in normal tissues. These results may suggest that down-regulation of pre-miR-488 reflects carcinogenesis rather than tumor progression of breast cancer.
Assessment of pre-miR-488 in blood cell fractions and plasma. The distribution of pre-miR-488 in blood was assessed by RT-qPCR. Whole blood is composed of polymorphonuclear leukocytes (PMN), mononuclear cells (MC), erythrocytes (RBC) and plasma. The expression level of pre-miR-488 in each fraction is shown in Figure 2D. The expression of pre-miR-488 in plasma was evaluated as follows: [putative plasma]=[whole blood]-[PMN]-[MC]-[RBC]. This approach was necessary because quantitative analysis of plasma RNA could not be performed because of the low levels present, although pre-miR-488 was detected in plasma by traditional RT-PCR (Figure 2E). The median relative ratios of RBC, MC, PMN and putative plasma to whole blood were 0.022, 0.000, 0.009 and 0.965, respectively. Pre-miR-488 expression was estimated to be much higher in plasma than in polymorphonuclear leukocytes, mononuclear cells, and erythrocytes. This finding indicates that pre-miR-488 is localized mainly in plasma in blood.
Correlation between amounts of pre-miR-488 and miR-488-5p in blood. We randomly chose 27 breast cancer cases from which we obtained high-quality RNA and examined the correlation between the amounts of pre-miR-488 and miR-488-5p in blood. Spearman's rank correlation co-efficient demonstrated a significantly positive correlation between pre-miR-488 expression and miR-488-5p levels in blood (r=0.457, p=0.017) (Figure 2F).
Prognostic significance of the amount of pre-miR-488 in blood. The prognostic significance of pre-miR-488 expression in blood was assessed using 330 IDC cases who underwent curative surgery for OS and RFS. Cases were divided into two groups as described in Materials and Methods. The cut-off value was 0.052. High pre-miR-488 expression was not significantly associated with poor OS (p=0.055), but with poor RFS (p=0.043) (Figure 3). Then, univariate and multivariate analyses for RFS were performed (Table I). Univariate analysis showed that lymphatic involvement, lymph node metastasis, PgR and high levels of pre-miR-488 were statistically significant prognostic factors for RFS (p=0.0005, 0.002, 0.028 and 0.032, respectively). Pre-miR-488 expression and clinicopathological factors such as tumor size, lymph node metastasis, PgR, and HER2 were included in the multivariate analysis. Notably, multivariate analysis demonstrated that pre-miR-488 expression in blood was an independent prognostic factor for RFS in patients with IDC (HR=3.19, p=0.030). Moreover, subgroup analyses for RFS were performed according to each subtype. As shown in Figure 3C, the high pre-miR-488 expression group exhibited poorer prognosis only in the HER2-positive group (p=0.039).
Clinicopathological significance of pre-miR-488 expression in blood. The relationship between clinicopathological factors and pre-miR-488 expression in blood was examined using the 330 patients with IDC. As shown in Table II, no significant difference was found between the high and low pre-miR-488 expression groups regarding clinicopathological factors such as age, tumor size, nuclear grade, venous involvement, lymphatic involvement, lymph node metastasis, the status of ER, PgR and HER2 and subtype.
Discussion
In this study of breast cancer, it was demonstrated that high expression of pre-miR-488 in blood was an independent biomarker that predicted recurrence of disease. Furthermore, pre-miR-488 expression was positively correlated with miR-488-5p levels in blood, which was shown to be a recurrence-related miR, by miR array analysis. To the best of our knowledge, this is the first study to explore pre-miR as a clinical biomarker in cancer patients and determine the clinical relevance of circulating pre-miR-488 in breast cancer.
Numerous reports have been published on the utility of miRs as diagnostic or prognostic biomarkers in malignancies. However, the expression analysis of mature miR requires high-quality RNA and specific primers for RT and PCR assessment. Moreover, an accurate internal control for mature miR is still controversial, although researchers have attempted to identify reliable and reproducible markers (6). These issues may lead to incorrect conclusions about the clinical significance of miR expression. On the other hand, expression analysis of pre-miR can be performed by conventional methods using universal materials. Focusing on pre-miR may allow laboratories to overcome the problems of miR expression analysis for use as biomarkers.
High expression of pre-miR-488 in blood was associated with a poor prognosis in breast cancer patients, especially in HER2-positive cases, although pre-miR-488 was down-regulated in the blood and tumor tissues of breast cancer patients compared with healthy controls. So far, miR-488 is reported to act as a tumor suppressor. It is down-regulated in a range of tumors, including pancreatic cancer (17), prostate cancer (18), hepatocellular carcinoma (19), ovarian cancer (20), gastric cancer (21), colorectal cancer (22), non-small cell lung cancer (23, 24), melanoma (25) and gastrointestinal stromal tumor (26) where it targets JAK1, AR, ADAM9, Six1, PAX6, PHF8, elF3a or PRKDC, respectively. Pre-miR-488 was also found to be down-regulated in breast tumor tissues compared with normal controls. This observation suggests that miR-488 may also act as a tumor suppressor in breast cancer tissue. Unexpectedly, a survival analysis using blood demonstrated that high expression of circulating pre-miR-488 was a poor prognostic biomarker in breast cancer.
Circulating pre-miR-488 was mainly located in the plasma. miRs are packaged inside exosomes (extracellular vesicles) that are secreted into the blood. Exosomal miRs are probably delivered to recipient cells where they regulate target gene expression. Pre-miRs are also enveloped by exosomes in body fluids and function as intercellular messengers (12, 13). We observed that the high expression of pre-miR-488 in blood was associated with a poor prognosis in breast cancer, especially in HER2-positive cases, although pre-miR-488 was down-regulated in the blood of breast cancer patients. These results imply that circulating pre-miR-488 plays an oncogenic role as an intercellular messenger, as does miR-203 in CRC as we reported previously (5). Circulating pre-miR-488 might affect the HER2 signal pathway because its high expression was associated with poor prognosis only in HER2-positive patients.
Collectively, miR-488 may have the potential to work as a tumor suppressor in tumor tissues, whereas as an oncogene in blood. Recent studies have focused on the contrasting roles of miRs in oncogenesis and tumor suppression (27). We speculate that pre-miR-488 may have dual roles in breast cancer development. Namely, it may function as a tumor suppressor during tumorigenesis and an enhancer during tumor progression as is the case with DNA damage response genes (28, 29). Furthermore, miR may have site-dependent functions that promote tumor development (5). miR localization should be taken into account when determining its clinical significance. Further examination of miR-488 to define its function and significance as a biomarker in various sites in breast cancer patients will be required for the proper clinical interpretation.
In summary, pre-miR-488 expression in blood could be a novel prognostic biomarker for predicting recurrence in breast cancer patients.
Acknowledgments
The Authors thank Dr. Tyler Lahusen for helpful comments and English proofreading. The Authors also thank K. Oda, M. Kasagi, S. Sakuma, N. Mishima and T. Kawano for their excellent technical assistance. This work was supported in part by the following grants and foundations: Japan Society for the Promotion of Science (JSPS) Grant-in-Aid for Science Research (Grant Numbers JP16K07177, JP16K10543, JP16K19197, JP17K16454, JP17K16521, JP17K10593 and JP17K19608); OITA Cancer Research Foundation; Daiwa Securities Health Foundation; Grant-in-Aid for Scientific Research on Innovative Areas (15H0912); Priority Issue on Post-K computer (hp170227) (hp170227, hp160219); JSPS KAKENHI (15H05707); Eli Lilly Japan K.K. Grant; Japanese Foundation for Multidisciplinary Treatment of Cancer.
Footnotes
Conflicts of Interest
The Authors have no conflicts of interest to disclose.
- Received June 6, 2018.
- Revision received June 14, 2018.
- Accepted June 15, 2018.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved