UTRN as a potential biomarker in breast cancer: a comprehensive bioinformatics and in vitro study

Utrophin (UTRN), known as a tumor suppressor, potentially regulates tumor development and the immune microenvironment. However, its impact on breast cancer's development and treatment remains unstudied. We conducted a thorough examination of UTRN using both bioinformatic and in vitro experiments in this study. We discovered UTRN expression decreased in breast cancer compared to standard samples. High UTRN expression correlated with better prognosis. Drug sensitivity tests and RT-qPCR assays revealed UTRN's pivotal role in tamoxifen resistance. Furthermore, the Kruskal-Wallis rank test indicated UTRN's potential as a valuable diagnostic biomarker for breast cancer and its utility in detecting T stage of breast cancer. Additionally, our results demonstrated UTRN's close association with immune cells, inhibitors, stimulators, receptors, and chemokines in breast cancer (BRCA). This research provides a novel perspective on UTRN's role in breast cancer's prognostic and therapeutic value. Low UTRN expression may contribute to tamoxifen resistance and a poor prognosis. Specifically, UTRN can improve clinical decision-making and raise the diagnosis accuracy of breast cancer.

Exploring Prognostic Immune Microenvironment-Related Genes in Head and Neck Squamous Cell Carcinoma from the TCGA Database

Purpose: Head and neck squamous cell carcinoma (HNSCC) has a high rate of local and distant metastases. In tumor tissues, the interaction between tumor cells and the tumor microenvironment (TME) is closely related to cancer development and prognosis. Therefore, screening for TME-related genes in HNSCC is crucial for understanding metastatic patterns. Methods: Our research relied mainly on a novel algorithm called Estimation of STromal and Immune cells in MAlignant Tumors using Expression data (ESTIMATE). Fragments Per Kilobase of exon model per Million mapped fragments (FPKM) data and HNSCC clinical data were obtained from the TCGA database, and the purity of HNSCC tissue and the features of stromal and immune cell infiltration were determined. Furthermore, differentially expressed genes (DEGs) were screened based on immune, stromal, and ESTIMATE scores, and their protein-protein interaction (PPI) networks and ClueGO functions were evaluated. Finally, the expression profiles of DEGs related to immunity in HNSCC were determined. Differential gene expression was verified in the highly invasive oral cancer cell lines (SCC-25, CAL-27, and FaDu) and oral cancer tissues. Results: Our analysis found that both the immune and ESTIMATE scores were significantly associated with the prognosis of HNSCC. Moreover, cross-validation using the Venn algorithm revealed that 433 genes were significantly upregulated, and 394 genes were significantly downregulated. All DEGs were associated with both ESTIMATE and immune scores. The enrichment of cytokine-cytokine receptor interactions and chemokine signaling pathways was observed using pathway enrichment analyses. We initially screened 25 genes after analyzing the key sub-networks of the PPI network. Survival analysis revealed the significance of CCR4, CXCR3, P2RY14, CCR2, CCR8, and CCL19 in relation to survival and their association with immune infiltration-related metastasis in HNSCC. Conclusions: The expression profiles of relevant TME-related genes were screened following stromal and immune cell scoring using ESTIMATE, and DEGs associated with survival were identified. These TME-related gene markers offer valuable utility as both prognostic indicators and markers denoting metastatic traits in HNSCC.

Identification of key immune cell-related genes involved in tumorigenesis and prognosis of cervical squamous cell carcinoma

The infiltration of immune cells can significantly affect the prognosis and immune therapy of patients with cervical squamous cell carcinoma (CSCC). This study aimed to explore key immune cell-related genes in the tumorigenesis and prognosis of CSCC. The module significantly related to immunity was screened by weighted gene co-expression network analysis (WGCNA) and ESTIMATE analysis, followed by correlation analysis with clinical traits. Key candidate genes were intersected with the protein-protein interaction (PPI) network genes for immune-related genes. The relationship between immune cell infiltration and key genes was analyzed. Tumor immune dysfunction and exclusion (TIDE) and immunophenoscore (IPS) predicted the response to immunotherapy in CSCC patients. Clinically, quantitative real-time polymerase chain reaction (qRT-PCR) and immunohistochemistry were manipulated for analyzing the changes in mRNA and protein expression of key genes in cancer. Western blot was conducted to assess the correlation between key genes and immune infiltration. The brown module was notably associated with the immune microenvironment of CSCC, from which three immune-related key genes (TYROBP, CCL5, and HLA-DRA) were obtained. High expression of these genes was significantly positively associated with the infiltration abundance of T cells, B cells, and other immune cells. High expression levels of three key genes were confirmed in para-cancer tissue and correlated with the abundance of immune cells. The high-expression group of key genes was more sensitive to immunotherapy. We provide a theoretical basis for searching for potential targets for effective treatment and diagnosis of CSCC and provide new ideas for developing novel immunotherapy strategies.

