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.     

Identification of immune-enhanced molecular subtype associated with BRCA1 mutations,immune checkpoints and clinical outcome in ovarian carcinoma

Ovarian carcinoma has the highest mortality among the malignant tumours in gynaecology, and new treatment strategies are urgently needed to improve the clinical status of ovarian carcinoma patients. The Cancer Genome Atlas (TCGA) cohort were performed to explore the immune function of the internal environment of tumours and its clinical correlation with ovarian carcinoma. Finally, four molecular subtypes were obtained based on the global immune-related genes. The correlation analysis and clinical characteristics showed that four subtypes were all significantly related to clinical stage; the immune scoring results indicated that most immune signatures were upregulated in C3 subtype, and the majority of tumour-infiltrating immune cells were upregulated in both C3 and C4 subtypes. Compared with other subtypes, C3 subtype had a higher BRCA1 mutation, higher expression of immune checkpoints, and optimal survival prognosis. These findings of the immunological microenvironment in tumours may provide new ideas for developing immunotherapeutic strategies for ovarian carcinoma.     

The joint role of methylation and immune-related lncRNAs in ovarian cancer: Defining molecular subtypes and developing prognostic signature

IntroductionComplex outcome of ovarian cancer (OC) stems from the tumor immune microenvironment (TIME) influenced by genetic and epigenetic factors. This study aimed to comprehensively explored the subclasses of OC through lncRNAs related to both N6-methyladenosine (m6A)/N1-methyladenosine (m1A)/N7-methylguanosine (m7G)/5-methylcytosine (m5C) in terms of epigenetic variability and immune molecules and develop a new set of risk predictive systems.Material and methodsThe lncRNA data of OC were collected from TCGA. Spearman correlation analysis on lncRNA data of OC with immune-related gene expression and with m6A/m5C/m1A/m7G were respectively conducted. The m6A/m5C/m1A/m7G-related m6A/m5C/m1A/m7G related immune lncRNA subtypes were identified on the basis of the prognostic lncRNAs. Heterogeneity among subtypes was evaluated by tumor mutation analysis, tumor microenvironment (TME) component analysis, response to immune checkpoint blocked (ICB) and chemotherapeutic drugs. A risk predictive system was developed based on the results of Cox regression analysis and random survival forest analysis of the differences between each specific cluster and other clusters.ResultsThree m6A/m5C/m1A/m7G-related immune lncRNA subtypes of OC showing distinct differences in prognosis, mutation pattern, TIME components, immunotherapy and chemotherapy response were identified. A set of risk predictive system consisting of 10 lncRNA for OC was developed, according to which the risk score of samples in each OC dataset was calculated and risk type was defined.ConclusionsThis study classified three m6A/m5C/m1A/m7G-related immune lncRNA subtypes with distinct heterogeneous mutation patterns, TME components, ICB therapy and immune response, and provided a set of risk predictive system consisted of 10 lncRNA for OC.     

Multi-Omics Analysis Elucidates The Immune And Intratumor Microbes Characteristics Of Ubiquitination Subtypes In Lung Adenocarcinoma

Ubiquitination modification is closely related to cancer and participates in the regulation of tumor microenvironment. However, the role of ubiquitination modification in the immune response and prognosis of lung adenocarcinoma has not been elucidated. This study aims to establish a disease classification associated with ubiquitination and reveal the landscape of intratumor microbes in patients with lung adenocarcinoma for the first time. A total of 1314 patients with lung adenocarcinoma in the GEO and TCGA databases were included in our study. We constructed a ubiquitination scoring model using WGCNA and constructed ubiquitination subtypes using unsupervised clustering, analyzed the clinical characteristics, immune characteristics, and intratumor microbes characteristics, and screened out the relevant gene signatures, which were verified by RT-qPCR in human cancer cells. The results showed that the high ubiquitination subtype had poor prognosis, low degree of immune infiltration, high index of tumor stemness, and poor effect of immunotherapy. The subtypes with lower ubiquitination scores have better prognosis, higher tumor microenvironment score and better immunotherapy effect. The C2 subtype has high level of immune infiltration, lower intratumor microbes diversity and abundance, and good prognosis. The C3 subtype has low level of immune infiltration, higher intratumor microbes diversity and abundance, and poor prognosis. The C1 subtype has characteristics between C2 and C3. In summary, this paper constructs a scoring system and several subtypes based on ubiquitination genes, and analyzed the characteristics, which can help provide new methods for clinical treatment.     

