Tumor-infiltrating immune cell score as an independent prognostic predictor for endometrial carcinoma: Insights from a comprehensive analysis of the immune landscape

BACKGROUND

Immune cells are crucial components in the tumor microenvironment and have a significant impact on the outcomes of patients.

AIMS

Here, we aimed to establish a prognostic score based on different types of tumor-infiltrating immune cells for Endometrial Carcinoma (EC).

METHODS AND RESULTS

We enrolled and analyzed 516 EC patients from The Cancer Genome Atlas. The relative abundance of 22 immune cells were estimated by using the CIBERSORTx algorithm. Cox regression was performed to identify potential prognostic immune cells, which were used to develop a Tumor-infiltrating Immune Cell Score (TICS). The prognostic and incremental value of TICS for overall survival were compared with traditional prognostic factors using the C-index and decision curves. Clustering analysis using all immune cells identified three immune landscape subtypes, which had weak correlation with survival. A TICS was constructed using CD8T cells, resting memory CD4 T cells, activated NK and activated DCs, and classified patients as low-, moderate- and high-risk subgroups. The low-risk subgroup had higher tumor mutation burden and activation of IL2/STAT5, IL2/STAT3 and IFN-gamma response pathways. Conversely, the high-risk subgroup was associated with DNA copy number variation, hypoxia and EMT process. The TICS subgroups significantly predicted overall survival, which was independent of patient age, tumor stage, grade and molecular classification. Moreover, we developed a nomogram incorporating TICS and clinicopathologic factors, which significantly improved the predictive accuracy compared to the clinicopathologic model alone.

CONCLUSION

The TICS is an effective and independent prognostic predictor for EC patients and may serve as a useful supplement to clinicopathological factors and molecular subtyping.

Human endogenous retroviruses of the HERV-K (HML-2) family are expressed in the brain of healthy individuals and modify the composition of the brain-infiltrating immune cells

Human endogenous retroviruses (HERVs) are remnants of ancient retroviral infections in the human genome. RNA expression of individual HERVs has frequently been observed in various pathologic conditions, but some activity can also be seen in healthy individuals, e.g. in the blood. To quantitate the basal expression levels of HERVs in the brain, we now used high-throughput sequencing-based metagenomic analysis to characterize the expression profiles of the HERV-K (HML-2) family proviruses in different brain regions of healthy brain tissue. To this end, RNA-seq data from the Genotype-Tissue Expression (GTEx) project was used. The GTEx project is a public resource to study tissue-specific gene expression and regulation, consisting of a large selection of sequenced samples from different tissues. The GTEx data used in this study consisted of 378 samples taken from 13 brain regions from 55 individuals. The data demonstrated that out of 99 intact proviruses in the family 58 were expressed, but the expression profiles were highly divergent and there were no significant differences in the expression profiles between the various anatomic regions of the brain. It is known that the brain contains a variety of infiltrating immune cells, which are probably of great importance both in the normal defense mechanisms as well as in the various pathogenic processes. Digital cytometry (CIBERSORTx) was used to quantify the proportions of the infiltrating immune cells in the same brain samples. Six most abundant (>5 % of the total population) cell types were observed to be CD4 memory resting T cells, M0 macrophages, plasma cells, CD8 T cells, CD4 memory activated T cells, and monocytes. Analysis of the correlations between the individual HERVs and infiltrating cell types indicated that a cluster of 6 HERVs had a notable correlation signature between T cell type infiltrating cell proportions and HERV RNA expression intensity. The correlations between inflammatory type infiltrating cells were negative or weak. Taken together, these data indicate that the expression of HERVs is associated with a T cell type immunity.

Defining trophoblast injury patterns in the transcriptomes of dysfunctional placentas

Trophoblast injury is central to clinically relevant placenta dysfunction. We hypothesized that the mRNA of primary human trophoblasts, exposed to distinct injuries in vitro, capture transcriptome patterns of placental biopsies obtained from common obstetrical syndromes. We deployed a CIBERSORTx deconvolution method to correlate trophoblastic RNAseq-based expression matrices with the transcriptome of omics-defined placental dysfunction patterns in vivo. We found distinct trophoblast injury patterns in placental biopsies from women with fetal growth restriction and a hypertensive disorder, or in biopsies clustered by their omics analysis. Our RNAseq data are useful for defining the contribution of trophoblast injuries to placental dysfunction syndromes.

