A novel ferroptosis-related microRNA signature with prognostic value in osteosarcoma

The induction of ferroptosis is suggested to be a potential therapeutic strategy for cancers. MicroRNAs (miRNAs) are reported to play an important role in cell death processes. This study aims to construct and validate a risk model based on ferroptosis-related miRNAs (FR_miRNAs) to predict prognosis and identify novel therapeutic targets for osteosarcoma. Data from the Therapeutically Applicable Research to Generate Effective Treatments database are used as the training cohort. A prognostic signature based on two FR_miRNAs (miR-635 and miR-593) is developed using univariate Cox regression, least absolute shrinkage and selection operator regression, and multivariate Cox regression analyses. The area under the curve values of the prognostic signature to predict the 1-year, 2-year, 3-year, and 5-year overall survival rates in patients with osteosarcoma are 0.782, 0.781, 0.722, and 0.777, respectively, indicating a good predictive ability. Based on the risk score, patients are divided into low-risk and high-risk groups. Patients with high-risk scores are associated with poor survival. The risk level is determined to be an independent prognostic factor. A nomogram is established for predicting prognosis. The expression levels of PRNP (miR-635-related ferroptosis-related gene (FRG); P=0.024) and HILPDA (miR-593-related FRG; P=0.025) are significantly different between the low-risk and high-risk groups. All results are validated in an external cohort (GSE39040). The results of the functional assay reveal that miR-635 mimics inhibit osteosarcoma (OS) cell proliferation and migration, whereas miR-593 overexpression exerts the opposite effect. In conclusion, miR-635 and miR-593 exert contrasting regulatory effects on OS cell proliferation and migration.

Development and validation of a novel necroptosis-related gene signature for predicting prognosis and therapeutic response in Ewing sarcoma

Ewing sarcoma (ES) is the second most common malignant bone tumor in children and has a poor prognosis due to early metastasis and easy recurrence. Necroptosis is a newly discovered cell death method, and its critical role in tumor immunity and therapy has attracted widespread attention. Thus, the emergence of necroptosis may provide bright prospects for the treatment of ES and deserves our further study. Here, based on the random forest algorithm, we identified 6 key necroptosis-related genes (NRGs) and used them to construct an NRG signature with excellent predictive performance. Subsequent analysis showed that NRGs were closely associated with ES tumor immunity, and the signature was also good at predicting immunotherapy and chemotherapy response. Next, a comprehensive analysis of key genes showed that RIPK1, JAK1, and CHMP7 were potential therapeutic targets. The Cancer Dependency Map (DepMap) results showed that CHMP7 is associated with ES cell growth, and the Gene Set Cancer Analysis (GSCALite) results revealed that the JAK1 mutation frequency was the highest. The expression of 3 genes was all negatively correlated with methylation and positively with copy number variation (CNV). Finally, an accurate nomogram was constructed with this signature and clinical traits. In short, this study constructed an accurate prognostic signature and identified 3 novel therapeutic targets against ES.

A novel signature incorporating lipid metabolism- and immune-related genes to predict the prognosis and immune landscape in hepatocellular carcinoma

Background

Liver hepatocellular carcinoma (LIHC) is a highly malignant tumor with high metastasis and recurrence rates. Due to the relation between lipid metabolism and the tumor immune microenvironment is constantly being elucidated, this work is carried out to produce a new prognostic gene signature that incorporates immune profiles and lipid metabolism of LIHC patients.

Methods

