Abstract

Aim. Biochemical markers, including microRNAs (miRs), may facilitate the diagnosis and prognosis of breast cancer. This study was aimed at assessing serum miR-155 expression in patients with breast cancer and receptors. Methods. This case-control study was conducted on 36 patients with breast cancer and 36 healthy individuals. After RNA extraction from the patient’s serum, cDNA was synthesized. The expression of miR-155 was measured using RT-qPCR. Demographic and histochemical data were extracted from patient documents. Data were analyzed using the Statistical Package for the Social Sciences (SPSS) software. Results. The mean age of subjects in breast cancer and control groups was and years, respectively. The serum miR-155 expression was higher in the cancer group () compared to the control group (). There was a significant relationship between serum miR-155 expression and the tumor grade (), tumor stage (), and tumor size () of the patients. However, no relationship between miR-155 expression and the presence of lymph node involvement (), HER2 (), Ki-67 (), progesterone receptor (), and estrogen receptors () was found. The ROC curve analysis showed that the AUC was 0.89 (77.78% sensitivity and 88.89% specificity), and the cutoff was 1.4 (Youden index: 0.6667) for detecting breast cancer. Conclusion. The findings of this study revealed that serum miR-155 may serve as a potential noninvasive molecular biomarker for breast cancer diagnosis and can help predict the grade of the disease.

1. Introduction

Breast cancer (BC) is the most prevalent cancer worldwide and also accounts for 22.8% of female cancers [1]. Breast cancer mortality was estimated to be 626679 in 2018 [2]. In Iran, breast cancer accounts for 76% of female cancers, with 8500 new cases each year [3]. The most important risk factors for breast cancer include female gender, age (30 years old and older) [4], positive family history for breast cancer [4], and familial genetic mutations, including mutations in the breast cancer A1 (BRCA1) and BRCA2 genes [5]. Furthermore, women with a history of breast cancer are more likely (20-25%) to develop microscopic cancer in the opposite breast [6]. A positive history for cancer in the endometrium, ovaries, or colon, as well as radiation therapy for Hodgkin’s lymphoma, was shown to increase the risk of breast cancer [6]. The gold standard for diagnosis of breast cancer is histopathology [6]. Several tumor markers have been suggested for the evaluation and management of breast cancer including estrogen and progesterone receptors (ER/PR), which are used for the assessment of susceptibility to hormone treatment, and human epidermal growth factor receptor 2 (HER2), which is used to assess the susceptibility to trastuzumab treatment [7].

Microribonucleic Acids (microRNAs) are a large subgroup of noncoding RNAs made up of 18-25 nucleotides [8]. MicroRNAs (miRs) regulate gene expression after transcription. The increased expression of some miRs, including miR-194 and miR-425, was shown in invading breast cancer cells [9]. One of the goals of this study is to investigate the role of miR-155 in women with breast cancer, but its primary goal is to find a biomarker for breast cancer diagnosis based on hormonal receptors. For the first time in this study, we investigated the role of miR-155 in contraceptive drugs and the number of pregnancies. We found a significant difference between the menarche age with the Ki-67 receptor and the tumor stage with contraceptive medication. It is also the first time to provide a diagnostic value of the tumor grade based on the receiver operating characteristic (ROC) curve as well as the Youden index in breast cancer patients.

2. Material and Methods

This case-control study was performed on 36 women with breast cancer (BC) who were referred to our Radiotherapy and Oncology Center from March 2017 to March 2018. The Medical Ethics code of the approved study protocol is IR.IAU.MSHD.REC.1396.83. Each subject, regardless of the allocated group, signed a written informed consent form before participation.

2.1. Study Population

