CD74 regulates complexity of tumor cell HLA class II peptidome in brain metastasis and is a positive prognostic marker for patient survival

Comprehensive Analysis Reveals Prognostic and Therapeutic Immunity-Related Biomarkers for Pediatric Metastatic Osteosarcoma

Background and Objectives: Osteosarcoma, the most prevalent malignant bone tumor in children and adolescents, presents a complex pathogenesis characterized by various genetic and epigenetic alterations. This study aims to identify key differentially expressed genes (DEGs) in pediatric osteosarcoma, with a focus on those influencing metastasis and patient survival. Materials and Methods: We utilized the GSE33382 dataset from the GEO database for a comprehensive bioinformatic analysis. This included a protein-protein interaction (PPI) network analysis, Cox regression, and Kaplan-Meier survival analysis to identify central DEGs associated with osteosarcoma metastasis and patient survival. Results: Our analysis identified 88 DEGs related to osteosarcoma metastasis. Among them, three survival-related central DEGs-C1QA, CD74, and HLA-DMA-were significantly linked to patient outcomes. Further correlation analysis established a strong relationship between these genes, tumor mutation burden (TMB), immune checkpoint gene expression, and overall survival. Notably, C1QA and CD74 exhibited higher expression in non-metastatic osteosarcoma cases, suggesting a potential role in disease progression. Conclusions: The identified DEGs, particularly C1QA, CD74, and HLA-DMA, may serve as critical biomarkers for pediatric osteosarcoma prognosis and potential targets for immunotherapy. These findings provide a deeper understanding of the molecular landscape of osteosarcoma and open new avenues for therapeutic intervention.

Identification of invasive subpopulations using spatial transcriptome analysis in thyroid follicular tumors

BACKGROUND

Follicular tumors include follicular thyroid adenomas and carcinomas; however, it is difficult to distinguish between the two when the cytology or biopsy material is obtained from a portion of the tumor. The presence or absence of invasion in the resected material is used to differentiate between adenomas and carcinomas, which often results in the unnecessary removal of the adenomas. If nodules that may be follicular thyroid carcinomas are identified preoperatively, active surveillance of other nodules as adenomas is possible, which reduces the risk of surgical complications and the expenses incurred during medical treatment. Therefore, we aimed to identify biomarkers in the invasive subpopulation of follicular tumor cells.

METHODS

We performed a spatial transcriptome analysis of a case of follicular thyroid carcinoma and examined the dynamics of CD74 expression in 36 cases.

RESULTS

We identified a subpopulation in a region close to the invasive area, and this subpopulation expressed high levels of CD74. Immunohistochemically, CD74 was highly expressed in the invasive and peripheral areas of the tumor.

CONCLUSIONS

Although high CD74 expression has been reported in papillary and anaplastic thyroid carcinomas, it has not been analyzed in follicular thyroid carcinomas. Furthermore, the heterogeneity of CD74 expression in thyroid tumors has not yet been reported. The CD74-positive subpopulation identified in this study may be useful in predicting invasion of follicular thyroid carcinomas.

Murine Bone Marrow Erythroid Cells Have Two Branches of Differentiation Defined by the Presence of CD45 and a Different Immune Transcriptome Than Fetal Liver Erythroid Cells

Mouse erythropoiesis is a multifaceted process involving the intricate interplay of proliferation, differentiation, and maturation of erythroid cells, leading to significant changes in their transcriptomic and proteomic profiles. While the immunoregulatory role of murine erythroid cells has been recognized historically, modern investigative techniques have been sparingly applied to decipher their functions. To address this gap, our study sought to comprehensively characterize mouse erythroid cells through contemporary transcriptomic and proteomic approaches. By evaluating CD71 and Ter-119 as sorting markers for murine erythroid cells and employing bulk NanoString transcriptomics, we discerned distinctive gene expression profiles between bone marrow and fetal liver-derived erythroid cells. Additionally, leveraging flow cytometry, we assessed the surface expression of CD44, CD45, CD71, and Ter-119 on normal and phenylhydrazine-induced hemolytic anemia mouse bone marrow and splenic erythroid cells. Key findings emerged: firstly, the utilization of CD71 for cell sorting yielded comparatively impure erythroid cell populations compared to Ter-119; secondly, discernible differences in immunoregulatory molecule expression were evident between erythroid cells from mouse bone marrow and fetal liver; thirdly, two discrete branches of mouse erythropoiesis were identified based on CD45 expression: CD45-negative and CD45-positive, which had been altered differently in response to phenylhydrazine. Our deductions underscore (1) Ter-119's superiority over CD71 as a murine erythroid cell sorting marker, (2) the potential of erythroid cells in murine antimicrobial immunity, and (3) the importance of investigating CD45-positive and CD45-negative murine erythroid cells separately and in further detail in future studies.