Characterization of the relationship between FLI1 and immune infiltrate level in tumour immune microenvironment for breast cancer

Breast cancer is the most common cancer and the leading cause of cancer death among women in the world. Tumour-infiltrating lymphocytes were defined as the white blood cells left in the vasculature and localized in tumours. Recently, tumour-infiltrating lymphocytes were found to be associated with good prognosis and response to immunotherapy in tumours. In this study, to examine the influence of FLI1 in immune system in breast cancer, we interrogated the relationship between the FLI1 expression levels with infiltration levels of 28 immune cell types. By splitting the breast cancer samples into high and low expression FLI1 subtypes, we found that the high expression FLI1 subtype was enriched in many immune cell types, and the up-regulated differentially expressed genes between them were enriched in immune system processes, immune-related KEGG pathways and biological processes. In addition, many important immune-related features were found to be positively correlated with the FLI1 expression level. Furthermore, we found that the FLI1 was correlated with the immune-related genes. Our findings may provide useful help for recognizing the relationship between tumour immune microenvironment and FLI1, and may unravel clinical outcomes and immunotherapy utility for FLI1 in breast cancer.     

A six-gene-based signature for breast cancer radiotherapy sensitivity estimation

Breast cancer (BRCA) represents the most common malignancy among women worldwide with high mortality. Radiotherapy is a prevalent therapeutic for BRCA that with heterogeneous effectiveness among patients. Here, we proposed to develop a gene expression-based signature for BRCA radiotherapy sensitivity estimation. Gene expression profiles of BRCA samples from the Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC) were obtained and used as training and independent testing dataset, respectively. Differential expression genes (DEGs) in BRCA samples compared with their paracancerous samples in the training set were identified by using the edgeR Bioconductor package. Univariate Cox regression analysis and LASSO Cox regression method were applied to screen optimal genes for constructing a radiotherapy sensitivity estimation signature. Nomogram combining independent prognostic factors was used to predict 1-, 3-, and 5-year OS of radiation-treated BRCA patients. Relative proportions of tumor infiltrating immune cells (TIICs) calculated by CIBERSORT and mRNA levels of key immune checkpoint receptors was adopted to explore the relation between the signature and tumor immune response. As a result, 603 DEGs were obtained in BRCA tumor samples, six of which were retained and used to construct the radiotherapy sensitivity prediction model. The signature was proved to be robust in both training and testing sets. In addition, the signature was closely related to the immune microenvironment of BRCA in the context of TIICs and immune checkpoint receptors’ mRNA levels. In conclusion, the present study obtained a radiotherapy sensitivity estimation signature for BRCA, which should shed new light in clinical and experimental research.     

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.     

Identification of autophagy-related genes signature predicts chemotherapeutic and immunotherapeutic efficiency in bladder cancer (BLCA)

Autophagy maintains cellular homeostasis by degrading and recycling cytoplasmic components under stress conditions, which is identified to be involved in tumorigenesis and now has been recognized as novel target in cancer treatment. In present study, we gathered total autophagy-related genes and established an autophagy-related genes signature (ATGRS) through LASSO cox regression analysis in BLCA. Kaplan-Meier survival and multivariate cox regression analyses both showed the ATGRS was a robust independent prognostic factor with high accuracy. Subsequently, integrated analyses indicated that ATGRS had a strong correlation with molecular subtypes, clinicopathological characteristics and somatic mutation alteration. Moreover, ATGRS was found to be positively correlated with the infiltration of immune cells in tumour microenvironment (TME) and immune checkpoint expression, indicating the potent role of autophagy by regulating the TME. In addition, ATGRS was proved to be efficient in predicting the clinical benefit of immune checkpoint inhibitors (ICIs) based immunotherapy and chemotherapy in BLCA. Furthermore, we observed abnormal expression levels of autophagy-related genes and found the different behaviour of ATGRS in pancancer by LASSO cox regression analysis. Therefore, construction of ATGRS in BLCA could help us to interpret the underlying mechanism of autophagy and sheds a light on the clinical application for a combination of autophagy modification with targeted immunotherapy and chemotherapy in BLCA.     