CD40LG as a Prognostic Molecular Marker Regulates Tumor Microenvironment Through Immune Process in Breast Cancer

Purpose

Breast cancer (BRCA) is the second most common malignancy in the world and the most common in women. Here, we utilized publicly available BRCA dataset to investigate potential prognosis-related genes through integrated bioinformatics analysis.

Materials and Methods

BRCA dataset was obtained from the Cancer Genome Atlas (TCGA) database. The ESTIMATE algorithm was used to calculate the ImmuneScores and StromalScores of the samples and then divided them into high- and low-score groups based on the median score. Common differentially expressed genes (DEGs) were identified through differential expression analysis. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed. The core prognostic genes were the intersection of hub genes from PPI network and prognostic genes from univariate Cox proportional hazard regression analysis. Finally, the CIBERSORT algorithm was used to calculate proportions of 22 tumor-infiltrating immune cells (TICs) in BRCA samples.

Results

A total of 486 DEGs were identified. These genes were mainly enriched in immune-related pathways. Crossover genes between the hub genes and the prognostic genes were CD2 and CD40LG. CD40LG was further investigated in this study. CD40LG was downregulated in BRCA samples compared with normal samples, and a lower CD40LG expression was associated with advanced tumor stages and a poor prognosis. CD40LG was shown to be involved in immune-related pathways of BRCA by Gene Set Enrichment Analysis. Finally, 14 TICs were found to have a close relationship with CD40LG.

Conclusion

CD40LG was found to be a core prognostic gene related to tumor microenvironment and deeply involved in immune-related pathways in BRCA. Our findings may provide new insights for exploring the molecular mechanisms of BRCA and developing new immunotherapies for the disease.

Identification of a novel tumour microenvironment-based prognostic biomarker in skin cutaneous melanoma

Skin cutaneous melanoma (SKCM) is one of the most destructive skin malignancies and has attracted worldwide attention. However, there is a lack of prognostic biomarkers, especially tumour microenvironment (TME)-based prognostic biomarkers. Therefore, there is an urgent need to investigate the TME in SKCM, as well as to identify efficient biomarkers for the diagnosis and treatment of SKCM patients. A comprehensive analysis was performed using SKCM samples from The Cancer Genome Atlas and normal samples from Genotype-Tissue Expression. TME scores were calculated using the ESTIMATE algorithm, and differential TME scores and differentially expressed prognostic genes were successively identified. We further identified more reliable prognostic genes via least absolute shrinkage and selection operator regression analysis and constructed a prognostic prediction model to predict overall survival. Receiver operating characteristic analysis was used to evaluate the diagnostic efficacy, and Cox regression analysis was applied to explore the relationship with clinicopathological characteristics. Finally, we identified a novel prognostic biomarker and conducted a functional enrichment analysis. After considering ESTIMATEScore and tumour purity as differential TME scores, we identified 34 differentially expressed prognostic genes. Using least absolute shrinkage and selection operator regression, we identified seven potential prognostic biomarkers (SLC13A5, RBM24, IGHV3OR16-15, PRSS35, SLC7A10, IGHV1-69D and IGHV2-26). Combined with receiver operating characteristic and regression analyses, we determined PRSS35 as a novel TME-based prognostic biomarker in SKCM, and functional analysis enriched immune-related cells, functions and signalling pathways. Our study indicated that PRSS35 could act as a potential prognostic biomarker in SKCM by investigating the TME, so as to provide new ideas and insights for the clinical diagnosis and treatment of SKCM.

Identification of novel subtypes based on ssGSEA in immune-related prognostic signature for tongue squamous cell carcinoma

Background

Tongue squamous cell carcinoma (TSCC) is characterized by aggressive invasion and poor prognosis. Currently, immune checkpoint inhibitors may prolong overall survival compared with conventional treatments. However, PD1/PDL1 remain inapplicable in predicting the prognosis of TSCC; thus, it is urgent to explore the genetic characteristics of TSCC.

Materials and methods

