Hypoxic niches attract and sequester tumor-associated macrophages and cytotoxic T cells and reprogram them for immunosuppression.

Glioblastoma (GBM), a highly lethal brain cancer, is notorious for immunosuppression, but the mechanisms remain unclear. Here, we documented a temporospatial patterning of tumor-associated myeloid cells (TAMs) corresponding to vascular changes during GBM progression. As tumor vessels transitioned from the initial dense regular network to later scant and engorged vasculature, TAMs shifted away from perivascular regions and trafficked to vascular-poor areas. This process was heavily influenced by the immunocompetence state of the host. Utilizing a sensitive fluorescent UnaG reporter to track tumor hypoxia, coupled with single-cell transcriptomics, we revealed that hypoxic niches attracted and sequestered TAMs and cytotoxic T lymphocytes (CTLs), where they were reprogrammed toward an immunosuppressive state. Mechanistically, we identified chemokine CCL8 and cytokine IL-1ß as two hypoxic-niche factors critical for TAM trafficking and co-evolution of hypoxic zones into pseudopalisading patterns. Therefore, perturbation of TAM patterning in hypoxic zones may improve tumor control.

Germline genetic contribution to the immune landscape of cancer.

Understanding the contribution of the host's genetic background to cancer immunity may lead to improved stratification for immunotherapy and to the identification of novel therapeutic targets. We investigated the effect of common and rare germline variants on 139 well-defined immune traits in ∼9000 cancer patients enrolled in TCGA. High heritability was observed for estimates of NK cell and T cell subset infiltration and for interferon signaling. Common variants of IFIH1, TMEM173 (STING1), and TMEM108 were associated with differential interferon signaling and variants mapping to RBL1 correlated with T cell subset abundance. Pathogenic or likely pathogenic variants in BRCA1 and in genes involved in telomere stabilization and Wnt-ß-catenin also acted as immune modulators. Our findings provide evidence for the impact of germline genetics on the composition and functional orientation of the tumor immune microenvironment. The curated datasets, variants, and genes identified provide a resource toward further understanding of tumor-immune interactions.

Subtyping and prognostic model construction based on vesicle-mediated transport-related genes in colorectal cancer.

BACKGROUND: Colorectal cancer (CRC) is impacted by various environmental and genetic variables. Dysregulation of vesicle-mediated transport-related genes (VMTRGs) has been observed in many malignancies, but their effect on prognosis in CRC remains unclear. METHODS: CRC samples were clustered into varying subtypes per differential expression of VMTRGs. R package was utilized to explore differences in survival, immune, and drug sensitivity among different disease subtypes. According to differentially expressed genes (DEGs) between subtypes, regression analysis was employed to build a riskscore model and identify independent prognostic factors. The model was validated through a Gene Expression Omnibus (GEO) dataset. Immune landscape, immunophenoscore (IPS), and Tumor Immune Dysfunction and Exclusion (TIDE) scores for different risk groups were calculated. RESULTS: Two subtypes of CRC were identified based on VMTRGs, which showed significant differences in survival rates, immune cell infiltration abundance, immune functional activation levels, and immune checkpoint expression levels. Cluster2 exhibited higher sensitivity to anti-tumor drugs such as Nilotinib, Cisplatin, and Oxaliplatin compared to Cluster1. DEGs were mainly enriched in biological processes such as epidermis development, epidermal cell differentiation, and receptor-ligand activity, and signaling pathways like pancreatic secretion. The constructed 13-gene riskscore model demonstrated good predictive ability for CRC patients' prognosis. Furthermore, differences in immune landscape, IPS, and TIDE scores were observed among different risk groups. CONCLUSION: This study successfully obtained two CRC subtypes with distinct survival statuses and immune levels based on differential expression of VMTRGs. A 13-gene risk model was constructed. The findings had important implications for prognosis and treatment of CRC.

Prognostic evaluation of the novel blueprint of DNA methylation sites by integrating bulk RNA-sequencing and methylation modification data in endometrial cancer.

