Increased MicroRNA-34b and -34c Predominantly Expressed in Stromal Tissues Is Associated with Poor Prognosis in Human Colon Cancer

Abstract

The microRNA-34 family (miR-34a, -34b and -34c) have been reported to be tumor suppressor microRNAs (miRNAs) that are regulated by the TP53 and DNA hypermethylation. However, the expression, regulation, and prognostic value of the miR-34 family have not been systematically studied in colon cancer. To elucidate the roles of miR-34 family in colon carcinogenesis, miR-34a/b/c were measured in tumors and adjacent noncancerous tissues from 159 American and 113 Chinese colon cancer patients using quantitative RT-PCR, and we examined associations between miR-34a/b/c expression with TNM staging, cancer-specific mortality, TP53 mutation status and Affymetrix microarray data. All miR-34 family members were significantly increased in colon tumors, counter to the proposed tumor suppressor role for these miRNAs. Increased miR-34b/c were observed in more advanced tumors in two independent cohorts and increased expression of miR-34b/c was associated with poor cancer-specific mortality. While the expression of miR-34 family was not associated with TP53 mutation status, TP53 transcriptional activity was associated with miR-34a/b/c expression that is consistent with the proposed regulation of miR-34a/b/c by TP53. To examine where the miR-34 family is expressed, the expression of miR-34 family was compared between epitheliums and stromal tissues using laser microdissection technique. The expression of miR-34b/c was increased significantly in stromal tissues, especially in cancer stroma, compared with epithelial tissue. In conclusion, increased miR-34b/c predominantly expressed in stromal tissues is associated with poor prognosis in colon cancer. MiR-34 may contribute to cancer-stromal interaction associated with colon cancer progression.

About the Authors:

Yukiharu Hiyoshi

Contributed equally to this work with: Yukiharu Hiyoshi, Aaron J. Schetter

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Aaron J. Schetter

Contributed equally to this work with: Yukiharu Hiyoshi, Aaron J. Schetter

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Hirokazu Okayama

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Kentaro Inamura

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Katsuhiro Anami

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Giang H. Nguyen

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Izumi Horikawa

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Jason E. Hawkes

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Elise D. Bowman

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Suet Yi Leung

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Curtis C. Harris

* E-mail: [email protected]

Affiliation: Laboratory of Human Carcinogenesis, Center for Cancer Research, National Cancer Institute, National Institutes of Health, Bethesda, Maryland

Introduction

Colorectal cancer is the third most common cancer and the fourth leading cause of cancer-related death worldwide. Approximately 1.2 million new cases and about 600,000 deaths estimated occur annually [1]. Further studies to discover risk factors, biomarkers and therapeutic targets for colorectal cancer may help reduce the burden of this cancer.

MicroRNAs (miRNAs) are non-coding, single-stranded RNAs that are 18–25 nucleotides long. Since their discovery, miRNAs have been found to down-regulate the expression of multiple genes by interacting with the 3’UTR of target mRNAs to signal their degradation and/or blocking their translation into proteins [2]. Furthermore, these molecules regulate genes involved in a variety of cellular processes such as apoptosis, differentiation, proliferation, stress response and metabolism, and their deregulation or dysfunction plays an important role in carcinogenesis and tumor progression [3–5]. Specific miRNA expression signatures have been found in many malignancies, including colorectal cancer [6–8].

The miR-34 family is largely considered to be a tumor suppressor miRNA [9]. This family consists of three members: miR-34a, generated from a transcriptional unit on the human chromosome 1p36, and miR-34b and miR-34c, which are generated by processing of a bicistronic transcript from chromosome 11q23. Both transcripts have been shown to be under the direct positive control of p53 [10–12]. Several converging evidence demonstrated that ectopic expression of miR-34 family induces apoptosis, senescence, cell cycle arrest and inhibits migration and invasion [9, 13]. Therefore, miR-34 family is thought to be important mediator of p53’s tumor-suppressive activities. Besides, the miR-34a and miR-34b/c genes display varying levels of DNA methylation in numerous types of tumors including colorectal cancer, and may therefore represent tumor suppressor genes [14–16]. Although previous reports have demonstrated that down-regulation of miR-34 family members was associated with poor prognosis in several types of malignancies [17–20], the prognostic value of miR-34 family in colon cancer has not been systematically studied.

To evaluate the issue, we set out to systematically study the expression of the miR-34 family in two large, independent cohorts of colon cancer patients. In addition, the association between miR-34 expression and TP53 mutation, global gene expression was examined to clarify possible mechanisms of miR-34 regulation. To our knowledge, this is the largest and systematic study analyzing miR-34 expression in human colon cancer.