Clinical and molecular characteristics of bladder urothelial carcinoma subtypes

Bladder urothelial carcinoma (BLCA) is a common malignancy with high heterogeneity. A reasonable molecular subtyping can facilitate biological study and personalized therapy of BLCA. In this study, unsupervised consensus clustering was used to acquire the molecular subtypes of BLCA based on messenger RNA (mRNA) and microRNA (miRNA) data. Gene signature markers and canonical signaling pathways were compared between different subtypes. The Database for Annotation, Visualization and Integrated Discovery (DAVID) was used for the functional annotation of overexpressed genes in different subtypes for Gene Ontology (GO) terms and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. Three molecular subtypes were identified including C1 (luminal-P53 like), C2 (luminal-other), and C3 (basal-immune-squamous). C2 was different from C1 and C3 in clinical characteristics, including younger, better prognosis, and a higher proportion of papillary, Asian, low-grade, early-stage, lymph node negative, and complete remission patients (P < 0.05). Three molecular subtypes also showed distinct mRNA and miRNA expression patterns. luminal and P53-like markers were highly expressed in subtype C1, luminal markers were highly expressed in subtype C2, and basal, EMT/claudin-low, immune and squamous-differentiation markers were highly expressed in subtype C3. In addition, highly expressed genes in C1 were involved in extracellular signal transduction and cell-cell interaction, highly expressed genes in C2 were associated with oxygen transport, energy and steroid metabolism, and highly expressed genes in C3 were related with inflammatory, immune, cytokine, and signal transduction. BLCA in different molecular subtypes showed different clinical and molecular characteristics and personalized therapy needed to be adopted according to the molecular subtypes.     

Comprehensive analysis reveals a prognostic and therapeutic biomarker CD3D in the breast carcinoma microenvironment

Breast carcinoma (BRCA) is the most common carcinoma among women worldwide. Despite the great progress achieved in early detection and treatment, morbidity and mortality rates remain high. In the present study, we make a systematic analysis of BRCA using TCGA database by applying CIBERSORT and ESTIMATE computational methods, uncovered CD3D as a prognostic biomarker by intersection analysis of univariate COX and protein–protein interaction (PPI). It revealed that high CD3D expression was strongly associated with poor survival of BRCA, based on The Cancer Genome Atlas (TCGA) database and online websites. Gene Set Enrichment Analysis (GSEA) revealed that the high CD3D expression group was mainly enriched for the immune-related pathways and the low CD3D expression group was mainly enriched for metabolic-related activities. Based on CIBERSORT analysis, the difference test and correlation test suggested that CD3D had a strong correlation with T cells, particularly CD8 + T cells, which indicated that CD3D up-regulation may increase T cell immune infiltration in the TME and induce antitumor immunity by activating T lymphocytes. Furthermore, the correlation analysis showed that CD3D expression had a strongly positive correlation with immune checkpoints, which indicating that the underlying mechanism involves CD3D mediated regulation of T cell functions in BRCA, and single cell RNA-seq analysis revealed that CD3D correlate with CD8 + T cells and it is itself highly expressed in CD8 + T cells. In summary, we identified a prognostic biomarker CD3D in BRCA, which was associated with lymphocyte infiltration, immune checkpoints and could be developed for innovative therapeutics of BRCA.     

Increased microglia/macrophage gene expression in a subset of adult and pediatric astrocytomas

