Exploration of a novel prognostic risk signatures and immune checkpoint molecules in endometrial carcinoma microenvironment

In this study, we devoted to investigate immune-related genes and tumor microenvironment (TME) in EC based on The Cancer Genome Atlas (TCGA) database. In total 799 up-regulated and 139 down-regulated immune-related and differentially expressed genes in EC were investigated for functional annotations and prognosis. By a conjoint Cox regression analysis, we built two risk models for OS and DFS, as well as the consistent nomograms. Immune-related pathways were revealed mostly enriched in the low-risk group. By further analyzing TME based on the risk signatures, the higher immune cell infiltration and activation, lower tumor purity and higher tumor mutational burden were found in low-risk group, which presented a better prognosis. Both the expression and immunophenoscore of immune checkpoints PD-1, CTLA4, PD-L1 and PD-L2 increased significantly in low-risk group. These findings may provide new ideas for novel biomarkers and immunotherapy targets in EC.     

Cancer testis antigen PRAME: An anti-cancer target with immunomodulatory potential

PReferentially expressed Antigen in Melanoma (PRAME) is a cancer testis antigen with restricted expression in somatic tissues and re-expression in poor prognostic solid tumours. PRAME has been extensively investigated as a target for immunotherapy, however, its role in modulating the anti-tumour immune response remains largely unknown. Here, we show that PRAME tumour expression is associated with worse survival in the TCGA breast cancer cohort, particularly in immune-unfavourable tumours. Using direct and indirect co-culture models, we found that PRAME overexpressing MDA-MB-468 breast cancer cells inhibit T cell activation and cytolytic potential, which could be partly restored by silencing of PRAME. Furthermore, silencing of PRAME reduced expression of several immune checkpoints and their ligands, including PD-1, LAG3, PD-L1, CD86, Gal-9 and VISTA. Interestingly, silencing of PRAME induced cancer cell killing to levels similar to anti-PD-L1 atezolizumab treatment. Comprehensive analysis of soluble inflammatory mediators and cancer cell expression of immune-related genes showed that PRAME tumour expression can suppress the expression and secretion of multiple pro-inflammatory cytokines, and mediators of T cell activation, differentiation and cytolysis. Together, our data indicate that targeting of PRAME offers a potential, novel dual therapeutic approach to specifically target tumour cells and regulate immune activation in the tumour microenvironment.     

Gene expression and DNA methylation analyses suggest that immune process-related ADCY6 is a prognostic factor of luminal-like breast cancer

Breast cancer is a malignant tumor that seriously threatens women's health, and luminal-like cancer subtypes account for the majority of the cases. The purpose of this study was to investigate the relationships among DNA methylation, gene expression profile, and the tumor-immune microenvironment of luminal-like breast cancer, and to identify the potential key genes that regulate immune cell infiltration in luminal-like breast cancer. The ESTIMATE algorithm was applied to calculate immune scores and stromal scores of patients with breast cancer. Kaplan-Meier curves were generated for survival analysis. The clusterProfile package was used for Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The protein-protein interaction (PPI) network was constructed using the STRING database and Cytoscape software. Correlations between ADCY6 expression and immune cell infiltration-related pathways were analyzed by gene set variation analysis. R software was used for the statistical analysis and figure generation. Disease-free survival was higher in the immune score-high group than it was in the immune score-low group, while the stromal score had no correlation with prognosis. There were 515 genes that differed in both gene expression and DNA methylation levels, and these genes were mainly enriched in immune process-related pathways. ADCY6 was enriched in module A of the PPI network. Patients with downregulation and hypermethylation of ADCY6 associated with a better prognosis. ADCY6 expression was negatively correlated with the activation of immune process-related signaling pathways, immune checkpoint receptors, and ligands, except for CLEC4G. DNA methylation was found to be involved in the regulation of the key cellular pathways of luminal-like breast cancer immune cell infiltration. Additionally, ADCY6 was identified as a prognostic factor involved in the DNA methylation-regulated immune processes in luminal-like breast cancer.     

TGFβ2 is a prognostic‐related biomarker and correlated with immune infiltrates in gastric cancer