Identification of hub genes associated with oxidative stress in heart failure and their correlation with immune infiltration using bioinformatics analysis

Both oxidative stress and the immune response are associated with heart failure (HF). In this study, our aim was to identify the hub genes associated with oxidative stress andimmune infiltration of HF by bioinformatics analysis and experimental verification. The expression profile of GSE36074 was obtained from the Gene Expression Omnibus (GEO) database. The differentially expressed genes (DEGs) were screened by GEO2R. The genes related to oxidative stress were extracted from GeneCards websites. Then, the functional enrichment analysis of oxidative stress-related DEGs (OSRDEGs) was performed using DAVID. In addition, we constructed a protein-protein interaction (PPI) network using the STRING database and screened for hub genes with Cytoscape software. We also used CIBERSORTx to analyze immune infiltration in mice heart tissues between the TAC and Sham groups and explored the correlation between immune cells and hub genes. Finally, the hub genes were carried out using reverse transcription quantitative PCR (RT-qPCR), immunohistochemistry (IHC) and western blot. A total of 136 OSRDEGs were found in GSE36074. Enrichment analysis revealed that these OSRDEGs were enriched in the mitochondrion, HIF-1, FoxO, MAPK and TNF signaling pathway. The five hub genes (Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1) were screened by the cytoHubba plugin. The correlation analysis between immune cells and hub genes showed that Mapk14 was positively correlated with Th2 Cells, while Hif1a and Hsp90ab1exhibited a negative correlation with Th2 Cells; Myc exhibited a negative correlation with Monocytes; whereas, Hsp90aa1 was negatively correlated with NK Resting. Finally, five hub genes were validated by RT-qPCR, IHC and western blot. Mapk14, Hif1a, Myc, Hsp90ab1, and Hsp90aa1 are hub genes of HF and may play a critical role in the oxidative stress of HF. This study may provide new targets for the treatment of HF, and the potential immunotherapies are worthy of further study.

Single-cell transcriptome analysis reveals heterogeneity and convergence of the tumor microenvironment in colorectal cancer

Introduction

Colorectal cancer (CRC) ranks second for mortality and third for morbidity among the most commonly diagnosed cancers worldwide. We aimed to investigate the heterogeneity and convergence of tumor microenvironment (TME) in CRC.

Methods

We analyzed the single-cell RNA sequencing data obtained from the Gene Expression Omnibus (GEO) database and identified 8 major cell types and 25 subgroups derived from tumor, para-tumor and peripheral blood.

Results

In this study, we found that there were significant differences in metabolic patterns, immunophenotypes and transcription factor (TF) regulatory patterns among different subgroups of each major cell type. However, subgroups manifested similar lipid metabolic patterns, immunosuppressive functions and TFs module at the end of the differentiation trajectory in CD8+ T cells, myeloid cells and Fibroblasts. Meanwhile, TFs regulated lipid metabolism and immunosuppressive ligand-receptor pairs were detected by tracing the differentiation trajectory. Based on the cell subgroup fractions calculated by CIBERSORTx and bulk RNA-sequencing data from The Cancer Genome Atlas (TCGA), we constructed an immune risk model and clinical risk model of CRC which presented excellent prognostic value.

Conclusion

This study identified that the differentiation was accompanied by remodeling of lipid metabolism and suppression of immune function, which suggest that lipid remodeling may be an important trigger of immunosuppression. More importantly, our work provides a new perspective for understanding the heterogeneity and convergence of the TME and will aid the development of prognosis and immunotherapies of CRC patients.

Bioinformatics Tools for Bulk Gene Expression Deconvolution in Diabetic Retinopathy

