Skip to main content

Main menu

  • Home
  • Current Issue
  • Archive
  • Info for
    • Authors
    • Subscribers
    • Advertisers
    • Editorial Board
  • Other Publications
    • In Vivo
    • Cancer Genomics & Proteomics
    • Cancer Diagnosis & Prognosis
  • More
    • IIAR
    • Conferences
    • 2008 Nobel Laureates
  • About Us
    • General Policy
    • Contact
  • Other Publications
    • Anticancer Research
    • In Vivo
    • Cancer Genomics & Proteomics

User menu

  • Register
  • Subscribe
  • My alerts
  • Log in
  • My Cart

Search

  • Advanced search
Anticancer Research
  • Other Publications
    • Anticancer Research
    • In Vivo
    • Cancer Genomics & Proteomics
  • Register
  • Subscribe
  • My alerts
  • Log in
  • My Cart
Anticancer Research

Advanced Search

  • Home
  • Current Issue
  • Archive
  • Info for
    • Authors
    • Subscribers
    • Advertisers
    • Editorial Board
  • Other Publications
    • In Vivo
    • Cancer Genomics & Proteomics
    • Cancer Diagnosis & Prognosis
  • More
    • IIAR
    • Conferences
    • 2008 Nobel Laureates
  • About Us
    • General Policy
    • Contact
  • Visit us on Facebook
  • Follow us on Linkedin
Research ArticleExperimental Studies

Extracellular Matrix–receptor Interaction Signaling Genes Associated with Inferior Breast Cancer Survival

MING-HSIN YEH, YAU-JIN TZENG, TING-YING FU, JUN-JIE YOU, HONG-TAI CHANG, LUO-PING GER and KUO-WANG TSAI
Anticancer Research August 2018, 38 (8) 4593-4605; DOI: https://doi.org/10.21873/anticanres.12764
MING-HSIN YEH
1Department of Surgery, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
2Department of Surgery, Chung Shan Medical University Hospital, Taichung, Taiwan, R.O.C.
3School of Medicine, Chung Shan Medical University, Taichung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
YAU-JIN TZENG
4Department of Medical Education and Research, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
5Institute of Biomedical Sciences, National Sun Yat-Sen University, Kaohsiung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
TING-YING FU
6Department of Pathology and Laboratory Medicine, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
JUN-JIE YOU
4Department of Medical Education and Research, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
HONG-TAI CHANG
1Department of Surgery, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
LUO-PING GER
4Department of Medical Education and Research, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
KUO-WANG TSAI
4Department of Medical Education and Research, Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan, R.O.C.
5Institute of Biomedical Sciences, National Sun Yat-Sen University, Kaohsiung, Taiwan, R.O.C.
7Department of Chemical Biology, National Pingtung University of Education, Pingtung, Taiwan, R.O.C.
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • For correspondence: kwtsai6733@gmail.com
  • Article
  • Figures & Data
  • Info & Metrics
  • PDF
Loading

Abstract

Background/Aim: Breast cancer is a common type of cancer in women, and metastasis frequently leads to therapy failure. Using next-generation sequencing (NGS), we aspired to identify the optimal differentially expressed genes (DEGs) for use as prognostic biomarkers for breast cancer. Materials and Methods: NGS was used to determine transcriptome profiles in breast cancer tissues and their corresponding adjacent normal tissues from three patients with breast cancer. Results: Herein, 15 DEGs (fold change >4 and <0.25) involved in extracellular matrix (ECM)–receptor interaction signaling were identified through NGS. Among them, our data indicated that high HMMR expression levels were correlated with a poor pathological stage (p<0.001) and large tumor size (p<0.001), whereas high COL6A6 and Reelin (RELN) expression levels were significantly correlated with an early pathological stage (COL6A6: p=0.003 and RELN: p<0.001). Multivariate analysis revealed that high HMMR and SDC1 expression levels were significantly correlated with poor overall survival (OS; HMMR: adjusted hazard ratio [aHR] 1.93, 95% confidence interval [CI]=1.10-3.41, p=0.023; SDC1: [aHR] 2.47, 95%CI=1.28-4.77, p=0.007) for breast cancer. Combined, the effects of HMMR and SDC1 showed a significant correlation with poor OS for patients with breast cancer (high expression for both HMMR and SDC1: [aHR] 3.29, 95%CI=1.52-7.12, p=0.003). Conclusion: These findings suggest that HMMR and SDC1 involved in the ECM–receptor interaction signaling pathway could act as effective independent prognostic biomarkers for breast ductal carcinoma.

  • Breast cancer
  • biomarkers
  • extracellular matrix

Breast cancer is the most common type of cancer in women worldwide. Although the incidence of breast cancer in some Western countries has recently fallen, its incidence and associated mortality rate continue to increase in Asia (1). Breast cancer is a heterogeneous disease with considerable variation in histological and biological features. Depending on the gene expression signature, breast cancers are classified into four subtypes, namely luminal A, luminal B, epidermal growth factor receptor 2 (HER2), and basal like (2, 3). In general, luminal A-type tumors have high estrogen receptor (ER) and progesterone receptor (PR) expression levels and no HER2 expression (ER+/PR+/HER2−), whereas luminal-B type tumors have not only ER and PR expression but also HER2 expression (ER+/PR+/HER2+). Moreover, HER2-type tumors exhibit high HER2 expression but no ER or PR expression (ER−/PR−/HER2+). By contrast, basal-like tumors do not have ER, PR, or HER2 expression (ER−/PR−/HER2−); thus, this type of breast cancer is also called triple-negative breast cancer (TNBC). Studies have revealed that patients with basal-like cancer have the least favorable prognosis and 5-year survival, whereas luminal A, luminal B, and HER2 have more favorable prognoses and a well-differentiated cellular morphology (3). The unfavorable prognosis of patients with breast cancer is strongly associated with metastatic events (4-6). Furthermore, metastasis is the major cause of cancer-related deaths, which has been a major challenge in treating breast cancer. Therefore, understanding the mechanism of metastasis and identifying favorable diagnostic and prognostic biomarkers would facilitate improving the survival rate of patients with breast cancer.