TGFβ2 is an essential regulator of immune cell functionality, but the mechanisms whereby it drives immune infiltration in gastric cancer remain uncertain. The Oncomine and Tumor Immunoassay Resource (TIMER) databases were used for assessing the expression of TGFβ2, after which TIMER was used to explore the relationship between TGFβ2 and tumour immune infiltration. Finally, we assessed how TGFβ2 expression correlated with the expression of a set of marker genes associated with immune infiltration using TIMER and GEPIA. We determined TGFβ2 expression to be significantly correlated with outcome in multiple types of cancer in the Cancer Genome Atlas (TCGA), with the effect being particularly pronounced in gastric cancer. Furthermore, elevated TGFβ2 expression was found to be significantly correlated with gastric cancer N staging, and with the expression of a variety of immune markers associated with particular immune cell subsets. These results indicate that TGFΒ2 is associated with patient outcome and tumour immune cell infiltration in multiple cancer types. This suggests that TGFβ2 is a key factor which governs immune cell recruitment to gastric cancer tumours, potentially playing a vital role in governing immune cell infiltration and thus representing a valuable prognostic biomarker in gastric cancer patients.     

Pan-cancer analysis of the CASP gene family in relation to survival,tumor-infiltrating immune cells and therapeutic targets

The cysteinyl aspartate protease (caspase, or CASP) gene family plays a significant role in programmed cell death, inflammation and immunity. However, the correlation between CASP family members and prognosis and tumor-infiltrating lymphocytes in different tumors has not been determined. We investigated the role of CASP genes in cancer prognosis and their relationship with clinicopathological parameters. We also evaluated the correlation between the expression of CASP family members and cancer immune infiltration and evaluated whether these molecules can be used as targets for immunotherapy. The CASP1/2/4/5/7/9 genes may represent prognostic factors and therapeutic targets for breast cancer, hepatocellular carcinoma and pancreatic cancer. Another finding is that the CASP1/4/5 genes help to regulate innate immunity and T cell immunity and may also have an important effect on tumor checkpoint inhibition. These findings may elucidate the roles played by CASP family members in cancer progression and identify strategies to promote collaborative activities in the context of immunotherapy.     

Immune-related long noncoding RNA signature for predicting survival and immune checkpoint blockade in hepatocellular carcinoma

Long noncoding RNAs (lncRNAs) show multiple functions, including immune response. Recently, the immune-related lncRNAs have been reported in some cancers. We first investigated the immune-related lncRNA signature as a potential target in hepatocellular carcinoma (HCC) survival. The training set (n = 368) and the independent external validation cohort (n = 115) were used. Immune genes and lncRNAs coexpression were constructed to identify immune-related lncRNAs. Cox regression analyses were perfumed to establish the immune-related lncRNA signature. Regulatory roles of this signature on cancer pathways and the immunologic features were investigated. The correlation between immune checkpoint inhibitors and this signature was examined. In this study, the immune-related lncRNA signature was identified in HCC, which could stratify patients into high- and low-risk groups. This immune-related lncRNA signature was correlated with disease progression and worse survival and was an independent prognostic biomarker. Our immune-related lncRNA signature was still a powerful tool in predicting survival in each stratum of age, gender, and tumor stage. This signature mediated cell cycle, glycolysis, DNA repair, mammalian target of rapamycin signaling, and immunologic characteristics (i.e., natural killer cells vs. Th1 cells down, etc). This signature was associated with immune cell infiltration (i.e., macrophages M0, Tregs, CD4 memory T cells, and macrophages M1, etc.,) and immune checkpoint blockade (ICB) immunotherapy-related molecules (i.e., PD-L1, PD-L2, and IDO1). Our findings suggested that the immune-related lncRNA signature had an important value for survival prediction and may have the potential to measure the response to ICB immunotherapy. This signature may guide the selection of the immunotherapy for HCC.     

Functional comparison of innate immune signaling pathways in primates

Humans respond differently than other primates to a large number of infections. Differences in susceptibility to infectious agents between humans and other primates are probably due to inter-species differences in immune response to infection. Consistent with that notion, genes involved in immunity-related processes are strongly enriched among recent targets of positive selection in primates, suggesting that immune responses evolve rapidly, yet providing only indirect evidence for possible inter-species functional differences. To directly compare immune responses among primates, we stimulated primary monocytes from humans, chimpanzees, and rhesus macaques with lipopolysaccharide (LPS) and studied the ensuing time-course regulatory responses. We find that, while the universal Toll-like receptor response is mostly conserved across primates, the regulatory response associated with viral infections is often lineage-specific, probably reflecting rapid host-virus mutual adaptation cycles. Additionally, human-specific immune responses are enriched for genes involved in apoptosis, as well as for genes associated with cancer and with susceptibility to infectious diseases or immune-related disorders. Finally, we find that chimpanzee-specific immune signaling pathways are enriched for HIV-interacting genes. Put together, our observations lend strong support to the notion that lineage-specific immune responses may help explain known inter-species differences in susceptibility to infectious diseases.     