Material and Methods

Clinical samples

Pairs of primary colon tumor and adjacent non-tumor tissues came from 159 patients recruited from the University of Maryland Medical Center or Baltimore Veterans Affairs Medical Center between 1993 and 2002, and from 113 patients recruited from Queen Mary Hospital in Hong Kong between 1991 and 2000. Detailed backgrounds for each tissue donor, including age, sex, clinical staging, tumor location and survival times from diagnosis have been collected. The final date of follow-up was December 31, 2008 for the American cohort and August 16, 2004, for the Chinese cohort. Tumor histopathology was classified according to the World Health Organization Classification of Tumor system.

Ethics Statement

This study was approved by the Institutional Review Board of the National Institutes of Health, the Institutional Review Board of the University of Hong Kong/Hospital Authority Hong Kong West Cluster, and the Institutional Review Board for Human Subject Research at the University of Maryland. All participants provided their written informed consent to participate in this study, and the consent procedure was approved by the Institutional Review Board.

RNA isolation and quantitative reverse transcriptase-PCR

RNA from frozen tissue samples and cell lines was extracted using standard TRIZOL (Invitrogen, Carlsbad, California) methods. RNA was reverse-transcribed using TaqMan MicroRNA Reverse Transcription Kit for miRNA assay (Applied Biosystems, Foster City, California). Quantitative-PCR of microRNAs was performed using TaqMan assays (Applied Biosystems, Foster City, California) according to the manufacturer’s instructions with the 7900 HT Real-Time RT-PCR System (Applied Biosystems, Foster City). All assays were performed in triplicate. All TaqMan probes were purchased from Applied Biosystems; hsa-miR-34a-5p (Assay ID: 000426), hsa-miR-34b-5p (ID: 000427), hsa-miR-34c-5p (ID: 000428), RNU66 (ID: 001002). All microRNAs were normalized to RNU66 using the -ΔCt method. MicroRNA expressions were further normalized against the average expression of nontumor for each microRNA for each cohort separately with the following equation, [-ΔCt(average nontumor)]–[-ΔCt(sample)]. This sets the average nontumor expression = 1, with the log2 nontumor expression = 0.

Affymetrix mRNA array and Ingenuity Pathway Analysis (IPA)

RNA from 56 tumors from the American cohort were analyzed on Affymetrix U133A mRNA microarrays. Array data was imported and RMA normalized using Partek Genomic Suites 6.6. Data is publically available at (GSE44861). The average expression of miR-34b and miR-34c was used to analyze the association of miR-34b/c expression and global mRNA expression. Ingenuity Pathway Analysis (IPA 8.0, Redwood city, CA) was performed to identify molecules and pathways potentially involved with the expression of miR-34 family.

Laser microdissection

For Laser microdissection (LMD), 5 pairs of frozen tumor and non-tumor tissues were randomly obtained from the American cohort. Seven micron frozen sections for each tissues were placed on the PEN membrane slides (Leica Microsystems). Sections were fixed in 75% ethanol for 1 min, stained with Nuclear fast red for 2min, and dehydrated in a graded series of ethanol and xylen. Slides were air-dried and then microdissected using the LMD6000 system (Leica Microsystems). Microdissected epithelial and stromal cells from tumor and non-tumor tissues were collected directly into Qiazol, and total RNA was isolated using the RNeasy Micro kit (Qiagen) according to the manufacturer’s instructions. qRT-PCR was performed as described earlier.

Statistical analysis

Relative expression quantitation of microRNA was calculated with RQ manager 1.2 (Applied Biosystems). Expression values were used to analyze differences in microRNA expression between tumors and adjacent non-tumor tissues for all qRT-PCR data using Graphpad Prism v5.0 (Graphpad Software Inc, San Diego, California). Wilcoxon matched-pairs tests were used to compare microRNA expression between tumors and paired nontumor tissues. All trend tests reported are Cuzick nonparamtric test for trend across ordered groups [21] and were performed using Stata 9.2 (StataCorp LP, College Station, Texas). Kaplan-Meier analysis was performed using Graphpad Prism v5.0. All p values were 2-sided and statistical significance was set at p<0.05.

Results