INTRODUCTION: Endometrial cancer (EC) is a prevalent malignancy affecting the female population, with an increasing incidence among younger age groups. DNA methylation, a common epigenetic modification, is well-established to play a key role in cancer progression. We suspected whether DNA methylation could be used as biomarkers for EC prognosis. METHODS: In the present study, we analyzed bulk RNA-sequencing data from 544 EC patients and DNA methylation data from 430 EC patients in the TCGA-UCEC cohort. We applied weighted correlation network analysis to select a key gene set associated with panoptosis. We conducted correlation analysis between transcriptomic data of the selected key genes and DNA methylation data to identify valuable DNA methylation sites. These sites were further screened by Cox regression and least absolute shrinkage and selection operator analysis. Immune microenvironment differences between high-risk and low-risk groups were assessed using single-sample gene set enrichment analysi, xCell and MCPcounter algorithms. RESULTS: Our results identified five DNA methylation sites (cg03906681, cg04549977, cg06029846, cg10043253 and cg15658376) with significant prognostic value in EC. We constructed a prognostic model using these sites, demonstrating satisfactory predictive performance. The low-risk group showed higher immune cell infiltration. Notably, methylation of site cg03906681 was negatively related to CD8 T cell infiltration, whereas cg04549977 exhibited positive correlations with immune infiltration, particularly in macrophages, activated B cells, dendritic cells and myeloid-derived suppressor cells. PD0325901_1060 was strongly correlated with risk scores, indicating a potential therapeutic response for high-risk EC patients. CONCLUSION: We have developed a robust DNA methylation-based prognostic model for EC, which holds promise for improving prognosis prediction and personalized treatment approaches. These findings may contribute to better management of EC patients, particularly in identifying those at higher risk who may benefit from tailored interventions.

Investigating the prognostic role of lncRNAs associated with disulfidptosis-related genes in clear cell renal cell carcinoma.

INTRODUCTION: Renal cell carcinoma (RCC) is a grave malignancy that poses a significant global health burden with over 400,000 new cases annually. Disulfidptosis, a newly discovered programmed cell death process, is linked to the actin cytoskeleton, which plays a vital role in maintaining cell shape and survival. The role of disulfidptosis is poorly depicted in the clear cell histologic variant of RCC (ccRCC). METHODS: Three sets of ccRCC cohorts, ICGC_RECA-EU (n = 91), GSE76207 (n = 32) and TCGA-KIRC (n = 607), were included in our study, the batch effect of which was removed using the "combat" function. Correlation was calculated using the "rcorr" function of the "Hmisc" package for Pearson analysis, which was visualized using the "pheatmap" package. Principal component analysis was performed by the "vegan" package, visualized using the "scatterplot3d" package. Long non-coding RNAs (lncRNAs) associated with disulfidptosis were screened out using least absolute shrinkage and selection operator (LASSO) and COX analysis. Tumor mutation, immune landscaping and immunotherapy prediction were performed for further characterization of two risk groups. RESULTS: A total of 1822 disulfidptosis-related lncRNAs was selected, among which 308 lncRNAs were found to be significantly associated with the clinical outcome of ccRCC patients. We retained 11 disulfidptosis-related lncRNAs, namely, AP000439.3, RP11-417E7.1, RP11-119D9.1, LINC01510, SNHG3, AC156455.1, RP11-291B21.2, EMX2OS, AC093850.2, HAGLR and RP11-389C8.2, through LASSO and COX analysis for prognosis model construction, which displayed satisfactory accuracy (area under the curve, AUC, values all above 0.6 in multiple cohorts) in stratification of ccRCC prognosis. A nomogram model was constructed by integrating clinical factors with risk score, which further enhanced the prediction efficacy (AUC values all above 0.7 in multiple cohorts). We found that patients of male gender, higher clinical stages and advanced pathological T stage were inclined to have higher risk score values. Dactinomycin_1911, Vinblastine_1004, Daporinad_1248 and Vinorelbine_2048 were identified as promising candidate drugs for treating ccRCC patients of higher risk score value. Moreover, patients of higher risk value were prone to be resistant to immunotherapy. CONCLUSION: We developed a prognosis predicting model based on 11 selected disulfidptosis-related lncRNAs, the efficacy of which was verified in different cohorts. Furthermore, we delineated an intricate portrait of tumor mutation, immune topography and pharmacosensitivity evaluations within disparate risk stratifications.

Construction of an immune predictive model and identification of TRIP6 as a prognostic marker and therapeutic target of CRC by integration of single-cell and bulk RNA-seq data.