Retinal neovascularization is one of the leading causes of vision loss and a hallmark of proliferative diabetic retinopathy (PDR). The immune system is observed to be involved in the pathogenesis of diabetic retinopathy (DR). The specific immune cell type that contributes to retinal neovascularization can be identified via a bioinformatics analysis of RNA sequencing (RNA-seq) data, known as deconvolution analysis. Previous study has identified the infiltration of macrophages in the retina of rats with hypoxia-induced retinal neovascularization and patients with PDR through a deconvolution algorithm, known as CIBERSORTx. Here, we describe the protocols of using CIBERSORTx to perform the deconvolution analysis and downstream analysis of RNA-seq data.

Role of immune cell infiltration and small molecule drugs in adhesive capsulitis: Novel exploration based on bioinformatics analyses

Background

Adhesive capsulitis (AC) is a type of arthritis that causes shoulder joint pain, stiffness, and limited mobility. The pathogenesis of AC is still controversial. This study aims to explore the role of immune related factors in the occurrence and development of AC.

Methods

The AC dataset was downloaded from Gene Expression Omnibus (GEO) data repository. Differentially expressed immune-related genes (DEIRGs) were obtained based on R package "DESeq2" and Immport database. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) were performed to explore the functional correlation of DEIRGs. MCC method and Least Absolute Shrinkage and Selection Operator (LASSO) regression were conducted to identify the hub genes. The immune cell infiltration in shoulder joint capsule between AC and control was evaluated by CIBERSORTx, and the relationship between hub genes and infiltrating immune cells was analyzed by Spearman's rank correlation. Finally, potential small molecule drugs for AC were screened by the Connectivity Map database (CMap) and further verified by molecular docking.

Results

A total of 137 DEIRGs and eight significantly different types of infiltrating immune cells (M0 macrophages, M1 macrophages, regulatory T cells, Tfh cells, monocytes, activated NK cells, memory resting CD4+T cells and resting dendritic cells) were screened between AC and control tissues. MMP9, FOS, SOCS3, and EGF were identified as potential targets for AC. MMP9 was negatively correlated with memory resting CD4+T cells and activated NK cells, but positively correlated with M0 macrophages. SOCS3 was positively correlated with M1 macrophages. FOS was positively correlated with M1 macrophages. EGF was positively correlated with monocytes. Additionally, dactolisib (ranked first) was identified as a potential small-molecule drug for the targeted therapy of AC.

Conclusions

This is the first study on immune cell infiltration analysis in AC, and these findings may provide a new idea for the diagnosis and treatment of AC.

Identification of hub genes correlated with tumor-associated M1-like macrophage infiltration in soft tissue sarcomas

Soft tissue sarcomas (STS) are a heterogeneous series of tumors that might result in severe disability and death. Tumor-associated M1-like macrophage infiltration plays a critical role in tumor development and progression. This study aimed at identifying the hub genes associated with M1-like macrophage infiltration in STS cells. First, the expression profiles from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were imported to calculate the level of M1-like macrophage infiltration by CIBERSORTx. Afterward, the Kaplan-Meier survival analysis was performed to evaluate the correlation between macrophage infiltration and prognosis. Then, weighted gene co-expression network analysis (WGCNA) and protein-protein interaction analysis of GEO data were applied to identify the key gene related to M1-like macrophage infiltration, followed by the functional analysis using TCGA cohort to validate downstream signaling associated with the gene. Finally, pan-cancer analysis was conducted to investigate the gene function in other types of tumors. We found LCK expression positively related to the M1-like macrophage infiltration level, and it positively regulated the expression level of genes regulated to macrophage polarization, and chemotaxis, including interferon-γ (INF-γ), interleukin-12 (IL12), tumor necrosis factor (TNF), PI3K, NF-κB, and CXCL9, 10, and 11. In summary, an 'LCK-INF-γ/IL-12-TNF/PI3K-NF-κB' axis might exist in STS cells that regulate M1-like macrophage infiltration.

An Integrative Analysis of Identified Schizophrenia-Associated Brain Cell Types and Gene Expression Changes