This study used two independent cohorts of colon cancer patients, one consisting of 159 cases recruited from Maryland and the second cohort of 113 cases recruited from Hong Kong, China (Table 1). The median follow-up times were 90.3 and 84.6 months for patients in the American and the Chinese cohorts, respectively. The cohorts were similar in TNM staging, tumor histology and cancer-specific mortality. The 5-year survival rate was 46.4% for the American cohort and 49.6% for the Chinese cohort and were not significantly different from one another (P = 0.36, Kaplan-Meier test). In addition to the racial, geographic and cultural differences between these 2 cohorts, the American cohort was considerably older (average 66.2 years vs 55.8 years) had a higher percentage of men (71% vs 50%). The Chinese cohort had a higher percentage of tumors located at distal colon (80% vs 20%). TP53 mutation evaluated by immunohistochemistry was detected in 50.6% of the American cohort and 51.8% of the Chinese cohort.

[Figure omitted. See PDF.]

Table 1. Characteristics of study populations and tumors.

https://doi.org/10.1371/journal.pone.0124899.t001

MiR-34a, -34b and -34c are increased in colon cancer, and increased miR-34b/c is associated with advanced tumors and poor prognosis

We performed TaqMan qRT-PCR of miR-34a, miR-34b and miR-34c in colon tumors and adjacent non-tumor tissues from both the American and the Chinese cohorts. Unexpectedly, The expression of these 3 miRNAs was significantly increased in tumors compared to non-tumor tissues in both cohorts (Fig 1A and 1B). These increases were relatively small, but highly significant. In the US cohort, the fold-change increase was 1.42, 1.27 and 1.25 for miR-34a, -34b and -34c respectively. In the Chinese cohorts, the fold change increase was 2.26, 1.64 and 1.40, respectively. There was no significant association between miR-34a expression and TNM staging (Fig 2A and 2B left panels), while increased miR-34b and -34c were associated with higher TNM staging in both cohorts (test for trend, P<0.05; Fig 2A middle and right panels, P<0.01; Fig 2B middle and right panels). MiR-34b and miR-34c are co-transcribed from a bicistronic transcript and their expression should be highly correlated [9]. In agreement with this, the expression of miR-34b and -34c were highly correlated in both tumors and non-tumor tissues from each cohort (Pearson’s correlation r>0.9, P<0.0001; S1A and S1B Fig). Thus, we used combined miR-34b/c expression into the following survival analysis. The American and the Chinese cohorts were combined to increase the statistical power for this analysis and expression levels for each miRNA were dichotomized based on median values. There was no significant association between miR-34a expression and prognosis (Fig 2C left panel). On the other hand, colon cancer patients with high expression of both miR-34b and -34c in tumor showed significant poor prognosis (log-rank P<0.05; Fig 2C right panel). When the survival analysis was performed on each cohort separately, the differences were not statistically significant (S2A and S2B Fig). In multivariate analysis, miR-34b/c was not significantly associated with survival when adjusted for TNM staging (data not shown).

[Figure omitted. See PDF.]

Fig 1. The miR-34 family members significantly increase in colon tumors compared to adjacent non-tumor tissues.

(A, B) Dot plots represent miR-34a/b/c threshold cycle values from TaqMan qRT-PCR normalized to U66 in American cohort (A) and Chinese cohort (B). Horizontal bars indicate median expression value. Wilcoxon matched-pairs test was used. MiR-34a was not detectable in one patient and this patient was excluded from all analyses.

https://doi.org/10.1371/journal.pone.0124899.g001

[Figure omitted. See PDF.]

Fig 2. The high expression of miR-34b/c is associated with advanced colon tumor and poor prognosis.

(A, B) Tissues were ordered from TNM stage I to IV tumors in the American cohort (A) and the Chinese cohort (B). Dot plots represent miR-34a/b/c threshold cycle values from TaqMan qRT-PCR normalized to U66. Horizontal bars indicate median expression value. The Cuzick nonparametric test for trend was used to evaluate trends. (C) Kaplan-Meier survival analysis of all stage cases in the American and Chinese combined cohorts stratified by median miR-34a expression and combined miR-34b/c expression. *Cases with low miR-34b and/or low miR-34c. **Cases with both high miR-34b and high miR-34c.

https://doi.org/10.1371/journal.pone.0124899.g002

The association between miR-34 expression and TP53 mutation

