Applications of single-cell and bulk RNA sequencing in onco-immunology

Mechanistic insights into Hippo-YAP pathway activation for enhanced DFU healing

With the increasing prevalence of diabetes, diabetic foot ulcers (DFUs) have become a global health challenge, significantly impacting patients' quality of life and placing a substantial burden on healthcare systems. Among various immune cell subsets, M2-polarized macrophages play a pivotal role in tissue repair and inflammation resolution. This study uses single-cell RNA sequencing (scRNA-seq) and bulk RNA sequencing to comprehensively investigate the role of the TFAP2A-LIFR-Hippo-YAP signaling axis in regulating macrophage M2 polarization and its critical function in DFU wound healing. Through scRNA-seq analysis, we identified nine major immune cell subsets in DFU samples, with macrophages emerging as key regulatory cells. In vitro experiments further confirmed that TFAP2A promotes macrophage M2 polarization (evidenced by increased expression of the M2 marker ARG1) and ameliorates endothelial dysfunction by enhancing tube formation, improving migration capacity, and upregulating relevant proteins such as VCAM-1. Moreover, TFAP2A serves as a central regulatory gene for macrophage function in DFU by upregulating LIFR expression and activating the Hippo-YAP signaling pathway, thereby inducing M2 polarization and mitigating endothelial dysfunction. Mouse model experiments further demonstrated that the TFAP2A-LIFR-Hippo-YAP signaling axis accelerates DFU wound healing through the induction of macrophage M2 polarization. This study unveils a novel immunoregulatory role of TFAP2A in DFU and provides a promising therapeutic target for the treatment of chronic diabetic wounds.NEW & NOTEWORTHY This study provides unprecedented insights into diabetic foot ulcer healing by demonstrating the novel immunoregulatory role of the TFAP2A-LIFR-Hippo-YAP signaling axis. Leveraging single-cell and bulk transcriptomics, we identify TFAP2A as a crucial regulator of macrophage M2 polarization, essential for wound healing and angiogenesis. These findings offer valuable mechanistic understanding and present TFAP2A as a promising therapeutic target for improving outcomes in chronic diabetic wounds.

Construction of a stromal cell-related prognostic signature based on a 101-combination machine learning framework for predicting prognosis and immunotherapy response in triple-negative breast cancer

Background

Triple-negative breast cancer (TNBC) is a highly aggressive subtype with limited therapeutic targets and poor immunotherapy outcomes. The tumor microenvironment (TME) plays a key role in cancer progression. Advances in single-cell transcriptomics have highlighted the impact of stromal cells on tumor progression, immune suppression, and immunotherapy. This study aims to identify stromal cell marker genes and develop a prognostic signature for predicting TNBC survival outcomes and immunotherapy response.

Methods

Single-cell RNA sequencing (scRNA-seq) datasets were retrieved from the Gene Expression Omnibus (GEO) database and annotated using known marker genes. Cell types preferentially distributed in TNBC were identified using odds ratios (OR). Bulk transcriptome data were analyzed using Weighted correlation network analysis (WGCNA) to identify myCAF-, VSMC-, and Pericyte-related genes (MVPRGs). A consensus MVP cell-related signature (MVPRS) was developed using 10 machine learning algorithms and 101 model combinations and validated in training and validation cohorts. Immune infiltration and immunotherapy response were assessed using CIBERSORT, ssGSEA, TIDE, IPS scores, and an independent cohort (GSE91061). FN1, a key gene in the model, was validated through qRT-PCR, immunohistochemistry, RNA interference, CCK-8 assay, apoptosis assay and wound-healing assay.

Results

In TNBC, three stromal cell subpopulations-myofibroblastic cancer-associated fibroblasts (myCAF), vascular smooth muscle cells (VSMCs), and pericytes-were enriched, exhibiting high interaction frequencies and strong associations with poor prognosis. A nine-gene prognostic model (MVPRS), developed from 23 prognostically significant genes among the 259 MVPRGs, demonstrated excellent predictive performance and was validated as an independent prognostic factor. A nomogram integrating MVPRS, age, stage, and tumor grade offered clinical utility. High-risk group showed reduced immune infiltration and increased activity in tumor-related pathways like ANGIOGENESIS and HYPOXIA, while low-risk groups responded better to immunotherapy based on TIDE and IPS scores. FN1, identified as a key oncogene, was highly expressed in TNBC tissues and cell lines, promoting proliferation and migration while inhibiting apoptosis.

Conclusion

This study reveals TNBC microenvironment heterogeneity and introduces a prognostic signature based on myCAF, VSMC, and Pericyte marker genes. MVPRS effectively predicts TNBC prognosis and immunotherapy response, providing guidance for personalized treatment. FN1 was validated as a key oncogene impacting TNBC progression and malignant phenotype, with potential as a therapeutic target.

Development and Validation of an Extracellular Matrix Gene Expression Signature for Prognostic Prediction in Patients with Uveal Melanoma

Uveal melanoma (UVM) is an aggressive cancer with a poor prognosis, particularly in metastatic cases. This study aimed to develop and validate a novel extracellular matrix (ECM) gene expression signature to predict prognosis and stratify patients by risk. ECM-related genes were identified and used to construct a prognostic model through Lasso-Cox regression analysis, leveraging RNA sequencing data from 80 UVM patients in The Cancer Genome Atlas (TCGA). The model was validated using an independent cohort of 63 UVM patients. Survival analyses, immune infiltration profiling, and functional enrichment analyses were conducted to evaluate the biological significance and clinical utility of the signature. The ECM signature stratified patients into high- and low-risk groups with significant differences in survival outcomes. High-risk patients showed elevated expression of MMP1 and MMP12, which are associated with ECM remodeling and immune modulation, alongside increased infiltration of immunosuppressive cells, such as M2 macrophages. Validation confirmed the prognostic value of the signature across cohorts. Functional analyses highlighted the involvement of ECM-related pathways, epithelial-mesenchymal transition, and immune system interactions in tumor progression. This ECM gene expression signature is a robust prognostic tool for UVM, offering insights into tumor biology and immune microenvironment interactions. It holds promise for improving patient stratification and guiding personalized therapeutic strategies. Further research is warranted to explore the functional roles of these genes in UVM progression.

Comprehensive bioinformatics analysis identified HMGB3 as a promising immunotherapy target for glioblastoma multiforme

OBJECTIVE

Glioblastoma multiforme (GBM) presents significant therapeutic challenges due to its heterogeneous tumorigenicity, drug resistance, and immunosuppression. Although several molecular markers have been developed, there still lack of sensitive molecular for accurately detection. Studying the mechanisms underlying the development of GBM and finding relevant prognostic biomarkers remains crucial.

METHODS

Single-cell RNA sequencing, bulk RNA-seq, and cancer immune cycle activities of GBM were used to assess the expression of different molecular related to GBM. Bioinformatics analyses were carried to evaluate the functional of the high mobility group protein B3 (HMGB3) in GBM.

RESULTS

HMGB3 was highly expressed in GBM tissues and influenced the interpatient and intratumoral transcriptomic heterogeneity as well as immunosuppression in GBM. HMGB3 also contributes to a no inflamed tumor microenvironment (TME) and has an inhibitory effect on tumor-associated immune cell infiltration. Besides, HMGB3 participated GBM chemotherapeutic sensitivity and negative correlation with 140 medicines.

CONCLUSION

HMGB3 as a heterogeneous and immunosuppressive molecule in the GBM TME, making it a potential target for precision therapy for GBM.

Deciphering IGH rearrangement complexity and detection strategies in acute lymphoblastic leukaemia

Acute lymphoblastic leukaemia is a highly heterogeneous malignancy characterised by various genomic alterations that influence disease progression and therapeutic outcomes. Gene fusions involving the immunoglobulin heavy chain gene represent a complex and diverse category. These fusions often result in enhancer hijacking, upregulation of partner proto-oncogenes and contribute to leukemogenesis. This review highlights the mechanisms underlying IGH gene fusions, the critical role they play in ALL pathogenesis, and current detection technologies.

Development of a MVI associated HCC prognostic model through single cell transcriptomic analysis and 101 machine learning algorithms

Hepatocellular carcinoma (HCC) is an exceedingly aggressive form of cancer that often carries a poor prognosis, especially when it is complicated by the presence of microvascular invasion (MVI). Identifying patients at high risk of MVI is crucial for personalized treatment strategies. Utilizing the single-cell RNA-sequencing dataset (GSE242889) of HCC, we identified malignant cell subtypes associated with microvascular invasion (MVI), in conjunction with the TCGA dataset, selected a set of MVI-related genes (MRGs). We developed an optimal prognostic model comprising 11 genes (NOP16, YIPF1, HMMR, NDC80, DYNLL1, CDC34, NLN, KHDRBS3, MED8, SLC35G2, RAB3B) based on MVI-related signature genes by integrating single-cell transcriptomic analysis with 101 machine learning algorithms. This model is meticulously crafted to forecast the prognosis of individuals afflicted with hepatocellular carcinoma (HCC). Additionally, we affirmed the predictive precision and superiority of our model through a meta-analysis against existing HCC models. Furthermore, we explored the differences between high- and low-risk groups through mutation and immune infiltration analyses. Lastly, we investigated immunotherapy responses and drug sensitivities between risk groups, providing novel therapeutic insights for liver cancer.

Single-cell RNA sequencing and machine learning provide candidate drugs against drug-tolerant persister cells in colorectal cancer