BACKGROUND: Investigations elucidating the complex immunological mechanisms involved in colorectal cancer (CRC) and accurately predicting patient outcomes via bulk RNA-Seq analysis have been notably limited. This study aimed to identify the immune status of CRC patients, construct a prognostic model, and identify prognostic signatures via bulk RNA sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq). METHODS: The scRNA-seq data of CRC were downloaded from Gene Expression Omnibus (GEO). The UCSC Xena database was used to obtain bulk RNA-seq data. Differentially expressed gene (DEG), functional enrichment, and random forest analyses were conducted in order to identify core genes associated with colorectal cancer (CRC) that were relevant to prognosis. A molecular immune prediction model was developed using logistic regression after screening features using the least absolute shrinkage and selection operator (LASSO). The differences in immune cell infiltration, mutation, chemotherapeutic drug sensitivity, cellular senescence, and communication between patients who were at high and low risk of CRC according to the predictive model were investigated. The prognostic genes that were closely associated with CRC were identified by random survival forest (RSF) analysis. The expression levels and clinical significance of the hub genes were analyzed in vitro. The LoVo cell line was employed to ascertain the biological role of thyroid hormone receptor-interacting protein 6 (TRIP6). RESULTS: A total of seven main cell subtypes were identified by scRNA-seq analysis. A molecular immune predictive model was constructed based on the risk scores. The risk score was significantly associated with OS, stage, mutation burden, immune cell infiltration, response to immunotherapy, key pathways, and cell-cell communication. The functions of the six hub genes were determined and further utilized to establish a regulatory network. Our findings unequivocally confirmed that TRIP6 upregulation was verified in the CRC samples. After knocking down TRIP6, cell proliferation, migration, and invasion of LoVo cells were inhibited, and apoptosis was promoted. CONCLUSIONS: The molecular predictive model reliably distinguished the immune status of CRC patients. We further revealed that TRIP6 may act as an oncogene in CRC, making it a promising candidate for targeted therapy and as a prognostic marker for CRC.

Anoikis-related signature predicts prognosis and characterizes immune landscape of ovarian cancer.

Ovarian cancer (OV) is the most lethal gynecological malignancy worldwide, with high recurrence rates. Anoikis, a newly-acknowledged form of programmed cell death, plays an essential role in cancer progression, though studies focused on prognostic patterns of anoikis in OV are still lacking. We filtered 32 potential anoikis-related genes (ARGs) among the 6406 differentially expressed genes (DEGs) between the 180 normal controls and 376 TCGA-OV samples. Through the LASSO-Cox analysis, a 2-gene prognostic signature, namely AKT2, and DAPK1, was finally distinguished. We then demonstrated the promising prognostic value of the signature through the K-M survival analysis and time-dependent ROC curves (p-value < 0.05). Moreover, based on the signature and clinical features, we constructed and validated a nomogram model for 1-year, 3-year, and 5-year overall survival, with reliable prognostic values in both TCGA-OV training cohort (p-value < 0.001) and ICGC-OV validation cohort (p-value = 0.030). We evaluated the tumor immune landscape through the CIBERSORT algorithm, which indicated the upregulation of resting Myeloid Dendritic Cells (DCs), memory B cells, and naïve B cells and high expression of key immune checkpoint molecules (CD274 and PDCD1LG2) in the high-risk group. Interestingly, the high-risk group exhibited better sensitivity toward immunotherapy and less sensitivity toward chemotherapies, including Cisplatin and Bleomycin. Especially, based on the IHC of tissue microarrays among 125 OV patients at our institution, we reported that aberrant upregulation of DAPK1 was related to poor prognosis. Conclusively, the anoikis-related signature was a promising tool to evaluate prognosis and predict therapy responses, thus assisting decision-making in the realm of OV precision medicine.

Screening of novel tumor-associated antigens for lung adenocarcinoma mRNA vaccine development based on pyroptosis phenotype genes.

This study aimed to identify new pyroptosis-associated tumor antigens for use in mRNA vaccines and the screening of sensitive LUAD populations suitable for vaccination. The association between tumor immune infiltrating cell abundance and potential tumor antigens was investigated and visualized using the analysis modules of gene expression, clinical outcomes, and somatic copy number variation. In addition, the pyroptosis-related genes (PRGs) were clustered, the relative pyroptosis subtypes (PSs) and gene modules were identified, and the prognostic value of the PSs was examined. The expression of key PRGs in two lung adenocarcinoma cell lines was verified by RT-qPCR. Four tumor pyroptosis-associated antigens, CARD8, NAIP, NLRP1, and NLRP3, were screened as potential candidates for LUAD mRNA vaccine development. In the construction of consensus clusters for PRGs, two PSs, PS1 and PS2, were classified, in which patients with PS1 LUAD had a better prognosis. In contrast, patients with PS2 LUAD may have better responsiveness to mRNA vaccine treatment. The key PRGs can be regarded as biomarkers to predict the LUAD prognosis and identify patients suitable for mRNA vaccines. The RT-qPCR results showed that the expression levels of CSMD3, LRP1B, MUC16 and TTN were significantly increased in the two lung adenocarcinoma cell lines, while the expression levels of CARD8, TP53 and ZFHX4 were significantly reduced. The antigens CARD8, NAIP, NLRP1, and NLRP3, which are associated with tumor pyroptosis, could be candidate molecules for LUAD mRNA vaccine development. Patients with PS2 LUAD may be suitable candidates for mRNA vaccine treatment.