Since miR-34a, miR-34b, and miR-34c were each systematically altered in colon tumors, we next attempted to determine possible reasons for this altered expression. Several studies indicate that p53 can bind to the promoter and activate the transcription of miR-34 family [10–12]. To address this, we examined whether miR-34 expression was associated with TP53 mutational status. TP53 mutational status can be examined either through immunohistochemistry or through direct sequencing. In the present study, the TP53 mutational status was assumed to be ‘wild-type’ or ‘mutant’ by immunohistochemical staining of p53 [22]. Contrary to our expectations, the expression of miR-34 family was not associated with TP53 mutational status in either the American and Chinese cohorts (Fig 3A and 3B). Additionally, TP53 mRNA expression was not associated with miR-34a (R = 0.005, p = 0.97) or miR-34b/c (R = 0.013, p = 0.92) expression based on microarray expression values of TP53 (S1 Table).

[Figure omitted. See PDF.]

Fig 3. TP53 mutation is not associated with the expression of miR-34 family.

(A, B) TP53 mutational status was evaluated by immunohistochemistry using Formalin-Fixed Paraffin-Embedded (FFPE) tissues. Dot plots represent miR-34a/b/c threshold cycle values of tumor tissue from TaqMan qRT-PCR normalized to U66 in American cohort (A) and Chinese cohort (B). Horizontal bars indicate median expression value. Mann-Whitney test was used. WT: wild-type p53, Mut: mutant p53.

https://doi.org/10.1371/journal.pone.0124899.g003

Ingenuity Pathway Analysis (IPA) predicts transcription regulators that potentially associated with miR-34 up-regulation

In order to investigate potential reasons for the altered expression of the miR-34 family, we used Affymetrix microarray data to identify genes that correlated with miR-34a and -34b/c and examined these genes with Ingenuity Pathway Analysis (IPA). MiR-34b/c correlated with 649 genes in Affymetrix microarray data of 56 tumors from the American cohort (S1 Table; P<0.001), and IPA was used to examine this set of genes. Increased expression of miR-34b/c was associated with active transcriptional activity of TP53, SMAD3, CBL, SNAI1, HTT, TWIST1 and CTNNB1 while it was associated with reduced transcriptional activity of MYCN, AHR, SMAD7, MYC, SPDEF and PPARG (Table 2, S2 Table). Similar analysis on miR-34a expression identified 385 genes that were correlated with miR-34a (S1 Table; P<0.001). Increased expression of miR-34a was associated with increased activity of 48 transcriptional regulators including TP53 and decreased activity of 17 transcriptional regulators (S3 Table).

[Figure omitted. See PDF.]

Table 2. List of Transcription Regulators that are potentially involved with miR-34b/c expression predicted by Ingenuity Pathway Analysis in 56 tumors.

https://doi.org/10.1371/journal.pone.0124899.t002

The comparison of miR-34 expression between colon epithelium and stroma with laser microdissection (LMD) technique

Altered expression of oncogenic miRNAs in colorectal cancer stromal tissue might be associated with cancer progression [23]. To examine where the miR-34 family is expressed in, we extracted RNA from cancer epithelium, cancer stroma, normal adjacent epithelium and normal adjacent stroma separately using laser microdissection for miR-34 expression analysis (S3 Fig). Five colon tumor and 5 adjacent non-tumor frozen tissue samples were obtained form the American cohort. MiR-34b and -34c expression levels in stroma were significantly higher as compared with that of epithelium in both cancer tissues and in normal adjacent tissues (Fig 4 middle and right panels). MiR-34a showed similar trend although the differences were not significant (Fig 4 left panel). We further examined the publicly available microarray dataset (GSE35602), in which miRNA expression profiles of micro-dissected colon cancer stromal and epithelial tissues were analyzed, demonstrating a clear trend of higher miR-34b and -34c expression in cancer stroma (S4 Fig). This microarray data are consistent with our own qRT-PCR analysis using LMD samples, indicating that miR-34b and -34c are expressed at higher levels in stomal tissues. In addition, we have performed H&E on 82 of the tissues from the American cohort to examine if there was an association with increased stromal area and any of the miR-34 family. There was no significant correlation between tumor/stroma ratio and expression levels of miR-34 (S5 Fig). These findings based on LMD analyses suggest that stromal cells serve the major source of miR-34b/c rather than epithelial cells. Hence, miR-34b/c expression detected in macro-dissected tumor samples (as shown in Fig 1) seems to mainly represent stromal expression of miR-34b/c.

[Figure omitted. See PDF.]

Fig 4. miR-34 family tends to be expressed predominantly in stromal tissue.

RNA samples from cancer epithelium, cancer stroma, normal adjacent epithelium and normal adjacent stroma were extracted separately using laser microdissection technique from 5 paired samples in the American cohort. Dot plots represent miR-34a/b/c threshold cycle values of tumor tissue from TaqMan qRT-PCR normalized to U48. Horizontal bars indicate median expression value. Un-paired t test was used.