Glioblastoma (GBM) is a highly malignant brain tumor with a dismal prognosis. Gene expression profiling of GBM has revealed clinically relevant tumor subtypes, and this provides exciting opportunities to better understand disease pathogenesis. Results from an increasing number of studies demonstrate a role for the immune response in cancer progression, yet it is unclear how the immune response differs across tumor subtypes and how it affects outcome. Utilizing gene expression data from The Cancer Genome Atlas Project and the Gene Expression Omnibus database, we demonstrate an enrichment of immune response-related gene expression in the mesenchymal subtype of adult GBM (n = 173) and pediatric high-grade gliomas (n = 53). In an independent cohort of pediatric astrocytomas (n = 24) from UCSF, we stratified tumors into subtypes and confirmed these findings. Using novel immune cell-specific gene signatures we demonstrate selective enrichment of microglia/macrophage-related genes in adult and pediatric GBM tumors of the mesenchymal subtype. Furthermore, immunostaining of adult GBM tumors showed significantly higher cell numbers of microglia/macrophages in mesenchymal versus non-mesenchymal tumors (p = 0.04). Interestingly, adult GBM tumors with the shortest survival had significant enrichment of microglia/macrophage-related genes but this was not true for pediatric GBMs. Consistent with an association with poor outcome, immune response-related genes were highly represented in an adult poor prognosis gene signature, with the expression of genes related to macrophage recruitment and activation being most strongly associated with survival (p<0.05) using CoxBoost multivariate modeling. Using a microglia/macrophage high gene signature derived from quantification of tumor-infiltrating cells in adult GBM, we identified enrichment of genes characteristic of CD4 T cells, granulocytes, and microglia/macrophages (n = 573). These studies support a role for the immune response, particularly the microglia/macrophage response, in the biology of an important subset of GBM. Identification of this subset may be important for future therapeutic stratification.     

Systematic analysis on multiple Gene Expression Omnibus data sets reveals fierce immune response in hepatitis B virus‐related acute liver failure

Acute liver failure (ALF) caused by hepatitis B virus (HBV) is common type of liver failure in the world, with high morbidity and mortality rates. However, the prevalence, genetic background and factors determining the development of HBV‐related ALF are rarely studied. In this study, we examined three Gene Expression Omnibus (GEO) data sets by bioinformatics analysis to identify differentially expressed genes (DEGs), key biological processes and pathways. Immune infiltration analysis showed high immune cells infiltration in HBV‐related ALF tissue. We then confirmed natural killer cells and macrophages infiltration in clinical samples by immunohistochemistry assay, implying these cells play a significant role in HBV‐ALF. We found 1277 genes were co‐up‐regulated and that 1082 genes were co‐down‐regulated in the 3 data sets. Inflammation‐related pathways were enriched in the co‐up‐regulated genes and synthetic metabolic pathways were enriched in the co‐down‐regulated genes. WGCNA also revealed a key module enriching in immune inflammation response and identified 10 hub genes, differentially expressed in an independent data set. In conclusion, we identified fierce immune inflammatory response to elucidate the immune‐driven mechanism of HBV‐ALF and 10 hub genes based on gene expression profiles.     

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.     

Functional characterization of a novel BRCA1-null ovarian cancer cell line in response to ionizing radiation

The breast and ovarian cancer susceptibility gene BRCA1 plays a major role in the DNA damage response pathway. The lack of well-characterized human BRCA1-null cell lines has limited the investigation of BRCA1 function, particularly with regard to its role in ovarian cancer. We propagated a novel BRCA1-null human ovarian cancer cell line UWB1.289 from a tumor of papillary serous histology, the most common form of ovarian carcinoma. UWB1.289 carries a germline BRCA1 mutation within exon 11 and has a deletion of the wild-type allele. UWB1.289 is estrogen and progesterone receptor negative and has an acquired somatic mutation in p53, similar to the commonly used BRCA1-null breast cancer cell line HCC1937. We used ionizing radiation to induce DNA damage in both UWB1.289 and in a stable UWB1.289 line in which wild-type BRCA1 was restored. We examined several responses to DNA damage in these cell lines, including sensitivity to radiation, cell cycle checkpoint function, and changes in gene expression using microarray analysis. We observed that UWB1.289 is sensitive to ionizing radiation and lacks cell cycle checkpoint functions that are a normal part of the DNA damage response. Restoration of wild-type BRCA1 function in these cells partially restores DNA damage responses. Expression array analysis not only supports this partial functional correction but also reveals interesting new information regarding BRCA1-positive regulation of the expression of claudin 6 and other metastasis-associated genes and negative regulation of multiple IFN-inducible genes.     