Immune cell patterns before and after neoadjuvant immune checkpoint blockade combined with chemoradiotherapy in locally advanced esophageal squamous cell carcinoma.

BACKGROUND: Neoadjuvant immune checkpoint blockade (ICB) combined with chemoradiotherapy offers high pathologic complete response (pCR) rate for patients with locally advanced esophageal squamous cell carcinomas (ESCC). But the dynamic tumor immune microenvironment modulated by such neoadjuvant therapy remains unclear. PATIENTS AND METHODS: A total of 41 patients with locally advanced ESCC were recruited. All patients received neoadjuvant toripalimab combined with concurrent chemoradiotherapy. Matched pre- and post-treatment tissues were obtained for fluorescent multiplex immunohistochemistry (mIHC) and IHC analyses. The densities and spatial distributions of immune cells were determined by HALO modules. The differences of immune cell patterns before and after neoadjuvant treatment were investigated. RESULTS: In the pre-treatment tissues, more stromal CD3 + FoxP3 + Tregs and CD86+/CD163 + macrophages were observed in patients with residual tumor existed in the resected lymph nodes (pN1), compared with patients with pCR. The majority of macrophages were distributed in close proximity to tumor nest in pN1 patients. In the post-treatment tissues, pCR patients had less CD86 + cell infiltration, whereas higher CD86 + cell density was significantly associated with higher tumor regression grades (TRG) in non-pCR patients. When comparing the paired pre- and post-treatment samples, heterogeneous therapy-associated immune cell patterns were found. Upon to the treatment, CD3 + T lymphocytes were slightly increased in pCR patients, but markedly decreased in non-pCR patients. In contrast, a noticeable increase and a less obvious decrease of CD86 + cell infiltration were respectively depicted in non-pCR and pCR patients. Furthermore, opposite trends of the treatment-induced alterations of CD8 + and CD15 + cell infiltrations were observed between pN0 and pN1 patients. CONCLUSIONS: Collectively, our data demonstrate a comprehensive picture of tumor immune landscape before and after neoadjuvant ICB combined with chemoradiotherapy in ESCC. The infiltration of CD86 + macrophage may serve as an unfavorable indicator for neoadjuvant toripalimab combined with chemoradiotherapy.

Identification of Lung Adenocarcinoma Subtypes Based on MHC-II Gene Expression Profile and Immunological Analysis.

INTRODUCTION: Major histocompatibility complex class II molecule (MHC-II) is pivotal in anti-tumor immunity, and targeting MHC-II in tumors may help improve patient survival. But function of MHC-II in the immunotherapy and prognosis of lung adenocarcinoma (LUAD) patients has not been thoroughly studied and reported. METHODS: We selected LUAD-related MHC-II genes from public databases based on previous literature reports. We identified different subtypes according to expression differences of these genes in different LUAD samples through cluster analysis. We used R package to conduct a series of analyses on different subtypes, exploring their survival differences, gene expression differences, pathway enrichment differences, and differences in immune characteristics and immune therapy. Finally, we screened potential drugs from the cMAP database. RESULTS: We identified two MHC-II-related LUAD subtypes. Our analyses presented that patients with cluster2 subtype showed better prognosis, higher immune scores, higher levels of immune cell infiltration and immune function activation. In addition, patients with this subtype had higher immunophenoscore, lower TIDE scores, and DEPTH scores. We also identified 10 small molecule drugs, such as lenalidomide, VX-745, and tyrphostin-AG-1295. CONCLUSION: Overall, MHC-II is not only a potential biomarker for accurately distinguishing LUAD subtypes but also a predictive factor for their survival. Our study offers novel insights into understanding of impact of MHC-II in LUAD and offers a new perspective for improving the accurate classification of LUAD patients and enhancing drug treatment.

A novel molecular classification based on efferocytosis-related genes for predicting clinical outcome and treatment response in acute myeloid leukemia.