Identification of immune microenvironment subtypes that predicted the prognosis of patients with ovarian cancer

Ovarian cancer (OC) is associated with high mortality rate. However, the correlation between immune microenvironment and prognosis of OC remains unclear. This study aimed to explore prognostic significance of OC tumour microenvironment. The OC data set was selected from the cancer genome atlas (TCGA), and 307 samples were collected. Hierarchical clustering was performed according to the expression of 756 genes. The immune and matrix scores of all immune subtypes were determined, and Kruskal-Wallis test was used to analyse the differences in the immune and matrix scores between OC samples with different immune subtypes. The model for predicting prognosis was constructed based on the expression of immune-related genes. TIDE platform was applied to predict the effect of immunotherapy on patients with OC of different immune subtypes. The 307 OC samples were classified into three immune subtypes A-C. Patients in subtype B had poorer prognosis and lower survival rate. The infiltration of helper T cells and macrophages in microenvironment indicated significant differences between immune subtypes. Enrichment analyses of immune cell molecular pathways showed that JAK–STAT3 pathway changed significantly in subtype B. Furthermore, predictive response to immunotherapy in subtype B was significantly higher than that in subtype A and C. Immune subtyping can be used as an independent predictor of the prognosis of OC patients, which may be related to the infiltration patterns of immune cells in tumour microenvironment. In addition, patients in immune subtype B have superior response to immunotherapy, suggesting that patients in subtype B are suitable for immunotherapy.     

A comprehensive analysis of IDO1 expression with tumour‐infiltrating immune cells and mutation burden in gynaecologic and breast cancers

Gynaecologic and breast cancers share some similarities at the molecular level. The aims of our study are to highlight the similarities and differences about IDO1, an important immune‐related gene in female cancers. The NGS data from TCGA of cervical squamous cell carcinoma (CESC), ovarian serous cystadenocarcinoma (OV), uterine corpus endometrial carcinoma (UCEC), uterine carcinosarcoma (UCS) and breast invasive carcinoma (BRCA) were analysed to identify molecular features, and clinically significant and potential therapeutic targets of IDO1. We found IDO1 was significantly up‐regulated in four gynaecologic cancers and breast cancer. According to breast cancer PAM50 classification scheme, IDO1 expression was higher in tumours of basal than other subtypes and showed better survival prognosis in BRCA and OV. Through immune infiltration analysis, we found a strong correlation between IDO1 and immune cell populations especially for dendritic cells and T cells. In addition, we investigated the association between IDO1 and tumour mutation burden (TMB) and found that IDO1 was significantly correlated with TMB in BRCA and CESC. GSVA revealed that hallmarks significantly correlated with IDO1 were involved in interferon gamma response, allograft rejection and inflammatory response. We also found PD‐L1 and LAG3 were highly positive related to IDO1 in gynaecologic cancers when comparing with their corresponding normal tissues. Our results indicated that IDO1 participated in anti‐tumour immune process and is correlated with mutation burden. These findings may expand our outlook of potential anti‐IDO1 treatments.     

The landscape and prognostic value of immune characteristics in uterine corpus endometrial cancer

In the present study, we explored the clinical and immunological characteristics of 575 uterine corpus endometrial carcinoma (UCEC) samples obtained from The Cancer Genome Atlas (TCGA) using the ESTIMATE and CIBERSORT algorithms. First, Kaplan–Meier and univariate Cox regression analyses indicated that the immune cell score was a prognostic factor for overall survival (OS) and recurrence-free survival (RFS). Multivariate Cox regression analysis further revealed that the immune cell score was an independent prognostic factor for UCEC patients. Second, we investigated the correlation between the infiltration levels of 22 types of immune cells and the immune score. Survival analysis based on the 22 immune cell types showed that higher levels of regulatory T cell, activated NK cell, and follicular helper T-cell infiltration were associated with longer OS, while higher levels of CD8+ T cell and naive B-cell infiltration were associated with longer RFS. Next, we performed differential expression and prognosis analyses on 1534 immune-related genes and selected five from 14 candidate genes to construct a prognostic prediction model. The area under the receiver-operating characteristic (ROC) curve (AUC) for 3- and 5-year survival were 0.711 and 0.728, respectively. Further validation using a stage I–II subgroup showed similar results, presenting AUC values for 3- and five-year survival of 0.677 and 0.692, respectively. Taken together, the present study provides not only a deeper understanding of the relationship between UCEC and the immune landscape but also guidance for the future development of UCEC immunotherapy.     

