Identification and Validation of Novel Metastasis-Related Immune Gene Signature in Breast Cancer

Overview

Journal Breast Cancer (Dove Med Press)

Publisher Dove Medical Press

Date 2024 Apr 18

PMID 38634039

Authors

Shen Ma

Ran Hao

Yi-Wei Lu

Hui-Po Wang

Jie Hu

Yi-Xin Qi

Affiliations

Soon will be listed here.

Abstract

Background: Distant metastasis remains the leading cause of death among patients with breast cancer (BRCA). The process of cancer metastasis involves multiple mechanisms, including compromised immune system. However, not all genes involved in immune function have been comprehensively identified.

Methods: Firstly 1623 BRCA samples, including transcriptome sequencing and clinical information, were acquired from Gene Expression Omnibus (GSE102818, GSE45255, GSE86166) and The Cancer Genome Atlas-BRCA (TCGA-BRCA) dataset. Subsequently, weighted gene co-expression network analysis (WGCNA) was performed using the GSE102818 dataset to identify the most relevant module to the metastasis of BRCA. Besides, ConsensusClusterPlus was applied to divide TCGA-BRCA patients into two subgroups (G1 and G2). In the meantime, the least absolute shrinkage and selection operator (LASSO) regression analysis was used to construct a metastasis-related immune genes (MRIGs)_score to predict the metastasis and progression of cancer. Importantly, the expression of vital genes was validated through reverse transcription quantitative polymerase chain reaction (RT-qPCR) and immunohistochemistry (IHC).

Results: The expression pattern of 76 MRIGs screened by WGCNA divided TCGA-BRCA patients into two subgroups (G1 and G2), and the prognosis of G1 group was worse. Also, G1 exhibited a higher mRNA expression level based on stemness index score and Tumor Immune Dysfunction and Exclusion score. In addition, higher MRIGs_score represented the higher probability of progression in BRCA patients. It was worth mentioning that the patients in the G1 group had a high MRIGs_score than those in the G2 group. Importantly, the results of RT-qPCR and IHC demonstrated that fasciculation and elongation protein zeta 1 (FEZ1) and insulin-like growth factor 2 receptor (IGF2R) were risk factors, while interleukin (IL)-1 receptor antagonist (IL1RN) was a protective factor.

Conclusion: Our study revealed a prognostic model composed of eight immune related genes that could predict the metastasis and progression of BRCA. Higher score represented higher metastasis probability. Besides, the consistency of key genes in BRCA tissue and bioinformatics analysis results from mRNA and protein levels was verified.

References

Zhang M, Liu Z, Aoshima K, Cai W, Sun H, Xu T . CECR2 drives breast cancer metastasis by promoting NF-κB signaling and macrophage-mediated immune suppression. Sci Transl Med. 2022; 14(630):eabf5473. PMC: 9003667. DOI: 10.1126/scitranslmed.abf5473. View

Teng M, Ngiow S, Ribas A, Smyth M . Classifying Cancers Based on T-cell Infiltration and PD-L1. Cancer Res. 2015; 75(11):2139-45. PMC: 4452411. DOI: 10.1158/0008-5472.CAN-15-0255. View

Wilkerson M, Hayes D . ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics. 2010; 26(12):1572-3. PMC: 2881355. DOI: 10.1093/bioinformatics/btq170. View

Xue Y, Li M, Hu J, Song Y, Guo W, Miao C . Ca2.2-NFAT2-USP43 axis promotes invadopodia formation and breast cancer metastasis through cortactin stabilization. Cell Death Dis. 2022; 13(9):812. PMC: 9500045. DOI: 10.1038/s41419-022-05174-0. View

Ashburner M, Ball C, Blake J, Botstein D, Butler H, Cherry J . Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000; 25(1):25-9. PMC: 3037419. DOI: 10.1038/75556. View

Zhou X, Zhang Y, Wang N . Systematic identification of key extracellular proteins as the potential biomarkers in lupus nephritis. Front Immunol. 2022; 13:915784. PMC: 9366080. DOI: 10.3389/fimmu.2022.915784. View