All patients with the documented diagnosis of breast cancer based on physical examination and imaging and laboratory assessments were included in the case group. Inclusion criteria for subjects in the case group included a documented diagnosis of breast cancer based on histopathology, age between 20 and 60 years old, and the possibility of obtaining blood samples from the subject. Any subject in the case group whose documentation for histopathological diagnosis of breast cancer was absent was excluded from the study. The control group subjects were randomly selected from healthy women who visited the center for a checkup. The inclusion criteria for control subjects were the absence of documented cancer and age between 20 and 60 years old. Exclusion criteria for the control group included a history of polycystic ovary syndrome and a history of any cancer in first-degree relatives.

The clinical data included grading of BC assessed by the Nottingham Grading System: well-differentiated (WD) tumor: grade I, moderately differentiated (MD) tumor: grade II, and poorly differentiated (PD) tumor: grade III. The well-differentiated tumor represented high homology to the normal terminal duct lobular unit, tubule formation (>75%), mild degree of nuclear pleomorphism, and low mitotic count. A moderately differentiated tumor (grade II) is characterized by tubule formation between 10% and 75%. A poorly differentiated tumor is characterized by a marked degree of cellular pleomorphism and frequent mitoses and no tubule formation (<10%) [10]. The type of involvement and stage of cancer, as well as the presence of HER2, PR, and ER, are detected in breast cancer cell biopsy.

2.2. Serum Preparation and RNA Extraction

A 5 ml blood sample was incubated at room temperature for 30 minutes and centrifugated at 4000 rpm for 10 minutes. In order to assess total RNA, 200 μl of the serum was extracted. Total RNA was extracted using the Norgen Biotek Plasma/Serum RNA Purification Mini RNA Kit (Ontario, Canada) Cat: 55000 according to the manufacturer’s instruction with modification.

2.3. cDNA Synthesis and qRT-PCR

Reverse transcription was performed using the BON-miR miRNA 1st-Strand cDNA Synthesis Kit (Bonyakhteh, Iran) cat: BON209001 based on the manufacturer’s instructions. The cDNA was synthesized using the thermocycler device for 10 minutes at 25°C, 60 minutes at 42°C, and 10 minutes at 70°C. The qRT-PCR was performed using the BON-miR QPCR (Bonyakhteh, Iran) cat: BON209002 kit. Primarily, 0.5 μl of forward primer (miR-155 or SNORD47), 0.5 μl of universal reverse primer, 6.5 μl of SYBR master mix, and 1 μl of the synthesized cDNA were mixed in the Eppendorf tube. The mixture was placed in the CFX96 Real-Time PCR Detection System (Bio-Rad) with the following temperature program: 1 cycle of 2 minutes at 95°C in the holding stage and 40 cycles of 5 seconds at 95°C and 30 seconds at 60°C in the cycling stage. The sequences of the forward and reverse primers are shown in Table 1.

2.4. Assessment of the Quantity of miR-155

The SNORD47 was used for the normalization based on the (Livak) 2-ΔΔCT method. A pooled healthy sample was prepared by mixing the vortexed samples of 36 healthy individuals. The pooled sample was used for a calibrator. was assessed for each sample based on the previously mentioned method, and the difference between of the sample and pooled sample () was calculated as follows:

Then, and the normalized value for each sample were calculated as follows:

2.5. Statistical Analysis

Data were assessed using the Statistical Package for the Social Sciences (SPSS) software (IBM Inc., Chicago, IL, USA) version 20. Graphs were created using GraphPad Prism 8.0 (GraphPad Software Inc., California). Data were checked for normality using the Kolmogorov-Smirnov test. Mean or median and standard deviation (SD) were used to present continuous variables, while frequency and percentage were used to present categorical variables. Comparison between groups was performed using the Student -test and one-way and two-way analysis of variance (ANOVA). A correlation between study parameters was assessed using the Pearson correlation coefficient. The receiver operating characteristic (ROC) curve and the area under the curve (AUC) were performed to assess the diagnostic value of miR-155 for the detection of BC and differentiation between grades, stages, lymph node metastasis, tumor size (T size), HER2, ER, PR, and Ki-67. The cutoff value for miR-155 for each diagnosis was calculated using the Youden index. Binary logistic regression and linear regression were performed to assess the relationship between study parameters and BC. The value lesser than 0.05 was considered statistically significant.