We utilized single‐sample gene set enrichment analysis (ssGSEA) to classify TSCC patients from the TCGA database into clusters with different immune cell infiltrations. ESTIMATE (immune‐related scores) and CIBERSORT (immune cell distribution) analyses were used to evaluate the immune landscape among clusters. GO, KEGG, and GSEA analyses were performed to analyze the different underlying molecular mechanisms in the clusters. Based on the immune characteristics, we applied the LASSO Cox regression to select hub genes and construct a prognostic risk model. Finally, we established an interactive network among these hub genes by using Cytoscape, and a pan‐cancer analysis to further verify and decipher the innate function of these genes.

Results

Using ssGSEA, we constructed three functional clusters with different overall survival and immune‐cell infiltration. ESTIMATE and CIBERSORT analyses revealed the different distributions of immune cells (T cells, B cells, and macrophages) with diverse immune‐related scores (ESTIMATE, immune, stromal, and tumor purity scores). Moreover, pathways including those of the interferon‐gamma response, hypoxia, and glycolysis of the different subtypes were investigated to elucidate their involvement in mediating the heterogeneous immune characteristics. Subsequently, after LASSO Cox regression, a signature of 15 immune‐related genes was established that is more prognostically effective than the TNM stage. Furthermore, three hub genes—PGK1, GPI, and RPE—were selected using Cytoscape evaluation and verified by immunohistochemistry. PGK1, the foremost regulator, was a comprehensively profiled pan‐cancer, and a PGK1‐based interactive network was established.

Conclusion

Our results suggest that immune‐related genes and clusters in TSCC have the potential to guide individualized treatments.

Characterization of Exosome-Related Gene Risk Model to Evaluate the Tumor Immune Microenvironment and Predict Prognosis in Triple-Negative Breast Cancer

Background

As a kind of small membrane vesicles, exosomes are secreted by most cell types from multivesicular endosomes, including tumor cells. The relationship between exosomes and immune response plays a vital role in the occurrence and development of tumors. Nevertheless, the interaction between exosomes and the microenvironment of tumors remains unclear. Therefore, we set out to study the influence of exosomes on the triple-negative breast cancer (TNBC) microenvironment.

Method

One hundred twenty-one exosome-related genes were downloaded from ExoBCD database, and IVL, CXCL13, and AP2S1 were final selected because of the association with TNBC prognosis. Based on the sum of the expression levels of these three genes, provided by The Cancer Genome Atlas (TCGA), and the regression coefficients, an exosome risk score model was established. With the median risk score value, the patients in the two databases were divided into high- and low-risk groups. R clusterProfiler package was employed to compare the different enrichment ways between the two groups. The ESTIMATE and CIBERSORT methods were employed to analyze ESTIMATE Score and immune cell infiltration. Finally, the correlation between the immune checkpoint-related gene expression levels and exosome-related risk was analyzed. The relationship between selected gene expression and drug sensitivity was also detected.

Results

Different risk groups exhibited distinct result of TNBC prognosis, with a higher survival rate in the low-risk group than in the high-risk group. The two groups were enriched by immune response and biological process pathways. A better overall survival (OS) was demonstrated in patients with high scores of immune and ESTIMATE rather than ones with low scores. Subsequently, we found that CD4⁺-activated memory T cells and M1 macrophages were both upregulated in the low-risk group, whereas M2 macrophages and activated mast cell were downregulated in the low-risk group in patients from the TCGA and GEO databases, respectively. Eventually, four genes previously proposed to be targets of immune checkpoint inhibitors were evaluated, resulting in the expression levels of CD274, CTLA4, LAG3, and TIM3 being higher in the low-risk group than high-risk group.

Conclusion

The results of our study suggest that exosome-related risk model was related to the prognosis and ratio of immune cell infiltration in patients with TNBC. This discovery may make contributions to improve immunotherapy for TNBC.