https://doi.org/10.1371/journal.pone.0124899.g004

Discussion

In the present study, we demonstrated that: 1) miR-34a/b/c were increased in human colon cancers compared to adjacent non-tumor tissues. 2) Increased expression of miR-34b/c was observed in more advanced tumors and associated with poor prognosis. 3) The expression of miR-34 family was not associated with TP53 mutational status although pathway analysis showed association between miR-34 expression and TP53 transcriptional activity. 4) The expression of miR-34b and -34c were increased significantly in stromal tissues compared with epitheliums.

In colorectal cancer, prior studies have focused on a tumor suppressive function for miR-34 family as shown in many kinds of malignancy [9]. In vitro experiments have demonstrated that ectopic expression of miR-34a into human colon cancer cell lines induced apoptosis, senescence and inhibited migration and invasion by targeting E2F, SIRT1 and Fra-1 [24–26]. Considering the tumor suppressive function of miR-34 as shown previously, we had expected that miR-34 family was suppressed in human colon cancer and decreased expression was associated with poor prognosis when we started the present study. Previously, Roy et al. demonstrated that the expression of miR-34a and miR-34c were down-regulated in colon cancer specimens compared to adjacent colonic mucosa in small number expression analysis (n = 10) [27]. Contrary to our expectation, all miR-34 family members were increased significantly in colon cancer tissues compared with adjacent normal tissues, and increased miR-34b/c was associated with poor prognosis. As a matter of fact, the present data is consistent with several recent reports analyzing miR-34 expression and its prognostic value. MiR-34 family members have been found to increase in colorectal cancer [28], gastric cancer [29], renal cell carcinoma [30, 31], lung squamous cell carcinoma [32] and pediatric leukemia [33]. In addition, Wang et al. demonstrated that the miR-34a expression was higher in rectal cancer having more advanced TNM stage [28], and Svoboda et al. demonstrated that miR-34b up-regulation was associated with poor prognosis in triple-negative (ER-, PgR-, HER2-) breast cancer patients [34]. The fact that miR-34 family is up-regulated in human cancer and its up-regulation is associated with poor prognosis, seems to be conflicted with tumor suppressive role of miR-34 demonstrated in previous reports.

On the basis of these findings, we analyzed epithelial and stromal expression of miR-34 separately using laser microdissection technique to examine if cancer cells express miR-34. Interestingly, all miR-34 family members, especially miR-34b and -34c, were increased in stromal tissues compared with epithelium. In addition, the expression levels of miR-34 seemed to be the highest in cancer stromal tissues. These findings suggest that colon cancer stroma is mainly responsible for increased miR-34b/c expression, which were associated with advanced TNM stage and poor prognosis. However, the reason for the up-regulation of miR-34 in stromal tissue is still unclear. Nishida et al. demonstrated that several miRNAs were significantly altered in cancer stroma, suggesting an important role of stromal miRNAs in cancer progression [23]. Stromal tissues consist of inflammatory cells, immunocompetent cells, endothelial cells and fibroblasts. One possible explanation of our findings is that miR-34 up-regulation in stromal tissues might reflect inflammation of tumor tissues, as it is well known that inflammation contributes to the development of cancer, including colon cancer [35]. Mathe et al. demonstrated that mice treated with C. parvum, to induce inflammation, showed significant increase of miR-21, miR-29b and miR-34a/b/c expression in the spleen [36]. Taken together, the fact that miR-34b/c was increased in cancer stromal tissues, especially in more advanced colon cancer, might be associated with inflammation contributing to cancer progression. In a recent study, Bu et al. have demonstrated that miR-34a regulated bimodal switch targets Notch in colon cancer stem cells [37]. The finding may demonstrate one more possibility that different level of miR-34 expression could be significant for miR-34-mediated interaction between human colon cancer and stroma.

In the present study, we demonstrated conflicting findings. As shown in Fig 3, the expression of miR-34 family was not associated with TP53 mutational status. However, Ingenuity Pathway Analysis demonstrated that Increased expression of miR-34a/b/c was associated with active transcriptional activity of TP53 (Table 2, S2 Table, S3 Table).

Although it has been well known that miR-34 was direct target of p53, a recent study demonstrated that mouse models that contain deletions for miR-34a, miR-34b and miR-34c still retain p53 function and the expression of these microRNAs is largely independent of p53 status [38]. In general, cancer stromal tissues are considered to harbor wild-type p53 regardless of p53 mutational status of cancer tissues. Considering the fact that miR-34 was expressed in stromal tissues predominantly, miR-34 detected in macro-dissected cancer tissues seems to be mainly regulated by cancer stromal tissues harboring wild-type p53. That may be a plausible explanation that the expression of miR-34 family was not associated with TP53 mutational status of cancer.