3. Results

A total of 72 subjects (36 BC patients and 36 controls) participated in this study. The mean age of the subjects in BC and control groups was and years, respectively. The mean of the body mass index (BMI) in the breast cancer patients and control group was and , respectively. The mean number of pregnancies in the BC patients and control group was and , respectively (Table 2).

Demographic characteristics of study subjects are presented in Tables 2 and 3 . There were no significant difference between the breast cancer and control groups in terms of age (), BMI (), number of pregnancies (), age of menarche (), history of abortion (), and contraceptive drug usage ().

Clinical characteristics of the breast cancer group are shown in Table 4. The most common tumor grade was MD (15, 41.7%), followed by PD (11, 30.6%). The most common cancer stage was stage II (17, 47.2%), followed by stage III (11, 30.6%). The most common types of receptors were HER2 negative (27, 75%), PR positive (19, 52.8%), ER positive (24, 66.7%), and (22, 61.1%).

The comparison of miR-155 expression between study groups is presented in Table 4. The expression of miR-155 in BC patients was times greater than that in the control group (). The miR-155 expression was significantly higher in all grades, stages and T sizes, and lymph node metastases as well as ER, PR, Ki-67, and HER2 categories compared to that of the control group () (Table 4 and Figure 1).

The miR-155 expression was significantly higher in WD, MD, and PD grades compared to that of controls (, , and , respectively). The miR-155 expression was significantly higher in grade III than in grade I () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in stage I, stage II, and stage III compared to that of the control group (, , and , respectively). However, there was not any significant difference between miR-155 in stage II compared to stages I and III. The multivariate analysis with BMI as a confounder revealed a considerable difference in terms of miR-155 expression and stage () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in the large tumor size (T size) compared to that of the control group (). The expression of miR-155 was significantly higher in T1, T2, and T3 compared to that of the control group (, , and , respectively), and also, there was a significant difference between T1, T2, and T3 in the patient group () (Table 4 and Figure 1).

The expression of miR-155 was significantly greater in lymph node involvement compared to that of the control group (). The miR-155 expression was higher in positive and negative lymph node metastases compared to that of the control group ( each group), but there was no significant difference between positive and negative lymph node involvement () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in the estrogen receptor (ER) compared to that of the control group (). The expression of miR-155 was significantly higher in ER+ and ER- compared to that of the control group ( each group). However, there was not any significant difference between miR-155 in ER+ compared to ER- () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in the progesterone receptor (PR) compared to that of the control group (). The expression of miR-155 was significantly higher in PR+ and PR- compared to that of the control group ( each). However, there was not any significant difference between miR-155 in PR+ compared to PR- () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in HER2 compared to that of the control group (). The expression of miR-155 was significantly higher in HER+ and HER- compared to that of the control group (). However, there was not any significant difference between miR-155 in HER+ compared to HER- () (Table 4 and Figure 1).

The expression of miR-155 was significantly higher in Ki-67 compared to that of the control group (). The expression of miR-155 was significantly higher in and compared to that of the control group ( each). However, there was not any significant difference between miR-155 in compared to () (Table 4 and Figure 1).

Two-way ANOVA results showed that age, BMI, number of pregnancies, age of menarche, contraceptive drug usage, and history of abortion had no significant effect on expression level (), and the difference between groups was due to BC for age, BMI, number of pregnancies, antipregnancy drugs, and abortion () (Table 5 and Figure 2).