Schizophrenia (SCZ) is a severe mental disorder that may result in hallucinations, delusions, and extremely disordered thinking. How each cell type in the brain contributes to SCZ occurrence is still unclear. Here, we leveraged the human dorsolateral prefrontal cortex bulk RNA-seq data, then used the RNA-seq deconvolution algorithm CIBERSORTx to generate SCZ brain single-cell RNA-seq data for a comprehensive analysis to understand SCZ-associated brain cell types and gene expression changes. Firstly, we observed that the proportions of brain cell types in SCZ differed from normal samples. Among these cell types, astrocyte, pericyte, and PAX6 cells were found to have a higher proportion in SCZ patients (astrocyte: SCZ = 0.163, control = 0.145, P.adj = 4.9 × 10^-4, effect size = 0.478; pericyte: SCZ = 0.057, control = 0.066, P.adj = 1.1 × 10^-4, effect size = 0.519; PAX6: SCZ = 0.014, control = 0.011, P.adj = 0.014, effect size = 0.377), while the L5/6_IT_CAR3 cells and LAMP5 cells are the exact opposite (L5/6_IT_Car3: SCZ = 0.102, control = 0.108, P.adj = 0.016, effect size = 0.369; LAMP5: SCZ = 0.057, control = 0.066, P.adj = 2.2 × 10^-6, effect size = 0.617). Next, we investigated gene expression in cell types and functional pathways in SCZ. We observed chemical synaptic transmission dysregulation in two types of GABAergic neurons (PVALB and LAMP5), and immune reaction involvement in GABAergic neurons (SST) and non-neuronal cell types (endothelial and oligodendrocyte). Furthermore, we observed that some differential expression genes from bulk RNA-seq displayed cell-type-specific abnormalities in the expression of molecules in SCZ. Finally, the cell types with the SCZ-related transcriptomic changes could be considered to belong to the same module since we observed two major similar coordinated transcriptomic changes across these cell types. Together, our results offer novel insights into cellular heterogeneity and the molecular mechanisms underlying SCZ.

The human blood transcriptome exhibits time-of-day-dependent response to hypoxia: Lessons from the highest city in the world

High altitude exposes humans to hypobaric hypoxia, which induces various physiological and molecular changes. Recent studies point toward interaction between circadian rhythms and the hypoxic response, yet their human relevance is lacking. Here, we examine the effect of different high altitudes in conjunction with time of day on human whole-blood transcriptome upon an expedition to the highest city in the world, La Rinconada, Peru, which is 5,100 m above sea level. We find that high altitude vastly affects the blood transcriptome and, unexpectedly, does not necessarily follow a monotonic response to altitude elevation. Importantly, we observe daily variance in gene expression, especially immune-related genes, which is largely altitude dependent. Moreover, using a digital cytometry approach, we estimate relative changes in abundance of different cell types and find that the response of several immune cell types is time- and altitude dependent. Taken together, our data provide evidence for interaction between the transcriptional response to hypoxia and the time of day in humans.

•

Low oxygen availability upon high altitude vastly affects human blood transcriptome

•

The transcriptomic changes upon altitude elevation are not necessarily monotonic

•

The daily variance in gene expression is dependent on altitude

•

The response of several immune cell types is time- and altitude dependent

Immune cell gene expression signatures in diffuse glioma are associated with IDH mutation status, patient outcome and malignant cell state, and highlight the importance of specific cell subsets in glioma biology

The tumor micro-environment (TME) plays an important role in various cancers, including gliomas. We estimated immune cell type-specific gene expression profiles in 3 large clinically annotated glioma datasets using CIBERSORTx and LM22/LM10 blood-based immune signatures and found that the proportions and estimated gene expression patterns of specific immune cells significantly varied according to IDH mutation status. When IDH-WT and IDH-MUT tumors were considered separately, cluster-of-cluster analyses of immune cell gene expression identified groups with distinct survival outcomes. We confirmed and extended these findings by applying a signature matrix derived from single-cell RNA-sequencing data derived from 19 glioma tumor samples to the bulk profiling data, validating findings from the LM22/LM10 results. To link immune cell signatures with outcomes in checkpoint therapy, we then showed a significant association of monocytic lineage cell gene expression clusters with patient survival and with mesenchymal gene expression scores. Integrating immune cell-based gene expression with previously described malignant cell states in glioma demonstrated that macrophage M0 abundance significantly correlated with mesenchymal state in IDH-WT gliomas, with evidence of a previously implicated role of the Oncostatin-M receptor and macrophages in the mesenchymal state. Among IDH-WT tumors that were enriched for the mesenchymal cell state, the estimated M0 macrophage expression signature coordinately also trended to a mesenchymal signature. We also examined IDH-MUT tumors stratified by 1p/19q status, showing that a mesenchymal gene expression signature the M0 macrophage fraction was enriched in IDH-MUT, non-codeleted tumors. Overall, these results highlight the biological and clinical significance of the immune cell environment related to IDH mutation status, patient prognosis and the mesenchymal state in diffuse gliomas.