SLAMF8, a potential new immune checkpoint molecule,is associated with the prognosis of colorectal cancer

Recently, immune checkpoint inhibitors (ICIs), such as programmed cell death 1 (PD-1) monoclonal antibodies (mAbs), have revolutionized the treatment of malignant tumors. Therefore, the number of studies aiming to screen and identify new immune checkpoint molecules for antitumor immunotherapy is increasing. Signaling lymphocytic activation molecule (SLAM) family members are mainly expressed by and regulate the functions of immune cells. Recent studies have shown that several SLAM family members are involved in the regulation of the tumor immune microenvironment and are promising targets for antitumor immunotherapy. Signaling lymphocytic activation molecule family member 8 (SLAMF8) is a type I cell surface glycoprotein and is encoded on chromosome 1q21. To further illustrate the clinical value of SLAMF8 in colorectal cancer (CRC), we retrospectively analyzed the relationship between SLAMF8 expression and the prognosis of CRC patients and the associations between SLAMF8 expression and the expression levels of other SLAM family members and other classic immune checkpoint molecules using The Cancer Genome Atlas (TCGA) data, RNA sequencing data, tissue immunohistochemistry staining, and systematic follow-up analysis. Here, high SLAMF8 expression was associated with poor overall survival (OS) in CRC. The mRNA expression level of SLAMF8 was positively correlated with the expression levels of multiple classic immune checkpoints and other SLAM family members. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis suggested that the pathways enriched in CRC tissues with high SLAMF8 expression were associated with the regulation of the tumor immune microenvironment.     

STAT1 is a key gene in a gene regulatory network related to immune phenotypes in bladder cancer: An integrative analysis of multi-omics data

The immunophenotype of bladder cancer plays a pivotal role in the prognosis of cancer, but the effect of different epigenetic factors on different immunophenotypes in bladder tumours remains unclear. This study used multi-omics data analysis to provide molecular basis support for different immune phenotypes. Unsupervised cluster analysis revealed distinct subclusters with higher (subcluster B2) or lower cytotoxic immune phenotypes (subcluster A1) related to PD-L1 and IFNG expression. Mutational landscape analyses showed that the mutation level of TP53 in subcluster B1 was highest than other subclusters, and subcluster B1 had a lower frequency of concurrent mutation than subcluster A2. A total of 2364 differentially expressed genes were identified between subclusters A2 and B1, and the main functions of the up-regulated genes in subcluster B1 were enriched in the activation of T cells and other related pathways. We found that STAT1 was a key gene in a gene regulatory network related to immune phenotypes in bladder cancer. Finally, we constructed a prognostic prediction model by LASSO Cox regression which could distinguish high-risk and low-risk cases significantly. In conclusion, the present study addressed a field synopsis between genetic and epigenetic events in immune phenotypes of bladder cancer.     

Immunogenomics Analysis Reveals that TP53 Mutations Inhibit Tumor Immunity in Gastric Cancer

Although immunotherapy continues to demonstrate efficacy in a variety of refractory cancers, currently, no any immunotherapeutic strategy is clinically used for gastric cancer (GC) except its microsatellite instable subtype. Thus, it is important to identify molecular biomarkers for predicting the responders to GC immunotherapy. TP53 mutations frequently occur in GC and are associated with unfavorable clinical outcomes in GC. We performed a comprehensive characterization of the associations between TP53 mutations and immune activities in GC based on two large-scale GC cancer genomics data. We compared expression and enrichment levels of 787 immune-related genes and 23 immune gene-sets among TP53-mutated GCs, TP53‐wildtype GCs, and normal tissue, and explored the correlations between p53-mediated pathways and immune activities in GC. Strikingly, almost all analyzed immune gene-sets were significantly downregulated in enrichment levels in TP53-mutated GCs compared to TP53‐wildtype GCs. These less active immune pathways and cell types in TP53-mutated GCs included 15 immune cell types and function, tumor-infiltrating lymphocytes, regulatory T cells, immune checkpoint, cytokine and cytokine receptor, human leukocyte antigen, pro‐inflammatory, and parainflammation. Moreover, we identified a number of p53-mediated pathways and proteins that were significantly associated with immune activities in GC. Furthermore, we demonstrated that the TP53 mutation itself could result in the depressed immune activities in GC and other cancer types. We revealed that chromosomal instability was an important mechanism for the depressed tumor immunity in TP53-mutated cancers. Finally, we showed that immune cell infiltration and immune activities were likely positively associated with survival prognosis in GC. Our findings suggest that p53 may play an important role in activating tumor immunity in GC and other cancer types and that the TP53 mutation status could be useful in stratifying cancer patients responsive to a certain immunotherapy.     

