A novel immune-related ceRNA network that predicts prognosis and immunotherapy response in lung adenocarcinoma

Wei-Jing Gong; Tao Zhou; San-Lan Wu; Yi-Fei Huang; Li-Ping Xiang; Jia-Qiang Xu; Yong Han; Yong-Ning Lv; Fang Zeng; Yu Zhang

doi:10.21037/atm-21-4151

A novel immune-related ceRNA network that predicts prognosis and immunotherapy response in lung adenocarcinoma

Ann Transl Med. 2021 Sep;9(18):1484. doi: 10.21037/atm-21-4151.

Authors

Wei-Jing Gong^{1

2}, Tao Zhou^{1

2}, San-Lan Wu^{1

2}, Yi-Fei Huang^{1

2}, Li-Ping Xiang^{1

2}, Jia-Qiang Xu^{1

2}, Yong Han^{1

2}, Yong-Ning Lv^{1

2}, Fang Zeng^{1

2}, Yu Zhang^{1

2}

Affiliations

¹ Department of Pharmacy, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China.
² Hubei Province Clinical Research Center for Precision Medicine for Critical Illness, Wuhan, China.

Abstract

Background: The tumor microenvironment plays an important role in the progression and malignancy of lung adenocarcinoma and affects the immunotherapy response. There is increasing evidence that long non-coding RNAs (lncRNAs) as competing endogenous RNAs (ceRNAs) have significant functions in the development and treatment response of various kinds of cancer. This study aimed to explore the association between immune-related lncRNA-microRNA (miRNA)-messenger RNA (mRNA)-ceRNA networks, and the prognosis of and immunotherapy response in lung adenocarcinoma.

Methods: RNA-sequencing (RNA-seq) and miRNA-seq data from The Cancer Genome Atlas (TCGA) were used to evaluate the infiltration of immune cells in lung adenocarcinoma samples by undertaking a single-sample gene set enrichment analysis (ssGSEA) to divide the cells into high and low immune cell infiltration groups. The differentially expressed mRNA (DEmRNA) was further analyzed by a weighted gene co-expression network analysis (WGCNA), search tool for recurring instances of neighboring genes (STRING), and Cytoscape to select hub genes. The ceRNA network was constructed using Cytoscape. Additionally, survival analyses were conducted to screen out prognostic candidate genes.

Results: Seven thousand five hundred and thirty-eight mRNAs, 540 lncRNAs, and 138 miRNAs were found to be differentially expressed between the high and low immune cell infiltration groups. The two DEmRNA modules most significantly associated with immune cell infiltration were further analyzed, and four clusters, including 179 DEmRNAs, were selected based on Molecular Complex Detection (MCODE) scores. The selected DEmRNAs in the four clusters were mainly enriched in pathways involved in regulating the immune response. Ultimately, a ceRNA network of SNHG6-hsa-miR-30e-5p-CYSLTR1 was identified as being associated with the prognosis of and immunotherapy response in lung adenocarcinoma.

Conclusions: The present study extends understandings of immune-related lncRNA-miRNA-mRNA-ceRNA networks and identifies novel targets and a regulatory pathway for anti-tumor immunotherapy.

Keywords: Competing endogenous RNA network (ceRNA network); immune infiltration; long non-coding RNA (lncRNA); lung adenocarcinoma; prognosis.