Diagnostic Value of Immune-Related Genes in Kawasaki Disease

Kawasaki disease (KD) is a systemic vasculitis that predominantly damages medium- and small-sized vessels, and mainly causes coronary artery lesions (CALs). The diagnostic criterion of KD mainly depends on clinical features, so children could be easily misdiagnosed and could suffer from CALs. Through analysis, a total of 14 immune-related DEGs were obtained, of which IL1B, ADM, PDGFC, and TGFA were identified as diagnostic markers of KD. Compared with the non-KD group, KD patients contained a higher proportion of naive B cells, activated memory CD4 T cells, gamma delta T cells, and neutrophils, while the proportions of memory B cells, CD8 T cells, activated memory CD4 T cells, and activated NK cells were relatively lower. In conclusion, immune-related genes can be used as diagnostic markers of KD, and the difference in immune cells between KD and non-KD might provide new insight into understanding the pathogenesis of KD.

Atlas of clinically distinct cell states and ecosystems across human solid tumors

Determining how cells vary with their local signaling environment and organize into distinct cellular communities is critical for understanding processes as diverse as development, aging, and cancer. Here we introduce EcoTyper, a machine learning framework for large-scale identification and validation of cell states and multicellular communities from bulk, single-cell, and spatially resolved gene expression data. When applied to 12 major cell lineages across 16 types of human carcinoma, EcoTyper identified 69 transcriptionally defined cell states. Most states were specific to neoplastic tissue, ubiquitous across tumor types, and significantly prognostic. By analyzing cell-state co-occurrence patterns, we discovered ten clinically distinct multicellular communities with unexpectedly strong conservation, including three with myeloid and stromal elements linked to adverse survival, one enriched in normal tissue, and two associated with early cancer development. This study elucidates fundamental units of cellular organization in human carcinoma and provides a framework for large-scale profiling of cellular ecosystems in any tissue.

The landscape of tumor cell states and ecosystems in diffuse large B cell lymphoma