Gene Expression Profile Analysis of Type 2 Diabetic Mouse Liver

Liver plays a key role in glucose metabolism and homeostasis, and impaired hepatic glucose metabolism contributes to the development of type 2 diabetes. However, the precise gene expression profile of diabetic liver and its association with diabetes and related diseases are yet to be further elucidated. In this study, we detected the gene expression profile by high-throughput sequencing in 9-week-old normal and type 2 diabetic db/db mouse liver. Totally 12132 genes were detected, and 2627 genes were significantly changed in diabetic mouse liver. Biological process analysis showed that the upregulated genes in diabetic mouse liver were mainly enriched in metabolic processes. Surprisingly, the downregulated genes in diabetic mouse liver were mainly enriched in immune-related processes, although all the altered genes were still mainly enriched in metabolic processes. Similarly, KEGG pathway analysis showed that metabolic pathways were the major pathways altered in diabetic mouse liver, and downregulated genes were enriched in immune and cancer pathways. Analysis of the key enzyme genes in fatty acid and glucose metabolism showed that some key enzyme genes were significantly increased and none of the detected key enzyme genes were decreased. In addition, FunDo analysis showed that liver cancer and hepatitis were most likely to be associated with diabetes. Taken together, this study provides the digital gene expression profile of diabetic mouse liver, and demonstrates the main diabetes-associated hepatic biological processes, pathways, key enzyme genes in fatty acid and glucose metabolism and potential hepatic diseases.     

Chromosome 9p copy number gains involving PD-L1 are associated with a specific proliferation and immune-modulating gene expression program active across major cancer types

Background

Inhibition of the PD-L1/PD-1 immune checkpoint axis represents one of the most promising approaches of immunotherapy for various cancer types. However, immune checkpoint inhibition is successful only in subpopulations of patients emphasizing the need for powerful biomarkers that adequately reflect the complex interaction between the tumor and the immune system. Recently, recurrent copy number gains (CNG) in chromosome 9p involving PD-L1 were detected in many cancer types including lung cancer, melanoma, bladder cancer, head and neck cancer, cervical cancer, soft tissue sarcoma, prostate cancer, gastric cancer, ovarian cancer, and triple-negative breast cancer.

Methods

Here, we applied functional genomics to analyze global mRNA expression changes associated with chromosome 9p gains. Using the TCGA data set, we identified a list of 75 genes that were strongly up-regulated in tumors with chromosome 9p gains across many cancer types.

Results

As expected, the gene set was enriched for chromosome 9p and in particular chromosome 9p24 (36 genes and 23 genes). Furthermore, we found enrichment of two expression programs derived from genes within and beyond 9p: one implicated in cell cycle regulation (22 genes) and the other implicated in modulation of the immune system (16 genes). Among these were specific cytokines and chemokines, e.g. CCL4, CCL8, CXCL10, CXCL11, other immunoregulatory genes such as IFN-G and IDO1 as well as highly expressed proliferation-related kinases and genes including PLK1, TTK, MELK and CDC20 that represent potential drug targets.

Conclusions

Collectively, these data shed light on mechanisms of immune escape and stimulation of proliferation in cancer with PD-L1 CNG and highlight additional vulnerabilities that may be therapeutically exploitable.

     

Analysis of PD-1 related immune transcriptional profile in different cancer types

Background

Programmed cell death 1 (PD-1) functions as an immune checkpoint in the process of anti-tumor immune response. The PD-1 blockade is now becoming a fundamental part in cancer immunotherapy. So it’s essential to elicit the PD-1 related immune process in different types of cancer.

Methods

The Cancer Genome Atlas was used to collect the RNA-seq data of 33 cancer types. The microenvironment cell populations-counter was used to analyze the immune cell infiltrates. KEGG and GO analysis were performed to investigate PD-1 associated biological process. Kaplan–Meier survival curves and Cox’s proportional hazards model were performed for prognostic value analysis.

Results

We demonstrated that PD-1 expression varied in different cancer types. The uveal melanoma had a low PD-1 expression and poor infiltrated with immune cells. But it showed the strong correlation of PD-1 with the most types of immune cells. The PD-1 demonstrated a robust relationship with other immunomodulators and showed its involvement in critical functions correlated with anti-tumor immune pathways. Survival analysis indicated the PD-1 expression suggested different prognosis in different cancer types.

Conclusions

Our investigations promote a better understanding of the PD-1 blockade and provide PD-1 related personized combined immunotherapy for different types of cancer patients.