Combined DNA Methylation and Transcriptomic Assessments to Determine a Prognostic Model for PD-1-Negative Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) has the highest incidence and mortality of any malignancy in the world. Immunotherapy has been a major breakthrough for HCC treatment, but immune checkpoint inhibitors (ICIs) are effective in only a small percentage of HCC patients. In the present study, we screened programmed cell death protein 1 (PD-1) -negative HCC samples, which are frequently resistant to ICIs, and identified their methylation and transcription characteristics through the assessment of differential gene methylation and gene expression. We also screened for potential targeted therapeutic drugs using the DrugBank database. Finally, we used a LASSO (least absolute shrinkage and selection operator) regression analysis to construct a prognostic model based on three differentially methylated and expressed genes (DMEGs). The results showed that ESTIMATE (Estimation of Stromal and Immune Cells in Malignant Tumors using Expression Data) scores for the tumor samples were significantly lower compared to normal sample ESTIMATE scores. In addition, we identified 31 DMEGs that were able to distinguish PD-1-negative samples from normal samples. A functional enrichment analysis showed that these genes were involved in a variety of tumor-related pathways and immune-related pathways, and the DrugBank screening identified potential therapeutic drugs. Finally, the prognostic model based on three DMEGs (UBD, CD5L, and CD213A2) demonstrated good predictive power for HCC prognosis and was verified using an independent cohort. The present study demonstrated the methylation characteristics of PD-1-negative HCC samples, identified several potential therapeutic drugs, and proposed a prognostic model based on UBD, CD5L, and CD213A2 methylation expression. In conclusion, this work provides an in-depth understanding of methylation in HCC samples that are not sensitive to ICIs.

Disparity of Hepatocellular Carcinoma in Tumor Microenvironment-Related Genes and Infiltrating Immune Cells between Asian and Non-Asian Populations

Hepatocellular carcinoma (HCC) is the most common cause of primary liver cancer deaths worldwide. The major risk factors for liver cancer development are cirrhosis, hepatitis B virus (HBV), hepatitis C virus (HCV) infection, and chronic alcohol abuse. HCC displays heterogeneity in terms of biology, etiology, and epidemiology. In Southeast Asia and Africa, chronic HBV infection is a major risk factor for HCC, whereas chronic HCV infection is a risk factor for HCC in western countries and Japan. Environmental and genetic conditions also play a role in the regional and temporal variations in the incidence of HCC. In this study, we used the ESTIMATE (ESTIMATE, Estimation of stromal and immune cells in malignant tumor tissues using expression data) algorithm and the CIBERSOFT tool to analyze gene expression profiles and infiltrating immune cells in HCC between Asian and non-Asian patients. The results showed that stromal and immune scores were dependent on overall survival (OS) in non-Asian patients but not in Asian patients. Kaplan-Meier survival analysis revealed four differentially expressed genes (DEGs) that were significantly associated with OS in non-Asian patients only. CIBERSORT (CIBERSORT, Cell type identification by estimating relative subsets of known RNA transcripts) analysis indicated that the composition of infiltrating immune cells was significantly different between Asian and non-Asian patients. By parsing the subclasses of HCC, the ability to predict prognosis and guide therapeutic targets for potentially actionable HCC may be improved.

Identification of Genes Related to Immune Infiltration in the Tumor Microenvironment of Cutaneous Melanoma

Cutaneous melanoma (CM) is the leading cause of skin cancer deaths and is typically diagnosed at an advanced stage, resulting in a poor prognosis. The tumor microenvironment (TME) plays a significant role in tumorigenesis and CM progression, but the dynamic regulation of immune and stromal components is not yet fully understood. In the present study, we quantified the ratio between immune and stromal components and the proportion of tumor-infiltrating immune cells (TICs), based on the ESTIMATE and CIBERSORT computational methods, in 471 cases of skin CM (SKCM) obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed genes (DEGs) were analyzed by univariate Cox regression analysis, least absolute shrinkage, and selection operator (LASSO) regression analysis, and multivariate Cox regression analysis to identify prognosis-related genes. The developed prognosis model contains ten genes, which are all vital for patient prognosis. The areas under the curve (AUC) values for the developed prognostic model at 1, 3, 5, and 10 years were 0.832, 0.831, 0.880, and 0.857 in the training dataset, respectively. The GSE54467 dataset was used as a validation set to determine the predictive ability of the prognostic signature. Protein–protein interaction (PPI) analysis and weighted gene co-expression network analysis (WGCNA) were used to verify “real” hub genes closely related to the TME. These hub genes were verified for differential expression by immunohistochemistry (IHC) analyses. In conclusion, this study might provide potential diagnostic and prognostic biomarkers for CM.

Complement C1q (C1qA, C1qB, and C1qC) May Be a Potential Prognostic Factor and an Index of Tumor Microenvironment Remodeling in Osteosarcoma