Biological heterogeneity in diffuse large B cell lymphoma (DLBCL) is partly driven by cell-of-origin subtypes and associated genomic lesions, but also by diverse cell types and cell states in the tumor microenvironment (TME). However, dissecting these cell states and their clinical relevance at scale remains challenging. Here, we implemented EcoTyper, a machine-learning framework integrating transcriptome deconvolution and single-cell RNA sequencing, to characterize clinically relevant DLBCL cell states and ecosystems. Using this approach, we identified five cell states of malignant B cells that vary in prognostic associations and differentiation status. We also identified striking variation in cell states for 12 other lineages comprising the TME and forming cell state interactions in stereotyped ecosystems. While cell-of-origin subtypes have distinct TME composition, DLBCL ecosystems capture clinical heterogeneity within existing subtypes and extend beyond cell-of-origin and genotypic classes. These results resolve the DLBCL microenvironment at systems-level resolution and identify opportunities for therapeutic targeting (https://ecotyper.stanford.edu/lymphoma).

IDO1 as a new immune biomarker for diabetic nephropathy and its correlation with immune cell infiltration

INTRODUCTION

Indoleamine 2,3-dioxygenase 1(IDO1) has complicated roles in immune-inflammatory response regulation, but its correlation with immune cell infiltration in diabetic nephropathy (DN) remains unknown.

METHODS

Gene expression data were extracted from the GEO database. Differentially expressed genes (DEGs) were identified and functional correlation analysis was performed. The immune hub gene was screened using Maximal Clique Centrality, and verified in DN model mice via western blotting, immunohistochemistry, and immunofluorescence analysis. CIBERSORTx was used to assign values to immune cell infiltration in DN and determine a correlation with the hub gene. The prognostic significance of the hub gene was then validated.

RESULTS

The 330 screened DEGs from the GEO dataset were most enriched in GO functions and KEGG pathways associated with immune inflammation. IDO1 was identified as a hub immune gene, with upregulated expression in DN model mice. IDO1 expression was positively correlated with M1 macrophages (R = 0.58, P < 0.001) and monocytes (R = 0.44, P = 0.049), and was negatively correlated with resting memory CD4 T cells (R = -0.51, P = 0.019). IDO1 expression was upregulated in peritoneal macrophages after high glucose stimulation, and inflammatory factor production was reversed by IDO1 inhibition. Higher IDO1 expression was associated with worse prognosis in DN patients via multivariate survival analysis (P < 0.001).

CONCLUSIONS

IDO1 was identified as a diagnostic and prognostic biomarker for DN and shown to play a vital role in immune cell infiltration in DN, ascertained using microarray data and CIBERSORTx for the first time.

LAMB1 Is Related to the T Stage and Indicates Poor Prognosis in Gastric Cancer

BACKGROUND AND OBJECTIVE

Gastric cancer (GC) is a common tumor malignancy with high incidence and poor prognosis. Laminin is an indispensable component of basement membrane and extracellular matrix, which is responsible for bridging the internal and external environment of cells and transmitting signals. This study mainly explored the association of the LAMB1 expression with clinicopathological characteristics and prognosis in gastric cancer.

METHODS

The expression data and clinical information of gastric cancer patients were downloaded from The Cancer Genome Atlas (TCGA) and Asian Cancer Research Group (ACRG). And we analyzed the relationship between LAMB1 expression and clinical characteristics through R. CIBERSORTx was used to calculate the absolute score of immune cells in gastric tumor tissues. Then COX proportional hazard models and Kaplan-Meier curves were performed to evaluate the role of LAMB1 and its influence on prognosis in gastric cancer patients. Finally, GO and KEGG analysis were applied for LAMB1-related genes in gastric cancer, and PPI network was constructed in Cytoscape software.

RESULTS

In the TCGA cohort, patients with gastric cancer frequently generated LAMB1 gene copy number variation, but had little effect on mRNA expression. Both in the TCGA and ACRG cohorts, the mRNA expression of LAMB1 in gastric cancer tissues was higher than it in normal tissues. All patients were divided into high expression group and low expression group according to the median expression level of LAMB1. The elevated expression group obviously had more advanced cases and higher infiltration levels of M2 macrophages. COX proportional hazard models and Kaplan-Meier curves revealed that patients with enhanced expression of LAMB1 have a worse prognosis. GO/KEGG analysis showed that LAMB1-related genes were enriched in PI3K-Akt signaling pathway, focal adhesion, ECM-receptor interaction, etc.

CONCLUSIONS

The high expression of LAMB1 in gastric cancer is related to the poor prognosis of patients, and it may be related to microenvironmental changes in tumors.

Identification of hub genes associated with neutrophils infiltration in colorectal cancer

Colorectal cancer (CRC) is the leading cause of cancer‐related mortality in the world. Accumulating evidence indicate that tumour infiltrating immune cells participated in cancer progression. Among them, tumour infiltrating neutrophils (TINs) are reported to play crucial role in various cancers. In this study, we used CIBERSORTx, a digital cytometry tool to evaluate the neutrophils infiltration in CRC based on gene expression data of CRC tissues from GSE39582 data set and The Cancer Genome Atlas data set (TCGA‐COAD and TCGA‐READ). Weighted gene co‐expression network analysis (WGCNA) was conducted in GSE39582 data set to identify hub genes associated with neutrophil infiltration. The association of hub gene and neutrophils was then validated in TCGA cohorts and an independent RJ cohort. Functional analysis was performed to investigate the molecular mechanisms of the interested hub gene. We found that neutrophil infiltration is elevated in CRC tissues, and it is related to a poorer prognosis. A total of 18 gene modules are identified by WGCNA in GSE39582 data set, among which lightcyan module is significantly correlated with neutrophils infiltration. Furthermore, Superoxide Dismutase 2 (SOD2) in lightcyan module was proved to correlated with neutrophils infiltration in various cancer types. In addition, SOD2 expression is highly associated with several chemokines, including CXCL8, a neutrophils‐related attractant, and functional analysis revealed that SOD2 is involved in neutrophils recruitment biological process. These results indicate that an ‘SOD2‐CXCL8‐neutrophil recruitment’ axis plays a potential role in colorectal cancer progression.

Identification of the tubulointerstitial infiltrating immune cell landscape and immune marker related molecular patterns in lupus nephritis using bioinformatics analysis

Background

Systemic lupus erythematosus (SLE) is a multisystem autoimmune disease that commonly affects the kidneys. Research into markers that can predict the prognosis of tubulointerstitial lupus nephritis (LN) has been impeded by the lack of well-designed studies.

Methods

In this study, we selected and merged 3 sets of renal biopsy tubulointerstitial data from GSE32591, GSE69438, and GSE127797, including 95 LN and 15 living healthy donors. CIBERSORTx was utilized for differentially infiltrating immune cell (DIIC) analysis. Weighted Gene Co-Expression network analysis (WGCNA) was employed to explore differentially expressed gene (DEG) related modules. Combined WGCNA hub genes and protein-protein interaction (PPI) validation was used for immune marker identification. Lastly, unsupervised clustering was carried out to validate the correlation between these markers and clinical characteristics.

Results

Our findings unveiled TYROBP, C1QB, LAPTM5, CTSS, PTPRC as the 5 immune markers, which were negatively correlated with glomerular filtration rate (GFR). Specifically, the expression levels of TYROBP and C1QB were significantly different between proliferative LN (PLN) and membranous LN (MLN). Unsupervised clustering could aggregate LN by these immune marker expression spectrums.

Conclusions

This study is the first to identify infiltrating immune cells and associated molecular patterns in the tubulointerstitium of LN by utilizing bioinformatics methods. These findings contribute to a better understanding of the mechanisms behind LN, and promote more precise diagnosis.

Evaluating Distribution and Prognostic Value of New Tumor-Infiltrating Lymphocytes in HCC Based on a scRNA-Seq Study With CIBERSORTx

Hepatocellular carcinoma (HCC) is a commonly diagnosed cancer with high mortality rates. The immune response plays an important role in the progression of HCC. Immunotherapies are becoming an increasingly promising tool for treating cancers. Advancements in scRNA-seq (single-cell RNA sequencing) have allowed us to identify new subsets in the immune microenvironment of HCC. Yet, distribution of these new cell types and their potential prognostic value in bulk samples from large cohorts remained unclear. This study aimed to investigate the tumor-infiltration and prognostic value of new cell subsets identified by a previous scRNA-seq study in a TCGA HCC cohort using CIBERSORTx, a machine learning method to estimate cell proportion and infer cell-type-specific gene expression profiles. We observed different distributions of tumor-infiltrating lymphocytes between tumor and normal cells. Among these, the CD4-GZMA cell subset showed association with prognosis (log-rank test, p < 0.05). We further analyzed CD4-GZMA cell specific gene expression with CIBERSORTx, and found 19 prognostic genes (univariable cox regression, p < 0.05). Finally, we applied Least absolute shrinkage and selection operator (LASSO) Cox regression to construct an immune risk score model and performed a prognostic assessment of our model in TCGA and ICGC cohorts. Taken together, the immune landscape in HCC bulk samples may be more complex than assumed, with heterogeneity and different tumor-infiltration relative to scRNA-seq results. Additionally, CD4-GZMA cells and their characteristics may yield therapeutic benefits in the immune treatment of HCC.

High Abundance of Intratumoral γδ T Cells Favors a Better Prognosis in Head and Neck Squamous Cell Carcinoma: A Bioinformatic Analysis

γδ T cells are a small subset of unconventional T cells that are enriched in the mucosal areas, and are responsible for pathogen clearance and maintaining integrity. However, the role of γδ T cells in head and neck squamous cell carcinoma (HNSCC) is largely unknown. Here, by using RNA-seq data from The Cancer Genome Atlas (TCGA), we discovered that HNSCC patients with higher levels of γδ T cells were positively associated with lower clinical stages and better overall survival, and high abundance of γδ T cells was positively correlated with CD8+/CD4+ T cell infiltration. Gene ontology and pathway analyses showed that genes associated with T cell activation, proliferation, effector functions, cytotoxicity, and chemokine production were enriched in the group with a higher γδ T cell abundance. Furthermore, we found that the abundance of γδ T cells was positively associated with the expression of the butyrophilin (BTN) family proteins BTN3A1/BTN3A2/BTN3A3 and BTN2A1, but only MICB, one of the ligands of NKG2D, was involved in the activation of γδ T cells, indicating that the BTN family proteins might be involved in the activation and proliferation of γδ T cells in the tumor microenvironment of HNSCC. Our results indicated that γδ T cells, along with their ligands, are promising targets in HNSCC with great prognostic values and treatment potentials.

Identification of Immune Cell Landscape and Construction of a Novel Diagnostic Nomogram for Crohn's Disease

Crohn’s disease (CD) has an increasing incidence and prevalence worldwide. The etiology of CD remains unclear and there is no gold standard for diagnosis. The dysregulated immune response and different infiltration status of immune cells are critical for CD pathogenesis; therefore, it is important to provide an overview of immune-cell alterations in CD and explore a novel method for auxiliary diagnosis. Here we analyzed microarray datasets from Gene Expression Omnibus (GEO), and an extended version of Cell-type Identification By Estimating Relative Subsets Of RNA Transcripts (CIBERSORTx) was utilized to estimate the fraction of 22 types of immune cells. Differentially expressed genes (DEGs) and a protein-protein interaction (PPI) network were identified, and we performed gene set enrichment analysis (GSEA) and gene set variation analysis (GSVA) to identify differentially regulated pathways in CD. Least absolute shrinkage and selection operator (LASSO) regression was conducted to filter features, and a diagnostic nomogram based on logistic regression was built and validated in an independent validation cohort. In the derivation cohort, we found a proportion of 17 immune-cell types to be significantly altered between CD and healthy controls and a total of 150 DEGs were identified, which were mostly related to the immune response. Among the 15 hub genes based on the PPI network, C-X-C chemokine ligand 8 (CXCL8) and interleukin-1B (IL-1B) showed the highest degree of interaction. Additionally, GSEA and GSVA identified five significantly enriched pathways, among which the nucleotide-binding oligomerization domain (NOD)-like receptor signaling pathway was critical in the CD development. Furthermore, six variables comprising of CXCL8, IL-1B, M1 macrophages, regulatory T cells, CD8⁺ T cells, and plasma cells were identified by LASSO regression and incorporated into a logistic regression model. The nomogram displayed a good prediction, with a 0.915 area under the receiver operating curve (AUC) and a C-index of 0.915 [95% confidence interval (CI): 0.875–0.955]. Similar results were found in the validation cohort, with an AUC of 0.884 and a 0.884 C-index (95% CI: 0.843–0.924). These results provide novel in silico insight into cellular and molecular characteristics of CD and potential biomarkers for diagnosis and targeted therapy.

Deconvolution of RNA-Seq Analysis of Hyperbaric Oxygen-Treated Mice Lungs Reveals Mesenchymal Cell Subtype Changes

Hyperbaric oxygen (HBO) is widely applied to treat several hypoxia-related diseases. Previous studies have focused on the immediate effect of HBO-exposure induced oxidative stress on the lungs, but knowledge regarding the chronic effects from repetitive HBO exposure is limited, especially at the gene expression level. We found that repetitive HBO exposure did not alter the morphology of murine lungs. However, by deconvolution of RNA-seq from those mice lungs using CIBERSORTx and the expression profile matrices of 8 mesenchymal cell subtypes obtained from bleomycin-treated mouse lungs, we identify several mesenchymal cell subtype changes. These include increases in Col13a1 matrix fibroblasts, mesenchymal progenitors and mesothelial cell populations and decreases in lipofibroblasts, endothelial and Pdgfrb high cell populations. Our data suggest that repetitive HBO exposure may affect biological processes in the lungs such as response to wounding, extracellular matrix, vasculature development and immune response.