The relationship between miR-155 expression and the tumor grade, tumor stage, T size, node metastasis, and tumor markers is presented in Tables 6 and 7. The binary logistic regression revealed that miR-155 expression and the grade and stage of the tumor were the predictors of BC (, , and ), respectively. It was shown that probability of BC increased by 6.15 times for every one-unit increase in mir-155 expression. Similarly, a one-unit increment in the grade and stage of the tumor was associated with 10.28 and 7.61 times increased risk of BC. Also, this analysis was performed for each parameter such as tumor grade and stage, T size, node metastasis, ER, PR, HER2, and Ki-67, compared to age, BMI, number of pregnancies, contraceptive drug usage, history of abortion, and age of menarche (Tables 6 and 7).

The linear regression analysis was performed between miR-155 and age, BMI, and number of pregnancies. In this study, it was found that miR-155 had no relationship with age and number of pregnancies ( and , respectively), while there was a significant relationship with BMI () (Figure 3).

The ROC curve was used to identify the sensitivity and specificity of the tumor grade (), tumor stage (), T size (), node metastases (), PR (), ER (), HER2 (), and Ki-67 () (Tables 8 and 9 and Figure 4).

Based on the ROC curve, the optimal cutoff in the expression of miR-155 for detecting BC was 1.40 (Youden index: 0.6667), which resulted in the sensitivity and specificity of 77.78% (95% CI: 61.92% to 88.28%) and 88.89% (95% CI: 74.69% to 95.59%), respectively (Table 8).

If the miR-155 expression was used as the biomarker for BC grades, at the Youden cutoff of 0.3626, it could identify low-grade (WD and MD) from high-grade (PD) BC with a sensitivity of 99.76% and specificity of 88.25% (Table 9 and Figure 5).

If the estrogen receptor (ER) was used as a biomarker for distinguishing BC in patients, expression of miR-155 at the Youden cutoff of 0.3397 could identify the healthy group from ER- and ER+ in BC with a sensitivity of 80.75% and specificity of 80.64% (Table 9).

4. Discussion

Although histochemistry is considered the standard method for the diagnosis of breast cancer, it still faces challenges in the preanalytical and analytical stages of cancer [11]. The preanalytical factors that may influence the histochemistry methods include delay in tissue fixation, type of chemical fixation, and duration of tissue fixation. At the same time, the analytical stage might be affected by the detection of cutoff and personal variations in the visual examination of the tissue samples [12]. Based on the challenges in the histochemistry methods, it was suggested to add biochemical parameters to assist the diagnosis of breast cancer [9].

The findings of this study revealed that the mean age of the subjects with breast cancer was 47.36 years which was similar to the finding previously reported in Iran [13] but lower than those in the previously reported studies conducted in other countries [14, 15]. The findings of this study, along with the results of the previous studies, indicate that the incidence of breast cancer in Iran is higher in younger ages than that in Western countries.

In this study, the expression of miR-155 in BC patients was times more than that in the control group. This finding was in line with the results of the previous studies that reported increased miR-155 expression by 2.62 to 8.8-fold in breast cancer patients compared to controls [16]. In the studies of Guo et al. [15], Sun et al. [17], and Zhang et al. [18], the increase in miR-155 expression in BC was 2.94 times, 2.62 times, and 2.87 times, respectively.

The findings of this study revealed that the highest expression of miR-155 was among grade 3 (PD) of breast cancer. Furthermore, the miR-155 expression in WD, MD, and PD was - (), - (), and - () fold higher than that of the healthy controls. This finding was in line with the results of a previous study that reported increased miR-155 expression with an increased tumor grade () [19].

The findings of this study revealed that the miR-155 expression increases with the increased stage of the tumor (). This finding was in line with the results of previous studies that reported increased miR-155 expression with an increased tumor stage ( and [19, 20]), and also, in previous studies, a significant relationship was observed between miR-155 expression and the stage of breast cancer tumor [15, 21]. However, Mar-Aguilar et al. [13] and Sun et al. [17] reported no significant relationship between miR-155 expression and the stage of breast cancer tumors ( and , respectively). In contrast to a previous study, the miR-155 expression was highest in stages II and III compared to stages I and IV in only one study [14].