The tumor microenvironment (TME) has important effects on the tumorigenesis and development of osteosarcoma (OS). However, the dynamic mechanism regulating TME immune and matrix components remains unclear. In this study, we collected quantitative data on the gene expression of 88 OS samples from The Cancer Genome Atlas (TCGA) database and downloaded relevant clinical cases of OS from the TARGET database. The proportions of tumor-infiltrating immune cells (TICs) and the numbers of immune and matrix components were determined by CIBERSORT and ESTIMATE calculation methods. Protein-protein interaction (PPI) network construction and Cox regression analysis were conducted to analyze differentially expressed genes (DEGs). The complement components C1qA, C1qB and C1qC were then determined to be predictive factors through univariate Cox analysis and PPI cross analysis. Further analysis found that the levels of C1qA, C1qB and C1qC expression were positively linked to OS patient survival time and negatively correlated with the clinicopathological feature percent necrosis at definitive surgery. The results of gene set enrichment analysis (GSEA) demonstrated that genes related to immune functions were significantly enriched in the high C1qA, C1qB and C1qC expression groups. Proportion analysis of TICs by CIBERSORT showed that the levels of C1qA, C1qB and C1qC expression were positively related to M1 and M2 macrophages and CD8+ cells and negatively correlated with M0 macrophages. These results further support the influence of the levels of C1qA, C1qB and C1qC expression on the immune activity of the TME. Therefore, C1qA, C1qB and C1qC may be potential indicators of remodeling in the OS TME, which is helpful to predict the prognosis of patients with OS and provide new ideas for immunotherapy for OS.

Identification of a tumor microenvironment-associated prognostic gene signature in bladder cancer by integrated bioinformatic analysis

Bladder cancer is a common malignancy in the urinary system. Stromal and immune cells in tumor microenvironments, including those in the bladder cancer microenvironment, can serve as prognostic markers. However, the complex processes of bladder cancer necessitate large-scale evaluation to better understand the underlying mechanisms and identify biomarkers for diagnosis and treatment. We used the Estimation of STromal and Immune cells in MAlignant Tumors using Expression data algorithm to assess the association between stromal and immune cell-related genes and overall survival of patients with bladder cancer. We also identified and evaluated differentially expressed genes between cancer and non-cancer tissues from The Cancer Genome Atlas. Patients were categorized into different prognosis groups according to their stromal/immune scores based on differential gene expression. In addition, the prognostic value of the differentially expressed genes was assessed in a separate validation cohort using the Gene Expression Omnibus microarray dataset GSE13507, which identified nine genes (TNC, CALD1, PALLD, TAGLN, TGFB1I1, HSPB6, RASL12, CPXM2, and CYR61) associated with overall survival. Multivariate regression analysis showed that three genes (TNC, CALD1, and PALLD) were possible independent prognostic markers for patients with bladder cancer. Multiple gene set enrichment analysis of individual genes showed strong correlations with stromal and immune interactions, indicating that these nine genes may be related to carcinogenesis, invasion, and metastasis of bladder cancer. These findings provide useful insight into the molecular mechanisms of bladder cancer development, and suggest candidate biomarkers for prognosis and treatment.

Prognostic value of LRRC4C in Colon and Gastric Cancers correlates with Tumour Microenvironment Immunity

In this study, we aimed to use ESTIMATE and CIBERSORT computational methods to analyse transcriptional information on COAD and STAD in TCGA. We downloaded transcriptome RNA-seq data of 446 patients with colon cancer from TCGA and estimated the amount of immune and stromal components in the COAD samples using CIBERSORT algorithms. We analysed differentially expressed genes in 446 TCGA samples and 585 Series GSE39582 samples, in high- and low-scoring groups, using Cox regression. The expression of LRRC4C, correlated positively with clinicopathological characteristics and negatively with the survival of patients with COAD. Single-gene survival analysis using Gene Expression Profiling Interactive Analysis 2.0 and Kaplan-Meier plotter revealed an association between high levels of LRRC4C expression and poor prognosis in patients with colon and gastric cancers. Gene set enrichment analysis of COAD and STAD samples indicated that genes in groups expressing high and low LRRC4C levels were mainly enriched in immune-related activities and metabolic pathways, respectively. Difference and correlation analyses of the relationship between LRRC4C expression and tumour-infiltrating immune cells, determined using CIBERSORT algorithms, revealed that monocytes, resting mast cells, and M2 macrophages were positively correlated with LRRC4C expression.

CCR7 Has Potential to Be a Prognosis Marker for Cervical Squamous Cell Carcinoma and an Index for Tumor Microenvironment Change

The tumor microenvironment (TME) has an essential role in the development of cervical squamous cell carcinoma (CSCC); however, the dynamic role of the stromal and immune cells is still unclear in TME. We downloaded data from The Cancer Genome Atlas (TCGA) database and applied ESTIMATE and CIBERSORT algorithms to measure the quantity of stromal and immune cells and the composition of tumor-infiltrating immune cell (TIC) in 253 CSCC cases. The protein-protein interaction (PPI) network and Cox regression analysis presented the differentially expressed genes (DEGs). Then, C-C chemokine receptor type 7 (CCR7) was screened out as a prognostic marker by the univariate Cox and intersection analysis of PPI. Further analysis showed a positive correlation between the expression of CCR7 and the survival of CSCC patients. The result of the Gene Set Enrichment Analysis (GSEA) of genes in the high CCR7 expression group displayed a predominant enrichment in immune-related pathways. An enrichment in metabolic activities was observed in the low CCR7 expression group. CIBERSORT analysis showed a positive correlation between Plasma cells, CD8⁺ T cells, and regulatory T cells and the CCR7 expression, suggesting that CCR7 might play a crucial role in maintaining the immunological dominance status for TME. Therefore, the expression level of CCR7 might help predict the survival of CSCC cases and be an index that the status of TME transitioned from immunological dominance to metabolic activation, which presented a new insight into the treatment of CSCC.

Identification of Prognostic Stromal-Immune Score-Based Genes in Hepatocellular Carcinoma Microenvironment

A growing amount of evidence has suggested the clinical importance of stromal and immune cells in the liver cancer microenvironment. However, reliable prognostic signatures based on assessments of stromal and immune components have not been well-established. This study aimed to identify stromal-immune score–based potential prognostic biomarkers for hepatocellular carcinoma. Stromal and immune scores were estimated from transcriptomic profiles of a liver cancer cohort from The Cancer Genome Atlas using the ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumors using Expression data) algorithm. Least absolute shrinkage and selection operator (LASSO) algorithm was applied to select prognostic genes. Favorable overall survivals and progression-free interval were found in patients with high stromal score and immune score, and 828 differentially expressed genes were identified. Functional enrichment analysis and protein–protein interaction networks further showed that these genes mainly participated in immune response, extracellular matrix, and cell adhesion. MMP9 (matrix metallopeptidase 9) was identified as a prognostic tumor microenvironment–associated gene by using LASSO and TIMER (Tumor IMmune Estimation Resource) algorithms and was found to be positively correlated with immunosuppressive molecules and drug response.

IGLL5 is correlated with tumor-infiltrating immune cells in clear cell renal cell carcinoma

Renal cell carcinomas (RCCs) account for about 90% of renal tumors, and their major histological subtype is ccRCC (clear cell RCC). Increasing evidence has indicated that the tumor microenvironment plays a significant role in the occurrence and development of ccRCC. In this study, we used ESTIMATE and CIBERSORT computational methods to calculate the proportion of immune and stromal components and the rate of TICs (tumor‐infiltrating immune cells) in 539 ccRCC samples from The Cancer Genome Atlas database. By examining the intersection of the differentially expressed genes obtained by the protein–protein interaction network and Cox regression analysis, we identified only one overlapping gene: IGLL5 (immunoglobulin lambda‐like polypeptide 5). We report that IGLL5 expression is correlated with TICs. Furthermore, our immunoinfiltration analyses revealed that three types of TIC are positively correlated with IGLL5 expression. IGLL5 may have potential as a prognostic biomarker of ccRCC.

Integrated Bioinformatical Analysis Identifies GIMAP4 as an Immune-Related Prognostic Biomarker Associated With Remodeling in Cervical Cancer Tumor Microenvironment

Tumor microenvironment (TME) is emerging as an essential part of cervical cancer (CC) tumorigenesis and development, becoming a hotspot of research these years. However, comprehending the specific composition of TME is still facing enormous challenges, especially the immune and stromal components. In this study, we downloaded the RNA-seq profiles and somatic mutation data of 309 CC cases from The Cancer Genome Atlas (TCGA) database, which were analyzed by integrative bioinformatical methods. Initially, ESTIMATE computational method was employed to calculate the amount of immune and stromal components. Then, based on the high- and low-immunity cohorts, we recognized the differentially expressed genes (DEGs) as well as the differentially mutated genes (DMGs). Additionally, we conducted an intersection analysis of DEGs and DMGs, ultimately determining an immune-related prognostic signature, GTPase, IMAP Family Member 4 (GIMAP4). Moreover, sequential analyses demonstrated that GIMAP4 was a protective factor in CC, positively correlated with the overall survival (OS) and negatively with distant metastasis. Besides, we utilized the Gene Set Enrichment Analysis (GSEA) to explore the enrichment-pathways in high and low-expression cohorts of GIMAP4. The results indicated that the genes of the high-expression cohort had a high enrichment in immune-related biological processes and metabolic activities in the low one. Furthermore, CIBERSORT analysis was applied to evaluate the proportion of tumor-infiltrating immune cells (TICs), illustrating that several activated TICs were strongly associated with GIMAP4 expression, which suggested that GIMAP4 had the potential to be an indicator for the immune state in TME of CC. Hence, GIMAP4 contributed to predicting the CC patients’ clinical outcomes, such as survival rate, distant metastasis and immunotherapy response. Moreover, GIMAP4 could serve as a promising biomarker for TME remodeling, suggesting the possible underlying mechanisms of tumorigenesis and CC progression, which may provide different therapeutic perceptions of CC, and therefore improve treatment.