Ngamcherdtrakul W, Yantasee W . siRNA therapeutics for breast cancer: recent efforts in targeting metastasis, drug resistance, and immune evasion. Transl Res. 2019; 214:105-120. PMC: 6848785. DOI: 10.1016/j.trsl.2019.08.005. View

Liu C, Jiang Z, Shao C, Zhang H, Yue H, Chen Z . [Preliminary study of the relationship between novel coronavirus pneumonia and liver function damage: a multicenter study]. Zhonghua Gan Zang Bing Za Zhi. 2020; 28(2):107-111. DOI: 10.3760/cma.j.issn.1007-3418.2020.02.003. View

Becht E, Giraldo N, Lacroix L, Buttard B, Elarouci N, Petitprez F . Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expression. Genome Biol. 2016; 17(1):218. PMC: 5073889. DOI: 10.1186/s13059-016-1070-5. View

10.

Coffelt S, Kersten K, Doornebal C, Weiden J, Vrijland K, Hau C . IL-17-producing γδ T cells and neutrophils conspire to promote breast cancer metastasis. Nature. 2015; 522(7556):345-348. PMC: 4475637. DOI: 10.1038/nature14282. View

11.

Chen M, Zeng G, Li J, Fu J, Long Y, Pan J . High infiltration of CD20 B lymphocytes in extranodal natural killer/T-cell lymphoma is associated with better prognosis. Br J Haematol. 2021; 191(5):e116-e120. DOI: 10.1111/bjh.17069. View

12.

Cao W, Chen H, Yu Y, Li N, Chen W . Changing profiles of cancer burden worldwide and in China: a secondary analysis of the global cancer statistics 2020. Chin Med J (Engl). 2021; 134(7):783-791. PMC: 8104205. DOI: 10.1097/CM9.0000000000001474. View

13.

Aran D, Hu Z, Butte A . xCell: digitally portraying the tissue cellular heterogeneity landscape. Genome Biol. 2017; 18(1):220. PMC: 5688663. DOI: 10.1186/s13059-017-1349-1. View

14.

Li T, Fu J, Zeng Z, Cohen D, Li J, Chen Q . TIMER2.0 for analysis of tumor-infiltrating immune cells. Nucleic Acids Res. 2020; 48(W1):W509-W514. PMC: 7319575. DOI: 10.1093/nar/gkaa407. View

15.

Racle J, de Jonge K, Baumgaertner P, Speiser D, Gfeller D . Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression data. Elife. 2017; 6. PMC: 5718706. DOI: 10.7554/eLife.26476. View

16.

Sturm G, Finotello F, Petitprez F, Zhang J, Baumbach J, Fridman W . Comprehensive evaluation of transcriptome-based cell-type quantification methods for immuno-oncology. Bioinformatics. 2019; 35(14):i436-i445. PMC: 6612828. DOI: 10.1093/bioinformatics/btz363. View

17.

Song W, Ren J, Xiang R, Kong C, Fu T . Identification of pyroptosis-related subtypes, the development of a prognosis model, and characterization of tumor microenvironment infiltration in colorectal cancer. Oncoimmunology. 2021; 10(1):1987636. PMC: 8526024. DOI: 10.1080/2162402X.2021.1987636. View

18.

Liu Y, Wang Y, Yang Y, Weng L, Wu Q, Zhang J . Emerging phagocytosis checkpoints in cancer immunotherapy. Signal Transduct Target Ther. 2023; 8(1):104. PMC: 9990587. DOI: 10.1038/s41392-023-01365-z. View

19.

Zhang H, Hughes C, Li W, He J, Surdez D, El-Naggar A . Proteomic Screens for Suppressors of Anoikis Identify IL1RAP as a Promising Surface Target in Ewing Sarcoma. Cancer Discov. 2021; 11(11):2884-2903. PMC: 8563374. DOI: 10.1158/2159-8290.CD-20-1690. View

20.

Wu T, Hu E, Xu S, Chen M, Guo P, Dai Z . clusterProfiler 4.0: A universal enrichment tool for interpreting omics data. Innovation (Camb). 2021; 2(3):100141. PMC: 8454663. DOI: 10.1016/j.xinn.2021.100141. View