Drug resistance often stems from drug-tolerant persister (DTP) cells in cancer. These cells arise from various lineages and exhibit complex dynamics. However, effectively targeting DTP cells remains challenging. We used single-cell RNA sequencing (scRNA-Seq) data and machine learning (ML) models to identify DTP cells in patient-derived organoids (PDOs) and computationally screened candidate drugs targeting these cells in familial adenomatous polyposis (FAP), associated with a high risk of colorectal cancer. Three PDOs (benign and malignant tumor organoids and a normal organoid) were evaluated using scRNA-Seq. ML models constructed based on public scRNA-Seq data classified DTP versus non-DTP cells. Candidate drugs for DTP cells in a malignant tumor organoid were identified from public drug sensitivity data. From FAP scRNA-Seq data, a specific TC1 cell cluster in tumor organoids was identified. The ML model identified up to 36 % of TC1 cells as DTP cells, a higher proportion than those for other clusters. A viability assay using a malignant tumor organoid demonstrated that YM-155 and THZ2 exert synergistic effects with trametinib. The constructed ML model is effective for DTP cell identification based on scRNA-Seq data for FAP and provides candidate treatments. This approach may improve DTP cell targeting in the treatment of colorectal and other cancers.

Single-Cell Sequencing: Genomic and Transcriptomic Approaches in Cancer Cell Biology

This article reviews the impact of single-cell sequencing (SCS) on cancer biology research. SCS has revolutionized our understanding of cancer and tumor heterogeneity, clonal evolution, and the complex interplay between cancer cells and tumor microenvironment. SCS provides high-resolution profiling of individual cells in genomic, transcriptomic, and epigenomic landscapes, facilitating the detection of rare mutations, the characterization of cellular diversity, and the integration of molecular data with phenotypic traits. The integration of SCS with multi-omics has provided a multidimensional view of cellular states and regulatory mechanisms in cancer, uncovering novel regulatory mechanisms and therapeutic targets. Advances in computational tools, artificial intelligence (AI), and machine learning have been crucial in interpreting the vast amounts of data generated, leading to the identification of new biomarkers and the development of predictive models for patient stratification. Furthermore, there have been emerging technologies such as spatial transcriptomics and in situ sequencing, which promise to further enhance our understanding of tumor microenvironment organization and cellular interactions. As SCS and its related technologies continue to advance, they are expected to drive significant advances in personalized cancer diagnostics, prognosis, and therapy, ultimately improving patient outcomes in the era of precision oncology.

Integration of single-cell and bulk RNA-sequencing data reveals the prognostic potential of epithelial gene markers for prostate cancer

Prognostic transcriptomic signatures for prostate cancer (PCa) often overlook the cellular origin of expression changes, an important consideration given the heterogeneity of the disorder. Current clinicopathological factors inadequately predict biochemical recurrence, a critical indicator guiding post-treatment strategies following radical prostatectomy. To address this, we conducted a meta-analysis of four large-scale PCa datasets and found 33 previously reported PCa-associated genes to be consistently up-regulated in prostate tumours. By analysing single-cell RNA-sequencing data, we found these genes predominantly as markers in epithelial cells. Subsequently, we applied 97 advanced machine-learning algorithms across five PCa cohorts and developed an 11-gene epithelial expression signature. This signature robustly predicted biochemical recurrence-free survival (BCRFS) and stratified patients into distinct risk categories, with high-risk patients showing worse survival and altered immune cell populations. The signature outperformed traditional clinical parameters in larger cohorts and was overall superior to published PCa signatures for BCRFS. By analysing peripheral blood data, four of our signature genes showed potential as biomarkers for radiation response in patients with localised cancer and effectively stratified castration-resistant patients for overall survival. In conclusion, this study developed a novel epithelial gene-expression signature that enhanced BCRFS prediction and enabled effective risk stratification compared to existing clinical- and gene-expression-derived prognostic tools. Furthermore, a set of genes from the signature demonstrated potential utility in peripheral blood, a tissue amenable to minimally invasive sampling in a primary care setting, offering significant prognostic value for PCa patients without requiring a tumour biopsy.

CCR7 in esophageal squamous cell carcinoma: an identification from single-cell and bulk transcriptome sequencing

BACKGROUND

Considering the involvement of CC-chemokine receptor type 7 (CCR7) in diverse tumors, the purpose of our current study is to reveal the specific mechanisms of CCR7 in esophageal squamous cell carcinoma (ESCC).

METHODS

We processed single-cell RNA sequencing data based on the Gene Expression Omnibus (GEO) database and utilizing Seurat software, and performed filtering and annotation to obtain different cell types. Genes related to gene ontology biological process and Hallmark were collected, and the enrichment of genes of interested in both single cell and TCGA-ESCA was quantified. Besides, genes related to CCR7 and KRAS signaling up were uploaded to construct the protein-protein interaction network. A series of cellular assays were incorporated to test the effects of CCR7 in ESCC cells.

RESULTS

28,281 cells were categorized into 4 non-immune cell classes (epithelial cells, smooth muscle cells, fibroblasts, endothelial cells) and 6 immune cell classes (mast cells, plasma B cells, B cells, neutrophils, macrophages, NK/T cells). CCR7 evidently expressed higher in epithelial cells of ESCA and was positively correlated with KRAS signaling up, inflammatory response and TNFA signaling via NFKB, with the most significant correlation witnessed between CCR7 and KRAS signaling up. Meanwhile, the positive correlation between KRAS signaling up and positive regulation of epithelial cell migration was observed. 12 common genes were related to both KRAS signaling up and positive regulation of epithelial cell migration, 3 of which (SOX9, FLT4 and HDAC9) were higher-expressed in M1 group of TCGA-ESCA. Besides, CCR7 as well as its ligands CCL19 and CCL21 was shown to express higher in ESCC cells, where increased level of immune response-related cytokines was seen. CCR7 knockdown diminished migration and proliferation as well as Macrophage M2 polarization.

CONCLUSION

CCR7 was highly-expressed in ESCC and positively correlated with KRAS signaling up, which may contribute to the migration of ESCC cells.

Advancements in single-cell RNA sequencing and spatial transcriptomics: transforming biomedical research

In recent years, significant advancements in biochemistry, materials science, engineering, and computer-aided testing have driven the development of high-throughput tools for profiling genetic information. Single-cell RNA sequencing (scRNA-seq) technologies have established themselves as key tools for dissecting genetic sequences at the level of single cells. These technologies reveal cellular diversity and allow for the exploration of cell states and transformations with exceptional resolution. Unlike bulk sequencing, which provides population-averaged data, scRNA-seq can detect cell subtypes or gene expression variations that would otherwise be overlooked. However, a key limitation of scRNA-seq is its inability to preserve spatial information about the RNA transcriptome, as the process requires tissue dissociation and cell isolation. Spatial transcriptomics is a pivotal advancement in medical biotechnology, facilitating the identification of molecules such as RNA in their original spatial context within tissue sections at the single-cell level. This capability offers a substantial advantage over traditional single-cell sequencing techniques. Spatial transcriptomics offers valuable insights into a wide range of biomedical fields, including neurology, embryology, cancer research, immunology, and histology. This review highlights single-cell sequencing approaches, recent technological developments, associated challenges, various techniques for expression data analysis, and their applications in disciplines such as cancer research, microbiology, neuroscience, reproductive biology, and immunology. It highlights the critical role of single-cell sequencing tools in characterizing the dynamic nature of individual cells.

Rewired chromatin structure and epigenetic gene dysregulation during HTLV-1 infection to leukemogenesis

Human T-cell leukemia virus type 1 (HTLV-1) broadly impacts host genes, affecting the infected cell population and inducing the development of a disease with a poor prognosis, adult T-cell leukemia-lymphoma (ATL). This study aimed to provide a comprehensive epigenomic characterization of the infected cell population and evaluated the transcriptome and chromatin structures of peripheral blood cells in HTLV-1-infected individuals using RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin sequencing (ATAC-seq). The infected cells showed significant changes in gene expression patterns from the polyclonal stage and before ATL onset while demonstrating similarities to tumor-forming ATL cells. These similarities were a result of large-scale open chromatin changes, supporting the independent early formation of epigenomic aberrations as an underlying mechanism for later clonal propagation. This study also demonstrated that HTLV-1 Tax directly affects the host chromatin structure, thereby developing fundamental epigenomic characteristics. Several Tax target genes, including the RASGRP3-ERK pathway, were recognized, indicating an impact on signaling pathways. This genome-wide variability in chromatin structural property is a novel feature of HTLV-1 infection and may contribute to pathogenic mechanisms. In addition, it has crucial implications for better understanding the impact of HTLV-1 on the host genome and identifying novel therapeutic targets.

Transcriptomic data integration and analysis revealing potential mechanisms of doxorubicin resistance in chondrosarcoma cells

Chondrosarcoma is a type of bone cancer that originates from cartilage cells. In clinical practice, surgical resection is the primary treatment for chondrosarcoma, but chemotherapy becomes essential for patients with metastasis or tumors in surgically inaccessible sites. However, drug resistance often leads to treatment failure. Tumor microenvironment proteins modulate intercellular communication, contributing to drug resistance. Doxorubicin (Dox) is a common chemotherapeutic agent. The present study aimed to establish Dox-resistant chondrosarcoma cells and compare their secretome with parental cells using antibody arrays. Results showed significantly heightened secretion of hepatocyte growth factor (HGF). Knockdown of both HGF and its receptor MET increased Dox sensitivity in chondrosarcoma cells. Treatment of chondrosarcoma cells with conditioned media (CM) from cells secreting high levels of HGF resulted in MET activation. Additionally, the expression levels of HGF and MET were significantly elevated in chondrosarcoma tissues compared to normal cartilage tissues, as confirmed by analysis of GEO database. RNA sequencing and Gene Set Enrichment Analysis (GSEA) elucidated the mechanism involving HGF. Additionally, genes with log fold change > 1 underwent bioinformatics analysis using the ShinyGO web server. The results from both GSEA and ShinyGO analyses corroborate each other, indicating the significance of HGF in cellular signal transduction, regulation of cell motility, developmental processes, immune-inflammatory responses, and functions related to blood and neural systems. In summary, highly secreted HGF can activate signaling pathways through its receptor MET, particularly Ras and Akt activation, enhancing drug resistance in chondrosarcoma cells. The present study may guide the development of novel therapeutic strategies targeting HGF, ultimately improving treatment outcomes and prognosis for malignant chondrosarcoma patients.

Epithelial and immune transcriptomic characteristics and possible regulatory mechanisms in asthma exacerbation: insights from integrated studies