The Immunopeptidome from a Genomic Perspective: Establishing the Noncanonical Landscape of MHC Class I-Associated Peptides

Tumor antigens can emerge through multiple mechanisms, including translation of noncoding genomic regions. This noncanonical category of tumor antigens has recently gained attention; however, our understanding of how they recur within and between cancer types is still in its infancy. Therefore, we developed a proteogenomic pipeline based on deep learning de novo mass spectrometry (MS) to enable the discovery of noncanonical MHC class I-associated peptides (ncMAP) from noncoding regions. Considering that the emergence of tumor antigens can also involve posttranslational modifications (PTM), we included an open search component in our pipeline. Leveraging the wealth of MS-based immunopeptidomics, we analyzed data from 26 MHC class I immunopeptidomic studies across 11 different cancer types. We validated the de novo identified ncMAPs, along with the most abundant PTMs, using spectral matching and controlled their FDR to 1%. The noncanonical presentation appeared to be 5 times enriched for the A03 HLA supertype, with a projected population coverage of 55%. The data reveal an atlas of 8,601 ncMAPs with varying levels of cancer selectivity and suggest 17 cancer-selective ncMAPs as attractive therapeutic targets according to a stringent cutoff. In summary, the combination of the open-source pipeline and the atlas of ncMAPs reported herein could facilitate the identification and screening of ncMAPs as targets for T-cell therapies or vaccine development.

SPASCER: spatial transcriptomics annotation at single-cell resolution

