Identification of significant genes and pathways associated with tenascin-C in cancer progression by bioinformatics analysis
Ilnaz Rahimmanesh1, Razieh Fatehi2, Hossein Khanahmad3
1 Applied Physiology Research Center, Cardiovascular Research Institute, Isfahan University of Medical Sciences; Department of Genetics and Molecular Biology, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran 2 Department of Genetics and Molecular Biology, School of Medicine, Isfahan University of Medical Sciences, Isfahan, Iran 3 Department of Genetics and Molecular Biology, School of Medicine, Isfahan University of Medical Sciences; Pediatric Inherited Diseases Research Center, Research Institute for Primordial Prevention of Non-Communicable Disease, Isfahan University of Medical Sciences, Isfahan, Iran
Date of Submission | 18-Aug-2020 |
Date of Acceptance | 07-Sep-2020 |
Date of Web Publication | 28-Feb-2022 |
Correspondence Address: Dr. Hossein Khanahmad Department of Genetics and Molecular Biology, School of Medicine, Isfahan University of Medical Sciences, Isfahan Iran
 Source of Support: None, Conflict of Interest: None  | Check |
DOI: 10.4103/abr.abr_201_20
Background: Tenascin-C (TNC) is a large glycoprotein of the extracellular matrix which associated with poor clinical outcomes in several malignancies. TNC over-expression is repeatedly observed in several cancer tissues and promotes several processes in tumor progression. Until quite recently, more needs to be known about the potential mechanisms of TNC as a key player in cancer progression and metastasis. Materials and Methods: In the present study, we performed a bioinformatics analysis of breast and colorectal cancer expression microarray data to survey TNC role and function with holistic view. Gene expression profiles were analyzed to identify differentially expressed genes (DEGs) between normal samples and cancer biopsy samples. The protein-protein interaction (PPI) networks of the DEGs with CluePedia plugin of Cytoscape software were constructed. Furthermore, after PPI network construction, gene-regulatory networks analysis was performed to predict long noncoding RNAs and microRNAs associated with TNC and cluster analysis was performed. Using the Clue gene ontology (GO) plugin of Cytoscape software, the GO and pathway enrichment analysis were performed. Results: PPI and DEGs-miRNA-lncRNA regulatory networks showed TNC is a significant node in a huge network, and one of the main gene with high centrality parameters. Furthermore, from the regulatory level perspective, TNC could be significantly impressed by miR-335-5p. GO analysis results showed that TNC was significantly enriched in cancer-related biological processes. Conclusions: It is important to identify the TNC underlying molecular mechanisms in cancer progression, which may be clinically useful for tumor-targeting strategies. Bioinformatics analysis provides an insight into the significant roles that TNC plays in cancer progression scenarios.
Keywords: Gene regulatory network, microarray analysis, protein interaction maps, Tenascin-C
How to cite this article: Rahimmanesh I, Fatehi R, Khanahmad H. Identification of significant genes and pathways associated with tenascin-C in cancer progression by bioinformatics analysis. Adv Biomed Res 2022;11:17 |
How to cite this URL: Rahimmanesh I, Fatehi R, Khanahmad H. Identification of significant genes and pathways associated with tenascin-C in cancer progression by bioinformatics analysis. Adv Biomed Res [serial online] 2022 [cited 2023 Jan 27];11:17. Available from: https://www.advbiores.net/text.asp?2022/11/1/17/338602 |
Introduction | |  |
The tumor microenvironment is highly heterogeneous meshwork, encompassing immune cells, fibroblasts, adipocytes, and endothelial cells, and a complex network termed the extracellular matrix (ECM).[1] The ECM's composition and organization are highly regulated and control cell behavior, but the evidence is increasing on the crucial role of the altered ECM components in providing oncogenic signals that promote cancer cell proliferation, invasion, new angiogenesis, and metastasis. Growing evidence showed that the composition and organization of ECM are underlying numerous changes in cancer when compared to normal tissue.[2],[3]
Tenascin-C (TNC) is a large hexameric glycoprotein that is located primarily in the ECM and interacts with cell surface receptors and other ECM proteins. TNC is a highly conserved ECM protein, without noticeable mutational hotspots, reflecting its importance for the maintenance of tissue homeostasis. TNC expression is tightly regulated, and TNC presents in stem cell niche regions and at sites of epithelial-mesenchymal interaction in adults and is frequently up-regulated in a variety of pathological conditions including chronic inflammation and cancer.[2],[4],[5] Its expression in several cancer tissues of breast, lung, colon, glioblastoma, liver, prostate, gastric cancer, pancreatic cancer, and gastrointestinal tract have been well documented in previous studies.[6],[7]
TNC expression in tumors is significantly higher than in normal tissue and correlated with distant metastasis and invasiveness, and tumor metastases can also express TNC.[7] Because of the abundance of TNC in tumor stroma and its association with cancer progression, invasion, angiogenesis, and metastasis, identification of the TNC-related pathways may be clinically useful for representing the cancer progression process and tumor-targeting.
Systems biology approaches have a holistic view of the cancer molecular mechanisms and in the present study, to identify TNC-related pathways in cancer, we performed a bioinformatics analysis on breast and colorectal cancer expression microarray data.
The differentially expressed genes (DEGs) in these datasets were analyzed and their interaction with TNC was investigated by protein-protein interaction (PPI) network analysis. In addition, TNC gene regulatory networks (GRN) were constructed to identify possible regulatory mechanisms at the transcriptional level. Moreover, possible functions of TNC and its related protein were evaluated by gene and pathway enrichment analyses. Therefore, we anticipated that the findings of the present study may provide further insight into TNC-related pathways in cancer pathogenesis and development at the molecular level.
Materials and Methods | |  |
Acquisition of microarray data
Gene expression profiles GSE37364, which was composed of 104 breast tumor samples and 17 normal breast tissue samples and GSE42568 included 36 samples from normal colonic biopsy and 55 samples related to different stages of colorectal cancer, were downloaded from the NCBI gene expression omnibus (GEO) database (http://www.ncbi.nlm.nih.gov/geo/). Both datasets were based on the Affymetrix GPL570 platform (Affymetrix Human Genome U133 Plus 2.0 Array). The reason for choosing these gene expression profiles was the appropriate number of cancers and control of human samples, also based on previous studies which represent the role of TNC in breast and colorectal cancer.
Data processing and DE genes determination
The quality control assessment of these datasets was assigned using prcomp function and “ggplot2” package[8] of R software.[9] The GEO2R, (http://www.ncbi.nlm.nih.gov/geo/geo2r/), the web tool of GEO, was used to identify genes that are differentially expressed across experimental conditions. The adjusted P values (adj. p) were applied to correct the false-positive results by default Benjamini-Hochberg false discovery rate method. The adjusted P ≤ 0.05 was considered as the cutoff value.
Protein-protein interaction network analysis
The PPI networks were constructed with CluePedia plugin V1.5.4[10] of Cytoscape software V3.6.2.[11] STRING v11 database (https://string-db.org) was used for retrieving interactions. The topological properties of the networks were calculated by the Network Analyzer tool based on undirected interactions.
Gene-regulatory networks and clustering analysis
The Enrichr database (http://amp.pharm.mssm. edu/Enrichr/) was used to predict long noncoding RNAs (lncRNAs) associated with DEGs. The lncHUB lncRNA Co-Expression based on Recount2 library was used. With the Network Analyst database,[12] experimentally validated miRNA-gene interactions were obtained from Tar Base and miRTarBase databases. The Cytoscape was utilized to build and visualize the GRN. Cytoscape merge tool was used for the construction of comprehensive intercoms gene regulators and genes obtained from PPI network analysis. Furthermore, the plug-in of molecular complex detection (MCODE) (Bader and Hogue, 2003) in Cytoscape software was applied to explore the significant modules in GRN networks and, with a graph-theoretic clustering algorithm, identifies densely connected regions and interconnected clusters in the networks. The advanced options set as degree cutoff = 2, K-Core = 2, and Node Score Cutoff = 0.2. Moreover, the function and pathway enrichment analysis was performed for genes in the modules.
Gene ontology and functional analysis
Function and pathway enrichment analysis was performed for DEGs in the modules using Clue gene ontology (GO) plugin V2.5.4,[13] a Cytoscape plug-in. For exploring functions and pathways related to DEGs, Kyoto Encyclopedia of Genes and Genomes (KEGG) and GO databases were used to identify their biological meaning. Enriched terms with >3 genes and adjusted P ≤ 0.05 was considered to be statistically significant.
Results | |  |
Identification of differentially expressed genes
Gene expression datasets GSE42568 and GSE37364 were downloaded from GEO datasets. DEGs were identified between normal tissue and cancerous samples based on the GEO2R analysis, using the adjusted P ≤ 0.05 and | logFC| ≥1 as the cutoff criteria.
A total of 9341 and 12240 DEGs were identified from GSE42568 and GSE37364 datasets, respectively, and Venn diagram results indicate that intersection included 6609 genes and TNC is common in both datasets [Figure 1]a and [Supplementary Table 1 [Additional file 1]] available on the website]. | Figure 1: (a) The intersection of the differentially expressed genes from GSE37364 and GSE42568 datasets. The intersection included 6609 genes. (b) The principle component analysis results indicate the acceptable quality of the GSE37364 and GSE42568 datasets
Click here to view |
In the quality assay step, the samples were segregated based on their states, indicate the acceptable quality of these datasets [Figure 1]b.
Protein-protein interaction network construction and analysis
Cytoscape was applied to evaluate the primary assay of the PPI interactions, protein interactome was investigated based on the STRING database. Activation, expression, inhibition, posttranslation modification, and binding interactions were retrieved.
Network topology was measured based on the graph theory concepts such as degree, betweenness, and closeness centrality. It is remarkable, TNC is one-fourth of the genes with high centrality. DEGs interacted directly with TNC in PPI networks are listed in [Table 1]. | Table 1: Protein-protein interaction networks analysis. Tenascin-C included the main node in protein-protein interaction networks, with high centrality parameters. GSE42568 and GSE37364 datasets, respectively
Click here to view |
Identification of Tenascin-C related regulatory network
After enrichment of genes at regulatory levels, GRNs were constructed and analyzed. Cluster analysis to determine significant modules was performed. Modules consider as protein complexes which are highly interconnected subnetworks. The results of module detection showed, In the GSE42568 dataset, TNC is clustered in MCODE Cluster 9 (score: 10.824) and in GSE37364 in MCODE Cluster 10 (score: 11.157). They are important clusters of the networks. Furthermore, the analytic results from lncRNAs and miRNAs prediction and DEGs-miRNA-lncRNA regulatory network construction indicated that TNC could be potentially modulated by hsa-mir-335-5p also by LINC01326, LINC00603, and TENM3-AS1 [Figure 2]a and [Figure 2]b. | Figure 2: Gene regulatory networks. Tenascin-C is a component of key modules from GSE42568 (a) and GSE37364 (b) datasets gene regulatory networks analysis, respectively
Click here to view |
Biological classification and pathway enrichment analysis
Aiming to evaluate the functions of identified DEGs, we uploaded DEGs clustered in modules to Clue GO plugin V2.5.4 to find significant GO categories and KEGG pathways. GO analysis and KEGG pathway analysis indicated that the TNC was significantly enriched in biological process, including cellular response to stress, focal adhesion, extracellular structure organization, PI3K-Akt signaling pathway, positive regulation of cell population proliferation, and regulation of cell adhesion [Table 2].
Discussion | |  |
Previously data demonstrated that TNC promotes several procedures in tumor progression and involved in a wide range of cellular processes including cell proliferation, migration, invasiveness, angiogenesis, and immunomodulation.[7],[14] TNC is one of the most highly expressed proteins in the tumor ECM which has poor prognostic value for patient survival in a variety of malignant tumors.
Despite advances in current molecular methods, the potential mechanism of TNC as a major player in cancer progression and metastasis was not well known.[1],[2] Therefore, it's meaningful to explore the TNC-related underlying molecular mechanisms in cancer progression procedure.
In the present study, we downloaded gene expression profiles from GEO to identify DEGs associated with breast cancer and colorectal cancer.
To explore the TNC protein interaction in these cancers, the PPI and GRN networks were constructed, and central parameters analyzed. The MCODE was applied to explore the significant modules in GRNs and TNC included in the main modules that were selected for further analysis.
Previous data reported that TNC induces epithelial-mesenchymal transition (EMT) in cancer tissues at the initial steps of cancer cell invasion, during which carcinoma cells lose their epithelial polarity and intercellular connections, which leads to cancer cell escape from the surrounding epithelium.[4],[5],[6] Data extracted from previous studies showed that binding of avb6 and avb1 integrins to TNC as their essential ligand induces EMT.[15] In addition, via bioinformatics analysis, we show some other integrin's family members are co-express with TNC, which included ITGA7, ITGB8, ITGAV, ITGBL1, ITGA5, and ITGA6. The co-expression of these integrins with TNC can involve in some biological processes such as ECM organization, extracellular structure organization, cell-matrix adhesion, blood vessel development, cell adhesion molecule binding, blood vessel morphogenesis, angiogenesis, integrin-mediated signaling pathway, and cell-cell adhesion. These pathways are significantly related to cancer metastasis, invasion, and angiogenesis.[7],[16],[17]
The versican (VCAN) a large aggregating chondroitin sulfate proteoglycan, is an important ECM component and associated with cell survival, growth, migration, invasion, angiogenesis, and metastasis in several cancers.[18] The PPI networks indicate the co-expression of TNC and VCAN both in breast and colorectal cancer.
ADAMTS, which known as adamalysins, is a matrix zinc-dependent metalloproteinases and plays an important role in the cell phenotype regulation and ECM remodeling through its role in cell adhesion and migration. ADAMTS can cleave TNC and generates TNC splice variants.[19]
As same as TNC, E-cadherin is a major contributor to cell-cell adhesion. TNC expression leads to the cytoplasmic delocalization of E-cadherin and β-catenin by Src-mediated focal adhesion kinase (FAK) signaling in combination with ανβ1 and ανβ6 integrins and promotes migratory behavior.[20]
The role of Neuropilin-1 (NRP-1), a multifunctional transmembrane protein that interacts with some of signaling receptors, in initiating tumorigenic pathways has been shown in several cancer types, including breast cancer.[21] TNC is associated with NRP-1-triggered cell migration and activation of the FAK/Akt signaling pathway. Upregulation of NRP-1 and TNC was also related to considerably increased Vimentin expression and downregulation of E-cadherin (mature form) and β-catenin that leads to EMT.[22],[23]
In addition, TNC stimulates endothelin receptor type A (EDNRA) expression, and EDNRA supports cell rounding by TNC through concomitant repression of RhoA and tumor suppressor-like molecule tropomyosin-1.[24]
TNC inhibits the interaction between fibronectin and syndecan-4, a transmembrane heparan sulfate proteoglycan, and as a result, activation of Rho family GTPases leads to the induction of actin-rich filopodia and stimulates cell invasion.[25]
The GO analysis results indicated that TNC was shown to be mostly involved in biological processes, including cellular response to stress, focal adhesion, extracellular structure organization, PI3K-Akt signaling pathway, positive regulation of cell population proliferation, and regulation of cell adhesion.
Therefore, studying these signaling pathways could facilitate the prediction of TNC-related cancer progression mechanisms.
The result from this study indicates that TNC can be regulated via some miRNA and lncRNA at the expression level. DEGs-miRNA-lncRNA regulatory network construction also showed that TNC could be significantly impressed by miR-335-5p. Previous data report that miR-335-5p playing a key role in TGF-β1-induced EMT, and its over-expression can negatively regulate cell migration and invasion in nonsmall-cell lung cancer.[26] In addition, a recent study showed that miR-335 could inhibit proliferation, invasion, and migration of the cancer cells, and arrest the cells in S phase.[27],[28] This notion is supported by a study which demonstrates that miR-335 suppresses breast cancer cell migration by negatively regulating the HGF/c-Met pathway.[29] Our bioinformatics analysis shows that miR-335 regulates TNC, and this interaction may have a key role in TNC-related EMT and TNC-induced cancer cell proliferation and invasion. Further analysis and comprehensive functional characterization of the TNC-related regulatory network can provide valuable insights into the molecular events which control TNC expression in tumor cells.
Obviously, this study has some limitations to consider. More gene expression profiles of most common cancer types must be included, and in addition, the experimental analysis must be conducted to ensure the accuracy and precision of the results obtained from this study. The evaluated level of TNC reported in previous studies in most cancer types, but the co-expression of TNC with different DEGs should be assessed.
Conclusion | |  |
In summary, based on the gene expression profile analysis, these screened genes and pathways contribute to the molecular mechanisms underlying the TNC-related cancer progression. These data may offer new clues for the identification of the potential targets, therapeutic agents, and biological mechanisms of breast and colorectal cancer after confirmation by further preclinical and clinical studies.
Financial support and sponsorship
This study was financially supported by grant No 199523 of the Biotechnology Development Council of the Islamic Republic of Iran.
Conflicts of interest
There are no conflicts of interest.
References | |  |
1. | Insua-Rodríguez J, Oskarsson T. The extracellular matrix in breast cancer. Adv Drug Deliv Rev 2016;97:41-55. |
2. | Jena MK, Janjanam J. Role of extracellular matrix in breast cancer development: A brief update. F1000Res 2018;7:274. |
3. | Walker C, Mojares E, Del Río Hernández A. Role of extracellular matrix in development and cancer progression. Int J Mol Sci 2018;19:3028. |
4. | Nagaharu K, Zhang X, Yoshida T, Katoh D, Hanamura N, Kozuka Y, et al. Tenascin C induces epithelial-mesenchymal transition-like change accompanied by SRC activation and focal adhesion kinase phosphorylation in human breast cancer cells. Am J Pathol 2011;178:754-63. |
5. | Sofat N, Robertson SD, Hermansson M, Jones J, Mitchell P, Wait R. Tenascin-C fragments are endogenous inducers of cartilage matrix degradation. Rheumatol Int 2012;32:2809-17. |
6. | Tsunoda T, Inada H, Kalembeyi I, Imanaka-Yoshida K, Sakakibara M, Okada R, et al. Involvement of large tenascin-C splice variants in breast cancer progression. Am J Pathol 2003;162:1857-67. |
7. | Yoshida T, Akatsuka T, Imanaka-Yoshida K. Tenascin-C and integrins in cancer. Cell Adh Migr 2015;9:96-104. |
8. | Wickham H. Elegant Graphics for Data Analysis (ggplot2). In: New York, NY: Springer-Verlag; 2009. |
9. | Core Team R. R: A Language and Environment for Statistical Computing. Vienna, Austria: R Foundation for Statistical Computing; 2013. |
10. | Bindea G, Galon J, Mlecnik B. Cluepedia cytoscape plugin: Pathway insights using integrated experimental and in silico data. Bioinformatics 2013;29:661-3. |
11. | Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res 2003;13:2498-504. |
12. | Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. Network analyst 3.0: A visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res 2019;47:W234-41. |
13. | Bindea G, Mlecnik B, Hackl H, Charoentong P, Tosolini M, Kirilovsky A, et al. Clue GO: A cytoscape plug-in to decipher functionally grouped gene ontology and pathway annotation networks. Bioinformatics 2009;25:1091-3. |
14. | Wawrzyniak D, Grabowska M, Głodowicz P, Kuczyński K, Kuczyńska B, Fedoruk-Wyszomirska A, et al. Down-regulation of tenascin-C inhibits breast cancer cells development by cell growth, migration, and adhesion impairment. PLoS One 2020;15:e0237889. |
15. | Katoh D, Nagaharu K, Shimojo N, Hanamura N, Yamashita M, Kozuka Y, et al. Binding of αvβ1 and αvβ6 integrins to tenascin-C induces epithelial-mesenchymal transition-like change of breast cancer cells. Oncogenesis 2013;2:e65. |
16. | Deng JL, Xu YH, Wang G. Identification of potential crucial genes and key pathways in breast cancer using bioinformatic analysis. Front Genet 2019;10:695. |
17. | Kubo Y, Ishikawa K, Mori K, Kobayashi Y, Nakama T, Arima M, et al. Periostin and tenascin-C interaction promotes angiogenesis in ischemic proliferative retinopathy. Sci Rep 2020;10:9299. |
18. | Mitsui Y, Shiina H, Kato T, Maekawa S, Hashimoto Y, Shiina M, et al. Versican promotes tumor progression, metastasis and predicts poor prognosis in renal carcinoma. Mol Cancer Res 2017;15:884-95. |
19. | Przemyslaw L, Boguslaw HA, Elzbieta S, Malgorzata SM. ADAM and ADAMTS family proteins and their role in the colorectal cancer etiopathogenesis. BMB Rep 2013;46:139-50. |
20. | Singhai R, Patil VW, Jaiswal SR, Patil SD, Tayade MB, Patil AV. E-Cadherin as a diagnostic biomarker in breast cancer. N Am J Med Sci 2011;3:227-33. |
21. | Chaudhary B, Khaled YS, Ammori BJ, Elkord E. Neuropilin 1: Function and therapeutic potential in cancer. Cancer Immunol Immunother 2014;63:81-99. |
22. | Chu W, Song X, Yang X, Ma L, Zhu J, He M, et al. Neuropilin-1 promotes epithelial-to-mesenchymal transition by stimulating nuclear factor-kappa B and is associated with poor prognosis in human oral squamous cell carcinoma. PLoS One 2014;9:e101931. |
23. | Naik A, Al-Yahyaee A, Abdullah N, Sam JE, Al-Zeheimi N, Yaish MW, et al. Neuropilin-1 promotes the oncogenic tenascin-C/integrin β3 pathway and modulates chemo resistance in breast cancer cells. BMC Cancer 2018;18:533. |
24. | Lange K, Kammerer M, Hegi ME, Grotegut S, Dittmann A, Huang W, et al. Endothelin receptor type B counteracts tenascin-C-induced endothelin receptor type A-dependent focal adhesion and actin stress fiber disorganization. Cancer Res 2007;67:6163-73. |
25. | Williams SA, Schwarzbauer JE. A shared mechanism of adhesion modulation for tenascin-C and fibulin-1. Mol Biol Cell 2009;20:1141-9. |
26. | Du W, Tang H, Lei Z, Zhu J, Zeng Y, Liu Z, et al. miR-335-5p inhibits TGF-β1-induced epithelial-mesenchymal transition in non-small cell lung cancer via ROCK1. Respir Res 2019;20:225. |
27. | Chen L, Song Y, Lu Y, Zheng W, Ma W, Zhang C. miR-335 inhibits cell proliferation, migration and invasion in HeLa cervical cancer cells. Int J Clin Exp Pathol 2016;9:10351-62. |
28. | Sandoval-Bórquez A, Polakovicova I, Carrasco-Véliz N, Lobos-González L, Riquelme I, Carrasco-Avino G, et al. MicroRNA-335-5p is a potential suppressor of metastasis and invasion in gastric cancer. Clin Epigenetics 2017;9:114. |
29. | Gao Y, Zeng F, Wu JY, Li HY, Fan JJ, Mai L, et al. MiR-335 inhibits migration of breast cancer cells through targeting oncoprotein c-Met. Tumour Biol 2015;36:2875-83. |
[Figure 1], [Figure 2]
[Table 1], [Table 2]
|