Background

Asthma exacerbation significantly contribute to disease mortality and result in heightened health care expenditures. This study was aimed at gaining important new insights into the heterogeneity of epithelial and immune cells and elucidating key regulatory genes involved in the pathogenesis of asthma exacerbation.

Methods

Functional enrichment, pseudotime, metabolism and cell-cell communication analyses of epithelial cells and immune cells in single-cell RNA sequencing (scRNA-seq) dataset were applied. Immune infiltration analysis was performed in bulk RNA sequencing (bulk RNA-seq) dataset. Key regulatory genes were obtained by taking the intersection of the differentially expressed genes (DEGs) between control and asthma group in epithelial cells, immune cells and bulk RNA-seq data. Asthma animal and in vitro cell line models were established to verify the key regulatory genes expression by employing quantitative reverse transcription polymerase chain reaction (qRT-PCR).

Results

ScRNA-seq analysis identified 7 epithelial subpopulations and 14 distinct immune cell types based on gene expression profiles. Further analysis demonstrated that these cells manifested high heterogeneity at the levels of functional variations, dynamics, communication patterns and metabolic changes. Notably, TMPRSS11A, TUBA1A, SCEL, ICAM4, TMPRSS11B, IGFBP2, CLC, NFAM1 and F13A1 were identified as key regulatory genes of asthma. The results of the qRT-PCR demonstrated that the 9 key regulatory genes were involved in asthma.

Conclusions

We systematically explored epithelial and immune characteristics in asthma exacerbation and identified 9 key regulatory genes underlying asthma occurrence and progression, which may be valuable for providing new insights into the cellular and molecular mechanisms driving asthma exacerbations.

Single-cell RNA-sequencing and genome-wide Mendelian randomisation along with abundant machine learning methods identify a novel B cells signature in gastric cancer

BACKGROUND

Gastric cancer (GC) has a poor prognosis, considerable cellular heterogeneity, and ranks fifth among malignant tumours. Understanding the tumour microenvironment (TME) and intra-tumor heterogeneity (ITH) may lead to the development of novel GC treatments.

METHODS

The single-cell RNA sequencing (scRNA-seq) dataset was obtained from the Gene Expression Omnibus (GEO) database, where diverse immune cells were isolated and re-annotated based on cell markers established in the original study to ascertain their individual characteristics. We conducted a weighted gene co-expression network analysis (WGCNA) to identify genes with a significant correlation to GC. Utilising bulk RNA sequencing data, we employed machine learning integration methods to train specific biomarkers for the development of novel diagnostic combinations. A two-sample Mendelian randomisation study was performed to investigate the causal effect of biomarkers on gastric cancer (GC). Ultimately, we utilised the DSigDB database to acquire associations between signature genes and pharmaceuticals.

RESULTS

The 18 genes that made up the signature were as follows: ZFAND2A, PBX4, RAMP2, NNMT, RNASE1, CD93, CDH5, NFKBIE, VWF, DAB2, FAAH2, VAT1, MRAS, TSPAN4, EPAS1, AFAP1L1, DNM3. Patients were categorised into high-risk and low-risk groups according to their risk scores. Individuals in the high-risk cohort exhibited a dismal outlook. The Mendelian randomisation study demonstrated that individuals with a genetic predisposition for elevated NFKBIE levels exhibited a heightened likelihood of acquiring GC. Molecular docking indicates that gemcitabine and chloropyramine may serve as effective therapeutics against NFKBIE.

CONCLUSIONS

We developed and validated a signature utilising scRNA-seq and bulk sequencing data from gastric cancer patients. NFKBIE may function as a novel biomarker and therapeutic target for GC.

Enhanced fibrotic potential of COL1A1hiNR4A1low fibroblasts in ischemic heart revealed by transcriptional dynamics heterogeneity analysis at both bulk and single-cell levels

Background

Fibroblasts in the fibrotic heart exhibit a heterogeneous biological behavior. The specific subsets of fibroblasts that contribute to progressive cardiac fibrosis remain unrevealed. Our aim is to identify the heart fibroblast (FB) subsets that most significantly promote fibrosis and the related critical genes as biomarkers for ischemic heart disease.

Methods

The single nuclei RNA sequencing (snRNA-seq) and bulk RNA sequencing datasets used in this study were obtained from the Gene Expression Omnibus (GEO). The activity of gene sets related to progressive fibrosis was quantified for each FB cluster using the AddmoleculeScore function. Differentially expressed genes (DEGs) for the specific cell cluster with the highest fibrotic transcription dynamics were identified and integrated with bulk RNA sequencing data for analysis. Multiple machine learning models were employed to identify the optimal gene panel for diagnosing ischemic heart disease (IHD) based on the intersected DEGs. The effectiveness and robustness of the gene-derived diagnostic tool were validated using two independent IHD cohorts.Subsequently, we validated the signature genes using a rat post-myocardial infarction heart failure model.

Results

We conducted an analysis on high-quality snRNA-seq data obtained from 3 IHD and 4 cardiac sarcoidosis heart samples, resulting in the identification of 16 FB clusters. Cluster2 exhibited the highest gene activity in terms of fibrosis-related transcriptome dynamics. The characteristic gene expression profile of this FB subset indicated a specific upregulation of COL1A1 and several pro-fibrotic factors, including CCDC102B, GUCY1A3, TEX41, NREP, TCAP, and WISP, while showing a downregulation of NR4A1, an endogenous inhibitor of the TGF-β pathway. Consequently, we designated this subgroup as COL1A1hiNR4A1low FB. Gene set enrichment analysis (GSEA) shows that the gene expression pattern of COL1A1hiNR4A1low FB was closer to pathways associated with cardiac fibrosis. Through machine learning, ten feature genes from COL1A1hiNR4A1low FB were selected to construct a diagnostic tool for IHD. The robustness of this new tool was validated using an independent cohort and heart failure rats.

Conclusion

COL1A1hiNR4A1low FB possess heightened capability in promoting cardiac fibrosis. Additionally, it offers molecular insights into the mechanisms underlying the regulation of the TGF-β pathway. Furthermore, the characteristic genes of COL1A1hiNR4A1 FB could serve as valuable tools for diagnosing of IHD.

T-cell receptor and B-cell receptor repertoires profiling in pleural tuberculosis

Background

Tuberculosis (TB) is a leading cause of death worldwide from a single infectious agent. In China the most common extra-pulmonary TB (EPTB) is pleural tuberculosis (PLTB). An important clinical feature of PLTB is that the lymphocytes associated with TB will accumulate in the pleural fluid. The adaptive immune repertoires play important roles in Mycobacterium tuberculosis (Mtb) infection.

Methods

In this study, 10 PLTB patients were enrolled, and their Peripheral Blood Mononuclear Cells(PBMCs) and Pleural Effusion Mononuclear Cells(PEMCs) were collected. After T cells were purified from PBMCs and PEMCs, high-throughput immunosequencing of the TCRβ chain (TRB), TCRγ chain(TRG), and B cell receptor(BCR) immunoglobulin heavy chain (IGH) were conducted on these samples.

Results

The TRB, TRG, and BCR IGH repertoires were characterized between the pleural effusion and blood in PLTB patients, and the shared clones were analyzed and collected. The binding activity of antibodies in plasma and pleural effusion to Mtb antigens was tested which indicates that different antibodies responses to Mtb antigens in plasma and pleural effusion in PLTB patients. Moreover, GLIPH2 was used to identify the specificity groups of TRB clusters and Mtb-specific TRB sequences were analyzed and collected by VJ mapping.

Conclusion

We characterize the adaptive immune repertoires and identify the shared clones and Mtb-specific clones in pleural effusion and blood in PLTB patients which can give important clues for TB diagnosis, treatment, and vaccine development.

Single cell landscape of potential mechanisms in primary and metastatic hepatocellular carcinoma

Hepatocellular carcinoma metastasis occurs mainly in the portal vein and lymph node metastasis, and is the main cause of patient death. However, the cellular origins and driving mechanisms of hepatocellular carcinoma metastasis have not been elucidated. In this study, we found that different tumor metastasis samples originated from different branches of the primary tumor subclone. A large number of sequence data confirmed the correlation between tumor metastasis index and metastasis and prognosis. Patients with a high index generally had a poor prognosis and abnormal metabolic function. Finally, we recommended a candidate therapy for each metastatic direction. This study explains the cellular origin and underlying mechanisms of hepatocellular carcinoma metastasis at the single-cell level and identifies drugs for its targeted therapy. It also provides a new research idea for the study of hepatocellular carcinoma metastasis and early identification of cancer metastasis from cellular resolution.

Identifying cancer prognosis genes through causal learning