Conclusions

In conclusion, we showed novel aspects of miR-34 family in human colon cancer. The demonstrated findings imply that increased miR-34 in stromal tissues may have important roles in colon cancer progression. Further analysis is needed to clarify the function and miR-34-mediated interaction between human colon cancer and stroma.

Supporting Information

[Figure omitted. See PDF.]

S1 Fig. The expression of miR-34b and -34c strictly parallels.

Correlation analysis of miR-34b/c expression in tumor and non-tumor tissue was performed. A, American cohort. B, Chinese cohort. Dot plots represent miR-34b/c threshold cycle values from TaqMan qRT-PCR normalized to U66. Pearson’s correlation test was used.

https://doi.org/10.1371/journal.pone.0124899.s001

(TIFF)

S2 Fig. Correlation between miR-34 expression and survival.

Kaplan-Meier survival analysis of all stage cases in the American cohort (A) and Chinese cohort (B) stratified by median miR-34a expression and combined miR-34b/c expression. *Cases with low miR-34b and/or low miR-34c. **Cases with both high miR-34b and high miR-34c.

https://doi.org/10.1371/journal.pone.0124899.s002

(TIFF)

S3 Fig. Laser microdissection technique.

RNA samples from cancer epithelium, cancer stroma, normal adjacent epithelium and normal adjacent stroma were extracted separately using laser microdissection technique from 5 colon tumors and 5 adjacent non-tumor tissues obtained from American cohort.

https://doi.org/10.1371/journal.pone.0124899.s003

(TIFF)

S4 Fig. The miR-34 family members increase in stroma compared to epithelium.

Dot plots represent miR-34a/b/c threshold cycle values from TaqMan qRT-PCR. Horizontal bars indicate median expression value.

https://doi.org/10.1371/journal.pone.0124899.s004

(TIFF)

S5 Fig. The miR-34 family does not associated with increased stromal area.

The H&E staining on 82 of the tumor tissues from the American cohort was performed, and the area ratio of tumor and stroma (TS-ratio) was determined. Dot plots represent miR-34a/b/c threshold cycle values from TaqMan qRT-PCR and TS-ratio. Pearson’s correlation test was used.

https://doi.org/10.1371/journal.pone.0124899.s005

(TIFF)

S1 Table. Correlation analysis between Affymetrix microarray data and miR-34 family expression in 56 colon tumors.

https://doi.org/10.1371/journal.pone.0124899.s006

(XLSX)

S2 Table. List of Transcription Regulators that are potentially involved with miR-34b/c expression predicted by Ingenuity Pathway Analysis in 56 colon tumors.

https://doi.org/10.1371/journal.pone.0124899.s007

(XLSX)

S3 Table. List of Transcription Regulators that are potentially involved with miR-34a expression predicted by Ingenuity Pathway Analysis in 56 colon tumors.

https://doi.org/10.1371/journal.pone.0124899.s008

(XLSX)

Acknowledgments

We thank Shelly Hoover for technical assistance for laser microdissection (LMD), and Teruhide Ishigame for helpful support and discussion.

Author Contributions

Conceived and designed the experiments: SYL CCH. Performed the experiments: AJS HO KI KA GHN JEH. Analyzed the data: IH. Contributed reagents/materials/analysis tools: EDB. Wrote the paper: YH.

Citation: Hiyoshi Y, Schetter AJ, Okayama H, Inamura K, Anami K, Nguyen GH, et al. (2015) Increased MicroRNA-34b and -34c Predominantly Expressed in Stromal Tissues Is Associated with Poor Prognosis in Human Colon Cancer. PLoS ONE 10(4): e0124899. https://doi.org/10.1371/journal.pone.0124899

References

1. Ferlay J, Shin HR, Bray F, Forman D, Mathers C, Parkin DM. Estimates of worldwide burden of cancer in 2008: GLOBOCAN 2008. Int J Cancer. 2010;127:2893–917. pmid:21351269

2. Wightman B, Ha I, Ruvkun G. Posttranscriptional regulation of the heterochronic gene lin-14 by lin-4 mediates temporal pattern formation in C. elegans. Cell. 1993;75:855–62. pmid:8252622

