Analysis of Differential Gene Expression and Core Canonical Pathways Involved in the Epithelial to Mesenchymal Transition of Triple Negative Breast Cancer Cells by Ingenuity Pathway Analysis

Elizabeth Cagle¹, Brent Lake¹, Anasua Banerjee¹, Jazmine Cuffee¹, Narendra Banerjee¹, Darla Gilmartin¹, Makaiyah Liverman¹, Shennel Brown¹, Erik Armstrong¹, Santanu Bhattacharya^2,3, Somiranjan Ghosh⁴, Tanmoy Mandal⁴, Hirendra Banerjee^1*

¹Department of Natural, Health and Human Sciences, Elizabeth City State University Campus of The University of North Carolina, Elizabeth, NC, USA.
²Department of Biochemistry and Molecular Biology, Mayo College of Medicine and Science, Jacksonville, FL, USA.
³Department of Physiology and Biomedical Engineering, Mayo College of Medicine and Science, Jacksonville, FL, USA.
⁴Departments of Pediatrics and Child Health, College of Medicine, Howard University, Washington, DC, USA.

DOI: 10.4236/cmb.2023.132002                

Abstract

Triple Negative Breast Cancer (TNBC) is a malignant form of cancer with very high mortality and morbidity. Epithelial to Mesenchymal Transition (EMT) is the most common pathophysiological change observed in cancer cells of epithelial origin that promotes metastasis, drug resistance and cancer stem cell formation. Since the information regarding differential gene expression in TNBC cells and cell signaling events leading to EMT is limited, this investigation was done by comparing transcriptomic data generated by RNA isolation and sequencing of a EMT model TNBC cell line in comparison to regular TNBC cells. RNA sequencing and Ingenuity Pathway Software Analysis (IPA) of the transcriptomic data revealed several upregulated and downregulated gene expressions along with novel core canonical pathways including Sirtuin signaling, Oxidative Phosphorylation and Mitochondrial dysfunction events involved in EMT changes of the TNBC cells.

Keywords

Triple Negative Breast Cancer, Epithelial to Mesenchymal Transition, Core Canonical Pathways

1. Introduction

Triple-negative breast cancer (TNBC) accounts for about 10% - 15% of all breast cancers [1] . The term triple-negative breast cancer refers to the fact that the cancer cells don’t have estrogen or progesterone receptors (ER or PR) and don’t make any or too much of the protein called human epidermal growth factor receptor-2 (HER2) [1] . These cancers tend to be more common in women younger than age 40, who are Black, or who have a BRCA1 mutation [2] . TNBC differs from other types of invasive breast cancer in that it tends to grow and spread faster, has fewer treatment options, and tends to have a worse prognosis [3] [4] . Triple-negative breast cancer (TNBC) is considered an aggressive cancer because it grows quickly, is more likely to have spread at the time it’s found and is more likely to come back after treatment than other types of breast cancer [1] [3] [4] . Triple-negative breast cancer has fewer treatment options than other types of invasive breast cancer. This is because the cancer cells do not have the estrogen or progesterone receptors or enough of the HER2 protein to make hormone therapy or targeted HER2 drugs work [1] . Hormone therapy and anti-HER2 drugs are not choices for women with triple-negative breast cancer, there is limited scope of chemotherapy; the overall survival rate of TNBC currently is around 77% (American Cancer Society) [4] .

Although epithelial-to-mesenchymal transition (EMT) and mesenchymal-to-epithelial transition (MET) have been implicated in the incidence of cancer metastasis and drug resistance, their impact on cancer progression and patient survival is not fully understood. During EMT, epithelial cells lose their polarity, as well as their cell-cell adhesions, and acquire the motile and invasive characteristics of mesenchymal cells [5] . Proteins such as vimentin (VIM) intermediate filament (IF) are generally upregulated when the cell is in the mesenchymal relative to the epithelial status [6] .

The VIM RFP reporter cell line (ATCC HTB-26MET) was created using CRISPR/Cas9 gene editing in the parental MDA-MB-231 breast adenocarcinoma cell line (ATCC HTB-26). HTB-26MET harbors a C-terminal red fluorescent protein (RFP) tag on the vimentin gene. This enables the tracking of the EMT status of cells in vitro by monitoring RFP expression. The integrity of the VIM RFP knock-in has been verified at the genomic, mRNA, and protein level for sequence and expression by scientists at American Type Culture Collection (ATCC, USA).