Codon optimization of the tat antigen of human immunodeficiency virus type 1 generates strong immune responses in mice following genetic immunization

DNA vaccines have been successful in eliciting potent immune responses in mice. Their efficiency, however, is restricted in larger animals. One reason for the limited performance of the DNA vaccines is the lack of molecular strategies to enhance immune responses. Additionally, genes directly cloned from pathogenic organisms may not be efficiently translated in a heterologous host expression system as a consequence of codon bias. To evaluate the influence of codon optimization on the immune response, we elected to use the Tat antigens of human immunodeficiency virus type 1 (HIV-1) (subtype C) and HIV-2, as these viral antigens are poorly immunogenic in natural infection and in experimental immunization and they are functionally important in viral infectivity and pathogenesis. Substituting codons that are optimally used in the mammalian system, we synthetically assembled Tat genes and compared them with the wild-type counterparts in two different mouse strains. Codon-optimized Tat genes induced qualitatively and quantitatively superior immune responses as measured in a T-cell proliferation assay, enzyme-linked immunospot assay, and chromium release assay. Importantly, while the wild-type genes promoted a mixed Th1-Th2-type cytokine profile, the codon-optimized genes induced a predominantly Th1 profile. Using a pepscan strategy, we mapped an immunodominant T-helper epitope to the core and basic domains of HIV-1 Tat. We also identified cross-clade immune responses between HIV-1 subtype B and C Tat proteins mapped to this T-helper epitope. Developing molecular strategies to optimize the immunogenicity of DNA vaccines is critical for inducing strong immune responses, especially to antigens like Tat. Our identification of a highly conserved T-helper epitope in the first exon of HIV-1 Tat of subtype C and the demonstration of a cross-clade immune response between subtypes B and C are important for a more rational design of an HIV vaccine.     

Identification of two molecular subtypes of dysregulated immune lncRNAs in ovarian cancer

Long non-coding RNA (lncRNA) has increasingly been identified as a key regulator in pathologies such as cancer. Multiple platforms were used for comprehensive analysis of ovarian cancer to identify molecular subgroups. However, lncRNA and its role in mapping the ovarian cancer subpopulation are still largely unknown. RNA-sequencing and clinical characteristics of ovarian cancer were acquired from The Cancer Genome Atlas database (TCGA). A total of 52 lncRNAs were identified as aberrant immune lncRNAs specific to ovarian cancer. We redefined two different molecular subtypes, C1(188) and C2(184 samples), in “iClusterPlus” R package, among which C2 grouped ovarian cancer samples have higher survival probability and longer median survival time (P <0.05) with activated IFN-gamma response, Wound Healing and Cytotoxic lymphocytes signal; 456 differentially expressed genes were acquired in C1 and C2 subtypes using limma (3.40.6) package, among which 419 were up-regulated and 37 were down-regulated, in TCGA dataset. Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) functional enrichment analysis revealed that these genes were actively involved in ECM-receptor interaction, PI3K-Akt signaling pathway interaction KEGG pathway. Compared with the existing immune subtype, the Cluster2 sample showed a substantial increase in the proportion of the existing C2 immune subtype, accounting for 81.37%, which was associated with good prognosis. Our C1 subtype contains only 56.49% of the existing immune C1 and C4, which also explains the poor prognosis of C1. Furthermore, 52 immune-related lncRNAs were used to divide the TCGA-endometrial cancer and cervical cancer samples into two categories, and C2 had a good prognosis. The differentially expressed genes were highly correlated with immune-cell-related pathways. Based on lncRNA, two molecular subtypes of ovarian cancer were identified and had significant prognostic differences and immunological characteristics.