Gene expression signatures are commonly used as biomarkers for predicting breast cancer progression and prognosis (7, 8). Next-generation sequencing (NGS) is a powerful, high-throughput technology for comprehensively profiling gene expression and enables for screening of aberrant genes and determining gene signatures for breast cancer (9, 10). The Cancer Genome Atlas (TCGA) is an online database that collects substantial NGS gene transcriptome data and clinical data related to breast cancer. These clinical data and information were de-identified and released for research use. Zhong et al. identified the six critical gene signatures by analyzing TCGA database. Furthermore, six-gene expression signatures can be used as independent prognostic markers in survival prediction of patients with breast cancer (11). Using TCGA database, Ong et al. successfully identified a powerful biomarker panel for predicting the prognosis progression of patients with head and neck cancer (12). Therefore, by analyzing these big data, researchers can perform prospective clinical trials and evaluate the impact of candidate genes on breast cancer (13).

In this study, the expression profiles of breast cancer tissue and adjacent normal tissue obtained from three patients were determined by using the NGS approach. Through pathway enrichment analysis, several differentially expressed genes (DEGs) involved in the extracellular matrix (ECM)–receptor interaction signaling pathway were identified. By further analyzing TCGA database, it was determined that HMMR and SDC1 deregulation could be an effective independent prognostic biomarker for breast cancer.

Materials and Methods

Clinical samples. Clinical samples were collected from patients with breast cancer who had signed an informed consent form and agreed to undergo a surgical procedure at Kaohsiung Veterans General Hospital. This study was approved by the Institutional Review Board (IRB) of Kaohsiung Veterans General Hospital, Kaohsiung, Taiwan (IRB No. VGHKS16-CT10-08). In this study, a total of six clinical samples, comprising three invasive ductal carcinoma (IDC) and three corresponding adjacent normal tissue samples, obtained from three patients with breast cancer were subjected to NGS.

Laser capture microdissection. Fresh-frozen breast cancer specimens were embedded in ornithine carbamyl transferase (Thermo Fisher Scientific Inc., Waltham, MA, USA) stored at −80°C and were then sectioned at a thickness of 10 μm and mounted onto MMI-Membrane Slides (Leica Microsystems, Wezlar, Germany). The frozen sections were fixed with ethanol for 30 sec and then stained with hematoxylin and eosin staining before dehydration (30 sec each in 95% and 100% ethanol). After air drying, the sections were laser microdissected using an MMI Microdissection Instrument (Olympus, Tokyo, Japan). Laser capture microdissection (LCM) caps were used to collect captured cells from the slices until the caps were filled with tumor cells obtained from a single sample from each patient.

Extraction of RNA from LCM-captured cells. The total RNA of LCM-captured clinical samples was extracted using the miRNeasy Micro Kit (Qiagen, Hilden, Germany) as directed in the instruction manual. Briefly, tissue samples were homogenized in 700 μl of QIAzoI Lysis Reagent at room temperature (15°C-25°C) for 5 min and mixed with 140 μl of chloroform to extract the DNA. The upper aqueous phase was transferred to a new collection tube, and 525 μl of 100% ethanol was added and mixed thoroughly through pipetting. In total, 700 μl of sample was added into a RNeasy MinElute spin column and then centrifuged at ≥8,000 × g for 15 sec. RNA integrity and quantity were evaluated using a Nanodrop 8000 platform (Thermo Fisher Scientific Inc.).

NGS process. Three high-purity N-T paired tissues were individually separated through LCM and were transferred to an RNA extraction kit (Qiagen) for RNA extraction. The library preparation and NGS processes were performed by a biotechnology company (Genomics, Taipei, Taiwan). Briefly, total RNA was treated with DNase I, and mRNA was further isolated using Oligo(dT) primer. Subsequently, mRNA was further fragmented into small templates for cDNA synthesis. A single nucleotide A (adenine) was then added to the terminal end of the fragments and ligated with adapters. Suitable fragments were selected for polymerase chain reaction (PCR) amplification. After the library preparation was completed, the library quality and concentration were examined using an Agilent 2100 Bioanalyzer and ABI StepOnePlus Real-Time PCR System. Subsequently, the library was sequenced using an Illumina HiSeq 4000 system.

Identification of DEGs. After NGS analysis, raw sequences were obtained from Illumina GA Pipeline software CASAVA Version 1.8 and were expected to generate more than 30 million reads per sample. After low-quality data were filtered, qualified reads were analyzed using TopHat/Cufflinks to estimate gene expression, which was calculated in fragments per kilobase of transcript per million mapped reads (FPKM). For differential expression analyses, CummeRbund was used to statistically analyze the gene expression profiles. The reference genome and gene annotations were retrieved from the Ensembl database. DEGs were screened with the following thresholds: false discovery rate (FDR)<0.001 and log2 (fold-change) >2 or <−2.

Functional annotation of DEGs. DEGs (tumor versus normal >4 and <0.25) were selected from NGS data and then mapped onto the Kyoto Encyclopedia of Genes and Genomes pathways using the Database for Annotation, Visualization, and Integrated Discovery (DAVID). Subsequently, the hypergeometric test was performed to identify significantly enriched pathways.

