Abstract
Background/Aim: P-Element-induced wimpy testis (PIWI) proteins in complex with PIWI-interacting RNA (piRNA) are involved in epigenetic regulation of gene expression in germline cells. Aberrant expression of piRNA and PIWI proteins have been identified in various types of tumour cells. The aim of this study was to evaluate the expression profiles of PIWI-like protein-1, -2 (PIWIL1 and PIWIL2), their immunohistochemical (IHC) characteristics in ductal breast cancer, and determine their correlation with clinicopathological parameters of this type of cancer. Materials and Methods: Material for IHC studies comprised of 101 invasive ductal carcinoma (IDC) cases and 31 mastopathy tissues. Frozen fragments of paired tissue specimens (tumour and adjacent non-malignant breast tissue) taken from 55 patients with IDC and 18 samples of mastopathy were used for molecular studies using real-time polymerase chain reaction (RT-PCR). Results: A statistically significantly higher level of PIWIL1 and PIWIL2 was found in IDC compared to mastopathy samples (p≤0.0001). Increased expression of PIWIL1 was correlated with increased PIWIL2 expression in breast cancer tissue. Surprisingly, PIWIL1 mRNA was detected only in cancer and mastopathy, but was not found in most normal breast tissues, although it is noteworthy that the PIWIL2 mRNA level was statistically significantly lower in mastopathy and IDC samples compared to normal breast tissues. Conclusion: Our results affirm the hypothesis that reactivation of PIWI expression in various caner types is crucial for cancer development.
Breast cancer is the most frequently diagnosed tumor in women and the leading cause of cancer-related death. Currently recommended classification of breast cancer is based on expression of estrogen receptor (ER), progesterone receptor (PR), human epidermal growth factor receptor 2 (HER2) and the proliferation-associated nuclear protein Ki67 to identify clinical supbgroups for appropriate therapy. Nevertheless, triple-negative breast cancer characterized by a lack of ER, PR and HER2 expression has a worse clinical outcome and an increased risk of recurrence due to lack of targeted therapies (1). Furthermore, heterogeneity of the tumour cells within the cancerous lesion has led to the concept of the existence of ‘cancer stem cells’ (CSCs or tumour-initiating cells), which present characteristics of stem cells, such as self-renewal capacity, fast proliferation and multilineage differentiation (2, 3). These unique properties of CSCs explain the failure of many antitumor therapies, which in most cases target rapidly dividing cells. Further extensive characterization of CSCs in individual tumour types is needed in order to identify specific markers for the heterogeneous population of CSCs and consequently be able to target CSCs selectively.
P-Element-induced wimpy testis (PIWI) proteins bind to a newly discovered class of non-coding small RNAs called PIWI-interacting RNAs (piRNA). Complexes of piRNA and PIWI proteins are predominantly expressed in germline cells, where they have been demonstrated to be involved in germline development, stem cell self-renewal and gametogenesis (4). The human PIWI proteins form a subfamily of the highly conserved Argonaute proteins and consist of four members: PIWIL1, PIWIL2, PIWIL3 and PIWIL4 (5). To date, a model of function for PIWI and piRNA in normal tissue and cancer has been proposed. piRNAs, along with abundant expression of PIWI, in germline stem cells regulate transposone silencing through DNA methylation during spermatogenesis (6-8). In this way, germline cells (expressing PIWI at a high level) develop normally into somatic tissues in which PIWI protein expression is absent. However, in cancer cells, high expression of PIWI and piRNAs results in aberrant DNA methylation, silencing of tumour-suppressor genes and an abnormal ‘stem-like’ state of cancer cells (6, 7).
Numerous data have shown the aberrant expression of PIWI genes/proteins in various cancer types, suggesting that both PIWI proteins and selected piRNAs are involved in tumorigenesis (4, 9). Preliminary studies suggest that over- and ectopic expression of PIWIL1 is associated with several different types of tumours (10, 11). Based mainly on immunohistochemical (IHC) studies, the increased expression of PIWIL1 has been detected in colon, breast, oesophageal, pancreatic, gastric and hepatocellular carcinoma and was correlated with histological grade of tumour, clinical stage and poorer clinical outcome of patients (12-17). Our recent study showed significantly higher mRNA level of PIWIL1 and decreased level of PIWIL2 mRNA in colorectal cancer tissues compared to corresponding non-cancerous samples (18). Recent studies suggest that deregulated expression of PIWI genes/proteins is linked to deregulated cell proliferation, altered apoptosis, genomic instability and invasion of tumours. PIWI proteins seem to be promising CSC markers due to the fact that their deregulated expression was also correlated with clinicopathological features and survival of cancer patients (19, 20).
The aim of our study was to evaluate the expression profiles of the PIWIL1 and PIWIL2 in invasive ductal carcinoma (IDC) samples in comparison to mastopathy tissues. We aimed to determine possible prognostic significance of PIWIL1 and PIWIL2 in breast cancer development and progression.
Materials and Methods
Clinical samples. The study was conducted on two separate patient cohorts: one for IHC and one for real-time polymerase chain reaction (PCR) studies. For IHC studies, paraffin blocks containing diffuse cystic mastopathy tissue derived from 31 patients, and IDC from 101 patients, treated at the Lower Silesian Oncology Center in Wrocław, Poland between years 1999-2005. All the tissue specimens utilized for the study were obtained before beginning treatment. The patient's clinicopathological data are listed in Table I. Paraffin sections stained with haematoxylin and eosin (H&E) were used to verify the diagnosis and malignancy grade of the tumors, according to World Health Organisation criteria (21). All patients with IDC were treated by mastectomy or quadrantectomy, and subsequent axillary lymph node resection. The mean patient age at diagnosis was 55.9 years. Adjuvant chemotherapy was applied to 79 patients with IDC. Only seven patients received neoadjuvant chemotherapy prior to primary tumor resection and 50 of the patients were treated with adjuvant radiotherapy.
The second patient cohort comprised 55 tissue samples of patients with IDC and 18 tissue mastopathy samples and were used for real-time polymerase chain reaction (RT-PCR) studies. Fresh frozen, paired tissue specimens (tumour IDC and adjacent non-malignant breast tissue samples) were obtained from patients with breast cancer operated at the Regional Specialist Hospital in Wroclaw, Research and Development Centre, between 2011 and 2014. Mastopathy samples and non-malignant breast tissue were used as control. The fresh tissue samples were collected in RNAlater™ Stabilization Solution (Ambion by Life Technologies, Carlsbad, CA, USA), frozen and stored at −80°C until use. All tumour and normal tissue samples were submitted for histopathological analysis. Detailed description of patients is presented in Table I.
Clinical and pathological data of studied patients with breast cancer.
The study protocol was approved by the Medical Ethics Committee of Regional Specialized Hospital, Research and Development Centre in Wroclaw (KB/ No. 10/2015). The study was conducted in accordance with the Helsinki Declaration of 1975 and its later amendments or comparable ethical standards.
Quantitative real-time RT-PCR assay. Total RNA was extracted from (30 mg) tumour, non-cancerous adjacent tissues and mastopathy tissues by a combination of TRI Reagent® (Ambion by Life Technologies) method with silica-based purification technology using NucleoSpin RNA isolation kit (Machery-Nagel GmbH & CO. Kg, Düren, Germany) according to the manufacturers' instructions. Briefly, disruption and homogenization steps with MagNA Lyser Green Beads were performed in 1 ml of TRI Reagent® (Ambion by Life Technologies) using the MagNA Lyser System (Roche Molecular Systems, Inc. Pleasanton, CA, USA) (3×30 s at 7000 rpm). The homogenized samples were incubated at room temperature for 5 min and then 0.2 ml of chloroform was added. The mixture was thoroughly shaken by vortexing for 15 s and incubated at room temperature for an additional 15 min. Phase separation was carried out by centrifugation at 4°C and 12,000 × g for 15 min. The aqueous phase was collected and RNA precipitation was performed by adding an equal volume of 70% ethanol. The RNA samples were further purified using spin columns according to the manufacturer's recommendation.
The concentration and purity of isolated RNA was measured in a NanoDrop2000 spectrophotometer (ThermoFischer Scientific, Wilmington, DE, USA). For the reverse transcription reaction, 250 ng of RNA was used. The reverse transcription reaction was performed using the High Capacity cDNA Reverse Transcription Kit (Applied Biosystem by Life Technologies) according to the procedure of the manufacturer. Diluted cDNA (4- fold in nuclease-free water) was used for the subsequent real time PCR using TaqMan Fast Universal PCR Master Mix and TaqMan Gene Expression Assays (Applied Biosystem by Life Technologies): Hs01041737_s1 specific for human PIWIL1, Hs01032720_m1 for PIWIL2 and Hs039299097_g1 for glyceraldehyde 3-phosphate dehydrogenase (GAPDH). All target genes and the housekeeping gene (GAPDH) were run simultaneously and reactions were carried out in triplicate on a 96-well optical reaction plate in 20 μl reaction volume, using the StepOnePlus Real Time PCR System (Applied Biosystem by Life Technologies). For quantification, the samples were normalized against the expression of GAPDH. Relative quantification (RQ) of mRNA for the examined genes was calculated using the Pfaffl's method (22). Fold-change cut-offs of >2 and <0.5 were used to identify up-regulated and down-regulated genes, respectively.
Immunohistochemistry. IHC was performed on 4 μm formalin-fixed paraffin-embedded sections of specimens. The slides were deparaffinized, rehydrated and antigen retrieval was carried out by boiling the sections in EnVision™ FLEX Target Retrieval Solution high pH (PIWIL1, PIWIL2) using a PTLink (Agilent Technologies, Dako, Santa Clara, CA, USA). All reactions were performed using an Autostainer Link 48 automated immunostainer (Agilent Technologies, Dako). Activity of endogenous peroxidase was blocked by 5 min exposure to EnVision™ FLEX Peroxidase-Blocking Reagent. Polyclonal goat anti-human PIWIL1 (sc-22685, 1:50 dilution; Santa Cruz Biotechnology, Inc., Dallas TX, USA), and polyclonal goat anti-human PIWIL2 (sc-67502, 1:100 dilution; Santa-Cruz Biotechnology, Inc.) were used as the primary antibodies. Primary antibodies were diluted in EnVision™ FLEX antibody diluent (Agilent Technologies, Dako). Tested sections were incubated with primary antibodies for 20 min at room temperature, followed by incubation with secondary antibody conjugated with horseradish peroxidase (EnVision™ FLEX/HRP; Agilent Technologies, Dako) for 20 min. The reactions were visualized using freshly prepared substrate for horseradish peroxidase (diaminobenzidine, EnVision™ FLEX Working Solution), with incubation for 10 min. Additionally, all slides were counterstained using EnVision™ FLEX Hematoxylin for 5 min. After dehydration in graded ethanol concentrations (70%, 96%, absolute) and in xylene, all slides were coversliped in SUB-X Mounting Medium (Agilent Technologies, Dako).
Histopathological examination and analysis of IHC reactions. All IHC preparations were examined using an Olympus BX-41 light microscope (Olympus, Shinjuku, Tokyo, Japan). The intensity of IHC reactions for PIWIL1 and PIWIL2 was estimated using the semi-quantitative immunoreactive (IRS) scale of Remmele and Stegner (23). IRS score represent the results of multiplication of the percentage of positively stained cells and the intensity of staining (values from 0-12). IRS of 8 or lower were defined as negative, values above 8 as positive.
Statistical analysis. Statistical analysis was performed with GraphPad Prism 5.0 (GraphPad Software, Inc., La Jolla, CA, USA). Shapiro–Wilk test was used to evaluate the normality assumption of examined groups. To compare the differences between the expression of examined markers in the all patient pairs of groups and clinicopathological data the unpaired t-test and Mann–Whitney test were used. To compare the differences between more than two groups, the Kruskal–Wallis and Dunn's multiple comparison tests were performed. Additionally, the Spearman correlation test was used to analyse the existing correlations. The Kaplan–Meier method was used to construct survival curves. In all analyses, the results were considered to be statistically significant at p<0.05.
Results
PIWIL1 and PIWIL2 expression by IHC and selected clinicopathological characteristics in patients with breast cancer. The protein expression level of PIWIL1 and PIWIL2 in sections of IDC and mastopathy samples was analyzed using IHC. In all cases, H&E staining was performed in order to determine tumour grade and evaluate the extent of necrosis in the tumours. Representative IHC staining of PIWIL1 and PIWIL2 in mastopathy and breast cancer samples is presented in Figure 1. All the studied samples demonstrated cytoplasmic localization of PIWIL1 and PIWIL2 in cancer cells.
Absent or weak PIWIL1 and PIWIL2 staining was noted in mastopathy tissues in contrast to the statistically significant strong expression in breast IDC samples (Figure 2). As presented in Table II, in analysed cancer cases, 75% demonstrated negative expression of PIWIL1 (IRS≤8), while 26.26% manifested positive expression (IRS>8). IHC reaction for PIWIL2 was negative (IRS<8) in 90.9% and positive in 11.11%. Moreover, significant differences for both PIWIL1 and PIWIL2 were noted using Mann–Whitney U-test between mastopathy cases and those of IDC (p<0.001) (Figure 2C and D).
By Spearman's rank correlation coefficient, high expression of PIWIL1 in IDC cases was statistically significant correlated with elevated expression of PIWIL2 (correlation coefficient r=0.7933; p<0.001) (Figure 3). However, no significant correlation was found between PIWIL1 and PIWIL2 expression in IDC and the expression of ER, PR and HER2 receptors, the size of primary tumours (pT), metastasis to lymph nodes (pN) and distant metastasis or the age of the patients. A positive association between PIWIL1 expression and tumour size (T) as well as PIWIL2 and distant metastasis on the border of statistical significance was observed. The correlations between the presence of PIWIL1 and PIWIL2 evaluated in IDC and clinicopathological parameters of the patients are summarized in Table II.
Immunohistochemical (IHC) identification of PIWI-like protein-1 (PIWIL1) (A, B) and PIWIL2 (C, D) in mastopathy (A,C) in comparison to invasive ductal breast cancer cases (B, D). Polyclonal goat anti-human PIWIL1 (N-17, sc-22685), polyclonal goat anti-human PIWIL2 (K-18, sc-67502) were used for PIWIL1 and PIWIL2 identification, respectively (Santa Cruz Biotechnology, Inc). Magnification, ×200.
In order to evaluate the prognostic value of PIWIL1 and PIWIL2 expression in IDC cases, the survival of patients depending on the expression level was analysed. However, no significant differences were noted in patient survival regarding the intensity of PIWIL1 and PIWIL2 expression (Figure 4A and B). The variables which affected the duration of patient survival without statistical significance were tumor size (T staging) (p=0.062) and tumor malignancy grade (p=0.65) (Figure 4C and D).
PIWIL1 and PIWIL2 mRNA expression level in IDC, adjacent non-malignant and mastopathy samples. In order to determine mRNA level of PIWIL1 and PIWIL2 in IDC, mastopathy and normal tissue adjacent to tumor, real-time PCR analysis was performed. As presented in Table III, in normal breast tissues, expression of PIWIL1 was mostly not detected (in 52 samples out of 55), while in mastopathy and IDC samples, the PIWIL1 expression was highly elevated. Significantly lower expression of PIWIL2 was detected in IDC tissues (RQ=0.33) and mastopathy (RQ=0.43) as compared to adjacent noncancerous samples (RQ=1.33) (Figure 5B, Table III). The data showed a statistically significant correlation between low PIWIL1 and PIWIL2 mRNA expression and increased tumour grade for breast IDC tissues (Figure 6A and B, respectively, Table IV).
Discussion
The complex of PIWI proteins and piRNAs regulate the gene expression at the epigenetic post-transcriptional level. In germline cells, PIWI proteins are essential for germline development, stem cell self-renewal and gametogenesis (24). Aberrant expression of different PIWI types on the mRNA and protein levels as well as different piRNAs have been reported in several types of tumour (9, 25). Interestingly, overexpression and ectopic expression of PIWIL1 has been detected in many types of cancer and in most cases an increased level of PIWIL1 was correlated with clinical and pathological parameters of tumour such as histological tumour grade and advanced clinical stage (4). Increased expression of PIWIL2 was identified among others in breast, cervical, gastric, ovarian, prostate and colorectal cancer (12, 26-30). Statistically significant increased level of PIWIL1 and decreased level of PIWIL2 was reported for colon and bladder cancer (18, 31). With regard to other members of the PIWIL protein family, increased expression of PIWIL4 was observed in renal carcinoma (32). Moreover, expression of four members of the PIWI protein family was changed in tumour tissue compared with non-tumorous adjacent tissue (33). Non-statistically significantly increased expression of PIWIL1 and PIWIL3 and statistically significant down-regulation of PIWIL2 and PIWIL4 in breast tumour tissues was described. However, among four PIWIL genes, only PIWIL3 and PIWIL4 were correlated with overall survival, and PIWIL3 alone with recurrence-free survival (33).
Association of PIWI-like protein-1 (PIWIL1) and PIWIL2 expression by immunoreactive score (IRS) with clinical and pathological data of the studied patients with breast cancer.
The expression level of PIWI-like protein-1 (PIWIL1) and PIWIL2 mRNA in invasive ductal carcinoma (IDC), adjacent non-malignant tissue and mastopathy samples. Data are the mean±standard deviation.
It has been proposed that reactivation of various PIWI proteins is crucial for cancer development and malignant progression, therefore it is very likely that PIWI proteins may be key markers for progression of various cancer types. However, the exact mechanism underlying the changes in the expression of PIWI proteins remains a subject of studies. Due to the fact that a limited number of studies demonstrated aberrant expression of different piRNAs in various cancer types, further studies are needed (33-37). Identification of piRNAs in blood or cancer tissues might be a valid method for identifying circulating or CSCs. Martinez and co-workers described somatic and malignant expression patterns of numerous piRNAs. Increased expression of numerous piRNAs were detected in tumour tissues and correlated with malignancy and clinical features of cancer (19). Although piRNA–PIWI complexes are found in somatic cells, their roles are not yet fully understood. Hashim and co-workers analysed piRNA patterns by direct small RNA sequencing in human breast cancer cells (38). They reported that PIWIL2 and PIWIL4, as well as other genes encoding components of the piRNA biogenesis pathway, are expressed in human breast cancer cells and that specific piRNAs are expressed at high levels both in cell lines (in vitro) and cancer tissues (in vivo). What is important is that a significant fraction of the detected piRNAs responde to oestrogen and transcription factor ERβ. Similar observations were made by Oner and co-workers (39), piRNAs may therefore be affected by hormones. The results of their studies showed that changes in the expression of piRNAs may be associated with hormone levels and dependent on the microenvironment of the tumour.
PIWI-like protein-1 (PIWIL1) (A, C) and PIWIL2 (B, D) expression in mastopathy compared with invasive ductal carcinoma (A, B) and according to tumour grade in patients with breast cancer (C, D). Data are presented as mean IHC score±standard deviation. Significantly different at ***p<0.001 and ****p<0.0001.
Correlation between expression of PIWI-like protein-1 (PIWIL1) and PIWIL2 in breast cancer samples by immunohistochemical score (r=0.7933, p<0.001 Spearman's correlation test).
Kaplan–Meier curves of overall survival for patients with invasive ductal breast carcinoma according to the expression of PIWI-like protein-1 (PIWIL1) (A), PIWIL2 (B), tumour grade (C) and histopathological grading (D). Cut-off points for the analysis were estimated based on the median expression.
Expression levels of PIWI-like protein-1 (PIWIL1) (A) and PIWIL2 (B) in invasive ductal carcinoma, adjacent non-malignant normal tissue and mastopathy samples. mRNA expression of PIWIL1 and PIWIL2 was measured by real-time reverse transcription-polymerase chain reaction assay. For quantification, the samples were normalized against the expression of glyceraldehyde 3-phosphate dehydrogenase (GAPDH) mRNA. Relative quantification (RQ) of mRNA for the examined genes was calculated using the method of Pfaffl (22). ***Significantly different at p<0.001. Data are presented as mean RQ± standard deviation.
In the present study, carried out on a well-characterized group of patients with IDC, the expression of PIWIL1 and PIWIL2 was determined with use of IHC and real-time PCR methods. To compare the level of expression of PIWIL1, and PIWIL2 in cancerous tissue, corresponding non-cancerous tissues, and mastopathy tissues were used. We found statistically increased expression of PIWIL1 and PIWIL2 in breast IDC samples in comparison to absence or weak staining of these proteins in mastopathy samples. To investigate whether gene expression directly affects the protein level of PIWIL1 and PIWIL2, mRNA expression was measured using the real-time PCR method. Interestingly, in normal breast tissue, expression of PIWIL1 in the vast majority of cases was not detected, while in paired cancerous samples was highly elevated. Simultaneously, significant lower expression of PIWIL2 was detected in IDC tissues and mastopathy compared to adjacent, non-cancerous tissues. These results are in accordance with our previous observations in colorectal cancer, where the mRNA levels of PIWIL1 in non-cancerous tissue was low or undetectable, while in paired cancerous tissue it was highly elevated. At the same time, the expression of PIWIL2 was significantly decreased (18). It is also important to emphasize that our study revealed negative correlation of decreased expression of PIWIL1 and PIWIL2 with increased grade of cancer. These findings are consistent with a previous study of Liu and co-workers demonstrating increased expression of PIWIL2 at the initiation stages of breast cancer (12), and suggest an important role of PIWIL2 for breast cancer development. Interestingly, PIWIL2 was observed to be ubiquitously and uniquely expressed in various stages of breast cancer and its expression pattern was associated with ER expression and proliferative marker Ki67 (12). Taken together, our study sheds new light on the molecular regulation of PIWIL in cancer development. The observed differences in protein expression may indicate reciprocal regulation between different piRNAs and PIWI genes in various cancer types. Reactivation of PIWI expression in cancer strongly suggests participation of these proteins in processes of tumour growth and differentiation (6, 7). The important issue awaiting further investigation is the epigenetic regulation in regard to the observed differences in PIWIL1 and PIWIL2 at the transcriptional and protein level. Based on current data, PIWI–piRNA complexes contribute to cancer development through aberrant DNA methylation resulting in genomic silencing and promoting a stem-like state of cancer cells (8). To conclude, aberrant expression of mRNA and protein level of PIWI proteins show a strong prognostic value.
The expression level of PIWI-like protein-1 (PIWIL1) (A) and PIWIL2 (B) according to tumor malignancy grade (G) in invasive ductal breast carcinoma. Expression of PIWIL1 and PIWIL2 was measured by real-time reverse transcription-polymerase chain reaction assay. Differences in mRNA expression level between individual groups according to grade were determined by the Mann–Whitney test. Significantly different at *p<0.05, **p<0.01 and ***p<0.001.
Association of PIWI-like protein-1 (PIWIL1) and PIWIL2 mRNA expression with clinical and pathological data of patients with invasive ductal breast cancer.
Acknowledgements
This publication was supported by the Leading National Research Center (KNOW, 2014-2018) of Wroclaw Center for Biotechnology.
- Received January 18, 2018.
- Revision received February 19, 2018.
- Accepted February 27, 2018.
- Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved











