Dataset Information

Interferon regulatory factor family influences tumor immunity and prognosis of patients with colorectal cancer.

ABSTRACT:

Background

Since interferon regulatory factor (IRF) family functions in immune response to viral infection, its role in colorectal cancer (CRC) has not been inspected before. This study tries to investigate members of IRF family using bioinformatics approaches in aspect of differential expressions, biological function, tumor immune infiltration and clinical prognostic value for patients with CRC.

Methods

Transcriptome profiles data, somatic mutations and clinical information of CRC were obtained from COAD/READ dataset of The Cancer Genome Atlas (TCGA) as a training set. Gene expression data (GSE17536 and GSE39582) were downloaded from the Gene Expression Omnibus as a validating set. A random forest algorithm was used to score the risk for every case. Analyzing gene and function enrichment, constructing protein-protein interaction and noncoding RNA network, identifying hub-gene, characterizing tumor immune infiltration, evaluating differences in tumor mutational burden (TMB) and sensitivity to chemotherapeutics or immunotherapy were performed by a series of online tools and R packages. Immunohistochemical (IHC) examinations were carried out validation in tissue samples.

Results

Principal-component analysis (PCA) suggested that the transcript expression levels of nine members of IRF family differed between normal colorectum and CRC. The risk score constructed by IRF family not only acted as an independent factor for predicting survival in CRC patients with different biological processes, signaling pathways and TMB, but also indicated different immunotherapy response with diverse immune and stromal cells infiltration. IRF3 and IRF7 were upregulated in CRC and suggested a shorter survival time in patients with CRC. Differentially expressed members of IRF family exhibited varying degrees of immune cell infiltration. IHC analysis showed a positive association between IRF3 and IRF7 expression and tumor-infiltrating immune cells, including CD4⁺ T cell and CD68⁺ macrophages.

Conclusions

On account of differential expression, IRF family members can help to predict both response to immunotherapy and clinical prognosis of patients with CRC. Our bioinformatic investigation not only gives a preliminary picture of the genetic features as well as tumor microenvironment, but it may provide a clue for further experimental exploration and verification on IRF family members in CRC.

SUBMITTER: Chen YJ

PROVIDER: S-EPMC8422700 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:ObjectiveThis study is aimed at investigating the regulating mechanisms of the interferon regulatory factor (IRF) family genes in head and neck squamous cell carcinoma.MethodsBased on the HNSC data in the 'The Cancer Genome Atlas (TCGA)' database, the expression pattern of IRF family genes was investigated. The association of IRFs family genes and survival outcomes were analyzed by Kaplan-Meier plotter web portal. The relation of IRF genes and tumor stages was evaluated by using stage plots and based on GEPIA portal. 50 genes interacting with IRFs were identified using the NetworkAnalyst's protein-protein interaction (PPI) network construction tool. The top 200 correlated genes with similar expression patterns in HNSC were obtained by the similar gene detection module of GEPIA. Furthermore, functional enrichment analysis was performed to determine the biological functions enriched by the interacting and correlated genes. The potential implication of IRFs in tumor immunity was investigated in terms of tumor-infiltrating immune cells, a pair of immune checkpoint genes (CD274 and PDCD1), and ESTIMATE-Stromal-Immune score.ResultsThe unpaired sample analysis shows that all of the IRF family genes were highly expressed in HNSC tumor samples compared to control samples. The survival analysis results showed that the overexpression of IRF1, IRF4, IRF5, IRF6, IRF8, and IRF9 was associated with better overall survival in HNSC, while the other IRFs genes (IRF2, IRF3. and IRF7) did not show prognostic values for overall survival outcome of HNSC. Four genes (STAT1, STAT2, FOXP3, and SPI1) were overlapping among 50 interacted genes in the PPI network and top 200 correlated genes identified by GEPIA. The 50 interacting genes in the PPI network and top 200 correlated genes were integrated into 246 genes. These 246 genes were found to be overrepresented in multiple KEGG pathways, e.g., Th17 cell differentiation, T cell receptor signaling pathway, cytokine-cytokine receptor interaction, natural killer (NK) cell-mediated cytotoxicity, FOXO signaling, PI3K-Akt signaling, and ErbB signaling. Most correlations between IRF gene members and TIICs were positive. The strongest positive correlation was identified between IRF8 and T cells (r = 0.849, p < 0.001). The majority of correlation between IRF family genes and ESTIMATE-Stromal-Immune score was found to be positive. The highest positive correlation was found to be between IRF8 and Immune score (r = 0.874, p = 1.09E - 158). Most correlations between IRFs and two immunoinhibitor genes (CD274 and PDCD1) were positive. IRF1 and PDCD1 were found to show the highest positive correlation (r = 0.764, p < 2.2e - 16).ConclusionsThe current analysis showed IRFs were differentially expressed in HNSC, indicated significant prognostic values, were involved in tumor immunity-related signaling pathways, and significantly regulated tumor-infiltrating immune cells. IRF family genes could be potential therapeutic biomarkers in targeting tumor immunity of head and neck cancer.