Prognostic Hub Genes in the Immune Microenvironment of Lung Adenocarcinoma by Estimation

BACKGROUND

The mortality of lung adenocarcinoma(LUAD) is high. Recent studies have found that the degree of immune infiltration and stromal cells in the tumour microenvironment or tumours makes a significant contribution to prognosis.

METHODS

During study, we screened differentially expressed genes (DEGs) of TCGA database for prognostic genes in LUAD immune microenvironment. Further, immune and stromal cells were quantified using ESTIMATE algorithm. To study the effects of immune and stromal cell-associated genes on the prognosis of LUAD, LUAD patients were divided into high and low groups according to their immune/ stromal scores. The obtained scores were found to be related to the phenotype and survival rate of LUAD patients. By selecting DEGs with high expression in immune and stromal cells, we performed functional enrichment analysis and found that most genes are associated with pathways of cancer, stimulus response and the MAPK signaling. The functions and enriched pathways of LUAD prognostic genes were shown by a protein-protein interaction (PPI) network. Nonetheless, an external database was used to validate the prognostic genes from the TCGA.

RESULTS

Prognostic genes were listed according to their expression position and protein function.

CONCLUSION

We provided a new targets for immunotherapy of LUAD, which further provides basic knowledge for future clinical research.

Microenvironment-related gene TNFSF13B predicts poor prognosis in kidney renal clear cell carcinoma

Background

Kidney renal clear cell carcinoma (KIRC) affects the genitourinary system. Although treatment of KIRC in early stages can be highly successful, this type of cancer is difficult to detect until later stages of disease that are less easily treatable. Previous studies have focused on tumor cells, but have ignored the contributions of the tumor microenvironment.

Methods

We analyzed KIRC gene expression data from The Cancer Genome Atlas with the ESTIMATE algorithm to identify differentially expressed genes. Through protein-protein interaction network analysis, we identified clusters and picked genes from the clusters for further analysis. Differential expression, Kaplan-Meier, and univariate Cox analyses were used to select prognostic biomarkers. Gene Set Enrichment Analysis (GSEA) and Tumor Immune Estimation Resource (TIMER) analysis were used to explore the immune characteristics of these genes as biomarkers.

Results

Through the ESTIMATE algorithm and other system biology tools, TNFSF13B was identified as a prognostic biomarker. TNFSF13B is highly expressed in tumors, and high expression of TNFSF13B leads to poor prognosis. Further GSEA and TIMER analysis revealed that the expression of TNFSF13B was related to the immune signaling pathway and lymphocyte infiltration. Our findings strongly suggest that TNFSF13B may be a potential biomarker or target related to the tumor microenvironment for KIRC.

Identification of novel genes in testicular cancer microenvironment based on ESTIMATE algorithm-derived immune scores

Testicular cancer is the most common solid malignancy among young men. We downloaded data of testicular cancer patients from The Cancer Genome Atlas database to find novel genes in the testicular cancer microenviroment based on ESTIMATE algorithm-derived immune scores. A total of 156 cases of testicular cancer were included in this study and 165 cases of normal testicular tissues were used. We divided the testicular cancer patients into high- and low-score groups based on their immune scores. We identified 1,226 differentially expressed genes (fold change > 2, false discovery rate < 0.05), including 688 downregulated genes and 538 upregulated genes, between these two groups. The top Gene Ontology terms were involved in the immune response-regulating cell surface receptor signaling pathway, immune response-activating cell surface receptor signaling pathway, external side of the plasma membrane, and receptor ligand activity. By performing the Kyoto Encyclopedia of Genes and Genomes analysis, we demonstrated that cAMP signaling pathway was highly enriched among these differentially expressed genes. High expression of LINC01564, LINC02208, ODAM, RNA5SP111, and RNU6-196P were found to be associated with poor overall survival. The expression of genes was further validated by the Human Protein Atlas and only ALB and IFNG were demonstrated to be differentially expressed between testis tissue and testicular cancer tissue.