Expression data from TCGA. From TCGA (https://cancer-genome.nih.gov) data set, 1215 pieces of transcriptome data were downloaded and the corresponding clinical data, including those related to 113 normal cases and 1,102 breast cancer cases. Among these transcriptome profiles, those of 753 cases were included in the clinical pathological features and survival analysis in this study (the expression data of patients were excluded if the patient was male or the tissue type was ductal carcinoma in situ [DCIS]).

Statistical analysis. The paired t-test, repeated-measures analysis of variance (ANOVA), and Wilcoxon matched-pairs signed-ranks test were used to compare transcriptome expression levels between the IDC tissues and the adjacent normal tissues. The Mann–Whitney U-test and Kruskal–Wallis one-way ANOVA were used for comparing transcriptome expression levels between different clinicopathologic groups in the tumor tissues. In addition, the log-rank test and Cox proportional hazards model were used to evaluate the effects of different transcriptome expression levels (high vs. low) on breast cancer patients' overall survival (OS), disease-free survival, and development of multiple primaries. The difference was considered significant at p<0.05.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table I.

Fifteen differentially expressed genes involved in extracellular matrix–receptor interaction are shown.

Results

Transcriptome profiles of breast cancer. The transcriptome profiles of three patients with breast cancer were obtained through NGS (Figure 1A). Briefly, the subtype of the three cases of breast cancer was determined to be triple-negative type, and the cell morphology was IDC. Three IDC tissues and their corresponding adjacent normal tissues were precisely collected from the three patients through LCM. After RNA extraction and library preparation, the transcriptome profiles of the six samples were determined using the Illumina HiSeq 4000 platform. More than 30 million clean reads were obtained in each library. After mapping the clean reads to the genome, the gene expression abundances were quantified and presented in FPKM. Total of 19,170 genes were identified on average from the sample ([FPKM]>1). Most of these genes were consistently expressed between the adjacent normal group and breast cancer group.

Identification of DEGs. After obtaining the transcriptome profiles, DEGs were identified with the following thresholds: FDR<0.001 and fold change (in FPKM values between tumor tissues and the corresponding adjacent normal tissues) >4 or <0.25. As shown in Figure 1A, our data revealed that 3219 genes were dysregulated (fold change >4 or <0.25) in tumor tissues compared with normal tissues obtained from patient #1; among these, 1,914 and 1,305 genes were up-regulated and down-regulated, respectively. Moreover, 5158 genes were dysregulated (fold change >4 or <0.25) in tumor tissues compared with normal tissues obtained from patient #2; among these, 3,256 and 1,902 genes were up-regulated and down-regulated, respectively (Figure 1A). Additionally, 3,824 genes were dysregulated (fold change >4 or <0.25) in tumor tissues compared with normal tissues obtained from patient #3; among these, 1,392 and 2,432 genes were up-regulated and down-regulated, respectively (Figure 1A). Combining these results revealed that 330 genes were synchronously up-regulated (fold change >4 or <0.25 and FDR<0.001) and that 401 genes were down-regulated (fold change >4 or <0.25) in breast tissues compared with the corresponding adjacent normal tissues (Figure 1B).

Analysis of putative biological function of DEGs through pathway enrichment analysis. The aforementioned data revealed that 330 and 401 genes were synchronously up-regulated and down-regulated in the three breast cancers, respectively. Subsequently, these DEGs (up-regulation: 330; down-regulation: 401) were subjected to pathway enrichment analysis using the DAVID database. As presented in Figure 1C, these DEGs were significantly enriched in 19 signaling pathways: the ECM–receptor interaction signaling pathway (p<0.001), malaria (p=0.013), staphylococcus aureus infection (p=0.002), lysosome (p=0.003), ABC transporters (p=0.003), focal adhesion (p=0.005), the PPAR signaling pathway (p=0.01), axon guidance (p=0.011), regulation of lipolysis in adipocytes (p=0.012), protein digestion and absorption (p=0.016), phagosome (p=0.019), hematopoietic cell lineage (p=0.036), biosynthesis of amino acids (p=0.047), the PI3K-Akt signaling pathway (p=0.058), proteoglycans in cancer (p=0.064), the Rap1 signaling pathway (p=0.088), the adipocytokine signaling pathway (p=0.092), biosynthesis of antibiotics (p=0.093), and the AMPK signaling pathway (p=0.095). Among them, the ECM–receptor interaction signaling pathway plays a crucial role in modulating breast cancer metastasis (14-16). Therefore, DEGs in the ECM–receptor interaction signaling pathway were prioritized for further examination in this study. As shown in Table I, eight genes (COL5A1, COL5A2, COL11A1, COMP, FN1, HMMR, SDC1, and SPP1) were significantly up-regulated in breast cancer and enriched in ECM–receptor interaction, whereas seven genes (CD36, COL6A6, ITGA7, ITGA10, LAMC3, RELN, and TNXB) were significantly down-regulated.

The expression levels of 15 candidate genes were examined in an independent cohort by analyzing TCGA database. A total of 1,215 transcriptome expression profiles and pieces of clinical pathological data were downloaded from TCGA database, comprising 113 adjacent normal tissues and 1,102 breast cancer tumors. As presented in Figure 2 and Table II, the expression levels of the eight genes, namely COL5A1, COL5A2, COL11A1, COMP, FN1, HMMR, SDC1, and SPP1, were significantly increased (fold change≥3 and p<0.001) in breast cancer tissues compared with the adjacent normal tissues. Furthermore, the expression levels of seven genes, CD36, COL6A6, ITGA7, ITGA10, LAMC3, RELN, and TNXB, were significantly decreased (fold change<0.5 and p<0.001) in breast cancer compared to adjacent normal tissues. In summary, our data revealed that the expression levels of all DEGs (15 out of 15) in breast cancer that were determined through TCGA database analysis were consistent with our NGS data.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Next-generation sequencing profiles of three N-T-paired breast cancers. (A) Flowchart of identification of differentially expressed genes through next-generation sequencing. Genes with differential expression in breast cancer compared to the adjacent normal samples were filtered with the following thresholds: fold change ≥4 or fold change <0.25 and FDR <0.001. The numbers of differentially expressed genes in the three patients with breast cancers are respectively showed in the figure panels. (B) Venn diagrams of the number of up-regulated and down-regulated gene candidates in the three N-T-paired breast cancers. (C) Protein coding genes with differential expression were subjected to pathway enrichment analysis (bottom panels).

Figure 2.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 2.

Expression levels of differentially expressed genes involved in the extracellular matrix–receptor interaction signaling pathway were examined by analyzing The Cancer Genome Atlas. (A) and (B) Expression data of eight up-regulated (COL5A1, COL5A2, COL11A1, COMP, FN1, HMMR, SDC1, SPP1) and seven down-regulated genes (CD36, COL6A6, ITGA7, ITGA10, LAMC3, RELN, TNXB) in 1,102 breast cancers compared to 113 corresponding adjacent normal tissues were obtained from The Cancer Genome Atlas. Boxplots of gene expression values (RPKM) for the up-regulated genes in breast cancer are shown.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table II.

Expression levels of differentially expressed genes analyzed in 1102 breast cancer cases (113 normal vs. 1102 tumor) using The Cancer Genome Atlas.

Evaluation of association between DEG expression and survival for breast cancer by Kaplan–Meier plotter. Overall, 15 DEGs involved in the ECM–receptor interaction signaling pathway were identified. To further evaluate the clinical impacts of these DEGs, an online survival analysis tool, the Kaplan–Meier plotter, was used to rapidly assess the effect of the 15 DEGs on breast cancer survival. This database collects the correlation between the survival curve and individual gene expression by analyzing the microarray data of 1809 patients (17). As shown in Table III, the expression levels of four genes, HMMR (hazard ratio [HR] 1.69, 95% confidence interval [CI]=1.36-2.1, p<0.001), SDC1 (HR=1.06, 95%CI=1.05-1.63, p<0.001), COL6A6 (HR=0.57, 95%CI=0.41-0.79, p<0.001), and RELN (HR=0.77, 95%CI=0.62-0.95, p=0.016), were significantly correlated with the OS curve of breast cancer (log-rank test p<0.05).

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table III.

Evaluation of impact of 15 differentially expressed genes on overall survival of breast cancer by using Kaplan–Meier plotter.

Evaluation of clinical impact of DEGs through TCGA database analysis. The clinical impacts of HMMR, SDC1, RELN, and COL6A6 were further examined by analyzing TCGA database. A total of 753 cases were used in this study to evaluate the correlation between survival and DEG expression (data on male individuals and DCIS were excluded). Using this independent cohort, we confirmed the clinical effects of these four genes in breast cancer. As presented in Table IV, high HMMR expression levels were significantly correlated with a poor American Joint Committee on Cancer (AJCC) pathological stage (p<0.001) and large tumor size (p<0.001), whereas no correlation was observed between clinical pathological features and SDC1 expression. Low COL6A6 expression levels were significantly associated with a high AJCC pathological stage (p=0.003) and large tumor size (p<0.001; Table V). High RELN expression levels were significantly correlated with the age of patients with breast cancer (<40, 40-59, vs. ≥60 years; p<0.001). Furthermore, low RELN expression levels were highly associated with a poor pathological stage (p<0.001; Table V). Notably, our data indicated that the expression levels of HMMR were significantly higher in IDC with the triple-negative subtype compared with those with the luminal A and luminal B subtypes (Table IV and Figure 3). Conversely, higher SDC1 expression levels were observed in breast cancer with Her+ compared with the luminal A, luminal B, and triple-negative subtypes (Table IV and Figure 3). In addition, the expression levels of COL6A6 and RELN did not correlate with any breast cancer subtypes (Table V and Figure 3). These results imply that the expression levels of HMMR and RELN might facilitate classifying breast cancer molecular subtypes.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table IV.

Correlation of HMMR and SDC1 expression with clinicopathological characteristics of patients with breast cancer.

HMMR and SDC1 expressions as independent survival biomarkers. The association of the expression levels of four genes, namely HMMR, SDC1, COL6A6, and RELN, with survival was also evaluated in patients with breast ductal carcinoma. Kaplan–Meier analysis revealed that high HMMR and SDC1 expression levels were significantly correlated with a poor OS curve (HMMR on OS: p=0.009 and SCD1 on OS: p=0.002; Figure 4A and B). A multivariate Cox regression model revealed a significant association between high HMMR and SDC1 expression levels with poor OS (HMMR: adjusted HR [aHR] 1.93, 95%CI=1.10-3.41, p=0.023; SDC1: aHR=2.47, 95%CI=1.28-4.77, p=0.007; Table VI).

In addition, high COL6A6 expression levels exhibited a borderline significant association with a satisfactory OS curve of breast cancer (p=0.049; Figure 4C), whereas those of RELN exhibited no significant association (p=0.219, Figure 4D). The multivariate Cox regression model revealed no significant association between COL6A6 and RELN expression and the OS of patients with breast cancer (COL6A6: aHR=0.57, 95%CI=0.29-1.11, p=0.097; RELA: aHR=0.76, 95%CI=0.42-1.40, p=0.385; Table VI). Furthermore, the effects of HMMR and SDC1 expression were combined to assess the correlation with survival in breast cancer. High HMMR and SDC1 expression levels were defined as a risk for increasing breast cancer deaths. As shown in Figure 5, high expression levels of both HMMR and SDC1 were associated with significantly shorter survival for patients with breast cancer (OS: p<0.001), even after adjustment for clinicopathologic factors (OS: aHR=3.29, 95%CI=1.52-7.12, p=0.003; Table VII).

Discussion

In the present study, 330 up-regulated and 401 down-regulated genes were identified in three TNBC samples through NGS. These DEGs were enriched in several cancer-related signaling pathways. Finally, we focused on ECM–receptor interaction signaling. The ECM comprises proteins and proteoglycans/glycosaminoglycans that provide structural support in human cells. Studies have shown that dysregulation of ECM components plays a crucial role in promoting cancer cell invasion and progression (14-16). HMMR are major components of the ECM and are involved in regulated cancer cell growth, migration, and invasion (16). In this study, our data reveal that HMMR overexpression was highly correlated with poor pathological stages and poor survival rates. Previous studies have revealed that HMMR is a multifunctional extracellular and intracellular protein and can regulate cell motility and the cell cycle (18). Stevens et al. reported that HMMR overexpression in primary lung adenocarcinoma was associated with poor prognosis (19). Akent'eva et al. showed that breast cancer cell RHAMM-selective peptides can induce apoptosis and necrosis of breast cancer cells, implying that HMMR-selective peptides have potential as a target therapy for cancer (20). Wang et al. reported that RHAMM overexpression was observed in breast cancer and that high RHAMM expression was significantly associated with poor prognostic value (21).

Figure 3.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 3.

Analysis of the expression levels of HMMR, SDC1, COL6A6, and RELN in specific subtypes of breast cancer. A total of 653 pieces of breast cancer data were downloaded from The Cancer Genome Atlas, comprising of 385 luminal A, 103 luminal B, 33 Her2 overexpression, and 132 TRBC, respectively. The expression levels of (A) HMMR, (B) SDC1, (C) COL6A6, and (D) RELN were evaluated with RPKM by analyzing The Cancer Genome Atlas.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table V.

Correlation of COL6A6 and RELN expression with clinicopathological characteristics of patients with breast cancer.

SDC1 belongs to the syndecan family and comprises cell surface heparan sulfate proteoglycans (HSPGs) (22). SDC1 has been demonstrated to play a critical role in regulating ECM architecture and facilitating breast cancer invasion and progression (23-25). Numerous studies have revealed that SDC1 was frequently down-regulated in most human cancers, including colon cancer, head and neck cancer, and gastric cancer (26-29). However, the expression levels of SDC1 have been reported to increase in breast cancer and pancreatic cancer (30, 31). A previous study reported that SDC1-overexpression was involved in promoting cell growth and cell invasion ability in breast cancer (32). Furthermore, SDC1 knockdown could interrupt colony and spheroid formation by regulating cancer stem cell–related genes in breast cancer (33). In this study, SDC1 was significantly up-regulated in breast cancer and high SDC1 expression was positively associated with poor survival of patients with breast cancer. Our results revealed that SDC1 might play an oncogenic role in regulating breast cancer progression. Previous studies have reported similar results that SDC1 was significantly up-regulated in breast cancer and that high SDC1 expression was correlated with poor prognostic value for breast cancer (31, 34, 35). Nguyen et al. reported that high SDC1 expression was significantly associated with a higher histological grade and had a negative effect on both OS and disease-free survival for patients with breast cancer (35). In addition, soluble SDC1 was reported to be detectable in sera, and SDC1 levels were increased in the sera of patients with breast cancer compared to healthy controls; moreover, high SDC1 levels were highly associated with breast tumor sizes (36). Overall, SDC1 could not only be a diagnostic biomarker, but also a prognostic biomarker for breast cancer.

Figure 4.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 4.

Analysis of the prognostic significance of HMMR, SDC1, COL6A6, and RELN in breast cancer. (A-D) A total of 743 overall survival curves were assessed depending on the HMMR, SDC1, COL6A6, and RELN expression levels in breast cancer tissues.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table VI.

Univariate and multivariate Cox regression analysis of gene expression for overall survival of patients with breast cancer patient

Figure 5.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 5.

Combination of HMMR and SDC1 expression levels to analyze their prognostic significance in breast cancer. Overall survival was compared according to the combined HMMR and SDC1 expression in breast cancer.

View this table:
  • View inline
  • View popup
  • Download powerpoint
Table VII.

Univariate and multivariate Cox's regression analysis of HMMR and SDC1 for overall survival of breast cancer patients.

In conclusion, the present study determined gene expression profiles and identified 15 DEGs involved in the ECM–receptor interaction signaling pathway. Using the TCGA database, further evaluation was executed and it was observed that two genes, HMMR and SDC1, can be used as an effective independent prognostic biomarkers in survival prediction for breast ductal carcinomas.

Acknowledgements

This work was supported by funding from the Ministry of Science Technology (MOST 106-2320-B-075B-004) and Kaohsiung Veterans General Hospital (VGHKS 105-072 and VGHKS105-135).

Footnotes

  • Conflicts of Interest

    The Authors declare that there are no conflicts of interest to disclose.

  • Received May 7, 2018.
  • Revision received June 13, 2018.
  • Accepted June 20, 2018.
  • Copyright© 2018, International Institute of Anticancer Research (Dr. George J. Delinasios), All rights reserved

References

  1. ↵
    1. Porter PL
    : Global trends in breast cancer incidence and mortality. Salud Publica Mex 51: s141-146, 2009.
    OpenUrlCrossRefPubMed
  2. ↵
    1. Alco G,
    2. Bozdogan A,
    3. Selamoglu D,
    4. Pilanci KN,
    5. Tuzlali S,
    6. Ordu C,
    7. Igdem S,
    8. Okkan S,
    9. Dincer M,
    10. Demir G,
    11. Ozmen V
    : Clinical and histopathological factors associated with Ki-67 expression in breast cancer patients. Oncol Lett 9: 1046-1054, 2015.
    OpenUrl
  3. ↵
    1. Goldhirsch A,
    2. Wood WC,
    3. Coates AS,
    4. Gelber RD,
    5. Thurlimann B,
    6. Senn HJ
    : Strategies for subtypes – dealing with the diversity of breast cancer: highlights of the St. Gallen International Expert Consensus on the Primary Therapy of Early Breast Cancer 2011. Ann Oncol 22: 1736-1747, 2011.
    OpenUrlCrossRefPubMed
  4. ↵
    1. Mehlen P,
    2. Puisieux A
    : Metastasis: a question of life or death. Nat Rev Cancer 6: 449-458, 2006.
    OpenUrlCrossRefPubMed
    1. Monteiro J,
    2. Fodde R
    : Cancer stemness and metastasis: therapeutic consequences and perspectives. Eur J Cancer 46: 1198-1203, 2010.
    OpenUrlCrossRefPubMed
  5. ↵
    1. Nguyen DX,
    2. Massague J
    : Genetic determinants of cancer metastasis. Nat Rev Genet 8: 341-352, 2007.
    OpenUrlCrossRefPubMed
  6. ↵
    1. van de Vijver MJ,
    2. He YD,
    3. van't Veer LJ,
    4. Dai H,
    5. Hart AA,
    6. Voskuil DW,
    7. Schreiber GJ,
    8. Peterse JL,
    9. Roberts C,
    10. Marton MJ,
    11. Parrish M,
    12. Atsma D,
    13. Witteveen A,
    14. Glas A,
    15. Delahaye L,
    16. van der Velde T,
    17. Bartelink H,
    18. Rodenhuis S,
    19. Rutgers ET,
    20. Friend SH,
    21. Bernards R
    : A gene-expression signature as a predictor of survival in breast cancer. N Engl J Med 347: 1999-2009, 2002.
    OpenUrlCrossRefPubMed
  7. ↵
    1. Paik S,
    2. Tang G,
    3. Shak S,
    4. Kim C,
    5. Baker J,
    6. Kim W,
    7. Cronin M,
    8. Baehner FL,
    9. Watson D,
    10. Bryant J,
    11. Costantino JP,
    12. Geyer CE Jr..,
    13. Wickerham DL,
    14. Wolmark N
    : Gene expression and benefit of chemotherapy in women with node-negative, estrogen receptor-positive breast cancer. J Clin Oncol 24: 3726-3734, 2006.
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Iyer MK,
    2. Chinnaiyan AM,
    3. Maher CA
    : ChimeraScan: a tool for identifying chimeric transcription in sequencing data. Bioinformatics 27: 2903-2904, 2011.
    OpenUrlCrossRefPubMed
  9. ↵
    1. Ma L,
    2. Nie L,
    3. Liu J,
    4. Zhang B,
    5. Song S,
    6. Sun M,
    7. Yang J,
    8. Yang Y,
    9. Fang X,
    10. Hu S,
    11. Zhao Y,
    12. Yu J
    : An RNA-seq-based gene expression profiling of radiation-induced tumorigenic mammary epithelial cells. Genomics Proteomics Bioinformatics 10: 326-335, 2012.
    OpenUrl
  10. ↵
    1. Zhong L,
    2. Lou G,
    3. Zhou X,
    4. Qin Y,
    5. Liu L,
    6. Jiang W
    : A six-long non-coding RNAs signature as a potential prognostic marker for survival prediction of ER-positive breast cancer patients. Oncotarget 8: 68861-67870, 2017.
    OpenUrl
  11. ↵
    1. Ong CJ,
    2. Shannon NB,
    3. Mueller S,
    4. Lek SM,
    5. Qiu X,
    6. Chong FT,
    7. Li K,
    8. Koh KKN,
    9. Tay GCA,
    10. Skanthakumar T,
    11. Hwang JSG,
    12. Lim TKH,
    13. Ang MK,
    14. Tan DSW,
    15. Tan NC,
    16. Tan HK,
    17. Soo KC,
    18. Iyer NG
    : A three gene immunohistochemical panel serves as an adjunct to clinical staging of patients with head and neck cancer. Oncotarget 8: 79556-79566, 2017.
    OpenUrl
  12. ↵
    1. Kim SY,
    2. Kawaguchi T,
    3. Yan L,
    4. Young J,
    5. Qi Q,
    6. Takabe K
    : Clinical relevance of microRNA expressions in breast cancer validated using the Cancer Genome Atlas (TCGA). Ann Surg Oncol 24: 2943-2949, 2017.
    OpenUrl
  13. ↵
    1. Huan JL,
    2. Gao X,
    3. Xing L,
    4. Qin XJ,
    5. Qian HX,
    6. Zhou Q,
    7. Zhu L
    : Screening for key genes associated with invasive ductal carcinoma of the breast via microarray data analysis. Genet Mol Res 13: 7919-7925, 2014.
    OpenUrl
    1. Schwertfeger KL,
    2. Cowman MK,
    3. Telmer PG,
    4. Turley EA,
    5. McCarthy JB
    : Hyaluronan, Inflammation, and Breast Cancer Progression. Front Immunol 6: 236, 2015.
    OpenUrlPubMed
  14. ↵
    1. Tolg C,
    2. McCarthy JB,
    3. Yazdani A,
    4. Turley EA
    : Hyaluronan and RHAMM in wound repair and the “cancerization” of stromal tissues. Biomed Res Int 2014: 103923, 2014.
    OpenUrl
  15. ↵
    1. Gyorffy B,
    2. Lanczky A,
    3. Eklund AC,
    4. Denkert C,
    5. Budczies J,
    6. Li Q,
    7. Szallasi Z
    : An online survival analysis tool to rapidly assess the effect of 22,277 genes on breast cancer prognosis using microarray data of 1,809 patients. Breast Cancer Res Treat 123: 725-731, 2010.
    OpenUrlCrossRefPubMed
  16. ↵
    1. Sohr S,
    2. Engeland K
    : RHAMM is differentially expressed in the cell cycle and down-regulated by the tumor suppressor p53. Cell Cycle 7: 3448-3460, 2008.
    OpenUrlCrossRefPubMed
  17. ↵
    1. Stevens LE,
    2. Cheung WKC,
    3. Adua SJ,
    4. Arnal-Estape A,
    5. Zhao M,
    6. Liu Z,
    7. Brewer K,
    8. Herbst RS,
    9. Nguyen DX
    : Extracellular matrix receptor expression in subtypes of lung adenocarcinoma potentiates outgrowth of micrometastases. Cancer Res 77: 1905-1917, 2017.
    OpenUrlAbstract/FREE Full Text
  18. ↵
    1. Akent'eva NP,
    2. Shushanov SS,
    3. Kotel'nikov AI
    : Effects of RHAMM/HMMR-Selective Peptides on Survival of Breast Cancer Cells. Bull Exp Biol Med 159: 658-661, 2015.
    OpenUrl
  19. ↵
    1. Wang C,
    2. Thor AD,
    3. Moore DH 2nd.,
    4. Zhao Y,
    5. Kerschmann R,
    6. Stern R,
    7. Watson PH,
    8. Turley EA
    : The overexpression of RHAMM, a hyaluronan-binding protein that regulates ras signaling, correlates with overexpression of mitogen-activated protein kinase and is a significant parameter in breast cancer progression. Clin Cancer Res 4: 567-576, 1998.
    OpenUrlAbstract
  20. ↵
    1. Choi Y,
    2. Chung H,
    3. Jung H,
    4. Couchman JR,
    5. Oh ES
    : Syndecans as cell surface receptors: Unique structure equates with functional diversity. Matrix Biol 30: 93-99, 2011.
    OpenUrlCrossRefPubMed
  21. ↵
    1. Yang N,
    2. Mosher R,
    3. Seo S,
    4. Beebe D,
    5. Friedl A
    : Syndecan-1 in breast cancer stroma fibroblasts regulates extracellular matrix fiber organization and carcinoma cell motility. Am J Pathol 178: 325-335, 2011.
    OpenUrlCrossRefPubMed
    1. Stanley MJ,
    2. Stanley MW,
    3. Sanderson RD,
    4. Zera R
    : Syndecan-1 expression is induced in the stroma of infiltrating breast carcinoma. Am J Clin Pathol 112: 377-383, 1999.
    OpenUrlCrossRefPubMed
  22. ↵
    1. Maeda T,
    2. Alexander CM,
    3. Friedl A
    : Induction of syndecan-1 expression in stromal fibroblasts promotes proliferation of human breast cancer cells. Cancer Res 64: 612-621, 2004.
    OpenUrlAbstract/FREE Full Text
  23. ↵
    1. Peretti T,
    2. Waisberg J,
    3. Mader AM,
    4. de Matos LL,
    5. da Costa RB,
    6. Conceicao GM,
    7. Lopes AC,
    8. Nader HB,
    9. Pinhal MA
    : Heparanase-2, syndecan-1, and extracellular matrix remodeling in colorectal carcinoma. Eur J Gastroenterol Hepatol 20: 756-765, 2008.
    OpenUrlCrossRefPubMed
    1. Inki P,
    2. Joensuu H,
    3. Grenman R,
    4. Klemi P,
    5. Jalkanen M
    : Association between syndecan-1 expression and clinical outcome in squamous cell carcinoma of the head and neck. Br J Cancer 70: 319-323, 1994.
    OpenUrlCrossRefPubMed
    1. Pap Z,
    2. Pavai Z,
    3. Denes L,
    4. Kovalszky I,
    5. Jung J
    : An immunohistochemical study of colon adenomas and carcinomas: E-cadherin, Syndecan-1, Ets-1. Pathol Oncol Res 15: 579-587, 2009.
    OpenUrlCrossRefPubMed
  24. ↵
    1. Hu XF,
    2. Yao J,
    3. Gao SG,
    4. Yang YT,
    5. Peng XQ,
    6. Feng XS
    : Midkine and syndecan1 levels correlate with the progression of malignant gastric cardiac adenocarcinoma. Mol Med Rep 10: 1409-1415, 2014.
    OpenUrl
  25. ↵
    1. Conejo JR,
    2. Kleeff J,
    3. Koliopanos A,
    4. Matsuda K,
    5. Zhu ZW,
    6. Goecke H,
    7. Bicheng N,
    8. Zimmermann A,
    9. Korc M,
    10. Friess H,
    11. Buchler MW
    : Syndecan-1 expression is up-regulated in pancreatic but not in other gastrointestinal cancers. Int J Cancer 88: 12-20, 2000.
    OpenUrlCrossRefPubMed
  26. ↵
    1. Cordone I,
    2. Masi S,
    3. Summa V,
    4. Carosi M,
    5. Vidiri A,
    6. Fabi A,
    7. Pasquale A,
    8. Conti L,
    9. Rosito I,
    10. Carapella CM,
    11. Villani V,
    12. Pace A
    : Overexpression of syndecan-1, MUC-1, and putative stem cell markers in breast cancer leptomeningeal metastasis: a cerebrospinal fluid flow cytometry study. Breast Cancer Res 19: 46, 2017.
    OpenUrl
  27. ↵
    1. Ibrahim SA,
    2. Yip GW,
    3. Stock C,
    4. Pan JW,
    5. Neubauer C,
    6. Poeter M,
    7. Pupjalis D,
    8. Koo CY,
    9. Kelsch R,
    10. Schule R,
    11. Rescher U,
    12. Kiesel L,
    13. Gotte M
    : Targeting of syndecan-1 by microRNA miR-10b promotes breast cancer cell motility and invasiveness via a Rho-GTPase- and E-cadherin-dependent mechanism. Int J Cancer 131: E884-896, 2012.
    OpenUrlCrossRefPubMed
  28. ↵
    1. Ibrahim SA,
    2. Gadalla R,
    3. El-Ghonaimy EA,
    4. Samir O,
    5. Mohamed HT,
    6. Hassan H,
    7. Greve B,
    8. El-Shinawi M,
    9. Mohamed MM,
    10. Gotte M
    : Syndecan-1 is a novel molecular marker for triple negative inflammatory breast cancer and modulates the cancer stem cell phenotype via the IL-6/STAT3, Notch and EGFR signaling pathways. Mol Cancer 16: 57, 2017.
    OpenUrl
  29. ↵
    1. Lim GH,
    2. Tan PH,
    3. Jara-Lazaro AR,
    4. Thike AA,
    5. Sim WC,
    6. Yap VB,
    7. Yip GW
    : Syndecan-1 is a potential biomarker for triple-positive breast carcinomas in Asian women with correlation to survival. Singapore Med J 55: 468-472, 2014.
    OpenUrl
  30. ↵
    1. Nguyen TL,
    2. Grizzle WE,
    3. Zhang K,
    4. Hameed O,
    5. Siegal GP,
    6. Wei S
    : Syndecan-1 overexpression is associated with nonluminal subtypes and poor prognosis in advanced breast cancer. Am J Clin Pathol 140: 468-474, 2013.
    OpenUrlCrossRef
  31. ↵
    1. Malek-Hosseini Z,
    2. Jelodar S,
    3. Talei A,
    4. Ghaderi A,
    5. Doroudchi M
    : Elevated Syndecan-1 levels in the sera of patients with breast cancer correlate with tumor size. Breast Cancer 24: 742-747, 2017.
    OpenUrl
PreviousNext
Back to top

In this issue

Anticancer Research: 38 (8)
Anticancer Research
Vol. 38, Issue 8
August 2018
  • Table of Contents
  • Table of Contents (PDF)
  • Index by author
  • Back Matter (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
Print
Download PDF
Article Alerts
Sign In to Email Alerts with your Email Address
Email Article

Thank you for your interest in spreading the word on Anticancer Research.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Extracellular Matrix–receptor Interaction Signaling Genes Associated with Inferior Breast Cancer Survival
(Your Name) has sent you a message from Anticancer Research
(Your Name) thought you would like to see the Anticancer Research web site.
CAPTCHA
This question is for testing whether or not you are a human visitor and to prevent automated spam submissions.
7 + 13 =
Solve this simple math problem and enter the result. E.g. for 1+3, enter 4.
Citation Tools
Extracellular Matrix–receptor Interaction Signaling Genes Associated with Inferior Breast Cancer Survival
MING-HSIN YEH, YAU-JIN TZENG, TING-YING FU, JUN-JIE YOU, HONG-TAI CHANG, LUO-PING GER, KUO-WANG TSAI
Anticancer Research Aug 2018, 38 (8) 4593-4605; DOI: 10.21873/anticanres.12764

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Reprints and Permissions
Share
Extracellular Matrix–receptor Interaction Signaling Genes Associated with Inferior Breast Cancer Survival
MING-HSIN YEH, YAU-JIN TZENG, TING-YING FU, JUN-JIE YOU, HONG-TAI CHANG, LUO-PING GER, KUO-WANG TSAI
Anticancer Research Aug 2018, 38 (8) 4593-4605; DOI: 10.21873/anticanres.12764
Reddit logo Twitter logo Facebook logo Mendeley logo
  • Tweet Widget
  • Facebook Like
  • Google Plus One

Jump to section

  • Article
    • Abstract
    • Materials and Methods
    • Results
    • Discussion
    • Acknowledgements
    • Footnotes
    • References
  • Figures & Data
  • Info & Metrics
  • PDF

Related Articles

  • No related articles found.
  • PubMed
  • Google Scholar

Cited By...

  • Galectin-3 Interacts with C/EBP{beta} and Upregulates Hyaluronan-Mediated Motility Receptor Expression in Gastric Cancer
  • Google Scholar

More in this TOC Section

  • Docosahexaenoic Acid Potentiates the Anticancer Effect of the Menadione/Ascorbate Redox Couple by Increasing Mitochondrial Superoxide and Accelerating ATP Depletion
  • Streptonigrin Mitigates Lung Cancer-induced Cachexia by Suppressing TCF4/TWIST1-induced PTHLH Expression
  • Atezolizumab Retains Cellular Binding to Programmed Death Ligand 1 Following Aerosolization via Mesh Nebulizer
Show more Experimental Studies

Similar Articles

Keywords

  • Breast cancer
  • biomarkers
  • extracellular matrix
Anticancer Research

© 2023 Anticancer Research

Powered by HighWire