Since EMT has been implicated for breast cancer metastasis, angiogenesis, drug resistance and eventually cancer stem cell formation [7] [8] [9] , we in this study investigated the differential gene expression and core canonical pathways involved in the EMT changes in a TNBC-EMT model Vimentin-RFP knock in cell line(ATCC HTB-26MET) in comparison to the same non EMT TNBC MDA-MB-231 breast adenocarcinoma cell line (ATCC HTB-26) by analyzing the transcriptomic data obtained by NGS RNA Sequencing and using the Ingenuity Pathway Analysis (IPA) Software as licensed by Qiagen Corporation, USA.

2. Materials and Methods

2.1. Cell Culture

MDA-MB-231-VIM-RFP (ATCC CRM-HTB-26) were obtained from the American Type Culture Collection (ATCC, Manassas, Virginia) and cultured in Eagle’s Minimum Essential Medium with the inclusion of 0.01 mg/mL of insulin and 10 µg/mL blasticidin, 10% FBS and antibiotics maintained in a 5% carbon dioxide incubator at 37˚C. The MDA-MB-231also purchased from ATCC was cultured in L-15 medium supplemented with 10% FBS and antibiotics and kept in an incubator at 37˚C.

2.2. Fluorescence Imaging

The breast cancer cell lines were split by Trypsinization and grown in 6-well cell culture plates, when confluent, the cells were photographed using an Olympus Fluorescence Microscope using the red filter.

2.3. RNA Isolation and Sequencing

All RNA isolation procedures were conducted according to the manufacturer’s protocol (Signosis LLC, Single Cell RT-PCR Assay Kit, Santa Clara, CA). Between 1000 and 10,000 cells, as confirmed using a cell counter (Denovix CellDrop Brightfield cell counter, Wilmington, DE), were isolated from cell culture and washed with 200 µL of ice cold 1X PBS. Ice-cold cell lysis buffer (50 µL) was added, and the solution was then snap-frozen at −80˚C for 5 min. Cells were incubated on ice for 10 min and centrifuged at 10,000 g for 2 min. Supernatant was transferred to a fresh nuclease-free microcentrifuge tube. DNAse I (1 µL) was added, and the sample was incubated at 37˚C in a water bath for 30 min and then inactivated at 75˚C for 10 min. Supernatant was then placed in ice and stored at −80˚C. The isolated RNA was sent to PrimBio Research Institute LLC (Exton, PA) for RNA transcriptome sequencing. Results obtained from RNA sequencing analysis yielded a fold change which was used to develop pathways using Ingenuity Pathway Analysis (IPA) Software licensed from Qiagen Corporation, USA.

2.4. Ingenuity Pathway Analysis Methodologies

IPA was used to organize and analyze data obtained from RNA sequences. RNA sequences were generated from Prim Bio Research Institute (Exton, PA) with cells cultured in this laboratory. Weight was assigned to certain gene products based on significance and fold change. Relevant pathways and import molecules emerged in the context of complex interrelated cellular processes [10] .

3. Results

A fluorescent image was taken of the MDA MB 231 breast cancer cells and Vimentin-RFP tagged knock in EMT MDA-MB-231 cells detecting the constitutively expressed Vimentin gene by red fluorescence. This image (Figure 1)

Figure 1. Fluorescent imaging of a non-transgenic breast cancer cell line (left) compared to the genetically engineered breast cancer cell (right). Red fluorescent demonstrates RFP labeled vimentin expression.

proves that the RFP tagged Vimentin gene which is upregulated in the EMT changes of breast cancer cells is correctly in frame inserted by the CRISPR-CAS9 technique in the genome of these TNBC cells and inducing EMT.

The Ingenuity Pathway Analysis software created by Qiagen was used to process the raw complete transcriptome data of the differential gene expression of the MDA-MB-231-VIM-RFP incomparison to MDA-MB-231 obtained by RNA Sequencing. A forecast model of potential upregulated and down regulated key canonical pathways that are important in cell life, cell morphology, and cellular functions of the cells was created by IPA based on the differential gene expression. Since there were several down and up regulated genes identified, we considered the top canonical pathways and taking that data ran a core analysis. The results deciphered three canonical pathways which were the most common overlapping pathway among all the different pathways detected analyzing the transcriptomic data. Figure 2 is based on a positive z-score, the orange hue indicates that the pathway is upregulated, and a negative z-score with blue color indicates that the pathway is down regulated. The gray hue indicates that the pathway is not activated, and the white color indicates that the activity is uncertain. The three core canonical pathways most frequently observed during the investigation are the Sirtuin Signaling, Oxidative Phosphorylation, and the Mitochondrial Dysfunction pathways. The Sirtuin Signaling pathway had the highest z-score, which helped to analyze the rest of the pathways.