In recent years, the explosive growth of spatial technologies has enabled the characterization of spatial heterogeneity of tissue architectures. Compared to traditional sequencing, spatial transcriptomics reserves the spatial information of each captured location and provides novel insights into diverse spatially related biological contexts. Even though two spatial transcriptomics databases exist, they provide limited analytical information. Information such as spatial heterogeneity of genes and cells, cell-cell communication activities in space, and the cell type compositions in the microenvironment are critical clues to unveil the mechanism of tumorigenesis and embryo differentiation. Therefore, we constructed a new spatial transcriptomics database, named SPASCER (https://ccsm.uth.edu/SPASCER), designed to help understand the heterogeneity of tissue organizations, region-specific microenvironment, and intercellular interactions across tissue architectures at multiple levels. SPASCER contains datasets from 43 studies, including 1082 sub-datasets from 16 organ types across four species. scRNA-seq was integrated to deconvolve/map spatial transcriptomics, and processed with spatial cell-cell interaction, gene pattern and pathway enrichment analysis. Cell-cell interactions and gene regulation network of scRNA-seq from matched spatial transcriptomics were performed as well. The application of SPASCER will provide new insights into tissue architecture and a solid foundation for the mechanistic understanding of many biological processes in healthy and diseased tissues.

Immune related biomarkers for cancer metastasis to the brain

Brain metastasis accounts for a large number of cancer-related deaths. The host immune system, involved at each step of the metastatic cascade, plays an important role in both the initiation of the brain metastasis and their treatment responses to various modalities, through either local and or systemic effect. However, few reliable immune biomarkers have been identified in predicting the development and the treatment outcome in patients with cancer brain metastasis. Here, we provide a focused perspective of immune related biomarkers for cancer metastasis to the brain and a thorough discussion of the potential utilization of specific biomarkers such as tumor mutation burden (TMB), genetic markers, circulating and tumor-infiltrating immune cells, cytokines, in predicting the brain disease progression and regression after therapeutic intervention. We hope to inspire the field to extend the research and establish practical guidelines for developing and validating immune related biomarkers to provide personalized treatment and improve treatment outcomes in patients with metastatic brain cancers.

The expression and methylation of PITX genes is associated with the prognosis of head and neck squamous cell carcinoma

Background: The PITX gene family, comprising PITX1, PITX2, and PITX3, is critical in organogenesis and has been evolutionary conserved in animals. PITX genes are associated with the advanced progression and poor prognosis of multiple cancers. However, the relationship between the PITX genes and head and neck squamous cell carcinoma (HNSC) has not been reported.

Methods: We used data from The Cancer Genome Atlas (TCGA) to analyze the association between PITX mRNA expression and clinicopathological parameters of patients with HNSC. The prognostic value of PITX genes was evaluated using the Kaplan-Meier plotter. Multivariate Cox analysis was used to screen out prognosis-associated genes to identify better prognostic indicators. The potential roles of PITX1 and PITX2 in HNSC prognosis were investigated using the protein-protein interaction (PPI) network, Gene Ontology (GO) analysis, and the Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis. The correlation between PITX1 and PITX2 expression or methylation and immune cell infiltration was evaluated using the tumor-immune system interaction database (TISIDB). MethSurv was used to identify DNA methylation and its effect on HNSC prognosis.

Results:PITX genes expression was correlated with different cancers. PITX1 and PITX2 expression was lower in the patients with HNSC. In HNSC, PITX1 expression was significantly related to the clinical stage, histologic grade, and N stage, while PITX2 expression was only significantly related to the histologic grade. The high expression of PITX3 was significantly related to the histologic grade, T stage, and N stage. Survival analysis revealed that PITX genes had prognostic value in HNSC, which was supported by multivariate Cox analysis. PPI network and enrichment analysis showed that the genes interacting with PITX1 and PITX2 belonged predominantly to signaling pathways associated with DNA binding and transcription. Of the CpG DNA methylation sites in PITX1 and PITX2, 28 and 22 were related to the prognosis of HNSC, respectively. Additionally, PITX1 and PITX2 expression and methylation was associated with tumor-infiltrating lymphocytes (TILs).

Conclusion: The PITX genes were differentially expressed in patients with HNSC, highlighting their essential role in DNA methylation and tumor-infiltrating immune cell regulation, as well as overall prognostic value in HNSC.

c-Rel-dependent monocytes are potent immune suppressor cells in cancer

Myeloid-derived suppressor cells (MDSCs) are a heterogeneous population of leukocytes that are important for tumorigenesis and tumor immunotherapy. They comprise up to 10% of leukocytes in the blood of tumor patients and their depletion may be required for successful tumor immunotherapy. However, the identity of MDSCs remains obscure, primarily due to their heterogeneity and lack of a known lineage-specific transcription factor specifying their differentiation. Using single-cell transcriptomics and gene knockout approaches, we now describe a subset of murine and human myeloid suppressor cells, named rel-dependent monocytes (rMos), which are programmed by the transcription factor c-Rel of the NF-κB family. Unlike MDSCs described previously, the c-Rel-dependent monocytes expressed a high amount of the proinflammatory cytokine IL-1β together with a low level of suppressive molecule arginase 1. Both in vitro and in tumor-bearing mice, these c-Rel⁺ IL-1β^hi Arg1^- monocytes promoted tumor growth by potently suppressing T cell function and showed a strong migratory phenotype, all of which were impaired by c-Rel deficiency or inhibition. Mechanistic studies revealed that c-Rel controlled the expression of monocyte signature genes through a unique transcriptional complex called the c-Rel enhanceosome, and IL-1β-CCL2 crosstalk between tumor cells and the rel-dependent monocytes maintained the suppressive tumor microenvironment. Thus, c-Rel specifies the development of a suppressive monocyte population and could be selectively targeted for treating cancer.

Integral Analyses of Competing Endogenous RNA Mechanisms and DNA Methylation Reveal Regulatory Mechanisms in Osteosarcoma

Background: Osteosarcoma (OS) is the most common primary malignant bone tumour in children and adolescents, with rapid growth, frequent metastasis, and a poor prognosis, but its pathogenesis has not been fully elucidated. Exploring the pathogenesis of OS is of great significance for improving diagnoses and finding new therapeutic targets.

Methods: Differentially expressed circRNAs (DECs), miRNAs (DEMs), methylated DNA sites (DMSs), and mRNAs (DEGs) were identified between OS and control cell lines. GSEA of DEGs and functional enrichment analysis of methylated DEGs were carried out to further identify potential biological processes. Online tools were used to predict the miRNA binding sites of DECs and the mRNA binding sites of DEMs, and then construct a circRNA-miRNA-mRNA network. Next, an analysis of the interaction between methylated DEGs was performed with a protein-protein interaction (PPI) network, and hub gene identification and survival analysis were carried out. The expression pattern of circRNA-miRNA-mRNA was validated by real-time PCR.

Results: GSEA and functional enrichment analysis indicated that DEGs and methylated DEGs are involved in important biological processes in cancer. Hsa_circ_0001753/has_miR_760/CD74 network was constructed and validated in cell lines. Low expression levels of CD74 are associated with poor overall survival times and show good diagnostic ability.

Conclusion: Methylated DEGs may be involved in the development of OS, and the hsa_circ_0001753/has_miR_760/CD74 network may serve as a target for the early diagnosis of and targeted therapy for OS.

Single-Molecule RNA Sequencing Reveals IFNγ-Induced Differential Expression of Immune Escape Genes in Merkel Cell Polyomavirus-Positive MCC Cell Lines

Merkel cell carcinoma (MCC) is a rare and highly aggressive cancer, which is mainly caused by genomic integration of the Merkel cell polyomavirus and subsequent expression of a truncated form of its large T antigen. The resulting primary tumor is known to be immunogenic and under constant pressure to escape immune surveillance. Because interferon gamma (IFNγ), a key player of immune response, is secreted by many immune effector cells and has been shown to exert both anti-tumoral and pro-tumoral effects, we studied the transcriptomic response of MCC cells to IFNγ. In particular, immune modulatory effects that may help the tumor evade immune surveillance were of high interest to our investigation. The effect of IFNγ treatment on the transcriptomic program of three MCC cell lines (WaGa, MKL-1, and MKL-2) was analyzed using single-molecule sequencing via the Oxford Nanopore platform. A significant differential expression of several genes was detected across all three cell lines. Subsequent pathway analysis and manual annotation showed a clear upregulation of genes involved in the immune escape of tumor due to IFNγ treatment. The analysis of selected genes on protein level underlined our sequencing results. These findings contribute to a better understanding of immune escape of MCC and may help in clinical treatment of MCC patients. Furthermore, we demonstrate that single-molecule sequencing can be used to assess characteristics of large eukaryotic transcriptomes and thus contribute to a broader access to sequencing data in the community due to its low cost of entry.

HLA Class II Histocompatibility Antigen γ Chain (CD74) Expression Is Associated with Immune Cell Infiltration and Favorable Outcome in Breast Cancer

Simple Summary

CD74 is a transmembrane protein normally expressed in immune cells, and aberrantly expressed in cancer cells. Although CD74 overexpression is mostly associated with hematologic malignancies, some studies have also reported CD74 expression in breast cancer especially associated to the triple negative subtype and metastatic breast cancer. The triple-negative breast cancer is generally more aggressive and with a poorer prognosis than the other subtypes. Immunotherapy holds great promise in clinical management of breast cancer, and CD74 may play a regulatory role in immune system responses. Our results showed that CD74 is associated with expression of programmed cell death ligand 1 (PD-L1), which in turn is involved in preventing anticancer immune responses. Overall, our results indicate that CD74 may be a therapeutic target for the treatment of breast cancer patients, in particular in triple negative breast cancer and metastatic breast cancers, where CD74 is commonly overexpressed.

Abstract

The triple-negative breast cancer (TNBC) subtype, defined as negative for ER, PgR, and HER2, is biologically more aggressive and with a poorer prognosis than the other subtypes, in part due to the lack of suitable targeted therapies. Consequently, identification of any potential novel therapeutic option, predictive and/or prognostic biomarker, or any other relevant information that may impact the clinical management of this group of patients is valuable. The HLA class II histocompatibility antigen γ chain, or cluster of differentiation 74 (CD74), has been associated with TNBCs, and poorer survival. However, discordant results have been reported for immunohistochemical studies of CD74 expression in breast cancer. Here we report validation studies for use of a novel CD74 antibody, UMAb231. We used this antibody to stain a TMA including 640 human breast cancer samples, and found no association with the TNBC subtype, but did find a positive correlation with outcome. We also found associations between CD74 expression and immune cell infiltration, and expression of programmed death ligand 1 (PD-L1). Given that CD74 may play a role in innate immune system responses and the potential of immunotherapy as a viable treatment strategy for TNBCs, CD74 expression may have predictive value for immune checkpoint therapies.

Therapeutic Effects of Dietary Soybean Genistein on Triple-Negative Breast Cancer via Regulation of Epigenetic Mechanisms

Consumption of dietary natural components such as genistein (GE) found in soy-rich sources is strongly associated with a lower risk of breast cancer. However, bioactive dietary component-based therapeutic strategies are largely understudied in breast cancer treatment. Our investigation sought to elucidate the potential mechanisms linking bioactive dietary GE to its breast cancer chemotherapeutic potential in a special subtype of aggressive breast cancer-triple-negative breast cancer (TNBC)-by utilizing two preclinical patient-derived xenograft (PDX) orthotopic mouse models: BCM-3204 and TM00091. Our study revealed that administration of GE resulted in a delay of tumor growth in both PDX models. With transcriptomics analyses in TNBC tumors isolated from BCM-3204 PDXs, we found that dietary soybean GE significantly influenced multiple tumor-regulated gene expressions. Further validation assessment of six candidate differentially expressed genes (DEGs)-Cd74, Lpl, Ifi44, Fzd9, Sat1 and Wwc1-demonstrated a similar trend at gene transcriptional and protein levels as observed in RNA-sequencing results. Mechanistically, GE treatment-induced Cd74 downregulation regulated the NF-κB/Bcl-xL/TAp63 signal pathway, which may contribute to soybean GE-mediated therapeutic effects on TNBC tumors. Additionally, our findings revealed that GE can modify expression levels of key epigenetic-associated genes such as DNA methyltransferases (Dnmt3b), ten-eleven translocation (Tet3) methylcytosine dioxygenases and histone deacetyltransferase (Hdac2), and their enzymatic activities as well as genomic DNA methylation and histone methylation (H3K9) levels. Collectively, our investigation shows high significance for potential development of a novel therapeutic approach by using bioactive soybean GE for TNBC patients who have few treatment options.

Epigenetic Regulation of Immunotherapy Response in Triple-Negative Breast Cancer

Triple-negative breast cancer (TNBC) is defined by the absence of estrogen receptor and progesterone receptor and human epidermal growth factor receptor 2 (HER2) overexpression. This malignancy, representing 15-20% of breast cancers, is a clinical challenge due to the lack of targeted treatments, higher intrinsic aggressiveness, and worse outcomes than other breast cancer subtypes. Immune checkpoint inhibitors have shown promising efficacy for early-stage and advanced TNBC, but this seems limited to a subgroup of patients. Understanding the underlying mechanisms that determine immunotherapy efficiency is essential to identifying which TNBC patients will respond to immunotherapy-based treatments and help to develop new therapeutic strategies. Emerging evidence supports that epigenetic alterations, including aberrant chromatin architecture conformation and the modulation of gene regulatory elements, are critical mechanisms for immune escape. These alterations are particularly interesting since they can be reverted through the inhibition of epigenetic regulators. For that reason, several recent studies suggest that the combination of epigenetic drugs and immunotherapeutic agents can boost anticancer immune responses. In this review, we focused on the contribution of epigenetics to the crosstalk between immune and cancer cells, its relevance on immunotherapy response in TNBC, and the potential benefits of combined treatments.

CD74 and CD44 Expression on CTCs in Cancer Patients with Brain Metastasis

Up to 40% of advance lung, melanoma and breast cancer patients suffer from brain metastases (BM) with increasing incidence. Here, we assessed whether circulating tumor cells (CTCs) in peripheral blood can serve as a disease surrogate, focusing on CD44 and CD74 expression as prognostic markers for BM. We show that a size-based microfluidic approach in combination with a semi-automated cell recognition system are well suited for CTC detection in BM patients and allow further characterization of tumor cells potentially derived from BM. CTCs were found in 50% (7/14) of breast cancer, 50% (9/18) of non-small cell lung cancer (NSCLC) and 36% (4/11) of melanoma patients. The next-generation sequencing (NGS) analysis of nine single CTCs from one breast cancer patient revealed three different CNV profile groups as well as a resistance causing ERS1 mutation. CD44 and CD74 were expressed on most CTCs and their expression was strongly correlated, whereas matched breast cancer BM tissues were much less frequently expressing CD44 and CD74 (negative in 46% and 54%, respectively). Thus, plasticity of CD44 and CD74 expression during trafficking of CTCs in the circulation might be the result of adaptation strategies.

CD74+ macrophages are associated with favorable prognosis and immune contexture in hepatocellular carcinoma

CD74 was initially thought to participate mainly in antigen presentation as an MHC class II chaperone. Recent studies have shown that CD74 plays an important role within the cell and throughout the immune system in a wide spectrum of neoplasms. However, the role of CD74 in hepatocellular carcinoma (HCC) remains elusive. In this study, HCC tissues from Zhongshan Hospital and data from The Cancer Genome Atlas (TCGA) were obtained and analyzed. Immunohistochemistry, flow cytometry, and single-cell RNA sequencing (scRNA-seq) were performed to detect the characteristics of CD74⁺ cells and explore their impact on the tumor microenvironment (TME) of HCC. Our data revealed that stromal CD74⁺ cell enrichment was associated with favorable prognosis in patients with HCC. CD74 was abundant in a large portion of HCC specimens and prominently distributed on stromal macrophages. scRNA-seq data also indicated that the pathways related to immune response were significantly upregulated in CD74⁺ macrophages. High infiltration of CD74⁺ macrophages was associated with increased infiltration of CD8⁺ cytotoxic T lymphocytes (CTLs) with enhanced effector functions in HCC. Besides, blocking CD74 weakened the antitumor activity and proliferation ability of CD8⁺ CTLs in HCC. Our findings highlight the critical role of CD74 in HCC. New drugs and antibodies targeting CD74 may be effective strategies for HCC therapy.

Identification of an immune prognostic 11-gene signature for lung adenocarcinoma

Background

The immunological tumour microenvironment (TME) has occupied a very important position in the beginning and progression of non-small cell lung cancer (NSCLC). Prognosis of lung adenocarcinoma (LUAD) remains poor for the local progression and widely metastases at the time of clinical diagnosis. Our objective is to identify a potential signature model to improve prognosis of LUAD.

Methods

With the aim to identify a novel immune prognostic signature associated with overall survival (OS), we analysed LUADs extracted from The Cancer Genome Atlas (TCGA). Immune scores and stromal scores of TCGA-LUAD were downloaded from Estimation of STromal and Immune cells in MAlignant Tumour tissues Expression using data (ESTIMATE). LASSO COX regression was applied to build the prediction model. Then, the prognostic gene signature was validated in the GSE68465 dataset.

Results

The data from TCGA datasets showed patients in stage I and stage II had higher stromal scores than patients in stage IV (P < 0.05), and for immune score patients in stage I were higher than patients in stage III and stage IV (P < 0.05). The improved overall survivals were observed in high stromal score and immune score groups. Patients in the high-risk group exhibited the inferior OS (P = 2.501e − 05). By validating the 397 LUAD patients from GSE68465, we observed a better OS in the low-risk group compared to the high-risk group, which is consistent with the results from the TCGA cohort. Nomogram results showed that practical and predicted survival coincided very well, especially for 3-year survival.

Conclusion

We obtained an 11 immune score related gene signature model as an independent element to effectively classify LUADs into different risk groups, which might provide a support for precision treatments. Moreover, immune score may play a potential valuable sole for estimating OS in LUADs.

Linking epigenetic signature and metabolic phenotype in IDH mutant and IDH wildtype diffuse glioma

AIMS

Changes in metabolism are known to contribute to tumour phenotypes. If and how metabolic alterations in brain tumours contribute to patient outcome is still poorly understood. Epigenetics impact metabolism and mitochondrial function. The aim of this study is a characterisation of metabolic features in molecular subgroups of isocitrate dehydrogenase mutant (IDHmut) and isocitrate dehydrogenase wildtype (IDHwt) gliomas.

METHODS

We employed DNA methylation pattern analyses with a special focus on metabolic genes, large-scale metabolism panel immunohistochemistry (IHC), qPCR-based determination of mitochondrial DNA copy number and immune cell content using IHC and deconvolution of DNA methylation data. We analysed molecularly characterised gliomas (n = 57) for in depth DNA methylation, a cohort of primary and recurrent gliomas (n = 22) for mitochondrial copy number and validated these results in a large glioma cohort (n = 293). Finally, we investigated the potential of metabolic markers in Bevacizumab (Bev)-treated gliomas (n = 29).

RESULTS

DNA methylation patterns of metabolic genes successfully distinguished the molecular subtypes of IDHmut and IDHwt gliomas. Promoter methylation of lactate dehydrogenase A negatively correlated with protein expression and was associated with IDHmut gliomas. Mitochondrial DNA copy number was increased in IDHmut tumours and did not change in recurrent tumours. Hierarchical clustering based on metabolism panel IHC revealed distinct subclasses of IDHmut and IDHwt gliomas with an impact on patient outcome. Further quantification of these markers allowed for the prediction of survival under anti-angiogenic therapy.

CONCLUSION

A mitochondrial signature was associated with increased survival in all analyses, which could indicate tumour subgroups with specific metabolic vulnerabilities.

Identification of novel prognosis-related genes in the endometrial cancer immune microenvironment

The incidence of endometrial cancer is increasing each year, and treatment effects are poor for patients with advanced and specific subtypes. Exploring immune infiltration-related factors in endometrial cancer can aid in the prognosis of patients and provide new immunotherapy targets. We downloaded immune metagene and functional data of patients with different subtypes of endometrial cancer from The Cancer Genome Atlas database and selected the lymphocyte-specific kinase (LCK) metagene as a representative genetic marker of the immune microenvironment in endometrial cancer. The results showed that LCK metagene expression is related to the prognosis of patients with endometrioid endometrial adenocarcinoma subtypes and highly correlated with the PTEN and PIK3CA mutational status. A search for LCK-related modules returned seven independent genetic predictors of survival in patients with endometrial cancer. The TIMER algorithm showed that the expression of these seven genes was positively correlated with the infiltration levels of six types of immune cells. The diagnostic value of these markers was validated using real-time quantitative PCR and immunohistochemical methods. Our results identified CD74, HLA-DRB5, CD52, HLA-DPB1 and HLA-DRB1 as possible valuable genetic markers for the diagnosis and prognosis of endometrial cancer and provided a theoretical basis for immunotherapy targets for its clinical treatment.

ZIAQ: a quantile regression method for differential expression analysis of single-cell RNA-seq data

MOTIVATION

Single-cell RNA sequencing (scRNA-seq) has enabled the simultaneous transcriptomic profiling of individual cells under different biological conditions. scRNA-seq data have two unique challenges that can affect the sensitivity and specificity of single-cell differential expression analysis: a large proportion of expressed genes with zero or low read counts ('dropout' events) and multimodal data distributions.

RESULTS

We have developed a zero-inflation-adjusted quantile (ZIAQ) algorithm, which is the first method to account for both dropout rates and complex scRNA-seq data distributions in the same model. ZIAQ demonstrates superior performance over several existing methods on simulated scRNA-seq datasets by finding more differentially expressed genes. When ZIAQ was applied to the comparison of neoplastic and non-neoplastic cells from a human glioblastoma dataset, the ranking of biologically relevant genes and pathways showed clear improvement over existing methods.

AVAILABILITY AND IMPLEMENTATION

ZIAQ is implemented in the R language and available at https://github.com/gefeizhang/ZIAQ.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Key genes with prognostic values in suppression of osteosarcoma metastasis using comprehensive analysis

BACKGROUND

Osteosarcoma is a primary malignant tumor originating from mesenchymal tissue, with a poor distant metastasis prognosis. The molecular mechanisms of osteosarcoma metastasis are extremely complicated.

METHODS

A public data series (GSE21257) was used to identify differentially expressed genes (DEGs) in osteosarcoma patients that did, or did not, develop metastases. Functional enrichment analysis, a protein-protein interaction network, and survival analysis of DEGs were performed. DEGs with a prognostic value were considered as candidate genes and their functional predictions, different expression in normal and malignant tissues, and immune infiltration were analyzed.

RESULTS

The DEGs were mainly enriched in the immune response. Three candidate genes (ALOX5AP, CD74, and FCGR2A) were found, all of which were expressed at higher levels in lungs and lymph nodes than in matched cancer tissues and were probably expressed in the microenvironment.

CONCLUSIONS

Candidate genes can help us understand the molecular mechanisms underlying osteosarcoma metastasis and provide targets for future research.

Surfactant Expression Defines an Inflamed Subtype of Lung Adenocarcinoma Brain Metastases that Correlates with Prolonged Survival

PURPOSE

To provide a better understanding of the interplay between the immune system and brain metastases to advance therapeutic options for this life-threatening disease.

EXPERIMENTAL DESIGN

Tumor-infiltrating lymphocytes (TIL) were quantified by semiautomated whole-slide analysis in brain metastases from 81 lung adenocarcinomas. Multi-color staining enabled phenotyping of TILs (CD3, CD8, and FOXP3) on a single-cell resolution. Molecular determinants of the extent of TILs in brain metastases were analyzed by transcriptomics in a subset of 63 patients. Findings in lung adenocarcinoma brain metastases were related to published multi-omic primary lung adenocarcinoma The Cancer Genome Atlas data (n = 230) and single-cell RNA-sequencing (scRNA-seq) data (n = 52,698).

RESULTS

TIL numbers within tumor islands was an independent prognostic marker in patients with lung adenocarcinoma brain metastases. Comparative transcriptomics revealed that expression of three surfactant metabolism-related genes (SFTPA1, SFTPB, and NAPSA) was closely associated with TIL numbers. Their expression was not only prognostic in brain metastasis but also in primary lung adenocarcinoma. Correlation with scRNA-seq data revealed that brain metastases with high expression of surfactant genes might originate from tumor cells resembling alveolar type 2 cells. Methylome-based estimation of immune cell fractions in primary lung adenocarcinoma confirmed a positive association between lymphocyte infiltration and surfactant expression. Tumors with a high surfactant expression displayed a transcriptomic profile of an inflammatory microenvironment.

CONCLUSIONS

The expression of surfactant metabolism-related genes (SFTPA1, SFTPB, and NAPSA) defines an inflamed subtype of lung adenocarcinoma brain metastases characterized by high abundance of TILs in close vicinity to tumor cells, a prolonged survival, and a tumor microenvironment which might be more accessible to immunotherapeutic approaches.

Grade II/III Glioma Microenvironment Mining and Its Prognostic Merit

OBJECTIVE

The tumor microenvironment greatly influences tumor formation, invasion, and progression. The ESTIMATE (Estimation of STromal and Immune cells in MAlignant Tumor tissues) algorithm quantifies stromal and immune components in a tumor, reflecting the tumor microenvironment. This study aimed to explore key prognostic genes in a grade II/III glioma microenvironment.

METHODS

We obtained stromal/immune scores for the Cancer Genome Atlas (TCGA) grade II/III glioma cohort from the online ESTIMATE portal. The associations of stromal/immune scores with clinicopathologic characteristics and overall survival of patients with grade II/III glioma were assessed by the Mann-Whitney U test and the Kaplan-Meier method, respectively. Functional enrichment analysis and protein-protein interaction network assessments were employed to analyze differentially expressed genes (DEGs). The top 7 genes with 5 or more edges in the protein-protein interaction network were selected. For validation, CGGA grade II/III glioma data were analyzed.

RESULTS

The results showed that elevated stromal/immune/ESTIMATE score was significantly associated with poor survival of patients with TCGA grade II/III glioma. Functional enrichment analysis showed that DEGs were associated with immune cell regulation, extracellular matrix, cytokine activation, and receptor binding. The selected DEGs (interleukin-10, beta-2 microglobulin, C-C motif chemokine ligand 5, cluster of differentiation 74, human leukocyte antigen-DRA, lymphocyte cytosolic protein 2, and myxovirus resistance protein 1) showed prognostic values in patients with grade II/III glioma of the TCGA and CGGA database.

CONCLUSIONS

Stromal/immune/ESTIMATE scores have prognostic values in patients with grade II/III glioma. The selected DEGs, including interleukin-10, beta-2 microglobulin, C-C motif chemokine ligand 5, cluster of differentiation 74, human leukocyte antigen-DRA, lymphocyte cytosolic protein 2, and myxovirus resistance protein 1, associated with tumor immunity and microenvironment, have prognostic values in grade II/III glioma. Further investigation of these genes could provide novel insights into the tumor microenvironment of glioma.

Expressions of HLA Class II Genes in Cutaneous Melanoma Were Associated with Clinical Outcome: Bioinformatics Approaches and Systematic Analysis of Public Microarray and RNA-Seq Datasets

Major histocompatibility complex (MHC) class II molecules, encoded by human leukocyte antigen (HLA) class II genes, play important roles in antigen presentation and initiation of immune responses. However, the correlation between HLA class II gene expression level and patient survival and disease progression in cutaneous melanoma is still under investigation. In the present study, we analyzed microarray and RNA-Seq data of cutaneous melanoma from The Cancer Genome Atlas (TCGA) using different bioinformatics tools. Survival analysis revealed higher expression level of HLA class II genes in cutaneous melanoma, especially HLA-DP and -DR, was significantly associated with better overall survival. Furthermore, the expressions of HLA class II genes were most closely associated with survival in cutaneous melanoma as compared with other cancer types. The expression of HLA class II co-expressed genes, which were found to associate with antigen processing, immune response, and inflammatory response, was also positively associated with overall survival in cutaneous melanoma. Therefore, the results indicated that increased HLA class II expression may contribute to enhanced anti-tumor immunity and related inflammatory response via presenting tumor antigens to the immune system. The expression pattern of HLA class II genes may serve as a prognostic biomarker and therapeutic targets in cutaneous melanoma.