We used the "DEseq2" R package and the "Venn" R package to identify differentially expressed genes related to lipid metabolism (LRDGs) in LIHC. Additionally, we performed unsupervised clustering of LIHC patients based on LRDGs to identify their subgroups and immuno-infiltration and Gene Ontology (GO) enrichment analysis on the subgroups. Next, we employed multivariate, LASSO and univariate Cox regression analyses to determine variables and to create a prognostic profile on the basis of immune- and lipid metabolism-related differential genes (IRDGs and LRDGs). We separated patients into low- and high-risk groups in accordance with the best cut-off value of risk score. We conducted Decision Curve Analysis (DCA), Receiver Operating Characteristic curve analysis as a function of time as well as Survival Analysis to evaluate this signature's prognostic value. We incorporated the clinical characteristics of patients into the risk model to obtain a nomogram prognostic model. GEO14520 and ICGC-LIRI JP datasets were employed to externally confirm the accuracy and robustness of signature. The gene set variation analysis (GSVA) and gene set enrichment analysis (GSEA) were applied for investigating the underlying mechanisms. Immune infiltration analysis was implemented to examine the differences in immune between both risk groups. Single-cell RNA sequencing (scRNA-SEQ) was utilized to characterize the genes that were involved in the distribution of signature and expression characteristics of different LIHC cell types. The patients' sensitivity in both risk groups to commonly used chemotherapeutic agents and semi-inhibitory concentrations (IC50) of the drugs was assessed using the GDSC database. On the basis of the differentially expressed genes (DEGs) in the two groups, the CMAP database was adopted for the prediction of potential small-molecule compounds. Small-molecule compounds were molecularly docked with prognostic markers. Lastly, we investigated the prognostic gene expression levels in normal and LIHC tissues with immunohistochemistry (IHC) and quantitative reverse transcription polymerase chain reaction(qRT-PCR).

Results

We built and verified a prognostic signature with seven genes that incorporated immune profiles and lipid metabolism. Patients were classified as low- and high-risk groups depending on their prognostic profiles. The overall survival (OS) was markedly lower in the high-risk group as compared to low-risk group. Time-dependent ROC curves more precisely predicted patients' survival at 1, 3 and 5 years; the area under the ROC curve was 0.81 (1 year), 0.75 (3 years) and 0.77 (5 years). The DCA curves showed the value of the prognostic genes in this signature for clinical applications. We included the patients' clinical characteristics in the risk model for both multivariate and univariate Cox regression analyses, and the findings revealed that the risk model represents an independent factor that influences OS in LIHC patients. With immune analysis, GSVA and GSEA, we identified that there are remarkable differences between the two risk groups in immune pathways, lipid metabolism, tumor development, immune cell infiltration and immune microenvironment, response to immunotherapy, and sensitivity to chemotherapy. Moreover, those with higher risk scores presented greater sensitivity to the chemotherapeutic agents. Experiments in vitro further elucidated the roles of SPP1 and FLT3 in the LIHC immune microenvironment. Furthermore, four small-molecule drugs that could target LIHC were screened. In vitro qRT-PCR , IHC revealed that the SPP1,KIF18A expressions were raised in LIHC in tumor samples, whereas FLT3,SOCS2 showed the opposite trend.

Conclusions

We developed and verified a new signature comprising immune- and lipid metabolism-associated markers and to assess the prognosis and the immune status of LIHC patients. This signature can be applied to survival prediction, individualized chemotherapy, and immunotherapeutic guidance for patients with liver cancer. This study also provides potential targeted therapeutics and novel ideas for the immune evasion and progression of LIHC.

Development and validation of a cuproptosis-associated prognostic model for diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is a highly heterogeneous disease. Therefore, more reliable biomarkers are required to better predict the prognosis of DLBCL. Cuproptosis is a novel identified form of programmed cell death (PCD) that is different from oxidative stress-related cell death (e.g., apoptosis, ferroptosis, and necroptosis) by Tsvetkov and colleagues in a recent study released in Science. Cuproptosis is copper-dependent PCD that is closely tied to mitochondrial metabolism. However, the prognostic value of cuproptosis-related genes (CRGs) in DLBCL remains to be further elucidated. In the present study, we systematically evaluated the molecular changes of CRGs in DLBCL and found them to be associated with prognosis. Subsequently, based on the expression profiles of CRGs, we characterized the heterogeneity of DLBCL by identifying two distinct subtypes using consensus clustering. Two isoforms exhibited different survival, biological functions, chemotherapeutic drug sensitivity, and immune microenvironment. After identifying differentially expressed genes (DEGs) between CRG clusters, we built a prognostic model with the Least absolute shrinkage and selection operator (LASSO) Cox regression analysis and validated its prognostic value by Cox regression analysis, Kaplan-Meier curves, and receiver operating characteristic (ROC) curves. In addition, the risk score can predict clinical characteristics, levels of immune cell infiltration, and prognosis. Furthermore, a nomogram incorporating clinical features and risk score was generated to optimize risk stratification and quantify risk assessment. Compared to the International Prognostic Index (IPI), the nomogram has demonstrated more accuracy in survival prediction. Furthermore, we validated the prognostic gene expression levels through external experiments. In conclusion, cuproptosis-related gene signature can serve as a potential prognostic predictor in DLBCL patients and may provide new insights into cancer therapeutic targets.