3. Esquela-Kerscher A, Slack FJ. Oncomirs—microRNAs with a role in cancer. Nat Rev Cancer. 2006;6:259–69. pmid:16557279

4. Brennecke J, Hipfner DR, Stark A, Russell RB, Cohen SM. bantam encodes a developmentally regulated microRNA that controls cell proliferation and regulates the proapoptotic gene hid in Drosophila. Cell. 2003;113:25–36. pmid:12679032

5. Liu W, Mao SY, Zhu WY. Impact of tiny miRNAs on cancers. World J Gastroenterol. 2007;13:497–502. pmid:17278213

6. Merkel O, Hamacher F, Laimer D, Sifft E, Trajanoski Z, Scheideler M, et al. Identification of differential and functionally active miRNAs in both anaplastic lymphoma kinase (ALK)+ and ALK- anaplastic large-cell lymphoma. Proc Natl Acad Sci USA. 2010;107:16228–33. pmid:20805506

7. Volinia S, Calin GA, Liu CG, Ambs S, Cimmino A, Petrocca F, et al. A microRNA expression signature of human solid tumors defines cancer gene targets. Proc Natl Acad Sci USA. 2006;103:2257–61. pmid:16461460

8. Cummins JM, He Y, Leary RJ, Pagliarini R, Diaz LA Jr., Sjoblom T, et al. The colorectal microRNAome. Proc Natl Acad Sci USA. 2006;103:3687–92. pmid:16505370

9. Hermeking H. The miR-34 family in cancer and apoptosis. Cell Death Differ. 2010;17:193–9. pmid:19461653

10. Chang TC, Wentzel EA, Kent OA, Ramachandran K, Mullendore M, Lee KH, et al. Transactivation of miR-34a by p53 broadly influences gene expression and promotes apoptosis. Mol Cell. 2007;26:745–52. pmid:17540599

11. He L, He X, Lim LP, de Stanchina E, Xuan Z, Liang Y, et al. A microRNA component of the p53 tumour suppressor network. Nature. 2007;447:1130–4. pmid:17554337

12. Bommer GT, Gerin I, Feng Y, Kaczorowski AJ, Kuick R, Love RE, et al. p53-mediated activation of miRNA34 candidate tumor-suppressor genes. Curr Biol. 2007;17:1298–307. pmid:17656095

13. Siemens H, Jackstadt R, Hunten S, Kaller M, Menssen A, Gotz U, et al. miR-34 and SNAIL form a double-negative feedback loop to regulate epithelial-mesenchymal transitions. Cell cycle. 2011;10:4256–71. pmid:22134354

14. Toyota M, Suzuki H, Sasaki Y, Maruyama R, Imai K, Shinomura Y, et al. Epigenetic silencing of microRNA-34b/c and B-cell translocation gene 4 is associated with CpG island methylation in colorectal cancer. Cancer Res. 2008;68:4123–32. pmid:18519671

15. Lodygin D, Tarasov V, Epanchintsev A, Berking C, Knyazeva T, Korner H, et al. Inactivation of miR-34a by aberrant CpG methylation in multiple types of cancer. Cell cycle. 2008;7:2591–600. pmid:18719384

16. Vogt M, Munding J, Gruner M, Liffers ST, Verdoodt B, Hauk J, et al. Frequent concomitant inactivation of miR-34a and miR-34b/c by CpG methylation in colorectal, pancreatic, mammary, ovarian, urothelial, and renal cell carcinomas and soft tissue sarcomas. Virchows Arch. 2011;458:313–22. pmid:21225432

17. Lee CH, Subramanian S, Beck AH, Espinosa I, Senz J, Zhu SX, et al. MicroRNA profiling of BRCA1/2 mutation-carrying and non-mutation-carrying high-grade serous carcinomas of ovary. PloS One. 2009;4:e7314. pmid:19798417

18. Hagman Z, Larne O, Edsjo A, Bjartell A, Ehrnstrom RA, Ulmert D, et al. miR-34c is downregulated in prostate cancer and exerts tumor suppressive functions. Int J Cancer. 2010;127:2768–76. pmid:21351256

19. Nakatani F, Ferracin M, Manara MC, Ventura S, Del Monaco V, Ferrari S, et al. miR-34a predicts survival of Ewing's sarcoma patients and directly influences cell chemo-sensitivity and malignancy. J Pathol. 2012;226:796–805. pmid:21960059

20. Jamieson NB, Morran DC, Morton JP, Ali A, Dickson EJ, Carter CR, et al. MicroRNA molecular profiles associated with diagnosis, clinicopathologic criteria, and overall survival in patients with resectable pancreatic ductal adenocarcinoma. Clin Cancer Res. 2012;18:534–45. pmid:22114136