These fundamental canonical pathways play a part in several processes and the development of diseases. There are a few routes that involve cancer genesis, proliferation, and cell survival in the setting of cancer progression, Figure 3 is a bubble chart indicating the predicted canonical pathway due to gene regulation in the comparison of EMT model MDA-MB-231 vs regular MDA breast cancer cells. The sizes of the bubbles state how many genes are associated with that specific pathway. The Sirtuin Signaling is involved in cellular immune response, cellular stress, injury, and many more cellular events as shown in Figure 4. Table 1 and Table 2 show the differentially expressed genes as identified from the transcriptomic data, whereas Table 3 and Table 4 describe the functionality of those differentially expressed genes in various disease process including cancer.

Figure 2. Based on IPA software analysis of RNA sequencing data. This figure shows the primary canonical pathways involved in EMT MDA-MB-231 vs regular breast cancer cell line. This graph displays two pathways that may be up- and down-regulated.

Figure 3. Based on the log fold change activation score of the canonical pathways of MDA-MB-231 vs. MDA-MB-231 EMT breast cancer cell line, this bubble chart displays the expected known canonical pathways gene regulation in regular MDA-MB-231 versus MDA-MB-231 EMT breast cancer cells in the context of function and disease development. The circles’ size reflects the number of genes involved in the pathway, although the color coding is the same as that of the canonical pathway bar graph.

4. Discussion

The main purpose of this investigation was to study the differential gene expression in the EMT model of MDA-MB-231-VIM-RFP comparison to regular

Figure 4. In the IPA data software the disease and function analysis of the RNA sequencing data from Primbio. Each one of these genes is known to be involved in cellular proliferation of breast cancer. These are the top three canonical pathway that are involved with cell proliferation with genes in common. The comparison of each gene and its known function is stated, in this case all the genes are decreasing cell proliferation (Table 3).

Table 1. Raw transcriptomic RNA sequencing data. The top four genes with the largest fold change (FC) upregulation in the EMT MDA-MB-231 vs. ordinary MDA breast cancer cell line are represented in this data. In breast cancer, CXCR4 exhibits the highest elevated fold change. Three out of the four identified as an oncogene which are CXCR4, SHBG, S100A1.

Table 2. RNA sequencing full transcriptome raw data. This information shows which four genes were down regulated by the greatest amount in EMT MDA-MB-231 breast cancer cells as compared to normal MDA breast cancer cells.

Table 3. Disease and Function analysis of upregulated genes associated with breast cancer proliferation in MDA-MB-232 vs regular MDA breast cancer cell line. In the comparison of Sirtuin Signaling, Oxidative Phosphorylation, and Mitochondrial Dysfunction pathway.

Table 4. Disease and Function analysis of down regulated genes associated with breast cancer proliferation in MDA-MB-232 vs regular MDA breast cancer cell line associated with the Sirtuin Signaling Pathway.

MDA-MB 231 breast cancer cell line. The IPA software played an important role in analyzing results, taking the large transcriptomic data obtained by RNA Sequencing. This detail analysis identified three main core canonical pathways, along with the differentially expressed genes and their role in cancer and other diseases. As evident from literature, these identified genes are highly associated with growth, EMT and metastasis of various cancers including breast cancer. We will outline their significant role in cancer development in the following section of the manuscript.

Regarding the upregulated genes, C-X-C chemokine receptor type 4 (CXCR4) is a G protein-coupled receptor (GPCR) that is expressed on the surface of various cell types, including immune cells, cancer cells, and nervous system cells [11] [12] . It binds to its ligand, CXCL12 (also known as stromal cell-derived factor 1 or SDF-1) and activates intracellular signaling pathways that regulate various cellular processes such as cell proliferation, survival, migration, and differentiation [13] . Some cancer therapies target CXCR4 to inhibit the growth and spread of cancer cells [14] . In cancer, CXCR4 has been found to be overexpressed in a variety of human tumors, including breast, ovarian, and lung cancer [14] . The binding of CXCL12 to CXCR4 on cancer cells promotes their growth, survival, and migration, which contributes to cancer progression and metastasis [11] [12] [13] [14] . Therefore, CXCR4 is considered a potential therapeutic target for cancer treatment. Some drugs that target CXCR4 have been developed for cancer treatment, such as plerixafor (AMD3100), which prevents the binding of CXCL12 to CXCR4 and thus inhibits the growth and spread of cancer cells [15] . These drugs have shown promising results in preclinical studies and are currently being evaluated in clinical trials for the treatment of various types of cancer [16] . Sex hormone-binding globulin (SHBG) is a protein, that binds to the sex hormones testosterone and estrogen in the bloodstream and is produced in the liver and regulates the levels of these hormones by binding to them and preventing them from interacting with target cells [17] [18] . Elevated levels of SHBG have been observed in certain types of cancer, including prostate cancer, breast cancer, and ovarian cancer [19] [20] [21] . Research suggests that SHBG may play a role in the progression of these cancers by regulating the levels of hormones that promote cell growth [19] [20] [21] .