The signature of pyroptosis-related gene prognostic and immune microenvironment in adrenocortical carcinoma

Adrenocortical carcinoma (ACC) has a low incidence but a poor prognosis. And ACC has complex clinical manifestations and limited treatment. Pyroptosis has a dual character and has both positive and negative effects on cancer. However, the role of pyroptosis-related genes (PRGs) in ACC and the impact on ACC progression remains unelucidated. This study performed systematic bioinformatics analysis and basic experimental validation to enable the establishment of prognostic models and demonstrate levels of immune infiltration. Pearson's correlation analysis was used to assess the association of PRGs with tumor immune infiltration, tumor mutation burden (TMB), microsatellite instability (MSI), and immune checkpoints. There 4 PRGs were upregulated, and 25 PRGs were downregulated in ACC. At the same time, we analyzed and reviewed the genetic mutation variation landscape of PRGs. Functional enrichment analysis was also performed to clarify the function of PRGs. Pyroptosis, the inflammatory response, the Toll-like receptor signaling pathway, and the NOD-like receptor signaling pathway are the functions and pathways mainly involved and exerted effects by these 33 PRGs. The results of the prognosis analysis revealed high expression of CASP3, CASP9, GSDMB, GSDMD, NLRC4, PRKACA, and SCAF11 caused a poor survival rate for ACC patients. The above seven PRGs were screened by the optimal λ value of LASSO Cox analysis, and the five selected genes (CASP3, CASP9, GSDMB, GSDMD, NLRC4) were involved in constructing a prognostic PRGs model which enables the overall survival in ACC patients can be predicted with moderate to high accuracy. Prognostic PRGs, especially CASP9, which is the independent factor of ACC prognosis, may be closely correlated with immune-cell infiltration, tumor mutation burden, microsatellite instability, and immune checkpoints. Quantitative Real-Time PCR (qRT-PCR), Western blot and immunohistochemical were performed to validate the mRNA expression levels of CASP9 in adjacent normal tissues and ACC tissues. According to the result of immune checkpoints analysis, NLRC4 and GSDMB may be identified as potential therapeutic targets. In conclusion, we established a prognostic model of PRG characteristics in ACC and analyzed the relationship between PRGs and immune infiltration. Through our study, it may be helpful to find the mechanism of pyroptosis in ACC.

Identification and validation of a novel cuproptosis-related gene signature in multiple myeloma