BACKGROUND: Previous studies have shown that macrophage-mediated efferocytosis is involved in immunosuppression in acute myeloid leukemia (AML). However, the regulatory role of efferocytosis in AML remains unclear and needs further elucidation. METHODS: We first identified the key efferocytosis-related genes (ERGs) based on the expression matrix. Efferocytosis-related molecular subtypes were obtained by consensus clustering algorithm. Differences in immune landscape and biological processes among molecular subtypes were further evaluated. The efferocytosis score model was constructed to quantify molecular subtypes and evaluate its value in prognosis prediction and treatment decision-making in AML. RESULTS: Three distinct efferocytosis-related molecular subtypes were identified and divided into immune activation, immune desert, and immunosuppression subtypes based on the characteristics of the immune landscape. We evaluated the differences in clinical and biological features among different molecular subtypes, and the construction of an efferocytosis score model can effectively quantify the subtypes. A low efferocytosis score is associated with immune activation and reduced mutation frequency, and patients have a better prognosis. A high efferocytosis score reflects immune exhaustion, increased activity of tumor marker pathways, and poor prognosis. The prognostic predictive value of the efferocytosis score model was confirmed in six AML cohorts. Patients exhibiting high efferocytosis scores may derive therapeutic benefits from anti-PD-1 immunotherapy, whereas those with low efferocytosis scores tend to exhibit greater sensitivity towards chemotherapy. Analysis of treatment data in ex vivo AML cells revealed a group of drugs with significant differences in sensitivity between different efferocytosis score groups. Finally, we validated model gene expression in a clinical cohort. CONCLUSIONS: This study reveals that efferocytosis plays a non-negligible role in shaping the diversity and complexity of the AML immune microenvironment. Assessing the individual efferocytosis-related molecular subtype in individuals will help to enhance our understanding of the characterization of the AML immune landscape and guide the establishment of more effective clinical treatment strategies.

Novel kinase 1 regulates CD8+T cells as a potential therapeutic mechanism for idiopathic pulmonary fibrosis.

Idiopathic pulmonary fibrosis (IPF) is a rare, chronic and progressively worsening lung disease that poses a significant threat to patient prognosis, with a mortality rate exceeding that of some common malignancies. Effective methods for early diagnosis and treatment remain for this condition are elusive. In our study, we used the GEO database to access second-generation sequencing data and associated clinical information from IPF patients. By utilizing bioinformatics techniques, we identified crucial disease-related genes and their biological functions, and characterized their expression patterns. Furthermore, we mapped out the immune landscape of IPF, which revealed potential roles for novel kinase 1 and CD8+T cells in disease progression and outcome. These findings can aid the development of new strategies for the clinical diagnosis and treatment of IPF.

Leveraging programmed cell death patterns to predict prognosis and therapeutic sensitivity in OSCC.

OBJECTIVES: Intricate associations between programmed cell death (PCD) and cancer development and treatment outcomes have been increasingly appreciated. Here, we integrated 12 PCD patterns to construct a novel biomarker, cell death index (CDI), for oral squamous cell carcinoma (OSCC) prognostication and therapeutic prediction. MATERIALS AND METHODS: Univariate Cox regression, Kaplan-Meier survival, and LASSO analyses were performed to construct the CDI. A nomogram combining CDI and selected clinicopathological parameters was established by multivariate Cox regression. The associations between CDI and immune landscape and therapeutic sensitivity were estimated. Single-cell RNA-seq data of OSCC was used to infer CDI genes in selected cell types and determine their expression along cell differentiation trajectory. RESULTS: Ten selected PCD genes derived a novel prognostic signature for OSCC. The predictive prognostic performance of CDI and nomogram was robust and superior across multiple independent patient cohorts. CDI was negatively associated with tumor-infiltrating immune cell abundance and immunotherapeutic outcomes. Moreover, scRNA-seq data reanalysis revealed that GSDMB, IL-1A, PRKAA2, and SFRP1 from this signature were primarily expressed in cancer cells and involved in cell differentiation. CONCLUSIONS: Our findings established CDI as a novel powerful predictor for prognosis and therapeutic response for OSCC and suggested its potential involvement in cancer cell differentiation.

Cuproptosis-Associated lncRNA Gene Signature Establishes New Prognostic Profile and Predicts Immunotherapy Response in Endometrial Carcinoma.