The findings of this study also revealed a significant link between the tumor size and the miR-155 expression (). It was in line with the study of Lu et al. [12]. Similar to our research, the most frequent tumors were 20 to 50 mm in size, and miR-155 expression was significantly higher in this size than in other sizes. However, there was no significant relationship between miR-155 expression and the tumor size in the studies of Sun et al. () and Chen et al. () [17, 19].

Increased expression of miR-155 was observed in the BC group, and there was a statistically significant relationship between the expression level and lymph node metastasis. Lymph node involvement was observed in 17 (47.2%) subjects, but there was no significant difference in miR-155 expression between lymph node involvement and noninvolvement in this study (). This finding was in line with Sun et al.’s study () [17] while this finding was in contrast with the previously reported relationship between miR-155 expression and lymph node invasion in previous studies which was confirmed by the studies of Chen et al. () [19], Zheng et al. () [20], Elshimy et al. () [22], and Amal Fawzy et al. [23].

In addition to myriad risk factors, most notably age, family history, and hormonal factors, some various behaviors and characteristics can be classified into breast cancer, including the histologic features of the malignant tumor grade, tumor stage, and indices measurable by immunohistochemistry, most commonly PR, ER, HER2, and Ki-67 [21, 24].

The result of the current study revealed no linkage between miR-155 expression and PR (), ER (), and HER2 () positivity. These results were partly in line with the findings of previous studies [12, 17]. There is controversy regarding the relationship between miR-155 expression and PR positivity. While similar findings were reported regarding the link between miR-155 expression and PR positivity in one study [17], no relationship was observed between miR-155 expression and PR positivity [12]. In this study, 25% of patients were HER2-positive and 75% HER2-negative, and 66.7% were ER-positive and 33.3% ER-negative. There was no significant relationship between miR-155 expression and HER2 () and ER () positivity. The results are corroborated by the studies of Lu et al. [12], Sun et al. (HER2 (), ER ()) [17], and Chen et al. [19] (ER (p =0.977), PR (p = 0.09)).

This study also failed to find a significant effect for Ki-67 on fold expression of miR-155 among BC patients (). This finding was in line with the previous study on 45 BC patients. Zheng et al. showed that upregulated miR-155 expression was associated with a higher proliferation index () () [20]. However, Bašová et al. [25] reported the link between miR-155 expression and in 134 patients ().

There was no significant difference in the expression of miR-155 between cancer and control groups in terms of age (). There was no any linkage between the expression of miR-155 in BC old compared to the healthy group with the same age () and also in BC compared to the healthy group (). This result was in line with the finding of Guo et al. [15]. They showed that there is no relationship between miR-155 in the BC and ≥45 years old (). Chen et al. [19] reported that they did not find any significant difference between miR-155 expression and age groups ().

This study revealed a significant effect for abortion on -fold expression of miR-155 in the BC group. In this study, a significant relationship was found between those who had a history of abortion and those who had no history (). This result was in line with the findings of previous studies. Guo et al. [15] reported that the history of abortion has a direct effect on upregulated miR-155 expression ().

To the best of our information, this study was the first paper that assessed the miR-155 expression in BC patients based on contraceptive drug usage (). The miR-155 expression in patients who had the background of using contraceptive drugs and in patients who had never use these drugs was, respectively, - () and - () fold higher than that in the healthy controls.

To the best of our knowledge, this study was the first study that assessed the miR-155 expression in BC patients based on their number of pregnancies. Although no significant difference in miR-155 fold expression and number of pregnancies (p = 0.266), the miR-155 expression in patients who had ≤4 parturitions and in patients who had >5 calving was, respectively, - (p <0.001) and - (p <0.001) fold higher than that in the healthy controls.