S100A1 is a multifaceted protein that plays a role in various physiological processes, particularly in muscle physiology and has been proposed as a potential therapeutic target for various pathological conditions related to calcium homeostasis disturbances, muscle and neuronal diseases [22] [23] . S100A1 has been found to be involved in the development and progression of cancer [24] . Studies have shown that S100A1 is overexpressed in various types of cancers, including breast cancer, melanoma, osteosarcoma and prostate cancer [24] [25] [26] [27] . The overexpression of S100A1 in cancer cells has been linked to several hallmarks of cancer, such as cell proliferation, survival, invasion, and metastasis [24] [25] . In breast cancer, S100A1 has been found to be overexpressed in invasive ductal carcinomas, which are the most common type of breast cancer [28] . Studies have shown that S100A1 promotes the proliferation and invasion of breast cancer cells and is associated with a poor prognosis in breast cancer patients [28] [29] . MON2 is a protein that belongs to the family of Mon1-Ccz1 complex, which is responsible for recruiting the small GTPase Rab9 to endosomes, which is an essential step in the transport of endosomal cargo to the Golgi [30] . There is limited research on the role of MON2 in cancer, but some studies have suggested that MON2 may play a role in the development and progression of certain types of cancer. One study has shown that MON2 is overexpressed in breast cancer cells and is associated with a poor prognosis in breast cancer patients. The study suggests that MON2 may promote the proliferation and invasion of breast cancer cells by regulating endosome-to-Golgi trafficking [31] .

The analysis of the down regulated genes showed guanine nucleotide-binding protein subunit beta-2-like 1 (GNB2L1) protein also known as receptor for activated protein kinase C1 (RACK1), that belongs to the family of G protein, is involved in cell proliferation, migration and chemoresistance [32] [33] . Short chain enoyl coenzyme A hydratase 1 (ECHS1) protein is an enzyme that belongs to the fatty acid metabolic pathway [34] . This gene has been shown to be involved in colon and breast cancers [35] [36] [37] . Y-box protein (YBX1) is a transcription factor that binds to specific DNA sequences called Y boxes and regulates gene expression. It has been found to be involved in several cellular processes including growth differentiation and stress response [38] . This gene has been shown to be involved in breast, colon and lung cancer and involved in chemoresistance and formation of cancer stem cells [39] [40] . Cystatin B (CSTB) is a cysteine protease inhibitor that is poorly expressed in lung and colon cancers with poor prognosis [41] .

Silent information regulation factor 1 (sirtuin Type 1, SIRT1), as a kind of NAD+ dependent class III histone deacetylation enzyme, has been found to be involved in tumor proliferation, invasion, and metastasis. The roles of SIRT1 in breast cancer is multifaceted depending on its substrate from upstream or downstream signaling pathway, overexpression of SIRT1 significantly promoted breast cancer growth both in vitro and in vivo, whereas knockdown of SIRT1 inhibited these phenotypes. Furthermore, SIRT1 was found to interact with Akt directly, consequently promoting the activity of Akt in breast cancer cells in vitro and positively correlating with expression of Akt, P-Akt, in breast cancer tissues in vivo [42] .

Mitochondria have been implicated in cell transformation since Otto Warburg considered “respiration damage” to be a pivotal feature of cancer cells. Numerous somatic mitochondrial DNA (mtDNA) mutations have been found in various types of neoplasms, including breast cancer. Studies have shown that TNBC cells have profound metabolic alterations characterized by decreased mitochondrial respiration and increased glycolysis. Due to their impaired mitochondrial function, TNBC cells are highly sensitive to glycolytic inhibition, suggesting that such metabolic intervention may be an effective therapeutic strategy for this subtype of breast cancer cells [43] . Oxidative phosphorylation (OXPHOS) is an active metabolic pathway in many cancers. RNA from pretreatment biopsies from patients with triple-negative breast cancer (TNBC) who received neoadjuvant chemotherapy demonstrated that the top canonical pathway associated with worse outcome was higher expression of OXPHOS signature [44] . Henceforth, our EMT model TNBC transcriptomic analyses by IPA have selected the above mentioned, three very important canonical pathways involved in TNBC cell signaling.