Background: Cuproptosis is a newly identified unique copper-triggered modality of mitochondrial cell death, distinct from known death mechanisms such as necroptosis, pyroptosis, and ferroptosis. Multiple myeloma (MM) is a hematologic neoplasm characterized by the malignant proliferation of plasma cells. In the development of MM, almost all patients undergo a relatively benign course from monoclonal gammopathy of undetermined significance (MGUS) to smoldering myeloma (SMM), which further progresses to active myeloma. However, the prognostic value of cuproptosis in MM remains unknown. Method: In this study, we systematically investigated the genetic variants, expression patterns, and prognostic value of cuproptosis-related genes (CRGs) in MM. CRG scores derived from the prognostic model were used to perform the risk stratification of MM patients. We then explored their differences in clinical characteristics and immune patterns and assessed their value in prognosis prediction and treatment response. Nomograms were also developed to improve predictive accuracy and clinical applicability. Finally, we collected MM cell lines and patient samples to validate marker gene expression by quantitative real-time PCR (qRT-PCR). Results: The evolution from MGUS and SMM to MM was also accompanied by differences in the CRG expression profile. Then, a well-performing cuproptosis-related risk model was developed to predict prognosis in MM and was validated in two external cohorts. The high-risk group exhibited higher clinical risk indicators. Cox regression analyses showed that the model was an independent prognostic predictor in MM. Patients in the high-risk group had significantly lower survival rates than those in the low-risk group (p < 0.001). Meanwhile, CRG scores were significantly correlated with immune infiltration, stemness index and immunotherapy sensitivity. We further revealed the close association between CRG scores and mitochondrial metabolism. Subsequently, the prediction nomogram showed good predictive power and calibration. Finally, the prognostic CRGs were further validated by qRT-PCR in vitro. Conclusion: CRGs were closely related to the immune pattern and self-renewal biology of cancer cells in MM. This prognostic model provided a new perspective for the risk stratification and treatment response prediction of MM patients.

A Novel Gene Signature Associated with Inflammatory Responses and Immune Status Assists in Prognosis and Intervention for Patients with HCC

BACKGROUND

Tumor growth depends on tumor cells and the tumor microenvironment, which are regulated by inflammation and immune responses. However, the roles of inflammation and immune status in hepatocellular carcinoma (HCC) remain unclear. The aim of this study was to evaluate the prognostic value of an inflammatory response- related gene signature associated with immune status, which may provide insight into new treatment options for HCC patients.

MATERIALS AND METHODS

Differentially expressed genes associated with inflammation were obtained from The Cancer Genome Atlas (TCGA), the Gene Expression Omnibus, and the Molecular Signatures Database. An inflammation-associated prognostic gene signature was constructed and validated using TCGA and the International Cancer Genome Consortium datasets, respectively, using LASSO Cox regression analysis. Log-rank was performed to compare the overall survival of low- and high-risk score cohorts. Immune cell infiltration and immune-related functions were analyzed using single-sample gene enrichment analysis. The structures of the drugs identified by the prognostic model were predicted using PubChem. The drugs sensitivity of bleomycin, simvastatin and zoledronate detected by CCK8 colorimetric assay. The mRNA levels of 7 genes in HCC after drug treatment analyzed via qRT-PCR.

RESULTS

Inflammation-associated genes, including ITGA5, MEP1A, P2RX4, RIPK2, SLC7A1 and SRI, were identified and found to be associated with the prognosis of HCC. We further found that the high-risk patients experienced poor prognosis, which was observed to be an independent and significant risk factor for prognosis. Moreover, we observed elevated expression levels in multiple immune cell types and immune function. Lastly, we validated that bleomycin, simvastatin and zoledronate could regulate these genes in HCC.

CONCLUSION

The inflammatory-response-associated gene signature could predict the prognosis and the immunological status of HCC patients. Additionally, bleomycin, simvastatin and zoledronate may represent potential drug candidates that could inhibit these genes. This may constitute a new approach for the treatment of HCC.

An inflammatory response-related gene signature associated with immune status and prognosis of acute myeloid leukemia

OBJECTIVE

To determine the prognostic significance of inflammatory response-associated genes in acute myeloid leukemia (AML).

METHODS

Transcriptomic profiles and related clinical information of AML patients were acquired from a public database. To establish a multi-gene prognosis signature, we performed least absolute shrinkage and selection operator Cox analysis for the TCGA cohort and evaluated the ICGC cohort for verification. Subsequently, Kaplan-Meier analysis was carried out to compare the overall survival (OS) rates between high- and low-risk groups. Biological function and single-sample gene set enrichment (ssGSEA) analyses were employed to investigate the association of risk score with immune status and the tumor microenvironment. Prognostic gene expression levels in AML samples and normal controls were confirmed by qRT-PCR and immunofluorescence.