Project description:BackgroundImmune related interferon regulatory factor 4 (IRF4) is a member of the IRF family, whereas the clinical significance and possible role of IRF4 in lung adenocarcinoma (LUAD) remains unclear. We aimed to investigate the role of IRF4 in predicting the prognosis of LUAD patients.MethodsUsing The Cancer Genome Atlas (TCGA) database and our immunohistochemical (IHC) cohort, we analyzed the correlation between IRF4 expression and clinical characteristics, and the prognostic value of IRF4 was also evaluated in LUAD. The potential biological functions of IRF4 in LUAD were analyzed by Gene Set Enrichment Analysis (GSEA). The relationship between IRF4 and immune cell infiltration were evaluated by TISIDB database and our own IHC cohort. In addition, an immune checkpoint inhibitor (ICI) treated cohort from Gene Expression Omnibus database was used to determine the role of IRF4 in LUAD patients with immunotherapy.ResultsWe found that either mRNA or protein expression level of IRF4 was significantly higher in LUAD than in normal tissues (P < 0.001). The elevate in IRF4 expression in LUAD was significantly associated with the earlier clinical stage (P = 0.002). Patients with LUAD and IRF4 high expression correlated with significant longer overall survival in both TCGA database (P < 0.05) and our IHC-cohort (P = 0.001). Our results also demonstrated that IRF4 could serve as an independent favorable prognostic factor in patients with LUAD. GSEA analysis indicated that high IRF4 expression group enriched with several immune-related pathways, such as B cell receptor signaling pathway, T cell receptor signaling pathway and cytokine-cytokine receptor interaction signaling pathway. In LUAD, IRF4 positively correlated with several different immune infiltrations including various B cells, CD8+ T cells and CD4+ T cells both in mRNA and protein levels. Additionally, we found that the expression of IRF4 was positively associated with PD-1 and PD-L1 mRNA expression levels, and IRF4 high expression predicted moderate better survival in LUAD with immunotherapy (P = 0.071).ConclusionsOur results suggested that IRF4 was associated with higher B cells and T cells infiltration levels and might be a favorable prognostic biomarker in LUAD patients, whereas the potential prognostic role of IRF4 in ICI-treated patients needed further exploration.

Project description:Background: Glioma is the most common primary tumor of the central nervous system. The conventional glioma treatment strategies include surgical excision and chemo- and radiation-therapy. Interferon Gamma (IFN-γ) is a soluble dimer cytokine involved in immune escape of gliomas. In this study, we sought to identify IFN-γ-related genes to construct a glioma prognostic model to guide its clinical treatment. Methods: RNA sequences and clinicopathological data were downloaded from The Cancer Genome Atlas (TCGA) and the China Glioma Genome Atlas (CGGA). Using univariate Cox analysis and the Least Absolute Shrinkage and Selection Operator (LASSO) regression algorithm, IFN-γ-related prognostic genes were selected to construct a risk scoring model, and analyze its correlation with the clinical features. A high-precision nomogram was drawn to predict prognosis, and its performance was evaluated using calibration curve. Finally, immune cell infiltration and immune checkpoint molecule expression were analyzed to explore the tumor microenvironment characteristics associated with the risk scoring model. Results: Four out of 198 IFN-γ-related genes were selected to construct a risk score model with good predictive performance. The expression of four IFN-γ-related genes in glioma tissues was significantly increased compared to normal brain tissue (p < 0.001). Based on ROC analysis, the risk score model accurately predicted the overall survival rate of glioma patients at 1 year (AUC: The Cancer Genome Atlas 0.89, CGGA 0.59), 3 years (AUC: TCGA 0.89, CGGA 0.68), and 5 years (AUC: TCGA 0.88, CGGA 0.70). Kaplan-Meier analysis showed that the overall survival rate of the high-risk group was significantly lower than that of the low-risk group (p < 0.0001). Moreover, high-risk scores were associated with wild-type IDH1, wild-type ATRX, and 1P/19Q non-co-deletion. The nomogram predicted the survival rate of glioma patients based on the risk score and multiple clinicopathological factors such as age, sex, pathological grade, and IDH Status, among others. Risk score and infiltrating immune cells including CD8 T-cell, resting CD4 memory T-cell, regulatory T-cell (Tregs), M2 macrophages, resting NK cells, activated mast cells, and neutrophils were positively correlated (p < 0.05). In addition, risk scores closely associated with expression of immune checkpoint molecules such as PD-1, PD-L1, CTLA-4, LAG-3, TIM-3, TIGIT, CD48, CD226, and CD96. Conclusion: Our risk score model reveals that IFN-γ -associated genes are an independent prognostic factor for predicting overall survival in glioma, which is closely associated with immune cell infiltration and immune checkpoint molecule expression. This model will be helpful in predicting the effectiveness of immunotherapy and survival rate in patients with glioma.

Project description:This manuscript presents the first extensive phylogenetics analysis of a key family of immune regulators, the interferon regulatory factor (IRF) family. The IRF family encodes transcription factors that play important roles in immune defense, stress responses, reproduction, development, and carcinogenesis. Several times during their evolution, the IRF genes have undergone expansion and diversification. These genes were also completely lost on two separate occasions in large groups of metazoans. The origin of the IRF family coincides with the appearance of multicellularity in animals. IRF genes are present in all principal metazoan groups, including sea sponges, placozoans, comb jellies, cnidarians, and bilaterians. Although the number of IRF family members does not exceed two in sponges and placozoans, this number reached five in cnidarians. At least four additional independent expansions lead up to 11 members in different groups of bilaterians. In contrast, the IRF genes either disappeared or mutated beyond recognition in roundworms and insects, the two groups that include most of the metazoan species. The IRF family separated very early into two branches ultimately leading to vertebrate IRF1 and IRF4 supergroups (SGs). Genes encoding the IRF-SGs are present in all bilaterians and cnidarians. The evolution of vertebrate IRF family members further proceeded with at least two additional steps. First, close to the appearance of the first vertebrate, the IRF family probably expanded to four family members, predecessors of the four vertebrate IRF groups (IRF1, 3, 4, 5 groups). In the second step, 10 vertebrate family members evolved from these four genes, likely as a result of the 2-fold duplication of the entire genome. Interestingly, the IRF family coevolved with the Rel/NF-kappaB family with which it shares some important evolutionary characteristics, including roles in defense responses, metazoan specificity, extensive diversification in vertebrates, and elimination of all family members in nematodes.