Our transcriptomic analysis of the CRISPR-CAS9 genome edited Vimentin-RFP knock in TNBC cell line deciphered several differentially regulated genes and pathways that are involved in the EMT of these highly malignant breast cancer cells which could be of use for both diagnostic, prognostic, and therapeutic targets for future drug design and development.

We conclude that there are differentially expressed genes in leading the Triple Negative Breast Cancer cells to Epithelial Mesenchymal Transition and increased malignancy.

Acknowledgement

The authors are grateful to Dr. Abedin of PRIMBIO Corporation for RNA Sequencing and data analysis. This research was supported by NIH Grant# T34-GM100831, NSF-NOYCE Graduate student training award and a US Department of Education Graduate student training award to Elizabeth City State University Campus of The University of North Carolina.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

[1]	Dass, S.A., Tan, K.L., Selva Rajan, R., Mokhtar, N.F., Mohd Adzmi, E.R., Wan Abdul Rahman, W.F., Al-Astani Tengku Din, T.A.D. and Balakrishnan, V. (2021) Triple Negative Breast Cancer: A Review of Present and Future Diagnostic Modalities. Medicina, 57, Article No. 62. https://doi.org/10.3390/medicina57010062
[2]	Pal, T., Bonner, D., Cragun, D., Monteiro, A.N.A., Phelan, C., Servais, L., Kim, J., Narod, S.A., Akbari, M.R. and Vadaparampil, S.T. (2015) A High Frequency of BRCA Mutations in Young Black Women with Breast Cancer Residing in Florida. Cancer, 121, 4173-4180. https://doi.org/10.1002/cncr.29645
[3]	Cruz-Lozano, M., González-González, A., Marchal, J.A., Munoz-Muela, E., Molina, M.P., Cara, F.E., Brown, A.M., Garcia-Rivas, G., Hernandez-Brenes, C., Lorente, J.A., Sanchez-Rovira, P., Chang, J.C. and Granados-Principal, S. (2019) Hydroxytyrosol Inhibits Cancer Stem Cells and the Metastatic Capacity of Triple-Negative Breast Cancer Cell Lines by the Simultaneous Targeting of Epithelial-to-Mesenchymal Transition, Wnt/β-Catenin and TGFβ Signaling Pathways. European Journal of Nutrition, 58, 3207-3219. https://doi.org/10.1007/s00394-018-1864-1
[4]	Dittmer, J. (2018) Breast Cancer Stem Cells: Features, Key Drivers and Treatment Options. Seminars in Cancer Biology, 53, 59-74. https://doi.org/10.1016/j.semcancer.2018.07.007
[5]	Chung, B.-M., Rotty, J.D. and Coulombe, P.A. (2013) Networking Galore: Intermediate Filaments and Cell Migration. Current Opinion in Cell Biology, 25, 600-612. https://doi.org/10.1016/j.ceb.2013.06.008
[6]	Strouhalova, K., Přechová, M., Gandalovičová, A., Brábek, J., Gregor, M. and Rosel, D. (2020) Vimentin Intermediate Filaments as Potential Target for Cancer Treatment. Cancers, 12, Article No. 184. https://doi.org/10.3390/cancers12010184
[7]	Strietz, J., Stepputtis, S.S., Follo, M., Bronsert, P., Stickeler, E. and Maurer, J. (2021) Human Primary Breast Cancer Stem Cells Are Characterized by Epithelial-Mesen-chymal Plasticity. International Journal of Molecular Sciences, 22, Article No. 1808. https://doi.org/10.3390/ijms22041808
[8]	Villard, A.V., Genna, A., Lambert, J., Volpert, M., Noel, A., Hollier, B., Polette, M., Vanwynsberghe, A.M. and Gilles, C. (2022) Regulation of Tissue Factor by CD44 Supports Coagulant Activity in Breast Tumor Cells. Cancers, 14, Article No. 3288. https://doi.org/10.3390/cancers14133288
[9]	Vitale, D., Kumar Katakam, S., Greve, B., Jang, B., Oh, E.-S., Alaniz, L. and Götte, M. (2019) Proteoglycans and Glycosaminoglycans as Regulators of Cancer Stem Cell Function and Therapeutic Resistance. The FEBS Journal, 286, 2870-2882. https://doi.org/10.1111/febs.14967
[10]	Balkwill, F. (2004) Cancer and the Chemokine Network. Nature Reviews Cancer, 4, 540-550. https://doi.org/10.1038/nrc1388
[11]	Lagerström, M.C. and Schiöth, H.B. (2008) Structural Diversity of G Protein-Coupled Receptors and Significance for Drug Discovery. Nature Reviews Drug Discovery, 7, 339-357. https://doi.org/10.1038/nrd2518
[12]	Murdoch, C. (2000) CXCR4: Chemokine Receptor Extraordinaire. Immunological Reviews, 177, 175-184. https://doi.org/10.1034/j.1600-065X.2000.17715.x
[13]	Teicher, B.A. and Fricker, S.P. (2010) CXCL12 (SDF-1)/CXCR4 Pathway in Cancer. Clinical Cancer Research, 16, 2927-2931. https://doi.org/10.1158/1078-0432.CCR-09-2329
[14]	Wang, J., Tannous, B.A., Poznansky, M.C. and Chen, H. (2020) CXCR4 Antagonist AMD3100 (Plerixafor): From an Impurity to a Therapeutic Agent. Pharmacological Research, 159, Article ID: 105010. https://doi.org/10.1016/j.phrs.2020.105010
[15]	Chen, I.X., Chauhan, V.P., Posada, J., Ng, M.R., Wu, M.W., Adstamongkonkul, P., Huang, P., Lindeman, N., Langer, R. and Jain, R.K. (2019) Blocking CXCR4 Alleviates Desmoplasia, Increases T-Lymphocyte Infiltration, and Improves Immunotherapy in Metastatic Breast Cancer. Proceedings of the National Academy of Sciences of the United States of America, 116, 4558-4566. https://doi.org/10.1073/pnas.1815515116
[16]	Simó, R., Sáez-López, C., Barbosa-Desongles, A., Hernández, C. and Selva, D.M. (2015) Novel Insights in SHBG Regulation and Clinical Implications. Trends in Endocrinology & Metabolism, 26, 376-383. https://doi.org/10.1016/j.tem.2015.05.001
[17]	Goldštajn, M.S., Toljan, K., Grgić, F., Jurković, I. and Baldani, D.P. (2016) Sex Hormone Binding Globulin (SHBG) as a Marker of Clinical Disorders. Collegium Antropologicum, 40, 211-218.
[18]	Nounu, A., Kar, S.P., Relton, C.L. and Richmond, R.C. (2022) Sex Steroid Hormones and Risk of Breast Cancer: A Two-Sample Mendelian Randomization Study. Breast Cancer Research, 24, Article No. 66. https://doi.org/10.1186/s13058-022-01553-9
[19]	Watts, E.L., Perez-Cornago, A., Fensom, G.K., Smith-Byrne, K., Noor, U., Andrews, C.D., Gunter, M.J., Holmes, M.V., Martin, R.M., Tsilidis, K.K., Albanes, D., Barricarte, A., Bueno-de-Mesquita, B., Chen, C., Cohn, B.A., Dimou, N.L., Ferrucci, L., Flicker, L., Freedman, N.D., Giles, G.G., Giovannucci, E.L., Goodman, G.E., Haiman, C.A., Hankey, G.J., Huang, J., Huang, W.Y., Hurwitz, L.M., Kaaks, R., Knekt, P., Kubo, T., Langseth, H., Laughlin, G., Le Marchand, L., Luostarinen, T., MacInnis, R.J., Mäenpää, H.O., Männistö, S., Metter, E.J., Mikami, K., Mucci, L.A., Olsen, A.W., Ozasa, K., Palli, D., Penney, K.L., Platz, E.A., Rissanen, H., Sawada, N., Schenk, J.M., Stattin, P., Tamakoshi, A., Thysell, E., Tsai, C.J., Tsugane, S., Vatten, L., Weiderpass, E., Weinstein, S.J., Wilkens, L.R., Yeap, B.B., Allen, N.E., Key, T.J. and Travis, R.C. (2022) Circulating Free Testosterone and Risk of Aggressive Prostate Cancer: Prospective and Mendelian Randomisation Analyses in International Consortia. International Journal of Cancer, 151, 1033-1046. https://doi.org/10.1002/ijc.34116
[20]	Varsavsky, M., Reyes-García, R., García-Martín, A., González-Ramírez, A.R., Avilés-Perez, M.D. and Muñoz-Torres, M. (2013) SHBG Levels Are Associated with Bone Loss and Vertebral Fractures in Patients with Prostate Cancer. Osteoporosis International, 24, 713-719. https://doi.org/10.1007/s00198-012-2012-z
[21]	Duarte-Costa, S., Castro-Ferreira, R., Neves, J.S. and Leite-Moreira, A.F. (2014) S100A1: A Major Player in Cardiovascular Performance. Physiological Research, 63, 669-681. https://doi.org/10.33549/physiolres.932712
[22]	Diaz-Romero, J. and Nesic, D. (2017) S100A1 and S100B: Calcium Sensors at the Cross-Roads of Multiple Chondrogenic Pathways. Journal of Cellular Physiology, 232, 1979-1987. https://doi.org/10.1002/jcp.25720
[23]	Wang, G., Li, H.-N., Cui, X.-Q., Xu, T., Dong, M.-L., Li, S.-Y. and Li, X.-R. (2021) S100A1 Is a Potential Biomarker for Papillary Thyroid Carcinoma Diagnosis and Prognosis. Journal of Cancer, 125, 760-5771. https://doi.org/10.7150/jca.51855
[24]	Xing, W., Zhao, Y., Lin, L., Zhao, Z., Yang, M., Wang, N., Cui, S., Bai, R. and Zhao, A. (2022) Overexpression of S100A1 in Osteosarcoma Inhibits Tumor Proliferation and Progression. Iranian Journal of Public Health, 51, 2773-2782. https://doi.org/10.18502/ijph.v51i12.11468
[25]	Young, B.D., Yu, W., Rodríguez, D.J.V., Varney, K.M., MacKerell Jr., A.D. and Weber, D.J. (2021) Specificity of Molecular Fragments Binding to S100B versus S100A1 as Identified by NMR and Site Identification by Ligand Competitive Saturation (SILCS). Molecules, 26, Article No. 381. https://doi.org/10.3390/molecules26020381
[26]	Eryilmaz, U., Demirci, B., Aksun, S., Boyacioglu, M., Akgullu, C., Ilgenli, T.F., Yalinkilinc, H.S. and Bilgen, M. (2015) S100A1 as a Potential Diagnostic Biomarker for Assessing Cardiotoxicity and Implications for the Chemotherapy of Certain Cancers. PLOS ONE, 10, e0145418. https://doi.org/10.1371/journal.pone.0145418
[27]	Zhang, S., Wang, Z., Liu, W., Lei, R., Shan, J., Li, L. and Wang, X. (2017) Distinct Prognostic Values of S100 mRNA Expression in Breast Cancer. Scientific Reports, 7, Article No. 39786. https://doi.org/10.1038/srep39786
[28]	Hatoum, D., Yagoub, D., Ahadi, A., Nassif, N.T. and McGowan, E.M. (2017) Annexin/S100A Protein Family Regulation through p14ARF-p53 Activation: A Role in Cell Survival and Predicting Treatment Outcomes in Breast Cancer. PLOS ONE, 12, e0169925. https://doi.org/10.1371/journal.pone.0169925
[29]	Zhao, S.-B., Dean, N., Gao, X.-D. and Fujita, M. (2020) MON2 Guides Wntless Transport to the Golgi through Recycling Endosomes. Cell Structure and Function, 45, 77-92. https://doi.org/10.1247/csf.20012
[30]	Li, Y., Yang, C., Liu, Z., Du, S., Can, S., Zhang, H., Zhang, L., Huang, X., Xiao, Z., Li, X., Fang, J., Qin, W., Sun, C., Wang, C., Chen, J. and Chen, H. (2022) Integrative Analysis of CRISPR Screening Data Uncovers New Opportunities for Optimizing Cancer Immunotherapy. Molecular Cancer, 21, Article No. 2. https://doi.org/10.1186/s12943-021-01462-z
[31]	Ron, D., Chen, C.H., Caldwell, J., Jamieson, L., Orr, E. and Mochly-Rosen, D. (1994) Cloning of an Intracellular Receptor for Protein Kinase C: A Homolog of the Beta Subunit of G Proteins. Proceedings of the National Academy of Sciences of the United States of America, 91, 839-843. https://doi.org/10.1073/pnas.91.3.839
[32]	Berns, H., Humar, R., Hengerer, B., Kiefer, F.N. and Battegay, E.J. (2000) RACK1 Is up-Regulated in Angiogenesis and Human Carcinomas. The FASEB Journal, 14, 2549-2558. https://doi.org/10.1096/fj.99-1038com
[33]	Janßen, U., Davis, E.M., Le Beau, M.M. and Stoffel, W. (1997) Human Mitochondrial Enoyl-CoA Hydratase Gene (ECHS1): Structural Organization and Assignment to Chromosome 10q26.2-q26.3. Genomics, 40, 470-475. https://doi.org/10.1006/geno.1996.4597
[34]	Liu, X., Feng, R. and Du, L. (2010) The Role of Enoyl-CoA Hydratase Short Chain 1 and Peroxiredoxin 3 in PP2-Induced Apoptosis in Human Breast Cancer MCF-7 Cells. FEBS Letters, 584, 3185-3192. https://doi.org/10.1016/j.febslet.2010.06.002
[35]	Zhao, Q.-M., Kuang, F., Wu, H. and Zhang, Y.-H. (2015) Attenuation of Enoyl Coenzyme A Hydratase 1 Expression in Colorectal Cancer Cells Using Small Interfering RNA Inhibits Cell Proliferation and Migration. Molecular Medicine Reports, 12, 470-474. https://doi.org/10.3892/mmr.2015.3418
[36]	Gong, X., Zhu, Y., Dong, J., Chen, J., You, J., Zheng, Q., Rao, Z., Mao, Q. and Jiang, J. (2013) Small Hepatitis B Surface Antigen Interacts with and Modulates Enoyl-Coenzyme A Hydratase Expression in Hepatoma Cells. Archives of Virology, 158, 1065-1070. https://doi.org/10.1007/s00705-012-1581-7
[37]	Uchiumi, T., Fotovati, A., Sasaguri, T., Shibahara, K., Shimada, T., Fukuda, T., Nakamura, T., Izumi, H., Tsuzuki, T., Kuwano, M. and Kohno, K. (2006) YB-1 Is Important for an Early Stage Embryonic Development: Neural Tube Formation and Cell Proliferation. Journal of Biological Chemistry, 281, 40440-40449. https://doi.org/10.1074/jbc.M605948200
[38]	Jürchott, K., Kuban, R.-J., Krech, T., Blüthgen, N., Stein, U., Walther, W., Friese, C., Kiełbasa, S.M., Ungethüm, U., Lund, P., Knösel, T., Kemmner, W., Morkel, M., Fritzmann, J., Schlag, P.M., Birchmeier, W., Krueger, T., Sperling, S., Sers, C., Royer, H.D., Herzel, H. and Schäfer, R. (2010) Identification of Y-Box Binding Protein 1 As a Core Regulator of MEK/ERK Pathway-Dependent Gene Signatures in Colorectal Cancer Cells. PLOS Genetics, 6, e1001231. https://doi.org/10.1371/journal.pgen.1001231
[39]	Matsumoto, K. and Bay, B.H. (2005) Significance of the Y-Box Proteins in Human Cancers. Journal of Molecular and Genetic Medicine, 1, 11-17. https://doi.org/10.4172/1747-0862.1000005
[40]	Ma, Y., Chen, Y. and Petersen, I. (2017) Expression and Epigenetic Regulation of Cystatin B in Lung Cancer and Colorectal Cancer. Pathology—Research and Practice, 213, 1568-1574. https://doi.org/10.1016/j.prp.2017.06.007
[41]	Jin, X., Wei, Y., Xu, F., Zhao, M., Dai, K., Shen, R., Yang, S. and Zhang, N. (2018) SIRT1 Promotes Formation of Breast Cancer through Modulating Akt Activity. Journal of Cancer, 9, 2012-2023. https://doi.org/10.7150/jca.24275
[42]	Pelicano, H., Zhang, W., Liu, J., Hammoudi, N., Dai, J., Xu, R.-H., Pusztai, L. and Huang, P. (2014) Mitochondrial Dysfunction in Some Triple-Negative Breast Cancer Cell Lines: Role of mTOR Pathway and Therapeutic Potential. Breast Cancer Research, 16, Article No. 434. https://doi.org/10.1186/s13058-014-0434-6
[43]	Evans, K.W., Yuca, E., Scott, S.S., Zhao, M., Paez Arango, N., Cruz Pico, C.X., Saridogan, T., Shariati, M., Class, C.A., Bristow, C.A., Vellano, C.P., Zheng, X., Gonzalez-Angulo, A.M., Su, X., Tapia, C., Chen, K., Akcakanat, A., Lim, B., Tripathy, D., Yap, T.A., Francesco, M.E.D., Draetta, G.F., Jones, P., Heffernan, T.P., Marszalek, J.R. and Meric-Bernstam, F. (2021) Oxidative Phosphorylation Is a Metabolic Vulnerability in Chemotherapy-Resistant Triple-Negative Breast Cancer. Cancer Research, 81, 5572-5581. https://doi.org/10.1158/0008-5472.CAN-20-3242
[44]	Krämer, A., Green, J., Pollard Jr., J. and Tugendreich, S. (2014) Causal Analysis Approaches in Ingenuity Pathway Analysis. Bioinformatics, 30, 523-530. https://doi.org/10.1093/bioinformatics/btt703