Accurate identification of causal genes for cancer prognosis is critical for estimating disease progression and guiding treatment interventions. In this study, we propose CPCG (Cancer Prognosis's Causal Gene), a two-stage framework identifying gene sets causally associated with patient prognosis across diverse cancer types using transcriptomic data. Initially, an ensemble approach models gene expression's impact on survival with parametric and semiparametric hazard models. Subsequently, an iterative conditional independence test combined with graph pruning is utilized to infer the causal skeleton, thereby pinpointing prognosis-related genes. Experiments on transcriptomic data from 18 cancer types sourced from The Cancer Genome Atlas Project demonstrate CPCG's effectiveness in predicting prognosis under four evaluation metrics. Validations on 24 additional datasets covering 12 cancer types from the Gene Expression Omnibus and the Chinese Glioma Genome Atlas Project further demonstrate CPCG's robustness and generalizability. CPCG identifies a concise but reliable set of genes, obviating the need for gene combination enumeration for survival time estimation. These genes are also proved closely linked to crucial biological processes in cancer. Moreover, CPCG constructs a stable causal skeleton and exhibits insensitivity to the order of data shuffling. Overall, CPCG is a powerful tool for extracting cancer prognostic biomarkers, offering interpretability, generalizability, and robustness. CPCG holds promise for facilitating targeted interventions in clinical treatment strategies.

Functional exploration and drug prediction on programmed cell death-related biomarkers in lung adenocarcinoma

Background

Our study aims to perform functional exploration and drug prediction of programmed cell death (PCD)-related biomarkers in lung adenocarcinoma (LUAD).

Methods

UCSC-Xena obtained LUAD-related genes. DESeq2 screened PCD-specific differentially expressed genes (DEGs), and these DEGs were intersected with genes identified by weighted gene co-expression network analysis (WGCNA) to pinpoint the key genes. KOBAS-i was used for enrichment analysis. String and GeneMania were used to construct protein interaction networks and gene-gene interaction networks, respectively. Using two machine learning algorithms to screen for key genes, and taking the intersection as biomarkers, validating via receiver operating characteristic (ROC) and in vitro experiments. Building a diagnostic model with a nomogram. Construct transcription factor (TF) regulatory network. CIBERSORT was used for immune infiltration analysis. Enrichr predicts targeted drugs and AutodockTools simulates molecular docking.

Results

120 hub genes related to PCD were identified, and an intersection of these genes with DEGs yielded 10 key genes, which were enriched in apoptosis-related pathways. Further machine learning screening of these genes led to the selection of 7 genes, among which 6 genes (FGR, LAPTM5, SIRPA, TLR4, ZEB2, and NLRC4) exhibited significant differences upon ROC validation, ultimately serving as biomarkers, in vitro experiments also confirmed. A nomogram demonstrated their excellent diagnostic performance. These six biomarkers are correlated with the infiltration status of most immune cells, suggesting that they affect LUAD through the immune system. TF regulation analysis identified the upstream miRNAs. Finally, drug prediction yielded three potential drugs: Lenvatinib, methadone, and trimethoprim.

Conclusion

PCD-related biomarkers in LUAD were explored, which may contribute to further understanding on PCD in LUAD.

The burgeoning spatial multi-omics in human gastrointestinal cancers

The development and progression of diseases in multicellular organisms unfold within the intricate three-dimensional body environment. Thus, to comprehensively understand the molecular mechanisms governing individual development and disease progression, precise acquisition of biological data, including genome, transcriptome, proteome, metabolome, and epigenome, with single-cell resolution and spatial information within the body's three-dimensional context, is essential. This foundational information serves as the basis for deciphering cellular and molecular mechanisms. Although single-cell multi-omics technology can provide biological information such as genome, transcriptome, proteome, metabolome, and epigenome with single-cell resolution, the sample preparation process leads to the loss of spatial information. Spatial multi-omics technology, however, facilitates the characterization of biological data, such as genome, transcriptome, proteome, metabolome, and epigenome in tissue samples, while retaining their spatial context. Consequently, these techniques significantly enhance our understanding of individual development and disease pathology. Currently, spatial multi-omics technology has played a vital role in elucidating various processes in tumor biology, including tumor occurrence, development, and metastasis, particularly in the realms of tumor immunity and the heterogeneity of the tumor microenvironment. Therefore, this article provides a comprehensive overview of spatial transcriptomics, spatial proteomics, and spatial metabolomics-related technologies and their application in research concerning esophageal cancer, gastric cancer, and colorectal cancer. The objective is to foster the research and implementation of spatial multi-omics technology in digestive tumor diseases. This review will provide new technical insights for molecular biology researchers.

Natural killer (NK) cells-related gene signature reveals the immune environment heterogeneity in hepatocellular carcinoma based on single cell analysis

The early diagnosis of liver cancer is crucial for the treatment and depends on the coordinated use of several test procedures. Early diagnosis is crucial for precision therapy in the treatment of the hepatocellular carcinoma (HCC). Therefore, in this study, the NK cell-related gene prediction model was used to provide the basis for precision therapy at the gene level and a novel basis for the treatment of patients with liver cancer. Natural killer (NK) cells have innate abilities to recognize and destroy tumor cells and thus play a crucial function as the "innate counterpart" of cytotoxic T cells. The natural killer (NK) cells is well recognized as a prospective approach for tumor immunotherapy in treating patients with HCC. In this research, we used publicly available databases to collect bioinformatics data of scRNA-seq and RNA-seq from HCC patients. To determine the NK cell-related genes (NKRGs)-based risk profile for HCC, we isolated T and natural killer (NK) cells and subjected them to analysis. Uniform Manifold Approximation and Projection plots were created to show the degree of expression of each marker gene and the distribution of distinct clusters. The connection between the immunotherapy response and the NKRGs-based signature was further analyzed, and the NKRGs-based signature was established. Eventually, a nomogram was developed using the model and clinical features to precisely predict the likelihood of survival. The prognosis of HCC can be accurately predicted using the NKRGs-based prognostic signature, and thorough characterization of the NKRGs signature of HCC may help to interpret the response of HCC to immunotherapy and propose a novel tumor treatment perspective.

Prediction of Prognosis and Immunotherapy Response with a Novel Natural Killer Cell Marker Genes Signature in Osteosarcoma

Background: Natural killer (NK) cells are characterized by their antitumor efficacy without previous sensitization, which have attracted attention in tumor immunotherapy. The heterogeneity of osteosarcoma (OS) has hindered therapeutic application of NK cell-based immunotherapy. The authors aimed to construct a novel NK cell-based signature to identify certain OS patients more responsive to immunotherapy. Materials and Methods: A total of eight publicly available datasets derived from patients with OS were enrolled in this study. Single-cell RNA sequencing data obtained from the Gene Expression Omnibus (GEO) database were analyzed to screen NK cell marker genes. Least Absolute Shrinkage and Selection Operator (LASSO) Cox regression analysis was used to construct an NK cell-based prognostic signature in the TARGET-OS dataset. The differences in immune cell infiltration, immune system-related metagenes, and immunotherapy response were evaluated among risk subgroups. Furthermore, this prognostic signature was experimentally validated by reverse transcription-quantitative real-time PCR (RT-qPCR). Results: With differentially expressed NK cell marker genes screened out, a five-gene NK cell-based prognostic signature was constructed. The prognostic predictive accuracy of the signature was validated through internal clinical subgroups and external GEO datasets. Low-risk OS patients contained higher abundances of infiltrated immune cells, especially CD8 T cells and naive CD4 T cells, indicating that T cell exhaustion states were present in the high-risk OS patients. As indicated from correlation analysis, immune system-related metagenes displayed a negative correlation with risk scores, suggesting the existence of immunosuppressive microenvironment in OS. In addition, based on responses to immune checkpoint inhibitor therapy in two immunotherapy datasets, the signature helped predict the response of OS patients to anti-programmed cell death protein 1 (PD-1) or anti-programmed cell death ligand 1 (PD-L1) therapy. RT-qPCR results demonstrated the roughly consistent relationship of these five gene expressions with predicting outcomes. Conclusions: The NK cell-based signature is likely to be available for the survival prediction and the evaluation of immunotherapy response of OS patients, which may shed light on subsequent immunotherapy choices for OS patients. In addition, the authors revealed a potential link between immunosuppressive microenvironment and OS.

Exploring the Molecular Mechanism of Skeletal Muscle Development in Ningxiang Pig by Weighted Gene Co-Expression Network Analysis

With the continuous improvement in living standards, people's demand for high-quality meat is increasing. Ningxiang pig has delicious meat of high nutritional value, and is loved by consumers. However, its slow growth and low meat yield seriously restrict its efficient utilization. Gene expression is the internal driving force of life activities, so in order to fundamentally improve its growth rate, it is key to explore the molecular mechanism of skeletal muscle development in Ningxiang pigs. In this paper, Ningxiang boars were selected in four growth stages (30 days: weaning period, 90 days: nursing period, 150 days: early fattening period, and 210 days: late fattening period), and the longissimus dorsi (LD) muscle was taken from three boars in each stage. The fatty acid content, amino acid content, muscle fiber diameter density and type of LD were detected by gas chromatography, acidolysis, hematoxylin eosin (HE) staining and immunofluorescence (IF) staining. After transcription sequencing, weighted gene co-expression network analysis (WGCNA) combined with the phenotype of the LD was used to explore the key genes and signaling pathways affecting muscle development. The results showed that 10 modules were identified by WGCNA, including 5 modules related to muscle development stage, module characteristics of muscle fiber density, 5 modules characteristic of muscle fiber diameter, and a module characteristic of palmitoleic acid (C16:1) and linoleic acid (C18:2n6C). Gene ontology (GO) enrichment analysis found that 52 transcripts relating to muscle development were enriched in these modules, including 44 known genes and 8 novel genes. The Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis showed that these genes were enriched in the auxin, estrogen and cyclic guanosine monophosphate-protein kinase G (cGMP-PKG) pathways. Twelve of these genes were transcription factors, there were interactions among 20 genes, and the interactions among 11 proteins in human, pig and mouse were stable. To sum up, through the integrated analysis of phenotype and transcriptome, this paper analyzed the key genes and possible regulatory networks of skeletal muscle development in Ningxiang pigs at various stages, to provide a reference for the in-depth study of skeletal muscle development.

Advances in Melanoma: From Genetic Insights to Therapeutic Innovations

Advances in melanoma research have unveiled critical insights into its genetic and molecular landscape, leading to significant therapeutic innovations. This review explores the intricate interplay between genetic alterations, such as mutations in BRAF, NRAS, and KIT, and melanoma pathogenesis. The MAPK and PI3K/Akt/mTOR signaling pathways are highlighted for their roles in tumor growth and resistance mechanisms. Additionally, this review delves into the impact of epigenetic modifications, including DNA methylation and histone changes, on melanoma progression. The tumor microenvironment, characterized by immune cells, stromal cells, and soluble factors, plays a pivotal role in modulating tumor behavior and treatment responses. Emerging technologies like single-cell sequencing, CRISPR-Cas9, and AI-driven diagnostics are transforming melanoma research, offering precise and personalized approaches to treatment. Immunotherapy, particularly immune checkpoint inhibitors and personalized mRNA vaccines, has revolutionized melanoma therapy by enhancing the body's immune response. Despite these advances, resistance mechanisms remain a challenge, underscoring the need for combined therapies and ongoing research to achieve durable therapeutic responses. This comprehensive overview aims to highlight the current state of melanoma research and the transformative impacts of these advancements on clinical practice.

Single-Cell RNA Sequencing and Its Applications in Pituitary Research

BACKGROUND

Mounting evidence underscores the significance of cellular diversity within the endocrine system and the intricate interplay between different cell types and tissues, essential for preserving physiological balance and influencing disease trajectories. The pituitary gland, a central player in the endocrine orchestra, exemplifies this complexity with its assortment of hormone-secreting and nonsecreting cells.

SUMMARY

The pituitary gland houses several types of cells responsible for hormone production, alongside nonsecretory cells like fibroblasts and endothelial cells, each playing a crucial role in the gland's function and regulatory mechanisms. Despite the acknowledged importance of these cellular interactions, the detailed mechanisms by which they contribute to pituitary gland physiology and pathology remain largely uncharted. The last decade has seen the emergence of groundbreaking technologies such as single-cell RNA sequencing, offering unprecedented insights into cellular heterogeneity and interactions. However, the application of this advanced tool in exploring the pituitary gland's complexities has been scant. This review provides an overview of this methodology, highlighting its strengths and limitations, and discusses future possibilities for employing it to deepen our understanding of the pituitary gland and its dysfunction in disease states.

KEY MESSAGE

Single-cell RNA sequencing technology offers an unprecedented means to study the heterogeneity and interactions of pituitary cells, though its application has been limited thus far. Further utilization of this tool will help uncover the complex physiological and pathological mechanisms of the pituitary, advancing research and treatment of pituitary diseases.

Development and evaluation of a disulfidoptosis-related lncRNA index for prognostication in clear cell renal cell carcinoma

Background

This study introduces a novel prognostic tool, the Disulfidoptosis-Related lncRNA Index (DRLI), integrating the molecular signatures of disulfidoptosis and long non-coding RNAs (lncRNAs) with the cellular heterogeneity of the tumor microenvironment, to predict clinical outcomes in patients with clear cell renal cell carcinoma (ccRCC).

Methods

We analyzed 530 tumor and 72 normal samples from The Cancer Genome Atlas (TCGA), employing k-means clustering based on disulfidoptosis-associated gene expression to stratify ccRCC samples into prognostic groups. lncRNAs correlated with disulfidoptosis were identified and used to construct the DRLI, which was validated by Kaplan-Meier and receiver operating characteristic curves. We utilized single-cell deconvolution analysis to estimate the proportion of immune cell types within the tumor microenvironment, while the ESTIMATE and TIDE algorithms were employed to assess immune infiltration and potential response to immunotherapy.

Results

The Disulfidoptosis-Related lncRNA Index (DRLI) effectively stratified ccRCC patients into high and low-risk groups, significantly impacting survival outcomes (P < 0.001). High-risk patients, marked by a unique lncRNA profile associated with disulfidoptosis, faced worse prognoses. Single-cell analysis revealed marked tumor microenvironment heterogeneity, especially in immune cell makeup, correlating with patient risk levels. In prognostic predictions, DRLI outperformed traditional clinical indicators, achieving AUC values of 0.779, 0.757, and 0.779 for 1-year, 3-year, and 5-year survival in the training set, and 0.746, 0.734, and 0.750 in the validation set. Notably, while the constructed nomogram showed exceptional predictive capability for short-term prognosis (AUC = 0.877), the DRLI displayed remarkable long-term predictive accuracy, with its AUC value reaching 0.823 for 10-year survival, closely approaching the nomogram's performance.

Conclusions

The study introduces the DRLI as a groundbreaking molecular stratification tool for ccRCC, enhancing prognostic precision and potentially guiding personalized treatment strategies. This advancement is particularly significant in the context of long-term survival predictions. Our findings also elucidate the complex interplay between disulfidoptosis, lncRNAs, and the immune microenvironment in ccRCC, offering a comprehensive perspective on its pathogenesis and progression. The DRLI and the nomogram together represent significant strides in ccRCC research, highlighting the importance of molecular-based assessments in predicting patient outcomes.

TimiGP-Response: the pan-cancer immune landscape associated with response to immunotherapy

Accumulating evidence suggests that the tumor immune microenvironment (TIME) significantly influences the response to immunotherapy, yet this complex relationship remains elusive. To address this issue, we developed TimiGP-Response (TIME Illustration based on Gene Pairing designed for immunotherapy Response), a computational framework leveraging single-cell and bulk transcriptomic data, along with response information, to construct cell-cell interaction networks associated with responders and estimate the role of immune cells in treatment response. This framework was showcased in triple-negative breast cancer treated with immune checkpoint inhibitors targeting the PD-1:PD-L1 interaction, and orthogonally validated with imaging mass cytometry. As a result, we identified CD8+ GZMB+ T cells associated with responders and its interaction with regulatory T cells emerged as a potential feature for selecting patients who may benefit from these therapies. Subsequently, we analyzed 3,410 patients with seven cancer types (melanoma, non-small cell lung cancer, renal cell carcinoma, metastatic urothelial carcinoma, hepatocellular carcinoma, breast cancer, and esophageal cancer) treated with various immunotherapies and combination therapies, as well as several chemo- and targeted therapies as controls. Using TimiGP-Response, we depicted the pan-cancer immune landscape associated with immunotherapy response at different resolutions. At the TIME level, CD8 T cells and CD4 memory T cells were associated with responders, while anti-inflammatory (M2) macrophages and mast cells were linked to non-responders across most cancer types and datasets. Given that T cells are the primary targets of these immunotherapies and our TIME analysis highlights their importance in response to treatment, we portrayed the pan-caner landscape on 40 T cell subtypes. Notably, CD8+ and CD4+ GZMK+ effector memory T cells emerged as crucial across all cancer types and treatments, while IL-17-producing CD8+ T cells were top candidates associated with immunotherapy non-responders. In summary, this study provides a computational method to study the association between TIME and response across the pan-cancer immune landscape, offering resources and insights into immune cell interactions and their impact on treatment efficacy.

Progress of single-cell RNA sequencing combined with spatial transcriptomics in tumour microenvironment and treatment of pancreatic cancer

In recent years, single-cell analyses have revealed the heterogeneity of the tumour microenvironment (TME) at the genomic, transcriptomic, and proteomic levels, further improving our understanding of the mechanisms of tumour development. Single-cell RNA sequencing (scRNA-seq) technology allow analysis of the transcriptome at the single-cell level and have unprecedented potential for exploration of the characteristics involved in tumour development and progression. These techniques allow analysis of transcript sequences at higher resolution, thereby increasing our understanding of the diversity of cells found in the tumour microenvironment and how these cells interact in complex tumour tissue. Although scRNA-seq has emerged as an important tool for studying the tumour microenvironment in recent years, it cannot be used to analyse spatial information for cells. In this regard, spatial transcriptomics (ST) approaches allow researchers to understand the functions of individual cells in complex multicellular organisms by understanding their physical location in tissue sections. In particular, in related research on tumour heterogeneity, ST is an excellent complementary approach to scRNA-seq, constituting a new method for further exploration of tumour heterogeneity, and this approach can also provide unprecedented insight into the development of treatments for pancreatic cancer (PC). In this review, based on the methods of scRNA-seq and ST analyses, research progress on the tumour microenvironment and treatment of pancreatic cancer is further explained.

Single-Cell Sequencing: High-Resolution Analysis of Cellular Heterogeneity in Autoimmune Diseases

Autoimmune diseases (AIDs) are complex in etiology and diverse in classification but clinically show similar symptoms such as joint pain and skin problems. As a result, the diagnosis is challenging, and usually, only broad treatments can be available. Consequently, the clinical responses in patients with different types of AIDs are unsatisfactory. Therefore, it is necessary to conduct more research to figure out the pathogenesis and therapeutic targets of AIDs. This requires research technologies with strong extraction and prediction capabilities. Single-cell sequencing technology analyses the genomic, epigenomic, or transcriptomic information at the single-cell level. It can define different cell types and states in greater detail, further revealing the molecular mechanisms that drive disease progression. These advantages enable cell biology research to achieve an unprecedented resolution and scale, bringing a whole new vision to life science research. In recent years, single-cell technology especially single-cell RNA sequencing (scRNA-seq) has been widely used in various disease research. In this paper, we present the innovations and applications of single-cell sequencing in the medical field and focus on the application contributing to the differential diagnosis and precise treatment of AIDs. Despite some limitations, single-cell sequencing has a wide range of applications in AIDs. We finally present a prospect for the development of single-cell sequencing. These ideas may provide some inspiration for subsequent research.

Pancreatic cancer environment: from patient-derived models to single-cell omics

Pancreatic cancer (PC) is a highly malignant cancer characterized by poor prognosis, high heterogeneity, and intricate heterocellular systems. Selecting an appropriate experimental model for studying its progression and treatment is crucial. Patient-derived models provide a more accurate representation of tumor heterogeneity and complexity compared to cell line-derived models. This review initially presents relevant patient-derived models, including patient-derived xenografts (PDXs), patient-derived organoids (PDOs), and patient-derived explants (PDEs), which are essential for studying cell communication and pancreatic cancer progression. We have emphasized the utilization of these models in comprehending intricate intercellular communication, drug responsiveness, mechanisms underlying tumor growth, expediting drug discovery, and enabling personalized medical approaches. Additionally, we have comprehensively summarized single-cell analyses of these models to enhance comprehension of intercellular communication among tumor cells, drug response mechanisms, and individual patient sensitivities.

Hypoxia-related THBD+ macrophages as a prognostic factor in glioma: Construction of a powerful risk model

Glioma is a prevalent malignant tumour characterized by hypoxia as a pivotal factor in its progression. This study aims to investigate the impact of the most severely hypoxic cell subpopulation in glioma. Our findings reveal that the THBD+ macrophage subpopulation is closely associated with hypoxia in glioma, exhibiting significantly higher infiltration in tumours compared to non-tumour tissues. Moreover, a high proportion of THBD+ cells correlates with poor prognosis in glioblastoma (GBM) patients. Notably, THBD+ macrophages exhibit hypoxic characteristics and epithelial-mesenchymal transition features. Silencing THBD expression leads to a notable reduction in the proliferation and metastasis of glioma cells. Furthermore, we developed a THBD+ macrophage-related risk signature (THBDMRS) through machine learning techniques. THBDMRS emerges as an independent prognostic factor for GBM patients with a substantial prognostic impact. By comparing THBDMRS with 119 established prognostic features, we demonstrate the superior prognostic performance of THBDMRS. Additionally, THBDMRS is associated with glioma metastasis and extracellular matrix remodelling. In conclusion, hypoxia-related THBD+ macrophages play a pivotal role in glioma pathogenesis, and THBDMRS emerges as a potent and promising prognostic tool for GBM, contributing to enhanced patient survival outcomes.

Analysis of cancer-associated fibroblasts related genes identifies COL11A1 associated with lung adenocarcinoma prognosis

BACKGROUND

The treatment of lung adenocarcinoma is difficult due to the limited therapeutic options. Cancer-associated fibroblasts play an important role in the development of cancers. This study aimed to identify a promising molecular target associated with cancer-associated fibroblasts for the treatment of lung adenocarcinoma.

METHODS

The Cancer Genome Atlas lung adenocarcinoma dataset was used to screen hub genes associated with cancer-associated fibroblasts via the EPIC algorithm and Weighted Gene Co-expression Network Analysis. Multiple databases were used together with our data to verify the differential expression and survival of COL11A1. Functional enrichment analysis and the single-cell TISCH database were used to elucidate the mechanisms underlying COL11A1 expression. The correlation between COL11A1 and immune checkpoint genes in human cancers was also evaluated.

RESULTS

Using the EPIC algorithm and Weighted Gene Co-expression Network Analysis, 13 hub genes associated with cancer-associated fibroblasts in lung adenocarcinoma were screened. Using the GEPIA database, Kaplan-Meier Plotter database, GSE72094, GSE75037, GSE32863, and our immunohistochemistry experiment data, we confirmed that COL11A1 overexpresses in lung adenocarcinoma and that high expression of COL11A1 is associated with a poor prognosis. COL11A1 has a genetic alteration frequency of 22% in patients with lung adenocarcinoma. COL11A1 is involved in the extracellular matrix activities of lung adenocarcinoma. Using the TISCH database, we found that COL11A1 is mainly expressed by cancer-associated fibroblasts in the tumor microenvironment rather than by lung adenocarcinoma cells. Finally, we found that COL11A1 is positively correlated with HAVCR2(TIM3), CD274 (PD-L1), CTLA4, and LAG3 in lung adenocarcinoma.

CONCLUSION

COL11A1 may be expressed and secreted by cancer-associated fibroblasts, and a high expression of COL11A1 may result in T cell exhaustion in the tumor microenvironment of lung adenocarcinoma. COL11A1 may serve as an attractive biomarker to provide new insights into cancer therapeutics.

Drug repositioning for immunotherapy in breast cancer using single-cell analysis

Immunomodulatory peptides, while exhibiting potential antimicrobial, antifungal, and/or antiviral properties, can play a role in stimulating or suppressing the immune system, especially in pathological conditions like breast cancer (BC). Thus, deregulation of these peptides may serve as an immunotherapeutic strategy to enhance the immune response. In this meta-analysis, we utilized single-cell RNA sequencing data and known therapeutic peptides to investigate the deregulation of these peptides in malignant versus normal human breast epithelial cells. We corroborated our findings at the chromatin level using ATAC-seq. Additionally, we assessed the protein levels in various BC cell lines. Moreover, our in-house drug repositioning approach was employed to identify potential drugs that could positively impact the relapse-free survival of BC patients. Considering significantly deregulated therapeutic peptides and their role in BC pathology, our approach aims to downregulate B2M and SLPI, while upregulating PIGR, DEFB1, LTF, CLU, S100A7, and SCGB2A1 in BC epithelial cells through our drug repositioning pipeline. Leveraging the LINCS L1000 database, we propose BRD-A06641369 for B2M downregulation and ST-4070043 and BRD-K97926541 for SLPI downregulation without negatively affecting the MHC complex as a significantly correlated pathway with these two genes. Furthermore, we have compiled a comprehensive list of drugs for the upregulation of other selected immunomodulatory peptides. Employing an immunotherapeutic approach by integrating our drug repositioning pipeline with single-cell analysis, we proposed potential drugs and drug targets to fortify the immune system against BC.

Molecular characterization of chronic cutaneous wounds reveals subregion- and wound type-specific differential gene expression

A limited understanding of the pathology underlying chronic wounds has hindered the development of effective diagnostic markers and pharmaceutical interventions. This study aimed to elucidate the molecular composition of various common chronic ulcer types to facilitate drug discovery strategies. We conducted a comprehensive analysis of leg ulcers (LUs), encompassing venous and arterial ulcers, foot ulcers (FUs), pressure ulcers (PUs), and compared them with surgical wound healing complications (WHCs). To explore the pathophysiological mechanisms and identify similarities or differences within wounds, we dissected wounds into distinct subregions, including the wound bed, border, and peri-wound areas, and compared them against intact skin. By correlating histopathology, RNA sequencing (RNA-Seq), and immunohistochemistry (IHC), we identified unique genes, pathways, and cell type abundance patterns in each wound type and subregion. These correlations aim to aid clinicians in selecting targeted treatment options and informing the design of future preclinical and clinical studies in wound healing. Notably, specific genes, such as PITX1 and UPP1, exhibited exclusive upregulation in LUs and FUs, potentially offering significant benefits to specialists in limb preservation and clinical treatment decisions. In contrast, comparisons between different wound subregions, regardless of wound type, revealed distinct expression profiles. The pleiotropic chemokine-like ligand GPR15L (C10orf99) and transmembrane serine proteases TMPRSS11A/D were significantly upregulated in wound border subregions. Interestingly, WHCs exhibited a nearly identical transcriptome to PUs, indicating clinical relevance. Histological examination revealed blood vessel occlusions with impaired angiogenesis in chronic wounds, alongside elevated expression of genes and immunoreactive markers related to blood vessel and lymphatic epithelial cells in wound bed subregions. Additionally, inflammatory and epithelial markers indicated heightened inflammatory responses in wound bed and border subregions and reduced wound bed epithelialization. In summary, chronic wounds from diverse anatomical sites share common aspects of wound pathophysiology but also exhibit distinct molecular differences. These unique molecular characteristics present promising opportunities for drug discovery and treatment, particularly for patients suffering from chronic wounds. The identified diagnostic markers hold the potential to enhance preclinical and clinical trials in the field of wound healing.

Prognostic Risk Models Using Epithelial Cells Identify β-Sitosterol as a Potential Therapeutic Target Against Esophageal Squamous Cell Carcinoma

Background

Esophageal squamous cell carcinoma (ESCC) is an aggressive and fatal malignancy that leads to epithelial cancer. The association between epithelial cell heterogeneity, prognosis, and immune response in this cancer remains uncertain. This study aimed to investigate epithelial cell heterogeneity in ESCC and develop a predictive risk model using the identified cell types.

Methods

Single-cell RNA sequencing (scRNA-seq) and differential ESCC gene data were accessed from the Gene Expression Omnibus. Functional enrichment analysis, inferCNV, cell development trajectories, and intercellular communication were analyzed following epithelial cell characterization. Differentially expressed ESCC (n = 773) and epithelial cell marker genes (n = 3407) were intersected to obtain core genes, and epithelial cell-related prognostic genes were identified. LASSO regression analysis was used to construct a prognostic model. The external dataset GSE53624 was used to further validate the stability of the model. Drug sensitivity predictions, and immune cell infiltration were analyzed. Molecular docking clarified the possible therapeutic role of β-sitosterol in ESCC. Finally, wound healing assay, cell colony, and transwell assay were constructed to detect the effects of the core gene PDLIM2 on ESCC cell proliferation, invasion, and migration.

Results

Eight cell clusters were identified, and epithelial cells were categorized into tumor and paratumor groups. The tumor group possessed more chromosomal variants than the paratumor group. Epithelial cells were associated with multiple cell types and significantly correlated with the Wnt, transforming growth factor, and epidermal growth factor signaling pathways. From 231 intersected genes, five core genes were screened for use in the risk model: CTSL, LAPTM4B, MYO10, NCF2, and PDLIM2. These genes may contribute to the cancerous transformation of normal esophageal epithelial cells and thereby act as biomarkers and potential therapeutic targets in patients with ESCC. β-Sitosterol furthermore displayed excellent docking potential with these genes. Meanwhile, further experiments demonstrated that the gene PDLIM2 plays a major role in the progression of oesophageal squamous carcinoma.

Conclusion

We successfully developed a risk model for the prognosis of ESCC based on epithelial cells that addresses the response of ESCC to immunotherapy and offers novel cancer treatment options.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we propose that RNA-seq should be considered a diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 196 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that-combined with sequence alignments and BLASTp-they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Transcriptomic signatures of chronic active antibody-mediated rejection deciphered by RNA sequencing of human kidney allografts

Natural killer (NK) cells mediate spontaneous cell-mediated cytotoxicity and antibody-dependent cell-mediated cytotoxicity. This dual functionality could enable their participation in chronic active antibody-mediated rejection (CA-ABMR). Earlier microarray profiling studies have not subcategorized antibody-mediated rejection into CA-ABMR and active-ABMR, and the gene expression pattern of CA-ABMR has not been compared with that of T cell-mediated rejection (TCMR). To fill these gaps, we RNA sequenced human kidney allograft biopsies categorized as CA-ABMR, active-ABMR, TCMR, or No Rejection (NR). Among the 15,910 genes identified in the biopsies, 60, 114, and 231 genes were uniquely overexpressed in CA-ABMR, TCMR, and active-ABMR, respectively; compared to NR, 50 genes were shared between CA-ABMR and active-ABMR, and 164 genes between CA-ABMR and TCMR. The overexpressed genes were annotated to NK cells and T cells in CA-ABMR and TCMR, and to neutrophils and monocytes in active-ABMR. The NK cell cytotoxicity and allograft rejection pathways were enriched in CA-ABMR. Genes encoding perforin, granzymes, and death receptor were overexpressed in CA-ABMR versus active-ABMR but not compared to TCMR. NK cell cytotoxicity pathway gene set variation analysis score was higher in CA-ABMR compared to active-ABMR but not in TCMR. Principal component analysis of the deconvolved immune cellular transcriptomes separated CA-ABMR and TCMR from active-ABMR and NR. Immunohistochemistry of kidney allograft biopsies validated a higher proportion of CD56+ NK cells in CA-ABMR than in active-ABMR. Thus, CA-ABMR was exemplified by the overexpression of the NK cell cytotoxicity pathway gene set and, surprisingly, molecularly more like TCMR than active-ABMR.

Highly Accurate Estimation of Cell Type Abundance in Bulk Tissues Based on Single-Cell Reference and Domain Adaptive Matching

Accurately identifies the cellular composition of complex tissues, which is critical for understanding disease pathogenesis, early diagnosis, and prevention. However, current methods for deconvoluting bulk RNA sequencing (RNA-seq) typically rely on matched single-cell RNA sequencing (scRNA-seq) as a reference, which can be limiting due to differences in sequencing distribution and the potential for invalid information from single-cell references. Hence, a novel computational method named SCROAM is introduced to address these challenges. SCROAM transforms scRNA-seq and bulk RNA-seq into a shared feature space, effectively eliminating distributional differences in the latent space. Subsequently, cell-type-specific expression matrices are generated from the scRNA-seq data, facilitating the precise identification of cell types within bulk tissues. The performance of SCROAM is assessed through benchmarking against simulated and real datasets, demonstrating its accuracy and robustness. To further validate SCROAM's performance, single-cell and bulk RNA-seq experiments are conducted on mouse spinal cord tissue, with SCROAM applied to identify cell types in bulk tissue. Results indicate that SCROAM is a highly effective tool for identifying similar cell types. An integrated analysis of liver cancer and primary glioblastoma is then performed. Overall, this research offers a novel perspective for delivering precise insights into disease pathogenesis and potential therapeutic strategies.

Single-cell RNA sequencing reveals immune microenvironment of small cell lung cancer-associated malignant pleural effusion

We used 10 × genomics single-cell transcriptome sequencing technology to reveal the tumor immune microenvironment characteristics of small cell lung cancer (SCLC) in a patient with malignant pleural effusion (MPE). A total of 8008 high-quality cells were finally obtained for subsequent bioinformatic analysis, which were divided into 10 cell clusters further identified as B cells, T cells, myeloid cells, NK cells, and cancer cells. Such SCLC related genes as NOTCH1, MYC, TSC22D1, SOX4, BLNK, YBX3, VIM, CD8A, CD8B, and KLF6 were expressed in different degrees during differentiation of T and B cells. Different ligands and receptors between T, B and tumor cells almost interact through MHC II, IL-16, galectin, and APP signaling pathway.

Identification of hub genes based on integrated analysis of single-cell and microarray transcriptome in patients with pulmonary arterial hypertension

BACKGROUND

Pulmonary arterial hypertension (PAH) is a devastating chronic cardiopulmonary disease without an effective therapeutic approach. The underlying molecular mechanism of PAH remains largely unexplored at single-cell resolution.

METHODS

Single-cell RNA sequencing (scRNA-seq) data from the Gene Expression Omnibus (GEO) database (GSE210248) was included and analyzed comprehensively. Additionally, microarray transcriptome data including 15 lung tissue from PAH patients and 11 normal samples (GSE113439) was also obtained. Seurat R package was applied to process scRNA-seq data. Uniform manifold approximation and projection (UMAP) was utilized for dimensionality reduction and cluster identification, and the SingleR package was performed for cell annotation. FindAllMarkers analysis and ClusterProfiler package were applied to identify differentially expressed genes (DEGs) for each cluster in GSE210248 and GSE113439, respectively. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genome (KEGG) were used for functional enrichment analysis of DEGs. Microenvironment Cell Populations counter (MCP counter) was applied to evaluate the immune cell infiltration. STRING was used to construct a protein-protein interaction (PPI) network of DEGs, followed by hub genes selection through Cytoscape software and Veen Diagram.

RESULTS

Nineteen thousand five hundred seventy-six cells from 3 donors and 21,896 cells from 3 PAH patients remained for subsequent analysis after filtration. A total of 42 cell clusters were identified through UMAP and annotated by the SingleR package. 10 cell clusters with the top 10 cell amounts were selected for consequent analysis. Compared with the control group, the proportion of adipocytes and fibroblasts was significantly reduced, while CD8+ T cells and macrophages were notably increased in the PAH group. MCP counter revealed decreased distribution of CD8+ T cells, cytotoxic lymphocytes, and NK cells, as well as increased infiltration of monocytic lineage in PAH lung samples. Among 997 DEGs in GSE113439, module 1 with 68 critical genes was screened out through the MCODE plug-in in Cytoscape software. The top 20 DEGs in each cluster of GSE210248 were filtered out by the Cytohubba plug-in using the MCC method. Eventually, WDR43 and GNL2 were found significantly increased in PAH and identified as the hub genes after overlapping these DEGs from GSE210248 and GSE113439.

CONCLUSION

WDR43 and GNL2 might provide novel insight into revealing the new molecular mechanisms and potential therapeutic targets for PAH.

Molecular subtyping based on immune cell marker genes predicts prognosis and therapeutic response in patients with lung adenocarcinoma

OBJECTIVE

Lung adenocarcinoma (LA) is one of the most common malignancies and is responsible for the greatest number of tumor-related deaths. Our research aimed to explore the molecular subtype signatures of LA to clarify the correlation among the immune microenvironment, clinical outcomes, and therapeutic response.

METHODS

The LA immune cell marker genes (LICMGs) identified by single-cell RNA sequencing (scRNA-seq) analysis were used to discriminate the molecular subtypes and homologous immune and metabolic traits of GSE72094 LA cases. In addition, the model-building genes were identified from 1441 LICMGs by Cox-regression analysis, and a LA immune difference score (LIDscore) was developed to quantify individual differences in each patient, thereby predicting prognosis and susceptibility to immunotherapy and chemotherapy of LA patients.

RESULTS

Patients of the GSE72094 cohort were divided into two distinct molecular subtypes based on LICMGs: immune activating subtype (Cluster-C1) and metabolically activating subtype (cluster-C2). The two molecular subtypes have distinct characteristics regarding prognosis, clinicopathology, genomics, immune microenvironment, and response to immunotherapy. Among the LICMGs, LGR4, GOLM1, CYP24A1, SFTPB, COL1A1, HLA-DQA1, MS4A7, PPARG, and IL7R were enrolled to construct a LIDscore model. Low-LIDscore patients had a higher survival rate due to abundant immune cell infiltration, activated immunity, and lower genetic variation, but probably the higher levels of Treg cells in the immune microenvironment lead to immune cell dysfunction and promote tumor immune escape, thus decreasing the responsiveness to immunotherapy compared with that of the high-LIDscore patients. Overall, high-LIDscore patients had a higher responsiveness to immunotherapy and a higher sensitivity to chemotherapy than the low-LIDscore group.

CONCLUSIONS

Molecular subtypes based on LICMGs provided a promising strategy for predicting patient prognosis, biological characteristics, and immune microenvironment features. In addition, they helped identify the patients most likely to benefit from immunotherapy and chemotherapy.

Single cell RNA sequencing – a valuable tool for cancer immunotherapy: a mini review

Single-cell RNA sequencing (scRNA-seq) technology has made great strides in research over the last decade. Data analysis has been aided by developments in bioinformatics tools and artificial intelligence, allowing biological and clinical researchers to get a deeper understanding of the different cell clusters and their dynamics within tumours. Combining conventional treatment modalities like chemotherapy and radiation with immunotherapy is a growing trend in cancer treatment. Hence, knowledge of the tumour microenvironment and the effect of each treatment modality on the TME, at a single cell level can provide treating clinicians with better clues for patient stratification and prognostication. With this knowledge, immunotherapy could become successful in treating a wide range of cancers, opening the path for the creation of even more effective treatment strategies. Despite the widespread availability of scRNA-seq technology, computational analysis and data interpretation are still challenges. Worldwide, such challenges are being addressed by various researchers, strengthening the contribution of this technology towards cancer elimination. In this mini-review, we primarily focus on the technique, its workflow, and the computational aspects of scRNA technology, along with an overview of the current challenges in the analysis and interpretation of the data generated.

Identification of cancer-associated fibroblasts signature for predicting the prognosis and immunotherapy response in hepatocellular carcinoma

Hepatocellular carcinoma (HCC) is one of the most common malignancies globally with poor prognosis. Cancer-associated fibroblasts (CAFs) play multiple functions in the regulation of tumorigenesis, metastasis and therapeutic resistance of cancer. The current study aimed to explore the role of CAFs-related genes in the prognosis and immunotherapy response in HCC. CAFs-related genes were identified by using single-cell RNA-sequencing analysis. Least absolute shrinkage and selection operator (LASSO) analysis was conducted to develop a CAFs-related prognostic signature (FRPS) in TCGA dataset and verified in ICGC, GSE14520 and GSE76427 cohorts. Several tools, including Tumor Immune Dysfunction and Exclusion (TIDE) score, immunophenoscore, and Tumor Mutation Burden (TMB) score were used to evaluate the value of FRPS in predicting immunotherapy benefits. The FRPS constructed based on 10 genes (RGS5, CNN3, PALLD, FLNA, KLHL23, MYC, NDRG2, SERPINE1, CD151 CALU) served as an independent risk factor and showed stable and powerful performance in predicting the overall survival rate of HCC patients with an AUCs of 0. 734, 0.727, and 0.717 in 2-, 3-, and 4-year ROC curve in TCGA cohort. Low risk score indicated a higher abundance of CD8+ T cells and NK, and lower abundance of Treg. Moreover, HCC patients with low risk score had a higher PD1&CTLA4 immunophenoscore, higher TMB score, and lower TIDE score. Moreover, high risk score indicated a lower IC50 value of 5-fluorouracil, camptothecin, cisplatin, docetaxel, gemcitabine, paclitaxel, afatinib, crizotinib, dasatinib, erlotinib, erlotinib, gefitinib, lapatinib, and osimertinib in HCC. Our study develops a novel FRPS HCC. The FRPS acts as a risk factor for the prognosis of HCC patients and it can predict the immunotherapy benefits of HCC patients.

Blood transcriptomics analysis offers insights into variant-specific immune response to SARS-CoV-2

Bulk RNA sequencing (RNA-seq) of blood is typically used for gene expression analysis in biomedical research but is still rarely used in clinical practice. In this study, we argue that RNA-seq should be considered a routine diagnostic tool, as it offers not only insights into aberrant gene expression and splicing but also delivers additional readouts on immune cell type composition as well as B-cell and T-cell receptor (BCR/TCR) repertoires. We demonstrate that RNA-seq offers vital insights into a patient's immune status via integrative analysis of RNA-seq data from patients infected with various SARS-CoV-2 variants (in total 240 samples with up to 200 million reads sequencing depth). We compare the results of computational cell-type deconvolution methods (e.g., MCP-counter, xCell, EPIC, quanTIseq) to complete blood count data, the current gold standard in clinical practice. We observe varying levels of lymphocyte depletion and significant differences in neutrophil levels between SARS-CoV-2 variants. Additionally, we identify B and T cell receptor (BCR/TCR) sequences using the tools MiXCR and TRUST4 to show that - combined with sequence alignments and pBLAST - they could be used to classify a patient's disease. Finally, we investigated the sequencing depth required for such analyses and concluded that 10 million reads per sample is sufficient. In conclusion, our study reveals that computational cell-type deconvolution and BCR/TCR methods using bulk RNA-seq analyses can supplement missing CBC data and offer insights into immune responses, disease severity, and pathogen-specific immunity, all achievable with a sequencing depth of 10 million reads per sample.

Deciphering the molecular and cellular atlas of immune cells in septic patients with different bacterial infections

BACKGROUND

Sepsis is a life-threatening organ dysfunction caused by abnormal immune responses to various, predominantly bacterial, infections. Different bacterial infections lead to substantial variation in disease manifestation and therapeutic strategies. However, the underlying cellular heterogeneity and mechanisms involved remain poorly understood.

METHODS

Multiple bulk transcriptome datasets from septic patients with 12 types of bacterial infections were integrated to identify signature genes for each infection. Signature genes were mapped onto an integrated large single-cell RNA (scRNA) dataset from septic patients, to identify subsets of cells associated with different sepsis types, and multiple omics datasets were combined to reveal the underlying molecular mechanisms. In addition, an scRNA dataset and spatial transcriptome data were used to identify signaling pathways in sepsis-related cells. Finally, molecular screening, optimization, and de novo design were conducted to identify potential targeted drugs and compounds.

RESULTS

We elucidated the cellular heterogeneity among septic patients with different bacterial infections. In Escherichia coli (E. coli) sepsis, 19 signature genes involved in epigenetic regulation and metabolism were identified, of which DRAM1 was demonstrated to promote autophagy and glycolysis in response to E. coli infection. DRAM1 upregulation was confirmed in an independent sepsis cohort. Further, we showed that DRAM1 could maintain survival of a pro-inflammatory monocyte subset, C10_ULK1, which induces systemic inflammation by interacting with other cell subsets via resistin and integrin signaling pathways in blood and kidney tissue, respectively. Finally, retapamulin was identified and optimized as a potential drug for treatment of E. coli sepsis targeting the signature gene, DRAM1, and inhibiting E. coli protein synthesis. Several other targeted drugs were also identified in other types of sepsis, including nystatin targeting C1QA in Neisseria sepsis and dalfopristin targeting CTSD in Streptococcus viridans sepsis.

CONCLUSION

Our study provides a comprehensive overview of the cellular heterogeneity and underlying mechanisms in septic patients with various bacterial infections, providing insights to inform development of stratified targeted therapies for sepsis.

Single-cell analysis reveals the immune heterogeneity and interactions in lungs undergoing hepatic ischemia-reperfusion

Hepatic ischemia-reperfusion IR (HIR) is an unavoidable pathophysiological process during liver transplantation, resulting in systematic sterile inflammation and remote organ injury. Acute lung injury (ALI) is a serious complication after liver transplantation with high postoperative morbidity and mortality. However, the underlying mechanism is still unclear. To assess the phenotype and plasticity of various cell types in the lung tissue microenvironment after HIR at the single-cell level, single-cell RNA sequencing (scRNA-seq) was performed using the lungs from HIR-induced mice. In our results, we identified 23 cell types in the lungs after HIR and found that this highly complex ecosystem was formed by subpopulations of bone marrow-derived cells that signaled each other and mediated inflammatory responses in different states and different intervals. We described the unique transcriptional profiles of lung cell clusters and discovered two novel cell subtypes (Tspo+Endothelial cells and Vcan+ monocytes), as well as the endothelial cell-immune cell and immune cell-T cell clusters interactome. In addition, we found that S100 calcium binding protein (S100a8/a9), specifically and highly expressed in immune cell clusters of lung tissues and exhibited detrimental effects. Finally, the cellular landscape of the lung tissues after HIR was established, highlighting the heterogeneity and cellular interactions between major immune cells in HIR-induced lungs. Our findings provided new insights into the mechanisms of HIR-induced ALI and offered potential therapeutic target to prevent ALI after liver transplantation.

Effect of age on the risk of immune-related adverse events in patients receiving immune checkpoint inhibitors

Identifying patients at increased risk of immune-related adverse events (irAEs) facilitates safe application of immune checkpoint inhibitors (ICIs). This retrospective study aimed to determine the effect of age on the risk of irAEs in patients receiving ICIs and to identify potential mechanisms underlying age-related irAE risk differences. We analyzed reports of FDA Adverse Event Reporting System from July 1, 2014, to September 30, 2021. The information component ratio (ICΔ) was used to compare the irAE risk between older adults (> 65 years) and younger adults (25-65 years), of which the 95% confidential interval lower limit (ICΔ025) exceeding zero indicated significantly increased risk. We found that older adults had a significantly higher overall irAE risk than younger adults (ICΔ025 0.38), which was observed in almost all organ systems. We further analyzed the correlation between age-related irAE risks and age-related transcriptional changes to identify potential genes and pathways underlying age-related irAE risk differences. We found that genes significantly correlated with ICΔ were enriched in processes including extracellular matrix organization, regulation of myeloid leukocyte mediated immunity, and regulation of c-Jun N-terminal kinase (JNK) cascade. In addition, single-cell RNA sequencing analysis confirmed that genes involved in collagen-containing extracellular matrix and JNK cascade were significantly upregulated in myeloid cells from ICI-associated colitis tissues compared with ICI-treated colon tissues without colitis. In conclusion, older adults receiving ICIs have higher irAE risks than younger adults. Upregulation of genes involved in JNK cascade and collagen-containing extracellular matrix in myeloid cells may contribute to increased irAE risks in older adults.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq reveals RNA N6-methyladenosine modification associated with prognosis and drug resistance in acute myeloid leukemia

Introduction

Acute myeloid leukemia (AML) is a type of blood cancer that is identified by the unrestricted growth of immature myeloid cells within the bone marrow. Despite therapeutic advances, AML prognosis remains highly variable, and there is a lack of biomarkers for customizing treatment. RNA N6-methyladenosine (m6A) modification is a reversible and dynamic process that plays a critical role in cancer progression and drug resistance.

Methods

To investigate the m6A modification patterns in AML and their potential clinical significance, we used the AUCell method to describe the m6A modification activity of cells in AML patients based on 23 m6A modification enzymes and further integrated with bulk RNA-seq data.

Results

We found that m6A modification was more effective in leukemic cells than in immune cells and induced significant changes in gene expression in leukemic cells rather than immune cells. Furthermore, network analysis revealed a correlation between transcription factor activation and the m6A modification status in leukemia cells, while active m6A-modified immune cells exhibited a higher interaction density in their gene regulatory networks. Hierarchical clustering based on m6A-related genes identified three distinct AML subtypes. The immune dysregulation subtype, characterized by RUNX1 mutation and KMT2A copy number variation, was associated with a worse prognosis and exhibited a specific gene expression pattern with high expression level of IGF2BP3 and FMR1, and low expression level of ELAVL1 and YTHDF2. Notably, patients with the immune dysregulation subtype were sensitive to immunotherapy and chemotherapy.

Discussion

Collectively, our findings suggest that m6A modification could be a potential therapeutic target for AML, and the identified subtypes could guide personalized therapy.

Application of single-cell sequencing to the research of tumor microenvironment

Single-cell sequencing is a technique for detecting and analyzing genomes, transcriptomes, and epigenomes at the single-cell level, which can detect cellular heterogeneity lost in conventional sequencing hybrid samples, and it has revolutionized our understanding of the genetic heterogeneity and complexity of tumor progression. Moreover, the tumor microenvironment (TME) plays a crucial role in the formation, development and response to treatment of tumors. The application of single-cell sequencing has ushered in a new age for the TME analysis, revealing not only the blueprint of the pan-cancer immune microenvironment, but also the heterogeneity and differentiation routes of immune cells, as well as predicting tumor prognosis. Thus, the combination of single-cell sequencing and the TME analysis provides a unique opportunity to unravel the molecular mechanisms underlying tumor development and progression. In this review, we summarize the recent advances in single-cell sequencing and the TME analysis, highlighting their potential applications in cancer research and clinical translation.