RESULTS

We identified a potential inflammatory response-related signature comprising 11 differentially expressed genes, including ACVR2A, CCL22, EBI3, EDN1, FFAR2, HRH1, ICOSLG, IL-10, INHBA, ITGB3, and LAMP3, and found that AML patients with high expression levels in the high-risk group had poor OS rates. Biological function analyses revealed that prognostic genes mainly participated in inflammation and immunity signaling pathways. Analyses of cancer-infiltrating immunocytes indicated that in high-risk patients, the immune suppressive microenvironment was significantly affected. The expression of the inflammation reaction-associated signature was found to be associated with susceptibility to chemotherapy. There was a significant difference in prognostic gene expression between AML and control tissues.

CONCLUSION

A novel inflammatory response-related signature was developed with 11 candidate genes to predict prognosis and immune status in AML patients.

A Multi-Omics Network of a Seven-Gene Prognostic Signature for Non-Small Cell Lung Cancer

There is an unmet clinical need to identify patients with early-stage non-small cell lung cancer (NSCLC) who are likely to develop recurrence and to predict their therapeutic responses. Our previous study developed a qRT-PCR-based seven-gene microfluidic assay to predict the recurrence risk and the clinical benefits of chemotherapy. This study showed it was feasible to apply this seven-gene panel in RNA sequencing profiles of The Cancer Genome Atlas (TCGA) NSCLC patients (n = 923) in randomly partitioned feasibility-training and validation sets (p < 0.05, Kaplan-Meier analysis). Using Boolean implication networks, DNA copy number variation-mediated transcriptional regulatory network of the seven-gene signature was identified in multiple NSCLC cohorts (n = 371). The multi-omics network genes, including PD-L1, were significantly correlated with immune infiltration and drug response to 10 commonly used drugs for treating NSCLC. ZNF71 protein expression was positively correlated with epithelial markers and was negatively correlated with mesenchymal markers in NSCLC cell lines in Western blots. PI3K was identified as a relevant pathway of proliferation networks involving ZNF71 and its isoforms formulated with CRISPR-Cas9 and RNA interference (RNAi) profiles. Based on the gene expression of the multi-omics network, repositioning drugs were identified for NSCLC treatment.

Immunogenomic Landscape and Immune-Related Gene-Based Prognostic Signature in Asian Gastric Cancer

Background

Asians have the highest incidence of gastric cancer (GC), and the prognosis of Asian GC is poor. Furthermore, the therapeutics for Asian GC is limited because of genetic heterogeneity and screening difficulty at the early stage. This study aimed to develop an immune-related gene (IRG)-based prognostic signature and to explore prognosis-related regulatory mechanism and therapeutic target for Asian GC.

Methods

To elucidate the prognostic value of IRGs in Asian GC, a comprehensive analysis of IRG expression profiles and overall survival times in 364 Asian GC patients from the Asian Cancer Research Group (ACRG) and The Cancer Genome Atlas (TCGA) databases was performed, and a novel prognostic index was established. To further explore regulatory prognosis mechanisms and therapeutic targets, a tumor immunogenomic landscape analysis, including stromal and immune subcomponents, cell types, panimmune gene sets, and immunomodulatory genes, was performed.

Result

Our analysis allowed the creation of an optimal risk assessment model, the Asian-specific IRG-based prognostic index (ASIRGPI), which showed a high accuracy in predicting survival in Asian GC. We also developed an ASIRGPI-based nomogram to predict the 3- and 5-year overall survival (OS) of Asian GC patients. The impact of the ASIRGPI on the worse prognosis of Asian GC was possibly related to the stromal component remodeling. Specifically, TGFβ gene sets were significantly associated with the ASIRGPI and worse prognosis. Immunomodulatory gene analysis further revealed that TGFβ1 and EDNRB may be the novel potential therapeutic targets for Asian GC.

Conclusions