The study also examined the association between miR-155 and menarche age. There was no significant association between miR-155 in patients under 13 years and over 13 years of menarche age (). The expression of miR-155 in patients younger than 13 years was 1.67 times higher than that in healthy (). Also, the expression of miR-155 was 1.75 times higher in BC old compared to the healthy group with the same age (). While this finding was in contrast with the previously reported relationship between miR-155 expression and menarche age, Guo et al. [15] showed that single-factor analysis of miR-155 expression among clinical pathologies indicated that miR-155 expression significantly differed among patients according to menarche age (). They also reported that subjects with a menarche age of <13 years, several artificial abortions, high BMI, and a family history of breast cancer had a relatively high miR-155 expression [15].

The current study found that the expression of miR-155 was significantly higher in the cancer group compared to controls in all BMI categories.

In Guo et al.’s study [15], menarche age under 13 and BMI over 24 kg/m2 were significantly associated with increased miR-155 expression, whereas in this study, there was only a relationship between BMI and expression level. There was an increase in BMI, although there was no statistically significant relationship between the patient groups (). Consistent with the previous study reported by Guo et al., the mean expression of miR-155 compared to that of the healthy group in terms of BMI less than 25 kg/m2, between 25 kg/m2 and 30 kg/m2, and more than 30 kg/m2 was (), (), and (), respectively () [15].

The ROC curve is a graphical presentation of screening properties to determine the best cutoff point. The AUC, sensitivity, and specificity of miR-155 were found to be 0.89, 77.78%, and 88.89%, respectively (), and the cutoff was 1.4 (Youden index: 0.6667). In a previous study, Mar-Aguilar et al. [13] reported that the AUC for miR-155 for the detection of BC was 0.99 (95% CI: 0.9866 to 1.0022), and the sensitivity and specificity of miR-155 were reported to be 94.40% and 100%, respectively, and the optimal cutoff was 7.92. In another study, the AUC, sensitivity, and specificity of the miR-155 for detecting BC were reported to be 0.879 (95% CI: 0.820-0.868), 84.2%, 88.1%, respectively, and the cutoff value was 1.24 [15]. Han et al. [14] reported that the AUC, sensitivity, and specificity of miR-155 for detecting BC were 0.749, 100%, and 51.02%, respectively, and the cutoff value was -1.17. Zhang et al. [18] showed that the AUC, sensitivity, and specificity of the miR-155 for detecting BC were 0.692 (95% CI: 0.625-0.754), 66.0%, and 68.9%, respectively, and the cutoff value was 0.321 (Youden index: 2.2). In another study, Sun et al. [17] reported that the AUC for miR-155 for the detection of BC was 0.801 (95% CI: 0.734 to 0.868), the sensitivity and specificity of miR-155 were reported to be 65.0% and 81.8%, respectively, and the optimal cutoff was 1.91.

The findings of this study also showed that miR-155 expression could be used in the differentiation of BC grades with a sensitivity of 81.82%, a specificity of 72%, and the cutoff of 1.71 (Youden index: 0.3626) (). To the best of our knowledge, no previous study assessed the sensitivity and specificity of miR-155 for differentiating between BC tumor grades.

One of the limitations of this study was the difficulty in obtaining consent from women to participate in the study along with the missing data in patient documents and the high cost of diagnostic kits, which resulted in the restriction of sampling due to the limited budget of the study. The findings of this study justify the need for further studies with a higher budget in the early detection of breast cancer by using biochemical markers.

In summary, the findings of this study indicate that the miR-155 expression can assist in diagnosis, prognosis, and TNM grading, including lymph node involvement and metastasis in breast cancer patients.

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflict of interest in the publication of this article.

Authors’ Contributions

Fatemeh Hosseini Mojahed and Amir Hossein Aalami have contributed equally to this study.

Acknowledgments

We would like to thank Dr. Hossein Abdeahad and the Mashhad University of Medical Sciences for all their technical support.

Supplementary Materials

Graphical abstract. (Supplementary Materials)