BTK Has Potential to Be a Prognostic Factor for Lung Adenocarcinoma and an Indicator for Tumor Microenvironment Remodeling: A Study Based on TCGA Data Mining

Tumor microenvironment (TME) plays a crucial role in the initiation and progression of lung adenocarcinoma (LUAD); however, there is still a challenge in understanding the dynamic modulation of the immune and stromal components in TME. In the presented study, we applied CIBERSORT and ESTIMATE computational methods to calculate the proportion of tumor-infiltrating immune cell (TIC) and the amount of immune and stromal components in 551 LUAD cases from The Cancer Genome Atlas (TCGA) database. The differentially expressed genes (DEGs) were analyzed by COX regression analysis and protein–protein interaction (PPI) network construction. Then, Bruton tyrosine kinase (BTK) was determined as a predictive factor by the intersection analysis of univariate COX and PPI. Further analysis revealed that BTK expression was negatively correlated with the clinical pathologic characteristics (clinical stage, distant metastasis) and positively correlated with the survival of LUAD patients. Gene Set Enrichment Analysis (GSEA) showed that the genes in the high-expression BTK group were mainly enriched in immune-related activities. In the low-expression BTK group, the genes were enriched in metabolic pathways. CIBERSORT analysis for the proportion of TICs revealed that B-cell memory and CD8+ T cells were positively correlated with BTK expression, suggesting that BTK might be responsible for the preservation of immune-dominant status for TME. Thus, the levels of BTK might be useful for outlining the prognosis of LUAD patients and especially be a clue that the status of TME transition from immune-dominant to metabolic activity, which offered an extra insight for therapeutics of LUAD.

Tumor-Infiltrating Leukocyte Composition and Prognostic Power in Hepatitis B- and Hepatitis C-Related Hepatocellular Carcinomas

Background: Tumor-infiltrating leukocytes (TILs) are immune cells surrounding tumor cells, and several studies have shown that TILs are potential survival predictors in different cancers. However, few studies have dissected the differences between hepatitis B- and hepatitis C-related hepatocellular carcinoma (HBV−HCC and HCV−HCC). Therefore, we aimed to determine whether the abundance and composition of TILs are potential predictors for survival outcomes in HCC and which TILs are the most significant predictors. Methods: Two bioinformatics algorithms, ESTIMATE and CIBERSORT, were utilized to analyze the gene expression profiles from 6 datasets, from which the abundance of corresponding TILs was inferred. The ESTIMATE algorithm examined the overall abundance of TILs, whereas the CIBERSORT algorithm reported the relative abundance of 22 different TILs. Both HBV−HCC and HCV−HCC were analyzed. Results: The results indicated that the total abundance of TILs was higher in non-tumor tissue regardless of the HCC type. Alternatively, the specific TILs associated with overall survival (OS) and recurrence-free survival (RFS) varied between subtypes. For example, in HBV−HCC, plasma cells (hazard ratio [HR] = 1.05; 95% CI 1.00–1.10; p = 0.034) and activated dendritic cells (HR = 1.08; 95% CI 1.01–1.17; p = 0.03) were significantly associated with OS, whereas in HCV−HCC, monocytes (HR = 1.21) were significantly associated with OS. Furthermore, for RFS, CD8+ T cells (HR = 0.98) and M0 macrophages (HR = 1.02) were potential biomarkers in HBV−HCC, whereas neutrophils (HR = 1.01) were an independent predictor in HCV−HCC. Lastly, in both HBV−HCC and HCV−HCC, CD8+ T cells (HR = 0.97) and activated dendritic cells (HR = 1.09) had a significant association with OS, while γ delta T cells (HR = 1.04), monocytes (HR = 1.05), M0 macrophages (HR = 1.04), M1 macrophages (HR = 1.02), and activated dendritic cells (HR = 1.15) were highly associated with RFS. Conclusions: These findings demonstrated that TILs are potential survival predictors in HCC and different kinds of TILs are observed according to the virus type. Therefore, further investigations are warranted to elucidate the role of TILs in HCC, which may improve immunotherapy outcomes.