As a tumor microenvironment-relevant gene set-based prognostic signature, the ASIRGPI model provides an effective approach for evaluating the prognosis of Asian GC and may even prolong OS by enabling the selection of individualized therapy with the novel targets.

Identification of a prognostic chemoresistance-related gene signature associated with immune microenvironment in breast cancer

Breast cancer is the most common form of cancer among women globally, and chemoresistance is a major challenge to disease treatment that is associated with a poor prognosis. This study was formulated to identify a reliable prognostic biosignature capable of predicting the survival of patients with chemoresistant breast cancer (CRBC) and evaluating the associated tumor immune microenvironment. Through a series of protein-protein interaction and weighted correlation network analyses, genes that were significantly associated with breast cancer chemoresistance were identified. Moreover, univariate Cox regression and lasso-penalized Cox regression analyses were employed to generate a prognostic model, and the prognostic utility of this model was then assessed using time-dependent receiver operating characteristic (ROC) and Kaplan-Meier survival curves. Finally, The CIBERSORT and ESTIMATE algorithms were additionally leveraged to assess relationships between the tumor immune microenvironment and patient prognostic signatures. Overall, a multigenic prognostic biosignature capable of predicting CRBC patient risk was successfully developed based on bioinformatics analysis and in vitro experiments. This biosignature was able to stratify CRBC patients into high- and low-risk subgroups. ROC curves also revealed that this biosignature achieved high diagnostic efficiency, and multivariate regression analyses indicated that this risk signature was an independent risk factor linked to CRBC patient outcomes. In addition, this signature was associated with the infiltration of the tumor microenvironment by multiple immune cell types. In conclusion, the chemoresistance-associated prognostic gene signature developed herein was able to effectively evaluate the prognosis of CRBC patients and to reflect the overall composition of the tumor immune microenvironment.

A Robust Prognostic Gene Signature Based on eRNAs-Driven Genes in Prostate Cancer

Prostate cancer (PCa) is the second most common malignancy in men, but its exact pathogenetic mechanisms remain unclear. This study explores the effect of enhancer RNAs (eRNAs) in PCa. Firstly, we screened eRNAs and eRNA -driven genes from The Cancer Genome Atlas (TCGA) database, which are related to the disease-free survival (DFS) of PCa patients;. screening methods included bootstrapping, Kaplan–Meier (KM) survival analysis, and Pearson correlation analysis. Then, a risk score model was established using multivariate Cox analysis, and the results were validated in three independent cohorts. Finally, we explored the function of eRNA-driven genes through enrichment analysis and analyzed drug sensitivity on datasets from the Genomics of Drug Sensitivity in Cancer database. We constructed and validated a robust prognostic gene signature involving three eRNA-driven genes namely MAPK15, ZNF467, and MC1R. Moreover, we evaluated the function of eRNA-driven genes associated with tumor microenvironment (TME) and tumor mutational burden (TMB), and identified remarkable differences in drug sensitivity between high- and low-risk groups. This study identified a prognostic gene signature, which provides new insights into the role of eRNAs and eRNA-driven genes while assisting clinicians to determine the prognosis and appropriate treatment options for patients with PCa.

In-Silico Evaluation of Genetic Alterations in Ovarian Carcinoma and Therapeutic Efficacy of NSC777201, as a Novel Multi-Target Agent for TTK, NEK2, and CDK1