21. Cuzick J. A Wilcoxon-type test for trend. Stat Med. 1985;4:87–90. pmid:3992076

22. Costa A, Marasca R, Valentinis B, Savarino M, Faranda A, Silvestrini R, et al. p53 gene point mutations in relation to p53 nuclear protein accumulation in colorectal cancers. J Pathol. 1995;176:45–53. pmid:7616356

23. Nishida N, Nagahara M, Sato T, Mimori K, Sudo T, Tanaka F, et al. Microarray analysis of colorectal cancer stromal tissue reveals upregulation of two oncogenic miRNA clusters. Clin Cancer Res. 2012;18:3054–70. pmid:22452939

24. Tazawa H, Tsuchiya N, Izumiya M, Nakagama H. Tumor-suppressive miR-34a induces senescence-like growth arrest through modulation of the E2F pathway in human colon cancer cells. Proc Natl Acad Sci USA. 2007;104:15472–7. pmid:17875987

25. Yamakuchi M, Ferlito M, Lowenstein CJ. miR-34a repression of SIRT1 regulates apoptosis. Proc Natl Acad Sci USA. 2008;105:13421–6. pmid:18755897

26. Wu J, Wu G, Lv L, Ren YF, Zhang XJ, Xue YF, et al. MicroRNA-34a inhibits migration and invasion of colon cancer cells via targeting to Fra-1. Carcinogenesis. 2012;33:519–28. pmid:22198213

27. Roy S, Levi E, Majumdar AP, Sarkar FH. Expression of miR-34 is lost in colon cancer which can be re-expressed by a novel agent CDF. J Hematol Oncol. 2012;5:58. pmid:22992310

28. Wang M, Zhang P, Li Y, Liu G, Zhou B, Zhan L, et al. The quantitative analysis by stem-loop real-time PCR revealed the microRNA-34a, microRNA-155 and microRNA-200c overexpression in human colorectal cancer. Med Oncol. 2012;29:3113–8. pmid:22562822

29. Katada T, Ishiguro H, Kuwabara Y, Kimura M, Mitui A, Mori Y, et al. microRNA expression profile in undifferentiated gastric cancer. Int J Oncol. 2009;34:537–42. pmid:19148490

30. Juan D, Alexe G, Antes T, Liu H, Madabhushi A, Delisi C, et al. Identification of a microRNA panel for clear-cell kidney cancer. Urology. 2010;75:835–41. pmid:20035975

31. Jung M, Mollenkopf HJ, Grimm C, Wagner I, Albrecht M, Waller T, et al. MicroRNA profiling of clear cell renal cell cancer identifies a robust signature to define renal malignancy. J Cell Mol Med. 2009;13:3918–28. pmid:19228262

32. Gao W, Shen H, Liu L, Xu J, Xu J, Shu Y. MiR-21 overexpression in human primary squamous cell lung carcinoma is associated with poor patient prognosis. J Cancer Res Clin Oncol. 2011;137:557–66. pmid:20508945

33. Zhang H, Luo XQ, Zhang P, Huang LB, Zheng YS, Wu J, et al. MicroRNA patterns associated with clinical prognostic parameters and CNS relapse prediction in pediatric acute leukemia. PloS One. 2009;4:e7826. pmid:19915715

34. Svoboda M, Sana J, Redova M, Navratil J, Palacova M, Fabian P, et al. MiR-34b is associated with clinical outcome in triple-negative breast cancer patients. Diagn Pathol. 2012;7:31. pmid:22439831

35. Mantovani A, Allavena P, Sica A, Balkwill F. Cancer-related inflammation. Nature. 2008;454:436–44. pmid:18650914

36. Mathe E, Nguyen GH, Funamizu N, He P, Moake M, Croce CM, et al. Inflammation regulates microRNA expression in cooperation with p53 and nitric oxide. Int J Cancer. 2012;131:760–5. pmid:22042537

37. Bu P, Chen KY, Chen JH, Wang L, Walters J, Shin YJ, et al. A microRNA miR-34a-regulated bimodal switch targets Notch in colon cancer stem cells. Cell Stem Cell. 2013;12:602–15. pmid:23642368

38. Concepcion CP, Han YC, Mu P, Bonetti C, Yao E, D'Andrea A, et al. Intact p53-dependent responses in miR-34-deficient mice. PLoS Genet 2012;8:e1002797 pmid:22844244