Abstract
Background: As microRNAs (miRNA) may play important roles in tumorigenesis by regulating the expression of proto-oncogenes or tumor suppressor genes, the present study analyzed single nucleotide polymorphisms (SNPs) located in miRNA and miRNA-binding sites of various genes and their impact on prognosis for 452 patients with early breast cancer. Materials and Methods: Three SNPs of miR-196a (rs3746444, rs11614913, and rs1044129) were selected using in silico analysis and genotyped using the Sequenom MassARRAY. Results: The median age of patients was 48 years, and 283 (62.6%) were estrogen and/or progesterone receptor (ER/PgR)-positive, 86 (19.0 %) had human epidermal growth factor receptor 2 (HER2)-overexpressing, and 77 (17.0%) had triple-negative early breast cancer. During the median follow-up of 6.9 years, 67 (14.8%) relapses and 55 (12.2%) deaths were recorded. Among the three polymorphisms, the C allele of miR-196a rs11614913T>C was significantly associated with worse disease-free (DFS) and distant DFS (DDFS) when adjusted for clinical and pathological parameters. In particular, the prognostic impact of rs11614913 was limited to the hormone receptor-expressing subtype, where the patients bearing the CC genotype showed worse survival in terms of DFS and DDFS compared with the patients with the TT or TC genotype as a recessive model (hazard ratio=2.610, p=0.003 for DFS; hazard ratio=2.730, p=0.013 for DDFS). Conclusion: The current study provides evidence that the miR-196a rs11614913T>C polymorphisms are possible prognostic biomarker for patients with hormone receptor-expressing early breast cancer.
Breast cancer is a common malignant tumor, affecting women at an increasing rate of incidence in many countries and is mostly diagnosed at an early stage via the widespread use of screening. Although several prognostic criteria have already been introduced to assist management after curative surgery for early breast cancer (EBC), the need for molecular markers has always been strongly suggested in order to discriminate individual variability and thus predict relapse or survival in patients with similar clinical status, especially when considering that adjuvant regimens containing more toxic chemotherapeutic agents, such as anthracylines, are acknowledged for their efficacy over survival in patients with EBC (1).
MicroRNAs (miRNAs) are a class of small, endogenous, non-coding RNAs able to regulate gene expression by translational repression or mRNA degradation of the target, thereby affecting critical functions in various physiological processes, ranging from cell proliferation to apoptosis (2, 3). Interestingly, most miRNAs are down-regulated in cancer cells and miRNA repression enhances cellular transformation and tumorigenesis, meaning that altering the function of miRNAs can affect tumor development and progression (4). Moreover, recent studies have demonstrated a relationship between the aberrant expression of miRNAs and cancer susceptibility, prognosis, and responsiveness to treatment (5-7). Most miRNAs bind to target sequences located within the 3’-untranslated region of mRNAs, resulting in the cleavage of the target mRNAs or repression of their translation (8). Thus, polymorphisms of miRNAs or within the miRNA-binding sites of the target genes that are implicated in cancer or that function as tumor suppressors or oncogenes could contribute to carcinogenesis or progression by altering the miRNA–mRNA interaction and thereby affecting the expression of the miRNA targets, as shown in previous studies (9-11).
Polymorphisms in miRNA genes and their target sites are rapidly being identified and investigated in human cancer as a novel class of variation. However, most studies have focused on how their de-regulation can lead to cancer in terms of risk of cancer development. Accordingly, this study investigated whether three target variants of miRNAs or miRNA binding sites, selected using web-based data, were associated with the prognosis for Korean patients with EBC who underwent curative surgery.
Materials and Methods
Patients' characteristics. Four hundred and fifty-two female patients who underwent surgery for EBC at Kyungpook National University Hospital (KNUH) between June 2000 and August 2008 were enrolled for evaluation in the current study. Patients with ductal carcinoma in situ, lobular carcinoma of the breast, or who underwent any type of neoadjuvant treatment prior to surgery were excluded. The patient data was obtained from the KNUH breast cancer registry and patient files. The tumors were classified and staged according to the WHO classification (12) and TNM staging system (13). The present study was approved by the local Ethics Committee of KNUH (No. 08-1008).
Selection of polymorphisms. This study selected three previously identified single nucleotide polymorphisms of miR-196a (SNPs; rs3746444, rs11614913, and rs1044129) that exhibited a potential association with breast cancer in previous studies and passed the selection criteria of a minor allelic frequency of more than 0.01 based on SNP databases (HapMap data, http://hapmap.ncbi.nlm.nih.gov).
Genotyping of polymorphisms. The genomic DNA of fresh-frozen breast tissue taken at the time of surgery was extracted using a Wizard genomic DNA purification kit (Promega, Madison, WI, USA). The three selected SNPs including two miRNA variants (rs3746444 and rs11614913) and one miRNA-binding site variant (rs1044129) were determined using the Sequenom MassARRAY (Sequenom Inc., San Diego, CA, USA) as described in detail in our previous publication (14) and the genotyping analysis was performed blindly as regards the subjects. The selected PCR-amplified DNA samples (n=2, for each genotype) were also examined using DNA sequencing to confirm the genotyping results.
Statistics. The SNP genotype was analyzed as a three-group categorial variable (referent model), and also grouped according to a dominant and recessive model. The Hardy–Weinberg equilibrium for each polymorphism was analyzed using a χ2-square test. For survival analysis, relapse was confirmed by biopsy, when possible, and categorized as local, regional, or distant; however, contralateral breast cancer during the follow-up period was not considered a relapse in this study. Disease-free survival (DFS) and overall survival (OS) was defined as the time from the date of surgery to the date of relapse or death from any cause or the date of the last follow-up, respectively. The cumulative incidence of relapse was defined as the time from the date of surgery to the date of the first event, where the curves were constructed based on the Kaplan–Meier method and analyzed using a log-rank test according to possible clinical (age, menopausal status, and use of adjuvant therapies) and histopathological risk factors [tumor size, number of involved lymph nodes, histological grade, and immunohistochemical expression of estrogen receptor (ER), progesterone receptor (PgR), and human epidermal growth factor receptor 2 (HER2)], and the genotype of the variants. The hazard ratio (HR) and 95% confidence interval (95% CI) for the genotypes of selected variants were calculated from a Cox regression analysis adjusted to age, stage, histological grade, and ER/PgR and HER2 status. In the multivariate analysis, the possible clinical and pathologic risk factors and SNPs significantly associated with survival in the adjusted univariate analysis were then analyzed as prognostic factors of relapse or survival for operated invasive ductal breast cancer. The differences in the continuous variables were compared using Student's t-test or an ANOVA test, while a χ2-test was used for the categorical variables. The statistical analyses were all performed using SPSS 15.0 (SPSS Inc., Chicago, IL, USA).
Results
Patients' characteristics and genotype frequency. For the 452 patients, median age was 48 years (range, 19 to 79 years) at the time of diagnosis, where 35.8% and 74.8% were node-positive and ER/PgR-expressing, respectively. The other basic clinical and pathological characteristics of the patients are listed in Table I. After curative surgery, 89.4% and 36.1% received adjuvant chemotherapy and radiotherapy, respectively, plus adjuvant hormonal treatment with tamoxifen or an aromatase inhibitor was also given if required, except for one patient. The median follow-up for the patients alive at the last follow-up was 7.3 years, at which point 67 (14.8%) had experienced relapse, including 17 locoregional and 50 distant relapses, where 12 distant relapses were identified after a locoregional relapse. In addition, 55 (12.2%) patients had died from breast cancer. The frequencies of each genotype are shown in Table II, and conformed to the Hardy–Weinberg equilibrium (p>0.05).
Survival analysis. At a median follow-up duration of 7.3 years, the estimated 5-and 10-year DFS, distant DFS, and OS were 87.8%, 91.0%, and 94.2%, and 81.0%, 87.0%, and 82.1%, respectively. Among the three target variants, miR-196a rs11614913 was significantly associated with DFS and distant DFS in the univariate analysis. Moreover, the multivariate analysis revealed that patients carrying the CC genotype for this polymorphism had a poor DFS and distant DFS when compared to patients with the T allele (DFS: HR=2.354, 95% CI=1.348-4.111, p=0.003; distant DFS: HR=2.399, 95% CI=1.236-4.656, p=0.010) regardless of patient and tumor characteristics (Figure 1 and Table III), yet there was no statistical association with OS (HR=0.942; p=0.865). In particular, the prognostic impact of miR-196a rs11614913 was only statistically associated with survival in patients with hormone-responsive breast cancers, based on the ER/PgR and HER2 immunohistochemistry (HR=2.610 and 2.730, p=0.003 and 0.013 for DFS and distant DFS, respectively; Figure 2 and Table IV). Moreover, no significant difference in the clinicopathological parameters was found according to the genotype or allele of miR-196a rs11614913.
Discussion
Although an association between the SNPs in protein-coding genes and the risk or prognosis of cancer has been introduced, only a few studies have been reported on the SNPs of miRNAs. Thus, based on several recent publications linking miRNA polymorphisms and susceptibility to breast cancer, the present study hypothesized that theses polymorphisms may be associated with the prognosis for early breast cancer. Among the three selected variants of miRNAs or miRNA-binding site previously identified as potential biomarkers, it was found that miR-196a variant rs11614913 may be a prognostic factor for estimating recurrence after curative surgery in patients with EBC.
miRNAs have already been shown to play an important role in a variety of carcinogenesis or tumor progression by altering the expression of oncogenes or tumor-suppressor genes. Previous studies have demonstrated an association between rs11614913 in miR-196a and cancer risk by alteration of the mature miR-196a expression and its binding to target mRNA (15-19). In particular, the C allele of rs11614913 was recently identified as a risk factor for breast cancer in Asians and Caucasians (16, 20), which is comparable with the present findings, although the impact of SNPs on cancer risk does not always correlate with their impact on prognosis. While the exact biological role of miR-196a and its variants has not yet been elucidated, Hoffman et al. demonstrated that the rs11614913 variant led to more efficient processing of the miRNA precursor to its mature form, this being more than two-times higher in breast cancer cells transfected with pre-miR-196a-C compared to cells transfected with pre-miR-196a-T (20). Furthermore, an Ingenuity Pathway Analysis network analysis revealed that several cancer-relevant transcripts known as prognostic indicators for poor clinical outcome in breast cancer including tumor protein-63 (TP63), inhibin beta-B (INHBB), the tumor suppressors growth arrest and DNA damage-inducible, gamma (GADD45G), and the S100 family of genes (A8 and A9) were significantly altered following the introduction of pre-miR-196a-C (20-23). Thus, when combined with epidemiological evidence, the rs11614913C allele appears to be associated with susceptibility to breast cancer in relation to increased mature miR-196a levels influencing genes with relevance for breast tumor development and progression.
Although there are no data regarding the prognostic role of this variant in breast cancer, consistent results associating the C allele of rs11614913 with an increased risk of cancer, including breast, gastric, head and neck, and non-small cell lung cancer, would seem to support this variant as a potential prognostic marker (16, 18, 20, 24, 25). Notably, miR-196a rs11614913 was only associated with survival in patients with the ER/PgR-positive and/or HER2-overexpressing breast cancer type. Therefore, rs11614913 could be a promising prognostic factor for patients with ER-positive EBC. Prognostic markers for EBC are crucial due to the ongoing controversies regarding adjuvant management. In particular, for patients with ER-positive EBC, the addition of chemotherapy remains contentious due to the use of highly toxic chemotherapeutic agents, such as anthracyclines. Therefore, the current findings are clearly important for this patient group with regard to determining the use of adjuvant chemotherapy after surgery. Meanwhile, several biological markers, such as Ki-67 (26, 27), and gene signatures (28-30) have recently been identified and introduced as important prognostic and predictive markers for hormone-responsive EBC, pending the validation results of ongoing clinical trials. However, the current study did not evaluate the association between rs11614913 and these recently discovered biological markers, thereby warranting further studies.
Notwithstanding, while SNPs are thought to be attractive biomarkers as they are stably-inherited, highly abundant and show diversity within and among populations, the application of individual SNPs is limited due to their penetrance and the difficulty involved in identifying their effects. Thus, until the present results are confirmed by replication studies with different populations, caution is warranted in terms of drawing definite conclusions from the current study. Furthermore, the current study found that rs11614913 CC was significantly associated with a poor DFS (HR=2.354; p=0.003) and distant DFS (HR=2.399; p=0.010), while no significant association was identified for OS (HR=0.942; p=0.865). This could be explained by the relatively short follow-up duration and thus small number of cases who died from breast cancer as disease in most patients was at an early stage (64.2% node-negative and 75.7% stage I or IIA). In addition, the diverse treatment modalities given after recurrence, including metastasectomy and active chemotherapy or hormone treatment, might have weakened the impact of this variant on OS, as previously described in various studies of recurrent or metastatic cancer.
In summary, rs11614913 is a possible prognostic factor for patients with ER/PgR-positive EBC after complete surgery, regardless of the clinical or pathological characteristics. However, since the exact mechanism and tumor specificity of the rs1044129 variant have not yet been defined and genetic polymorphisms often vary between different ethnic groups, the present findings need to be confirmed in further studies with other populations of patients with breast cancer to clarify the association between these polymorphisms and the prognosis of breast cancer. Moreover, the present findings warrant a further study to explore rs11614913 as a predictive marker of treatment, such as chemotherapy and hormonal therapy for EBC or as a novel therapeutic target for the management of breast cancer.
Acknowledgements
This work was supported by a Korea Science and Engineering Foundation (KOSEF) grant funded by the Korean Government (MEST) (No. 2009-0091573).
Footnotes
-
↵* These Authors equally contributed as first authors.
- Received February 14, 2014.
- Revision received April 5, 2014.
- Accepted April 7, 2014.
- Copyright© 2014 International Institute of Anticancer Research (Dr. John G. Delinassios), All rights reserved