Ovarian cancer is often detected at the advanced stages at the time of initial diagnosis. Early-stage diagnosis is difficult due to its asymptomatic nature, where less than 30% of 5-year survival has been noticed. The underlying molecular events associated with the disease’s pathogenesis have yet to be fully elucidated. Thus, the identification of prognostic biomarkers as well as developing novel therapeutic agents for targeting these markers become relevant. Herein, we identified 264 differentially expressed genes (DEGs) common in four ovarian cancer datasets (GSE14407, GSE18520, GSE26712, GSE54388), respectively. We constructed a protein-protein interaction (PPI) interaction network with the overexpressed genes (72 genes) and performed gene enrichment analysis. In the PPI networks, three proteins; TTK Protein Kinase (TTK), NIMA Related Kinase 2 (NEK2), and cyclin-dependent kinase (CDK1) with higher node degrees were further evaluated as therapeutic targets for our novel multi-target small molecule NSC777201. We found that the upregulated DEGs were enriched in KEGG and gene ontologies associated with ovarian cancer progression, female gamete association, otic vesicle development, regulation of chromosome segregation, and therapeutic failure. In addition to the PPI network, ingenuity pathway analysis also implicate TTK, NEK2, and CDK1 in the elevated salvage pyrimidine and pyridoxal pathways in ovarian cancer. The TTK, NEK2, and CDK1 are over-expressed, demonstrating a high frequency of genetic alterations, and are associated with poor prognosis of ovarian cancer cohorts. Interestingly, NSC777201 demonstrated anti-proliferative and cytotoxic activities (GI₅₀ = 1.6 µM~1.82 µM and TGI₅₀ = 3.5 µM~3.63 µM) against the NCI panels of ovarian cancer cell lines and exhibited a robust interaction with stronger affinities for TTK, NEK2, and CDK1, than do the standard drug, paclitaxel. NSC777201 displayed desirable properties of a drug-like candidate and thus could be considered as a novel small molecule for treating ovarian carcinoma.

Immune Landscape of the Tumor Microenvironment Identifies Prognostic Gene Signature CD4/CD68/CSF1R in Osteosarcoma

Background: Osteosarcoma (OSA), the most common primary bone malignancy in children and adolescents, is prone to metastases and unfavorable prognosis. Owing to its strong genomic heterogeneity, traditional chemotherapy, or targeted immunotherapy has not effectively improved the related overall survival for decades. Since the landscape of the OSA tumor immune microenvironment is scarcely known, despite it playing a crucial role in predicting clinical outcomes and therapeutic efficacies, we aimed to elucidate its molecular characteristics.

Methods: The immune signature of 101 OSA samples was explored using transcriptome profiling and clinical characteristics retrieved from the Therapeutically Applicable Research to Generate Effective Treatments (TARGET) program. Correlations between the prognostic immune markers and their clinical chemotherapy responses were assessed and verified based on 45 OSA primary tumors.

Findings: We identified the heterogeneity underlying tumor immune signature in OSA, and found CD4+ T cells and macrophage markers CD4/IFNGR2/CD68 to be feasible prognostic factors, exerting significantly positive correlation with each other. Specifically, CSF1R, which plays an essential role in the regulation of proliferation and differentiation of macrophages, was found to be a specific signature associated with CD4/CD68, with improved OSA clinical outcomes.

Interpretation: The immune landscape based on CD4/CD68/CSF1R gene signatures showed considerable promise for prognostic and therapeutic stratification in OSA patients. A specific immune signature for OSA, abundantly consisting of Th1-polarized CD4+ T cells and CSF1R-related CD68+ macrophages, may improve the predictive efficacy of chemotherapy and improve prognosis in patients with OSA.

Inflammation and lymphocyte infiltration are associated with shorter survival in patients with high-grade glioma

Glioma represents a serious health burden in terms of morbidity and mortality. The prognostic significance of the lymphoid and myeloid infiltrates in glioma is not clearly determined. Moreover, the characterization of different leukocyte subsets in the tumor microenvironment relies mainly on immunohistochemistry observations, and data about their association with prognosis are contradictory. Here, we performed acomprehensive study of both the tumor-infiltrating and circulating immune compartments of patients with high-grade glioma. Nineteen tumor biopsies and 30 PBMC samples were analyzed by RNA sequencing. Validation was performed on The Cancer Genome Atlas (TCGA) RNA sequencing data from glioma and on additional 39 tumor biopsies analyzed by flow cytometry. We identified prognostic tumor and peripheral immune signatures, which associate increased inflammation, immune infiltration and activation with shorter overall survival in high-grade glioma patients. Importantly, we confirmed our observations by flow cytometry analysis and validated the tumor-signature using the TCGA dataset. In addition, both tumor genotype and grade associated with the degree of glioma immune infiltration. Unlike in the majority of cancers, lymphocyte infiltration at the tumor site is anegative prognostic factor in glioma, suggesting the ambivalent pro-tumorigenic role of immune responses in glioma.

