- 1Department of Breast Surgery, Guangxi Medical University Cancer Hospital, Nanning, China
- 2Fujian Maternity and Child Health Hospital, College of Clinical Medicine for Obstetrics & Gynecology and Pediatrics, Fujian Medical University, Fujian, China
- 3Department of Breast Surgery, Harbin Medical University Cancer Hospital, Harbin, China
- 4Department of Pathology, Dunedin School of Medicine, University of Otago, Dunedin, New Zealand
Background: Breast cancer (BC) is one of the most common malignancies in women, and long non-coding RNAs (lncRNAs) are key regulators of its development. T cells can recognize and kill cancer cells, and CD4+ T conventional (Tconv) cells are the main orchestrators of cancer immune function. However, research on CD4+ Tconv-related lncRNAs (CD4TLAs) prognostic signature in patients with BC is still lacking.
Method: A TCGA database and a GEO database were used to collect the BC patients. Through LASSO Cox regression analysis CD4TLAs-related prognostic models were further constructed, and risk scores (RS) were generated and developed a nomogram based on CD4TLAs. The accuracy of this model was validated in randomized cohorts and different clinical subgroups. Gene set enrichment analysis (GSEA) was used to explore potential signature-based functions. The role of RS has been further explored in the tumor microenvironment (TME), immunotherapy, and chemotherapy.
Result: A prognostic model based on 16 CD4TLAs was identified. High-RS was significantly associated with a poorer prognosis. RS was shown to be an independent prognostic indicator in BC patients. The low-RS group had a significant expression of immune infiltrating cells and significantly enriched immune-related functional pathways. In addition, the results of immunotherapy prediction indicated that patients with low-RS were more sensitive to immunotherapy.
Conclusions: Our signature has potential predictive value for BC prognosis and immunotherapy response. The findings of this work have greatly increased our understanding of CD4TLA in BC.
Background
Breast cancer (BC) is the leading cause of cancer-related deaths and one of the most common cancer in women around the world (1, 2). In 2020, the International Agency for Research on Cancer of the World Health Organization announced that new breast cancer diagnoses surpassed lung cancer diagnoses for the first time. Research on the prevention and treatment of breast cancer remains a crucial concern around the globe. BC is currently treated primarily by surgery and or by combining multiple adjuvant therapies, such as chemotherapy, immunotherapy, targeted therapy, and radiation therapy (3–5). Although early detection and treatment of BC has improved, patient outcomes still remain unfavorable. Therefore, there is a need to intensify the research for novel high-efficiency molecular biomarkers for the treatment of BC patients, and the development of new prognostic signatures must be intensified.
Long non-coding RNAs (lncRNAs) are those with over 200 nucleotides in length (6, 7). In addition to coding for proteins, they participate directly or indirectly in the regulation of numerous life activities. Recent studies suggest that lncRNAs can play a role in innate and adaptive immunity in cancer by mediating the functional status of immune cells, pathways, and genes (8, 9). Findings of Dongdi et al. revealed that lncRNA HOTAIR could affect cell growth, migration, invasion, and apoptosis in BC through the miR-20a-5p/HMGA2 axis (10). LncRNA PCAT7 can promote the malignant progression of BC by regulating the ErbB/PI3K/Akt pathway (11). The lncRNA HAL plays an influential role in regulating the BC microenvironment as well as promoting BC stemness (12). These studies suggest that lncRNAs can be valuable biomarkers and therapeutic targets for treating BC. Recently, progress in immunology have led to the identification of new mechanisms of immune-related disease. Cellular immune responses in tumors are often attributed to CD8+ T cells. However, CD4+ T cells also play a key role. CD4+ T cells can recognize tumor antigens, supporting the initiation and expansion of CD8+ T cells in lymphoid tissues (13–15). Matthew et al. showed that CD4+ T cells at the tumor site and throughout the body are a key component of immune-mediated tumor rejection (16). As an indispensable part of the immune composition, CD4+ T conventional (Tconv) cells play a crucial role in BC. Nevertheless, the application of CD4+ Tconv-related lncRNA (CD4TLA) combination in prognosis prediction of BC patients has not been established.
Researchers have been able to study BC in greater depth thanks to the recent development of single-cell sequencing technology. The leading goal of this study was to construct a CD4TLA signature that could be used to predict the prognosis and immunotherapy of BC patients, and possibly insight into clinical treatment options. Based on single-cell data, we screened CD4+ Tconv-related genes in BC and, using bioinformatics, constructed a prognostic signature based on CD4TLAs for BC. This study examined the clinical characteristics, prognostic value, and tumor microenvironment (TME) of the obtained risk score (RS) in BC. Furthermore, we tested the efficacy of RS in BC immunotherapy and chemotherapy.
Materials and Methods
Breast Cancer Data Source and Preprocessing
The gene expression data and clinical annotations for BC samples (TCGA-BRCA, n=1109) were obtained from The Cancer Genome Atlas (TCGA, https://portal.gdc.cancer.gov/). Patients without survival information were excluded. Fragments per kilobase million (FPKM) were transformed into transcripts per million (TPM) values in the TCGA-BRCA cohort. Data on BC single cells RNA-seq (GSE110686 and GSE11472) were derived from GEO databases (https://www.ncbi.nlm.nih.gov/geo/) (17, 18).
Identifying CD4+ Tconv-Related Differential Expressed lncRNAs Breast Cancer
Tumor Immune Single-cell Hub (TISCH) is a scRNA-seq database which focus on TME (19). To identify CD4TLAs, we first extracted CD4+ Tconv-related genes from BC single-cell data through the TISCH database. The correlation between lncRNAs and CD4+ Tconv-related genes was calculated using the Pearson correlation. The lncRNA with an absolute |correlation coefficient| >0.5 and a P < 0.001 was considered as a CD4TLA. Through the R package “limma”, CD4TLAs that are significantly different (P < 0.05) between normal tissue and tumor tissue were identified (20, 21).
Construction of a CD4+ Tconv-Related lncRNA Risk Scores Model
To screen lncRNAs associated with OS in BC patients, we performed Univariate Cox regression analysis on the significantly different CD4TLAs. These were then subjected to LASSO regression analysis to identify more meaningful prognostic lncRNAs. Finally, 16 CD4TLAs significantly associated with OS were identified, and the risk score (RS) of each patient was calculated based on the expression level of CD4TLAs and the Cox regression coefficient. Patients were divided into high- and low-risk subgroups by median risk score. The ROC curve was used to check the accuracy of the model. In addition, we validated the model by randomly assigning BC patients to validation group 1 and validation group 2. Univariate and multivariate Cox regression analyses examined whether RS could be an independent prognostic factor in BC patients.
Functional Enrichment Analysis
The Metascape database is a resource for gene annotation and analysis (22). In this study, the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses of CD4+ Tconv-related genes were conducted using the Metascape database. We performed GSVA enrichment analysis using “GSVA” R package in order to explore the differences in biological processes between high- and low-risk groups. The R package “cluster profile” was used for functional annotation and the gene set file (c2.cp.kegg.v7.2.symbols.gmt) was obtained from the MSigDB database (https://www.gsea-msigdb.org).
Tumor Microenvironment Analysis
The “estimate” package was used to estimate the composition of the immune stroma in the BC tumor microenvironment (TME) and to derive immune scores, stroma scores, and estimate scores. The ssGSEA algorithm was used to quantify differences in immune cell infiltration subsets and immunological functional enrichment between high- and low-RS subgroups.
Evaluation of the Immunotherapy Response and Chemotherapy Based on RS
To further evaluate the potential application of RS in clinical treatment, we predicted the response of high- and low-RS subgroups to immunotherapy by the TIDE algorithm (23). Immunotherapy data from patients were obtained from the Cancer Immuneome Database (TCIA), and we further assessed Immunophenoscore (IPS) in high- and low-RS subgroups to predict patient response to immunotherapy (24). In addition, in order to explore the sensitivity of high- and low-RS patients to different chemotherapeutic drugs, we predicted the IC50 of chemotherapeutic drugs by the “pRRophetic” package.
Statistical Analysis
All statistical analyses were carried out utilizing the R version 4.1.2 and SPSS 23.0. The student’s t-test (unpaired, two-tailed) was used to evaluate the differences between the two independent groups. One-way analysis of variance (ANOVA) and Kruskal–Wallis test were used as parametric and non-parametric methods, respectively, for data from more than two groups. The R packages “survival” and “survminer” were used for survival analysis. Volcano and heatmaps were drawn by the “ggplots” package. P < 0.05 was deemed to be significant.
Results
Single-Cell Data Analysis of BC Patients
The workflow of this study is shown in Figure 1. Based on the TISCH database, we analyzed the GSE110686 and GSE114727 datasets, and compared and visualized the major cell types in their microenvironments (Figures 2A, 3A). We found that CD4+ Tconv was the largest subpopulation of cells (Figures 2B, 3B). Furthermore, we analyzed functional differences between different cell types by GSEA enrichment (Figures 2C, D, 3C, D). Findings revealed that CD4+ Tconv was significantly enriched in carbohydrate metabolism, signal transduction, translation, signaling molecules and interaction, endocrine and metabolic disease, immune disease, cardiovascular disease, immune system, and cell motility-related pathways. These results suggest that CD4+ Tconv has an important role in BC and deserves further study.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g002.jpg)
Figure 2 Breast cancer single-cell data analysis based on the GSE110686 dataset. (A) UMAP plots with cells colored by cell type are displayed. (B) The pie plot shows the cell number distribution of each cell type. Heatmap showing enriched up- (C) or down-regulated (D) pathways identified based on differential genes in each cell type.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g003.jpg)
Figure 3 Breast cancer single-cell data analysis based on the GSE114727 dataset. (A) UMAP plots with cells colored by cell type are displayed. (B) The pie plot shows the cell number distribution of each cell type. Heatmap showing enriched up- (C) or down-regulated (D) pathways identified based on differential genes in each cell type.
Identification and Screening of CD4+ Tconv -Related Genes and lncRNAs in BC
The CD4+ Tconv-related genes were first extracted from the GSE110686 and GSE114727 datasets and then intersected, yielding 192 CD4+ Tconv-related genes (Supplementary Figure 1; Supplementary Table 1). The results of functional enrichment showed that they were mainly enriched in immune-related pathways, such as regulation of T cell activation, Cytokine Signaling in the Immune system, negative regulation of immune system process, regulation of immune effector process (Figure 4A). To further explore the relationship between enriched pathways, we present it as a network graph, where terms with similarity > 0.3 are connected by edges (Figure 4B). In addition, we further performed protein-protein interaction enrichment analysis, which is mainly related to the regulation of cytokine signaling in the Immune system, cell activation, regulation of T cell activation (Figure 4C). The MCODE network is shown in Figure 4D. A CD4+ Tconv-related gene lncRNA (CD4TLA) coexpression network was constructed to identify CD4TLAs (Figure 5A). They were then subjected to differential analysis, and finally, 244 CD4TLAs that were significantly different between tumor and normal groups were identified (Supplementary Figures 2A, B; Supplementary Table 2).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g004.jpg)
Figure 4 Functional enrichment of CD4+ Tconv-related genes and visualization of interactome analysis results. (A) Metascape enrichment analysis for the CD4+ Tconv-related genes. (B) Metascape enrichment network visualization showing the intra-cluster and inter-cluster similarities of enriched terms. (C) Metascape visualization of the interactome network formed by CD4+ Tconv-related genes candidates, where the MCODE complexes are colored according to their identities. (D) Seven MCODE complexes automatically identified in Metascape, colored by their identities. Their functional labels are generated based on the top-three functional enriched terms.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g005.jpg)
Figure 5 Development and validation of CD4TLAs-Related Risk Signature. (A) Co-expression network of lncRNAs and CD4+ Tconv-related genes. (B) Cross‐validation for tuning parameter selection in the lasso regression. (C) Validation was performed for tuning parameter selection through the least absolute shrinkage and selection operator (LASSO) regression model for overall survival (OS). (D) The heatmap shows the expression distribution of A. (E) CD4TLAs-related prognostic model constructed in total BC patients. (F) Validation cohort 1. (G) Validation cohort 2. From top to bottom: Kaplan-Meier survival analysis of high- and low-RS subgroups; risk score between high- and low-risk groups; patient survival time; plots of the AUC for time-dependent ROC performance.
Development and Validation of Prognostic Signatures Associated With CD4 TLAs in BC
From single-cell data from BC patients, we obtained CD4+ Tconv-related genes and identified CD4TLAs that were significantly dysregulated in BC. To explore the role of these CD4TLAs in BC, we further analyzed them based on TCGA-BRCA data. By using univariate COX analysis on the significantly different CD4TLAs, we screened lncRNAs significantly associated with prognosis. Further, we constructed a BC prognostic model by suing the LASSO algorithm and generated a risk score (RS). We finally identified risk signatures of 16 lncRNAs (Figures 5B, C). The risk score was calculated as follows: RS = (-0.145 * AC137630.3 exp.) + (-0.041 * SIDT1-AS1 exp.) + (-0.510 * AL133467.1 exp.) + (-0.595 * LIPE-AS1 exp.) + (0.369 * AC012213.3 exp.) + (-0.036 * LRRC8C-DT exp.) + (-0.113 * USP30-AS1 exp.) + (-0.064 * AL691482.3 exp.) + (0.130 * YTHDF3-AS1 exp.) + (0.606 * AC000067.1 exp.) + (-0.028 * LINC02446 exp.) + (-0.027 * U62317.1 exp.) + (-0.012 * MIAT exp.) + (-0.033 * SLC12A5-AS1 exp.) + (0.038 * LINC01235 exp.) + (1.089 * AL138789.1 exp.). The heatmap shows the expression levels of 16 lncRNAs in BC, and we found that most of the lncRNAs were significantly highly expressed in BC, except LIPE-AS1, LINC01235, AL133467.1, LRRC8C-DT (Figure 5D). According to the median RS, BC patients were divided into low-risk and high-risk groups. We found that patients in the high-RS group had a poorer prognosis (Figure 5E). To assess the predictive ability of the model, we constructed time-dependent ROC curves, and the results showed that the CD4 TLAs-related prognostic model had high accuracy in BC patients (Figure 5E; AUC at 1, 3, and 5 years were 0.742, 0.751, and 0.723, respectively).
TCGA-BC patients were divided in a 1:1 ratio into a validation cohort 1 and a validation cohort 2 to further test the accuracy of the model. RS were calculated by the same algorithm, noting that high-RS was significantly associated with poorer prognosis in both validation cohort 1 (Figure 5F) and validation cohort 2 (Figure 5G). The results of AUC further confirmed the stability of the model (Figure 5F, AUC at 1, 3, and 5 years were 0.675, 0.646, and 0.651, respectively; Figure 5G, AUC at 1, 3, and 5 years were 0.701, 0.696, and 0.691, respectively). Furthermore, by univariate and multivariate COX analysis, we found that RS was an independent prognostic factor in BC patients (Supplementary Figure 3).
Identification of the Role of RS in Clinical Subgroups
We analyzed RS in different clinical characteristics in order to further develop the application value and test the accuracy of RS (Age, Stage, T, M, N). The results showed that there were no significant differences in RS distribution among the different clinical subgroups (Figures 6A–E). Furthermore, it is worth noting that RS also had a prognostic role in most clinical subgroups and that high RS was significantly associated with poorer prognosis. To improve the clinical utility of the CD4TLAs model, we constructed a predictive nomogram model incorporating clinicopathological features and RS in the TCGA dataset. The accuracy of the nomogram was checked with a calibration curve (Supplementary Figure 4). The results showed that RS was the main predictor (Figure 7).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g006.jpg)
Figure 6 Relationship between risk scores and clinical characteristics. Relationship between risk score and Age (A), Stage (B), T (C), M (D), N (E). From left to right are the distribution level exploration and survival analysis of RS in different clinical subgroups.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g007.jpg)
Figure 7 A nomogram was established integrating the RS, Age, Stage, T, M, N, for predicting 1-, 3-, and 5-year survival outcomes of BC. To calculate survival probability, identify patient values on each axis and then for each draw a vertical line upward to the points axis. Add the points for all variables and locate this sum on the total points line. * P < 0.05, ** P < 0.01, and *** P < 0.001.
Analysis of the Relationship Between RS and Tumor Microenvironment
The TME has been confirmed to play a key role in the occurrence and development of BC. Exploring the relationship between RS and tumor microenvironment is helpful in the development of new functions of RS. We evaluated immune scores (Figure 8A), stroma scores (Figure 8B), and estimate scores (Figure 8C) in BC based on the estimate algorithm, and we found that high-RS groups tended to have lower scores. Furthermore, we explored the relationship between TME scores and prognosis, and results revealed that a high-immune score had a significant survival advantage (Figures 8D–F). As a further analysis, we used the ssGSEA to examine the distribution of immune cell infiltration and the enrichment of immune-related functional pathways in high- and low-RS subgroups (Figures 8G, H). In the low-RS group, immune infiltrating cells and immune-related pathways were significantly enriched. Another important finding was that RS was negatively correlated with most immune infiltrating cells (eg T cell follicular helper, T cell CD8+, B cell memory, Macrophage M1, NK cell activated, T cell regulatory, T cell gamma delta, T cell CD4+ memory activated), except Macrophage M0, Macrophage M2 (Figure 8I).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g008.jpg)
Figure 8 The relationship between RS and tumor microenvironment. Comparison of the immune score (A), stromal score (B), estimate score (C) in high- and low-RS groups. Relationship between immune score (D), stromal score (E), estimated score (F), and BC prognosis. Enrichment analysis of immune cell infiltration (G) and immune-related pathways (H) in high- and low-RS groups. (I) Correlation analysis of immune infiltrating cells and RS. *** P < 0.001.
Evaluation of the Immunotherapy Response and Chemotherapy Based on RS
CD4TLAs-associated RS has been shown to play a very important role in the prognosis and TME of BC patients. Immunotherapy has been widely used in the clinical treatment of tumor patients, therefore, we further explored the role of RS in immunotherapy. Based on our analysis of immune checkpoint expression levels, we found that the low-RS group exhibited significantly higher expression levels of all immune checkpoints except CD276 (Figure 9A). In addition, we evaluated the response to immunotherapy in the high and low RS groups by the TIDE algorithm, and the results showed that the low RS group had a lower TIDE score and was more sensitive to immunotherapy (Figure 9B). We obtained the IPS data of BC patients through the TCIA database, and further evaluated the role of RS in IPS immunotherapy (Figures 9C–F). The results showed that the median line value of the low-RS group was significantly higher than that of the high-RS group, which once again proved that the low-RS group was more sensitive to immunotherapy. These findings can assist us in developing potential functions of RS as well as providing clinical assistance. GSEA enrichment analysis clarified the functional pathways that were significantly enriched in the high and low RS groups, and the results showed that they were mainly associated with cell cycle and cell receptor regulators (Figure 9G). We also evaluated chemotherapeutic drug responses in patients with high- and low-RS by the R package “pRRophetic”, where the half-maximal inhibitory concentration (IC50) was estimated by ridge regression and the prediction accuracy (Figures 10A–P). The results showed that most chemotherapeutic drugs had higher IC50 scores in the high-RS group, except for Bicalutamide (Figure 10O) and Imatinib (Figure 10P). Based on the above results, we believe that RS may help guide the clinical treatment of BC.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g009.jpg)
Figure 9 Application of RS in Immunotherapy. (A) Expression levels of immune checkpoints in high- and low-RS groups. (B) Prediction of immunotherapy response based on the TIDE algorithm. IPS score for immunotherapy. (C) CTLA4− PD1+. (D) CTLA4− PD1−. (E) CTLA4+ PD1−. (F) CTLA4+ PD1+. (G) GSEA enrichment analysis results. * P < 0.05, ** P < 0.01, and *** P < 0.001.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/880769/fimmu-13-880769-HTML/image_m/fimmu-13-880769-g010.jpg)
Figure 10 (A–P) Analysis of drug sensitivity in high- and low-RS groups. The distribution of IC50 score. The abscissa represents the high- and low-RS groups, and the ordinate represents the distribution of the IC50 score.
Discussion
In the clinical treatment of various tumors, immunotherapy has been widely used. Advances in therapy to block inhibitory T cells checkpoint molecules such as CTLA-4 or PD-1/PD-L1 have shown clinical success, but only in a subset of BC patients (25–31). Adaptive T cell responses are very important for controlling tumor growth through the production of inflammatory cytokines and direct cytolytic targeting. CD4+ Tconv is also one of the major contributors to the antitumor response and can enhance patient immunity through the dendritic cell (DC) licensing and stimulation of pro-inflammatory myeloid cell programs (32–35). CD4+ T cells have been shown to play a crucial role in tumor progression and clinical treatment of BC patients (36–38). Therefore, studying CD4+ Tconv contributes to our understanding of anti-tumor and the development of novel immunotherapeutic approaches. More and more studies have shown that lncRNAs are involved in the regulation of autoimmunity and inflammation-related processes. However, research on CD4TLAs prognostic signature in patients with BC is still lacking.
As a first step in this study, we identified the major cell subsets of BC by single-cell data analysis, and we found that CD4+ Tconv was the largest cell subset and was significantly associated with multiple signaling pathways in BC. CD4+ Tconv is therefore believed to play a very important role in the development and occurrence of BC. The lncRNA network in Figure 5A shows the regulatory relationship between CD4+ Tconv-related genes and their targeted lncRNAs (CD4TLAs). LncRNAs are fundamental components of cellular processes and are involved in the occurrence and development of BC (39). CD4TLAs are thought to have important regulatory roles in BC and are associated with CD4+ Tconv cells. We screened and validated CD4TLAs and identified 244 CD4TLAs between tumor and normal groups. Then lncRNAs significantly associated with prognosis were screened out by univariate COX analysis, and the LASSO algorithm was used to construct a BC prognosis model and generate RS. We finally identified risk signatures for 16 lncRNAs. The accuracy and stability of the model were validated on TCGA-BRAC data, validation cohort 1, validation cohort 2, and different clinical subgroups. Having good accuracy and applicability, our model shows good performance. The results of univariate and multivariate COX analysis further validated RS as an independent prognostic factor in BC patients. The results of RS-TME relationship analysis showed that immune infiltrating cells were significantly enriched in the low-RS group, and Immune-related functional pathways are more active. We also found that RS was negatively correlated with the majority of immune cells in the TME. We believe that RS can characterize and represent TME to some extent. In addition, we further evaluated the application of RS in clinical treatment, and we found that patients in the low-RS group were more sensitive to immunotherapy and had lower IC50 scores for most chemotherapeutic drugs. This study suggest that the prognostic risk model of 16 CD4TLAs may be a reliable independent prognostic factor for breast cancer in clinical practice and can provide potential help for clinical treatment.
In recent years, precision genomic medicine has been focused on identifying specific prognostic indicators of survival based on large medical datasets with clinical outcomes. The development of prognostic risk models in BC has been increasing, and lncRNAs in BC have also attracted attention, such as ferroptosis-related lncRNAs, immune-related lncRNAs, stemness-related lncRNAs, autophagy-related lncRNA, etc (40–47). However, the role of CD4TLAs in BC prognosis remains unclear. To that end, we developed a CD4TLAs prognostic risk model to assist clinicians in determining personalized treatment strategies and to inform CD4TLAs research.
There are also some limitations to this study. First, this is a public database-based study that still needs further validation in vivo and in vitro, especially for clinical applications. Second, future studies should identify the regulatory T cells (Treg cells) apart from conventional T cells and Treg-related lncRNA signature and its relationship with prognosis and immunotherapy response in BC. This would provide a more complete picture. Finally, to further validate our bioinformatics predictions, in-depth studies of the 16 CD4TLAs, including functional experiments and molecular mechanisms, are required.
Conclusion
In conclusion, we identified a novel BC prognostic signature consisting of 16 CD4TLAs. With further prospective validation, these signatures may become a therapeutic target for BC and inform the personalized treatment of BC patients.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Materials.
Author Contributions
SN and QH designed study, analyzed data, and wrote the manuscript. QH, JW, YP, KQ, and LL analyzed data and contributed in writing the manuscript. All authors contributed to the article and approved the submitted version.
Funding
The work was supported by National Natural Science Foundation of China grants (82002779, 81760478); and Guangxi Provincial Natural Science Foundation of China grants (2019GXNSFAA245083); and Youth Science Foundation of Guangxi Medical University (GXMUYSF202224).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fimmu.2022.880769/full#supplementary-material
Supplementary Figure 1 | Cross-shared CD4+ Tconv-related genes in the GSE110686 and GSE114727 datasets.
Supplementary Figure 2 | Identification of CD4 TLAs that were significantly different between tumor and normal groups. (A) Volcano map. (B) Heatmap.
Supplementary Figure 3 | RS was an independent prognostic factor in BC patients. Univariate (A) and Multivariate (B) COX analysis.
Supplementary Figure 4 | The accuracy of the nomogram was checked with a calibration curve.
Abbreviations
TCGA, The Cancer Genome Atlas; GEO, Gene-Expression Omnibus; Tconv, T conventional; CD4TLA, CD4+ Tconv-related lncRNA; FPKM, Fragments per Kilobase Million; TPM, Transcripts per Million; RS, risk score; TME, Tumor microenvironment; OS, Overall survival; ROC, receiver operating characteristic; TME, Tumor microenvironment; GO, Gene ontology; KEGG, Kyoto Encyclopedia of Genes and Genomes; DEGs, differentially expressed genes.
References
1. Barzaman K, Karami J, Zarei Z, Hosseinzadeh A, Kazemi MH, Moradi-Kalbolandi S, et al. Breast Cancer: Biology, Biomarkers, and Treatments. Int Immunopharmacol (2020) 84:106535. doi: 10.1016/j.intimp.2020.106535
2. Liang Y, Zhang H, Song X, Yang Q. Metastatic Heterogeneity of Breast Cancer: Molecular Mechanism and Potential Therapeutic Targets. Semin Cancer Biol (2020) 60:14–27. doi: 10.1016/j.semcancer.2019.08.012
3. Thorat MA, Balasubramanian R. Breast Cancer Prevention in High-Risk Women. Best Pract Res Clin Obstet Gynaecol (2020) 65:18–31. doi: 10.1016/j.bpobgyn.2019.11.006
4. Britt KL, Cuzick J, Phillips KA. Key Steps for Effective Breast Cancer Prevention. Nat Rev Cancer (2020) 20(8):417–36. doi: 10.1038/s41568-020-0266-x
5. Braden A M, Stankowski R V, Engel J M, A Onitilo A. Breast Cancer Biomarkers: Risk Assessment, Diagnosis, Prognosis, Prediction of Treatment Efficacy and Toxicity, and Recurrence. Curr Pharm Des (2014) 20(30):4879–98. doi: 10.2174/1381612819666131125145517
6. Aalijahan H, Ghorbian S. Long Non-Coding RNAs and Cervical Cancer. Exp Mol Pathol (2019) 106:7–16. doi: 10.1016/j.yexmp.2018.11.010
7. Fang Y, Fullwood MJ. Roles, Functions, and Mechanisms of Long Non-Coding RNAs in Cancer. Genomics Proteomics Bioinf (2016) 14(1):42–54. doi: 10.1016/j.gpb.2015.09.006
8. Chandra Gupta S, Nandan Tripathi Y. Potential of Long non-Coding RNAs in Cancer Patients: From Biomarkers to Therapeutic Targets. Int J Cancer (2017) 140(9):1955–67. doi: 10.1002/ijc.30546
9. Ashrafizaveh S, Ashrafizadeh M, Zarrabi A, Husmandi K, Zabolian A, Shahinozzaman M, et al. Long Non-Coding RNAs in the Doxorubicin Resistance of Cancer Cells. Cancer Lett (2021) 508:104–14. doi: 10.1016/j.canlet.2021.03.018
10. Zhao W, Geng D, Li S, Chen Z, Sun M. LncRNA HOTAIR Influences Cell Growth, Migration, Invasion, and Apoptosis via the miR-20a-5p/HMGA2 Axis in Breast Cancer. Cancer Med (2018) 7(3):842–55. doi: 10.1002/cam4.1353
11. Zhou J, Zhang S, Luo M. LncRNA PCAT7 Promotes the Malignant Progression of Breast Cancer by Regulating ErbB/PI3K/Akt Pathway. Future Oncol (2021) 17(6):701–10. doi: 10.2217/fon-2020-0273
12. García-Venzor A, Mandujano-Tinoco EA, Lizarraga F, Zampedri C, Krötzsch E, Salgado RM, et al. Microenvironment-Regulated lncRNA-HAL is Able to Promote Stemness in Breast Cancer Cells. Biochim Biophys Acta Mol Cell Res (2019) 1866(12):118523. doi: 10.1016/j.bbamcr.2019.118523
13. Liu LL, Zhu SG, Jiang XY, Ren J, Lin Y, Zhang NN, et al. LncRNA Expression in CD4+ T Cells in Neurosyphilis Patients. Front Cell Infect Microbiol (2017) 7:461. doi: 10.3389/fcimb.2017.00461
14. Fu J, Shi H, Wang B, Zhan T, Shao Y, Ye L, et al. LncRNA PVT1 Links Myc to Glycolytic Metabolism Upon CD4+ T Cell Activation and Sjögren's Syndrome-Like Autoimmune Response. J Autoimmun (2020) 107:102358. doi: 10.1016/j.jaut.2019.102358
15. DeRogatis JM, Viramontes KM, Neubert EN, Tinoco R. PSGL-1 Immune Checkpoint Inhibition for CD4+ T Cell Cancer Immunotherapy. Front Immunol (2021) 12:636238. doi: 10.3389/fimmu.2021.636238
16. Spitzer MH, Carmi Y, Reticker-Flynn NE, Kwek SS, Madhireddy D, Martins MM, et al. Systemic Immunity Is Required for Effective Cancer Immunotherapy. Cell (2017) 168(3):487–502.e15. doi: 10.1016/j.cell.2016.12.022
17. Azizi E, Carr AJ, Plitas G, Cornish AE, Konopacki C, Prabhakaran S, et al. Single-Cell Map of Diverse Immune Phenotypes in the Breast Tumor Microenvironment. Cell (2018) 174(5):1293–1308.e36. doi: 10.1016/j.cell.2018.05.060
18. Savas P, Virassamy B, Ye C, Salim A, Mintoff CP, Caramia F, et al. Single-Cell Profiling of Breast Cancer T Cells Reveals a Tissue-Resident Memory Subset Associated With Improved Prognosis. Nat Med (2018) 24(7):986–93. doi: 10.1038/s41591-018-0078-7
19. Sun D, Wang J, Han Y, Dong X, Ge J, Zheng R, et al. TISCH: A Comprehensive Web Resource Enabling Interactive Single-Cell Transcriptome Visualization of Tumor Microenvironment. Nucleic Acids Res (2021) 49(D1):D1420–30. doi: 10.1093/nar/gkaa1020
20. Xu M, Li Y, Li W, Zhao Q, Zhang Q, Le K, et al. Immune and Stroma Related Genes in Breast Cancer: A Comprehensive Analysis of Tumor Microenvironment Based on the Cancer Genome Atlas (TCGA) Database. Front Med (Lausanne) (2020) 7:64. doi: 10.3389/fmed.2020.00064
21. Ning S, Li H, Qiao K, Wang Q, Shen M, Kang Y, et al. Identification of Long-Term Survival-Associated Gene in Breast Cancer. Aging (Albany NY) (2020) 12(20):20332–49. doi: 10.18632/aging.103807
22. Zhou Y, Zhou B, Pache L, Chang M, Khodabakhshi AH, Tanaseichuk O, et al. Metascape Provides a Biologist-Oriented Resource for the Analysis of Systems-Level Datasets. Nat Commun (2019) 10(1):1523. doi: 10.1038/s41467-019-09234-6
23. Fu J, Li K, Zhang W, Wan C, Zhang J, Jiang P, et al. Large-Scale Public Data Reuse to Model Immunotherapy Response and Resistance. Genome Med (2020) 12(1):21. doi: 10.1186/s13073-020-0721-z
24. Van Allen EM, Miao D, Schilling B, Shukla SA, Blank C, Zimmer L, et al. Genomic Correlates of Response to CTLA-4 Blockade in Metastatic Melanoma. Science (2015) 350(6257):207–211. pii: aaf8264. doi: 10.1126/science.aaf8264
25. de Melo Gagliato D, Buzaid AC, Perez-Garcia J, Cortes J. Immunotherapy in Breast Cancer: Current Practice and Clinical Challenges. BioDrugs (2020) 34(5):611–23. doi: 10.1007/s40259-020-00436-9
26. Lyons TG. Targeted Therapies for Triple-Negative Breast Cancer. Curr Treat Options Oncol (2019) 20(11):82. doi: 10.1007/s11864-019-0682-x
27. Li CW, Lim SO, Chung EM, Kim YS, Park AH, Yao J, et al. Eradication of Triple-Negative Breast Cancer Cells by Targeting Glycosylated PD-L1. Cancer Cell (2018) 33(2):187–201.e10. doi: 10.1016/j.ccell.2018.01.009
28. Liu L, Wang Y, Miao L, Liu Q, Musetti S, Li J, et al. Combination Immunotherapy of MUC1 mRNA Nano-Vaccine and CTLA-4 Blockade Effectively Inhibits Growth of Triple Negative Breast Cancer. Mol Ther (2018) 26(1):45–55. doi: 10.1016/j.ymthe.2017.10.020
29. Jang BS, Han W, Kim IA. Tumor Mutation Burden, Immune Checkpoint Crosstalk and Radiosensitivity in Single-Cell RNA Sequencing Data of Breast Cancer. Radiother Oncol (2020) 142:202–9. doi: 10.1016/j.radonc.2019.11.003
30. Peng Z, Su P, Yang Y, Yao X, Zhang Y, Jin F, et al. Identification of CTLA-4 Associated With Tumor Microenvironment and Competing Interactions in Triple Negative Breast Cancer by Co-Expression Network Analysis. J Cancer (2020) 11(21):6365–75. doi: 10.7150/jca.46301
31. Zhang Y, Li Z, Chen M, Chen H, Zhong Q, Liang L, et al. lncRNA TCL6 Correlates With Immune Cell Infiltration and Indicates Worse Survival in Breast Cancer. Breast Cancer (2020) 27(4):573–85. doi: 10.1007/s12282-020-01048-5
32. Behrens G, Li M, Smith CM, Belz GT, Mintern J, Carbone FR, et al. Helper T Cells, Dendritic Cells and CTL Immunity. Immunol Cell Biol (2004) 82(1):84–90. doi: 10.1111/j.1440-1711.2004.01211.x
33. Phiwpan K, Zhou X. MicroRNAs in Regulatory T Cells. Cancer Lett (2018) 423:80–5. doi: 10.1016/j.canlet.2018.03.012
34. Corthay A, Skovseth DK, Lundin KU, Røsjø E, Omholt H, Hofgaard PO, et al. Primary Antitumor Immune Response Mediated by CD4+ T Cells. Immunity (2005) 22(3):371–83. doi: 10.1016/j.immuni.2005.02.003
35. Corthay A, Lundin KU, Lorvik KB, Hofgaard PO, Bogen B. Secretion of Tumor-Specific Antigen by Myeloma Cells is Required for Cancer Immunosurveillance by CD4+ T Cells. Cancer Res (2009) 69(14):5901–7. doi: 10.1158/0008-5472.CAN-08-4816
36. Gu-Trantien C, Loi S, Garaud S, Equeter C, Libin M, de Wind A, et al. CD4+ Follicular Helper T Cell Infiltration Predicts Breast Cancer Survival. J Clin Invest (2013) 123(7):2873–92. doi: 10.1172/JCI67428
37. Song PN, Mansur A, Dugger KJ, Davis TR, Howard G, Yankeelov TE, et al. CD4 T-Cell Immune Stimulation of HER2 + breast Cancer Cells Alters Response to Trastuzumab In Vitro. Cancer Cell Int (2020) 20(1):544. doi: 10.1186/s12935-020-01625-w
38. Li R, Chibbar R, Xiang J. Novel EXO-T Vaccine Using Polyclonal CD4+ T Cells Armed With HER2-Specific Exosomes for HER2-Positive Breast Cancer. Onco Targets Ther (2018) 11:7089–93. doi: 10.2147/OTT.S184898
39. Crudele F, Bianchi N, Reali E, Galasso M, Agnoletto C, Volinia S. The Network of non-Coding RNAs and Their Molecular Targets in Breast Cancer. Mol Cancer (2020) 19(1):61. doi: 10.1186/s12943-020-01181-x
40. Zhang K, Ping L, Du T, Liang G, Huang Y, Li Z, et al. A Ferroptosis-Related lncRNAs Signature Predicts Prognosis and Immune Microenvironment for Breast Cancer. Front Mol Biosci (2021) 8:678877. doi: 10.3389/fmolb.2021.678877
41. Ma W, Zhao F, Yu X, Guan S, Suo H, Tao Z, et al. Immune-Related lncRNAs as Predictors of Survival in Breast Cancer: A Prognostic Signature. J Transl Med (2020) 18(1):442. doi: 10.1186/s12967-020-02522-6
42. Li X, Li Y, Yu X, Jin F. Identification and Validation of Stemness-Related lncRNA Prognostic Signature for Breast Cancer. J Transl Med (2020) 18(1):331. doi: 10.1186/s12967-020-02497-4
43. Shen Y, Peng X, Shen C. Identification and Validation of Immune-Related lncRNA Prognostic Signature for Breast Cancer. Genomics (2020) 112(3):2640–6. doi: 10.1016/j.ygeno.2020.02.015
44. Xiao R, Yang M, Tan Y, Ding R, Li D. Identification of Five Immune-Related lncRNAs Predicting Survival and Tumor Microenvironment Characteristics in Breast Cancer. Comput Math Methods Med (2021) 2021:6676692. doi: 10.1155/2021/6676692
45. Li X, Jin F, Li Y. A Novel Autophagy-Related lncRNA Prognostic Risk Model for Breast Cancer. J Cell Mol Med (2021) 25(1):4–14. doi: 10.1111/jcmm.15980
46. Zhang Y, Li Z, Chen M, Chen H, Zhong Q, Liang L, et al. Identification of a New Eight-Long Noncoding RNA Molecular Signature for Breast Cancer Survival Prediction. DNA Cell Biol (2019) 38(12):1529–39. doi: 10.1089/dna.2019.5059
Keywords: long non-coding RNA, CD4+ conventional T cells, breast cancer, The Cancer Genome Atlas, prognostic signature
Citation: Ning S, Wu J, Pan Y, Qiao K, Li L and Huang Q (2022) Identification of CD4+ Conventional T Cells-Related lncRNA Signature to Improve the Prediction of Prognosis and Immunotherapy Response in Breast Cancer. Front. Immunol. 13:880769. doi: 10.3389/fimmu.2022.880769
Received: 21 February 2022; Accepted: 12 April 2022;
Published: 04 May 2022.
Edited by:
Murugaiyan Gopal, Harvard Medical School, United StatesReviewed by:
Bhalchandra Mirlekar, University of North Carolina at Chapel Hill, United StatesEva Reali, University of Milano-Bicocca, Italy
Copyright © 2022 Ning, Wu, Pan, Qiao, Li and Huang. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lei Li, wshililei@buaa.edu.cn; Qinghua Huang, huangqinghua@stu.gxmu.edu.cn
†These authors have contributed equally to this work