Uterine corpus endometrial carcinoma (UCEC), a prevalent kind of cancerous tumor in female reproductive system that has a dismal prognosis in women worldwide. Given the very limited studies of cuproptosis-related lncRNAs (CRLs) in UCEC. Our purpose was to construct a prognostic profile based on CRLs and explore its assess prognostic value in UCEC victims and its correlation with the immunological microenvironment. METHODS: 554 UCEC tumor samples and 23 normal samples' RNA-seq statistics and clinical details were compiled from data in the TCGA database. CRLs were obtained using Pearson correlation analysis. Using LASSO Cox regression, multivariate Cox regression, and univariate Cox regression analysis, six CRLs are confirmed to develop a risk prediction model at last.We identified two main molecular subtypes and observed that multilayer CRLs modifications were related to patient clinicopathological features, prognosis, and tumor microenvironment (TME) cell infiltration characteristics, and then we verified the prognostic hallmark of UCEC and examined its immunological landscape.Finally, using qRT-PCR, model hub genes' expression patterns were confirmed. RESULTS: A unique CRL signature was established by the combination of six differently expressed CRLs that were highly linked with the prognosis of UCEC patients. According to their CRLs signatures, the patients were divided into two groups: the low-risk and the high-risk groups. Compared to individuals at high risk, patients at low risk had higher survival rates (p < 0.001). Additionally, Cox regression reveals that the profiles of lncRNAs linked to cuproptosis may independently predict prognosis in UCEC patients. The 1-, 3-, and 5-year risks' respective receiver operating characteristics (ROC) exhibited AUC values of 0.778, 0.810, and 0.854. Likewise, the signature could predict survival in different groups based on factors like stage, age, and grade, among others. Further investigation revealed differences between the different risk score groups in terms of drug sensitivity,immune cell infiltration,tumor mutation burden (TMB) score and microsatellite instability (MSI) score. Compared to the group of high risk, the low-risk group had greater rates of TMB and MSI. Results from qRT-PCR revealed that in UCEC vs normal tissues, AC026202.2, NRAV, AC079466.2, and AC090617.5 were upregulated,while LINC01545 and AL450384.1 were downregulated. CONCLUSIONS: Our research clarified the relationship between CRLs signature and the immunological profile and prognosis of UCEC.This signature will establish the framework for future investigations into the endometrial cancer CRLs mechanism as well as the exploitation of new diagnostic tools and new therapeutic.

Identification and Analysis of Subtypes of Liver Cancer Based on Genes Related to E3 Ubiquitin Ligases and Deubiquitinating Enzymes.

Recent studies have reported a correlation between ubiquitination or deubiquitination and cancer development. But mechanisms underlying the roles of genes associated with E3 ubiquitin ligases and deubiquitinating enzymes (DUB) in liver cancer remain to be explored. We analyzed and screened differentially expressed genes related to E3 ubiquitin ligases and DUB in liver cancer on the basis of public databases. Cluster analysis was utilized to classify liver cancer samples into different subtypes. Survival analysis, immune analysis, and pathway enrichment analysis were performed on the subtypes. We constructed a protein-protein interaction network using STRING to screen hub genes. Finally, we used the Connectivity Map (CMap) database to predict targeted small molecules. The results show that a total of 139 differentially expressed E3/DUB genes in liver cancer were screened. Then, liver cancer was classified into two subtypes, cluster 1 and cluster 2, based on E3-related and DUB-related genes. Patients in cluster 1 had higher survival rates and immune levels than those in cluster 2. Four hub genes (RPSA, RPS5, RPL30, and RPL8) significantly affecting the survival of the two subtypes of liver cancer patients were identified based on cluster 1 and cluster 2. Finally, the CMap database predicted that small-molecule drugs including probenecid, dexamethasone, and etomidate may improve the prognosis of liver cancer patients. These findings may offer a reference for risk stratification studies and drug development in liver cancer.

Coordinating single-cell and bulk RNA-seq in deciphering the intratumoral immune landscape and prognostic stratification of prostate cancer patients.

INTRODUCTION: Prostate cancer is a common cancer among male population. The aberrant expression of histone modifiers has been identified as a potential driving force in numerous cancer types. However, the mechanism of histone modifiers in the development of prostate cancer remains unknown. METHODS: Expression profiles and clinical data were obtained from GSE70769, GSE46602, and GSE67980. Seruat R package was utilized to calculate the gene set enrichment of the histone modification pathway and obtain the Histone score. Least absolute shrinkage and selection operator (LASSO) and Cox regression analyses were employed to identify marker genes with prognostic value. Kaplan-Meier survival analysis was conducted to assess the efficacy of the prognostic model. In addition, microenvironment cell populations counter (MCPcounter), single-sample gene set enrichment analysis (ssGSEA), and xCell algorithms were employed for immune infiltration analysis. Drug sensitivity prediction was performed using oncoPredict R package. RESULTS: We screened differentially expressed genes (DEGs) between Histone-high score (Histone-H) and Histone-low score (Histone-L) groups, which were enriched in RNA splicing and DNA-binding transcription factor binding pathways. We retained four prognostic marker genes, including TACC3, YWHAH, TAF1C and TTLL5. The risk model showed significant efficacy in stratification of the prognosis of prostate cancer patients in both internal and external cohorts (p < .0001 and p = .032, respectively). In addition, prognostic gene YWHAH was infiltrated in abundance of fibroblasts and highly correlated with Entinostat_1593 drug sensitivity score and the value of risk score. CONCLUSION: We innovatively developed a histone modification-related prognostic model with high prognostic potency and identified YWHAH as possible diagnostic and therapeutic biomarkers for prostate cancer. It provides novel insights to address prostate cancer and enhance clinical outcomes, thereby opening up a new avenue for customized treatment alternatives.