Characterization of a Myeloid Activation Signature that Correlates with Survival in Melanoma Patients

Understanding the cellular interactions within the tumor microenvironment (TME) of melanoma paved the way for novel therapeutic modalities, such as T cell-targeted immune checkpoint inhibitors (ICI). However, only a limited fraction of patients benefits from such therapeutic modalities, highlighting the need for novel predictive and prognostic biomarkers. As myeloid cells orchestrate the tumor-specific immune response and influence the efficacy of ICI, assessing their activation state within the TME is of clinical relevance. Here, we characterized a myeloid activation (MA) signature, comprising the three genes Cxcl11, Gbp1, and Ido1, from gene expression data of human myeloid cells stimulated with poly(I:C) or cGAMP. This MA signature positively correlated to overall survival in melanoma. In addition, increased expression of the MA signature was observed in melanoma patients responding to ICI (anti-PD-1), as compared to non-responders. Furthermore, the MA signature was validated in the murine B16F10 melanoma model where it was induced and associated with decreased tumor growth upon intratumoral administration of poly(I:C) and cGAMP. Finally, we were able to visualize co-expression of the MA signature genes in myeloid cells of human melanoma tissues using RNAscope in situ hybridization. In conclusion, the MA signature indicates the activation state of myeloid cells and represents a prognostic biomarker for the overall survival in melanoma patients.

Transposon insertional mutagenesis in mice identifies human breast cancer susceptibility genes and signatures for stratification

Despite concerted efforts to identify causal genes that drive breast cancer (BC) initiation and progression, we have yet to establish robust signatures to stratify patient risk. Here we used in vivo transposon-based forward genetic screening to identify potentially relevant BC driver genes. Integrating this approach with survival prediction analysis, we identified six gene pairs that could prognose human BC subtypes into high-, intermediate-, and low-risk groups with high confidence and reproducibility. Furthermore, we identified susceptibility gene sets for basal and claudin-low subtypes (21 and 16 genes, respectively) that stratify patients into three relative risk subgroups. These signatures offer valuable prognostic insight into the genetic basis of BC and allow further exploration of the interconnectedness of BC driver genes during disease progression.

Robust prognostic gene signatures and therapeutic targets are difficult to derive from expression profiling because of the significant heterogeneity within breast cancer (BC) subtypes. Here, we performed forward genetic screening in mice using Sleeping Beauty transposon mutagenesis to identify candidate BC driver genes in an unbiased manner, using a stabilized N-terminal truncated β-catenin gene as a sensitizer. We identified 134 mouse susceptibility genes from 129 common insertion sites within 34 mammary tumors. Of these, 126 genes were orthologous to protein-coding genes in the human genome (hereafter, human BC susceptibility genes, hBCSGs), 70% of which are previously reported cancer-associated genes, and ∼16% are known BC suppressor genes. Network analysis revealed a gene hub consisting of E1A binding protein P300 (EP300), CD44 molecule (CD44), neurofibromin (NF1) and phosphatase and tensin homolog (PTEN), which are linked to a significant number of mutated hBCSGs. From our survival prediction analysis of the expression of human BC genes in 2,333 BC cases, we isolated a six-gene-pair classifier that stratifies BC patients with high confidence into prognostically distinct low-, moderate-, and high-risk subgroups. Furthermore, we proposed prognostic classifiers identifying three basal and three claudin-low tumor subgroups. Intriguingly, our hBCSGs are mostly unrelated to cell cycle/mitosis genes and are distinct from the prognostic signatures currently used for stratifying BC patients. Our findings illustrate the strength and validity of integrating functional mutagenesis screens in mice with human cancer transcriptomic data to identify highly prognostic BC subtyping biomarkers.