Manufacture and Characterization of Good Manufacturing Practice-Compliant SARS-COV-2 Cytotoxic T Lymphocytes.

BACKGROUND: Coronavirus disease 2019 (COVID-19) is caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). SARS-CoV-2 virus-specific cytotoxic T-cell lymphocytes (vCTLs) could provide a promising modality in COVID-19 treatment. We aimed to screen, manufacture, and characterize SARS-CoV-2-vCTLs generated from convalescent COVID-19 donors using the CliniMACS Cytokine Capture System (CCS). METHODS: Donor screening was done by stimulation of convalescent COVID-19 donor peripheral blood mononuclear cells with viral peptides and identification of interferonγ (IFN-γ)+ CD4 and CD8 T cells using flow cytometry. Clinical-grade SARS-CoV-2-vCTLs were manufactured using the CliniMACS CCS. The enriched SARS-CoV-2-vCTLs were characterized by T-cell receptor sequencing, mass cytometry, and transcriptome analysis. RESULTS: Of the convalescent donor blood samples, 93% passed the screening criteria for clinical manufacture. Three validation runs resulted in enriched T cells that were 79% (standard error of the mean 21%) IFN-γ+ T cells. SARS-CoV-2-vCTLs displayed a highly diverse T-cell receptor repertoire with enhancement of both memory CD8 and CD4 T cells, especially in CD8 TEM, CD4 TCM, and CD4 TEMRA cell subsets. SARS-CoV-2-vCTLs were polyfunctional with increased gene expression in T-cell function, interleukin, pathogen defense, and tumor necrosis factor superfamily pathways. CONCLUSIONS: Highly functional SARS-CoV-2-vCTLs can be rapidly generated by direct cytokine enrichment (12 hours) from convalescent donors. CLINICAL TRIALS REGISTRATION: NCT04896606.

The prognostic value and immune landscaps of m6A/m5C-related lncRNAs signature in the low grade glioma.

BACKGROUND: N6-methyladenosine (m6A) and 5-methylcytosine (m5C) are the main RNA methylation modifications involved in the oncogenesis of cancer. However, it remains obscure whether m6A/m5C-related long non-coding RNAs (lncRNAs) affect the development and progression of low grade gliomas (LGG). METHODS: We summarized 926 LGG tumor samples with RNA-seq data and clinical information from The Cancer Genome Atlas and Chinese Glioma Genome Atlas. 105 normal brain samples with RNA-seq data from the Genotype Tissue Expression project were collected for control. We obtained a molecular classification cluster from the expression pattern of sreened lncRNAs. The least absolute shrinkage and selection operator Cox regression was employed to construct a m6A/m5C-related lncRNAs prognostic signature of LGG. In vitro experiments were employed to validate the biological functions of lncRNAs in our risk model. RESULTS: The expression pattern of 14 sreened highly correlated lncRNAs could cluster samples into two groups, in which various clinicopathological features and the tumor immune microenvironment were significantly distinct. The survival time of cluster 1 was significantly reduced compared with cluster 2. This prognostic signature is based on 8 m6A/m5C-related lncRNAs (GDNF-AS1, HOXA-AS3, LINC00346, LINC00664, LINC00665, MIR155HG, NEAT1, RHPN1-AS1). Patients in the high-risk group harbored shorter survival times. Immunity microenvironment analysis showed B cells, CD4 + T cells, macrophages, and myeloid-derived DC cells were significantly increased in the high-risk group. Patients in high-risk group had the worse overall survival time regardless of followed TMZ therapy or radiotherapy. All observed results from the TCGA-LGG cohort could be validated in CGGA cohort. Afterwards, LINC00664 was found to promote cell viability, invasion and migration ability of glioma cells in vitro. CONCLUSION: Our study elucidated a prognostic prediction model of LGG by 8 m6A/m5C methylated lncRNAs and a critical lncRNA regulation function involved in LGG progression. High-risk patients have shorter survival times and a pro-tumor immune microenvironment.

Demystifying the impact of prenatal tobacco exposure on the placental immune microenvironment: Avoiding the tragedy of mending the fold after death.

Prenatal tobacco exposure (PTE) correlates significantly with a surge in adverse pregnancy outcomes, yet its pathological mechanisms remain partially unexplored. This study aims to meticulously examine the repercussions of PTE on placental immune landscapes, employing a coordinated research methodology encompassing bioinformatics, machine learning and animal studies. Concurrently, it aims to screen biomarkers and potential compounds that could sensitively indicate and mitigate placental immune disorders. In the course of this research, two gene expression omnibus (GEO) microarrays, namely GSE27272 and GSE7434, were included. Gene set enrichment analysis (GSEA) and immune enrichment investigations on differentially expressed genes (DEGs) indicated that PTE might perturb numerous innate or adaptive immune-related biological processes. A cohort of 52 immune-associated DEGs was acquired by cross-referencing the DEGs with gene sets derived from the ImmPort database. A protein-protein interaction (PPI) network was subsequently established, from which 10 hub genes were extracted using the maximal clique centrality (MCC) algorithm (JUN, NPY, SST, FLT4, FGF13, HBEGF, NR0B2, AREG, NR1I2, SEMA5B). Moreover, we substantiated the elevated affinity of tobacco reproductive toxicants, specifically nicotine and nitrosamine, with hub genes through molecular docking (JUN, FGF13 and NR1I2). This suggested that these genes could potentially serve as crucial loci for tobacco's influence on the placental immune microenvironment. To further elucidate the immune microenvironment landscape, consistent clustering analysis was conducted, yielding three subtypes, where the abundance of follicular helper T cells (p < 0.05) in subtype A, M2 macrophages (p < 0.01), neutrophils (p < 0.05) in subtype B and CD8+ T cells (p < 0.05), resting NK cells (p < 0.05), M2 macrophages (p < 0.05) in subtype C were significantly different from the control group. Additionally, three pivotal modules, designated as red, blue and green, were identified, each bearing a close association with differentially infiltrated immunocytes, as discerned by the weighted gene co-expression network analysis (WGCNA). Functional enrichment analysis was subsequently conducted on these modules. To further probe into the mechanisms by which immune-associated DEGs are implicated in intercellular communication, 20 genes serving as ligands or receptors and connected to differentially infiltrating immunocytes were isolated. Employing a variety of machine learning techniques, including one-way logistic regression, LASSO regression, random forest and artificial neural networks, we screened 11 signature genes from the intersection of immune-associated DEGs and secretory protein-encoding genes derived from the Human Protein Atlas. Notably, CCL18 and IFNA4 emerged as prospective peripheral blood markers capable of identifying PTE-induced immune disorders. These markers demonstrated impressive predictive power, as indicated by the area under the curve (AUC) of 0.713 (0.548-0.857) and 0.780 (0.618-0.914), respectively. Furthermore, we predicted 34 potential compounds, including cyclosporine, oestrogen and so on, which may engage with hub genes and attenuate immune disorders instigated by PTE. The diagnostic performance of these biomarkers, alongside the interventional effect of cyclosporine, was further corroborated in animal studies via ELISA, Western blot and immunofluorescence assays. In summary, this study identifies a disturbance in the placental immune landscape, a secondary effect of PTE, which may underlie multiple pregnancy complications. Importantly, our research contributes to the noninvasive and timely detection of PTE-induced placental immune disorders, while also offering innovative therapeutic strategies for their treatment.

Single-cell RNA sequencing reveals peripheral blood mononuclear immune cell landscape associated with operational tolerance in a kidney transplant recipient.

Operational tolerance (OT) after kidney transplantation is defined as stable graft acceptance without the need for immunosuppression therapy. However, it is not clear which cellular and molecular pathways are driving tolerance in these patients. In this first-of-its-kind pilot study, we assessed the immune landscape associated with OT using single-cell analyses. Peripheral mononuclear cells from a kidney transplant recipient with OT (Tol), 2 healthy individuals (HC), and a kidney transplant recipient with normal kidney function on standard-of-care immunosuppression (SOC) were evaluated. The immune landscape of the Tol was drastically different from that of SOC and emerged closer to the profile of HC. TCL1A+ naive B cells and LSGAL1+ regulatory T cells (Tregs) were in higher proportions in Tol. We were unable to identify the Treg subcluster in SOC. The ligand-receptor analysis in HC and Tol identified interactions between B cells, and Tregs that enhance the proliferation and suppressive function of Tregs. SOC reported the highest proportion of activated B cells with more cells in the G2M phase. Our single-cell RNA sequencing study identified the mediators of tolerance; however, it emphasizes the requirement of similar investigations on a larger cohort to reaffirm the role of immune cells in tolerance.