Identifying phenotype-associated subpopulations by integrating bulk and single-cell sequencing data

Heterogeneity of mature oligodendrocytes in the central nervous system

Mature oligodendrocytes form myelin sheaths that are crucial for the insulation of axons and efficient signal transmission in the central nervous system. Recent evidence has challenged the classical view of the functionally static mature oligodendrocyte and revealed a gamut of dynamic functions such as the ability to modulate neuronal circuitry and provide metabolic support to axons. Despite the recognition of potential heterogeneity in mature oligodendrocyte function, a comprehensive summary of mature oligodendrocyte diversity is lacking. We delve into early 20 th -century studies by Robertson and Río-Hortega that laid the foundation for the modern identification of regional and morphological heterogeneity in mature oligodendrocytes. Indeed, recent morphologic and functional studies call into question the long-assumed homogeneity of mature oligodendrocyte function through the identification of distinct subtypes with varying myelination preferences. Furthermore, modern molecular investigations, employing techniques such as single cell/nucleus RNA sequencing, consistently unveil at least six mature oligodendrocyte subpopulations in the human central nervous system that are highly transcriptomically diverse and vary with central nervous system region. Age and disease related mature oligodendrocyte variation denotes the impact of pathological conditions such as multiple sclerosis, Alzheimer's disease, and psychiatric disorders. Nevertheless, caution is warranted when subclassifying mature oligodendrocytes because of the simplification needed to make conclusions about cell identity from temporally confined investigations. Future studies leveraging advanced techniques like spatial transcriptomics and single-cell proteomics promise a more nuanced understanding of mature oligodendrocyte heterogeneity. Such research avenues that precisely evaluate mature oligodendrocyte heterogeneity with care to understand the mitigating influence of species, sex, central nervous system region, age, and disease, hold promise for the development of therapeutic interventions targeting varied central nervous system pathology.

Single-cell multi-omics in the study of digestive system cancers

Digestive system cancers are prevalent diseases with a high mortality rate, posing a significant threat to public health and economic burden. The diagnosis and treatment of digestive system cancer confront conventional cancer problems, such as tumor heterogeneity and drug resistance. Single-cell sequencing (SCS) emerged at times required and has developed from single-cell RNA-seq (scRNA-seq) to the single-cell multi-omics era represented by single-cell spatial transcriptomics (ST). This article comprehensively reviews the advances of single-cell omics technology in the study of digestive system tumors. While analyzing and summarizing the research cases, vital details on the sequencing platform, sample information, sampling method, and key findings are provided. Meanwhile, we summarize the commonly used SCS platforms and their features, as well as the advantages of multi-omics technologies in combination. Finally, the development trends and prospects of the application of single-cell multi-omics technology in digestive system cancer research are prospected.

The molecular subtypes of autoimmune diseases

Autoimmune diseases (ADs) are characterized by their complexity and a wide range of clinical differences. Despite patients presenting with similar symptoms and disease patterns, their reactions to treatments may vary. The current approach of personalized medicine, which relies on molecular data, is seen as an effective method to address the variability in these diseases. This review examined the pathologic classification of ADs, such as multiple sclerosis and lupus nephritis, over time. Acknowledging the limitations inherent in pathologic classification, the focus shifted to molecular classification to achieve a deeper insight into disease heterogeneity. The study outlined the established methods and findings from the molecular classification of ADs, categorizing systemic lupus erythematosus (SLE) into four subtypes, inflammatory bowel disease (IBD) into two, rheumatoid arthritis (RA) into three, and multiple sclerosis (MS) into a single subtype. It was observed that the high inflammation subtype of IBD, the RA inflammation subtype, and the MS "inflammation & EGF" subtype share similarities. These subtypes all display a consistent pattern of inflammation that is primarily driven by the activation of the JAK-STAT pathway, with the effective drugs being those that target this signaling pathway. Additionally, by identifying markers that are uniquely associated with the various subtypes within the same disease, the study was able to describe the differences between subtypes in detail. The findings are expected to contribute to the development of personalized treatment plans for patients and establish a strong basis for tailored approaches to treating autoimmune diseases.

Integrating bulk and single-cell sequencing reveals metastasis-associated CAFs in head and neck squamous cell carcinoma

AIMS

Cancer-associated fibroblasts (CAFs) have been shown to promote the metastasis of head and neck squamous cell carcinoma (HNSCC), but the underlying mechanisms remain unclear. The purpose of this study is to identify gene in CAFs that are associated with metastasis and to preliminarily validate its impact on the metastasis of HNSCC.

MATERIALS AND METHODS

Scissor analysis was utilized to process single-cell and bulk RNA sequencing datasets, identifying genes associated with the metastasis of HNSCC through differential gene expression analysis. A risk model was constructed using LASSO regression analysis. Quantitative real time-PCR and Western blot were employed to measure mRNA and protein expressions, respectively. Multiplex immunohistochemistry (mIHC) was used to assess protein expression in CAFs. siRNA was utilized to achieve gene knockdown. CCK-8 and Transwell assays were employed to evaluate the biological characteristics of HNSCC cells.

KEY FINDINGS

Compare to the nonmetastatic primary CAFs (nmCAFs), tissue inhibitors of metalloproteinase-1 (TIMP1) was founded to be overexpressed in the cells and tissues of metastatic primary CAFs (mCAFs). Knocking down TIMP1 in CAFs can signifucantly inhibit the proliferation, invasion, and migration of HNSCC cells.

SIGNIFICANCE

CAFs facilitate HNSCC cell metastasis by upregulating TIMP1, highlighting TIMP1 as a potential therapeutic target in HNSCC metastasis management.

Tumor-associated macrophage clusters linked to immunotherapy in a pan-cancer census

Transcriptional heterogeneity of tumor-associated macrophages (TAMs) has been investigated in individual cancers, but the extent to which these states transcend tumor types and represent a general feature of cancer remains unclear. We performed pan-cancer single-cell RNA sequencing analysis across nine cancer types and identified distinct monocyte/TAM composition patterns. Using spatial analysis from clinical study tissues, we assessed TAM functions in shaping the tumor microenvironment (TME) and influencing immunotherapy. Two specific TAM clusters (pro-inflammatory and pro-tumor) and four TME subtypes showed distinct immunological features, genomic profiles, immunotherapy responses, and cancer prognosis. Pro-inflammatory TAMs resided in immune-enriched niches with exhausted CD8+ T cells, while pro-tumor TAMs were restricted to niches associated with a T-cell-excluded phenotype and hypoxia. We developed a machine learning model to predict immune checkpoint blockade response by integrating TAMs and clinical data. Our study comprehensively characterizes the common features of TAMs and highlights their interaction with the TME.

Gas6-Axl signal promotes indoor VOCs exposure-induced pulmonary fibrosis via pulmonary microvascular endothelial cells-fibroblasts cross-talk

Volatile organic compounds (VOCs) as environmental pollutants were associated with respiratory diseases. Pulmonary fibrosis (PF) was characterized by an increase of extracellular matrix, leading to deterioration of lung function. The adverse effects on lung and the potential mechanism underlying VOCs induced PF had not been elucidated clearly. In this study, the indoor VOCs exposure mouse model along with an ex vivo biosensor assay was established. Based on scRNA-seq analysis, the adverse effects on lung and potential molecular mechanism were studied. Herein, the results showed that VOCs exposure from indoor decoration contributed to decreased lung function and facilitated pulmonary fibrosis in mice. Then, the whole lung cell atlas after VOCs exposure and the heterogeneity of fibroblasts were revealed. We explored the molecular interactions among various pulmonary cells, suggesting that endothelial cells contributed to fibroblasts activation in response to VOCs exposure. Mechanistically, pulmonary microvascular endothelial cells (MPVECs) secreted Gas6 after VOCs-induced PANoptosis phenotype, bound to the Axl in fibroblasts, and then activated fibroblasts. Moreover, Atf3 as the key gene negatively regulated PANoptosis phenotype to ameliorate fibrosis induced by VOCs exposure. These novel findings provided a new perspective about MPVECs could serve as the initiating factor of PF induced by VOCs exposure.

Comprehensive analysis of single cell and bulk RNA sequencing reveals the heterogeneity of melanoma tumor microenvironment and predicts the response of immunotherapy

BACKGROUND

Tumor microenvironment (TME) heterogeneity is an important factor affecting the treatment response of immune checkpoint inhibitors (ICI). However, the TME heterogeneity of melanoma is still widely characterized.

METHODS

We downloaded the single-cell sequencing data sets of two melanoma patients from the GEO database, and used the "Scissor" algorithm and the "BayesPrism" algorithm to comprehensively analyze the characteristics of microenvironment cells based on single-cell and bulk RNA-seq data. The prediction model of immunotherapy response was constructed by machine learning and verified in three cohorts of GEO database.

RESULTS

We identified seven cell types. In the Scissor+ subtype cell population, the top three were T cells, B cells and melanoma cells. In the Scissor- subtype, there are more macrophages. By quantifying the characteristics of TME, significant differences in B cells between responders and non-responders were observed. The higher the proportion of B cells, the better the prognosis. At the same time, macrophages in the non-responsive group increased significantly. Finally, nine gene features for predicting ICI response were constructed, and their predictive performance was superior in three external validation groups.

CONCLUSION

Our study revealed the heterogeneity of melanoma TME and found a new predictive biomarker, which provided theoretical support and new insights for precise immunotherapy of melanoma patients.

Data acquisition approaches for single cell proteomics

Single-cell proteomics (SCP) aims to characterize the proteome of individual cells, providing insights into complex biological systems. It reveals subtle differences in distinct cellular populations that bulk proteome analysis may overlook, which is essential for understanding disease mechanisms and developing targeted therapies. Mass spectrometry (MS) methods in SCP allow the identification and quantification of thousands of proteins from individual cells. Two major challenges in SCP are the limited material in single-cell samples necessitating highly sensitive analytical techniques and the efficient processing of samples, as each biological sample requires thousands of single cell measurements. This review discusses MS advancements to mitigate these challenges using data-dependent acquisition (DDA) and data-independent acquisition (DIA). Additionally, we examine the use of short liquid chromatography gradients and sample multiplexing methods that increase the sample throughput and scalability of SCP experiments. We believe these methods will pave the way for improving our understanding of cellular heterogeneity and its implications for systems biology.

A pro-metastatic tRNA fragment drives aldolase A oligomerization to enhance aerobic glycolysis in lung adenocarcinoma

Despite being the leading cause of lung cancer-related deaths, the underlying molecular mechanisms driving metastasis progression are still not fully understood. Transfer RNA-derived fragments (tRFs) have been implicated in various biological processes in cancer. However, the role of tRFs in lung adenocarcinoma (LUAD) remains unclear. Our study identified a tRF, tRF-Val-CAC-024, associated with the high-risk component of LUAD, through validation using 3 cohorts. Our findings demonstrated that tRF-Val-CAC-024 acts as an oncogene in LUAD. Mechanistically, tRF-Val-CAC-024 was revealed to bind to aldolase A (ALDOA) dependent on Q125/E224 and promote the oligomerization of ALDOA, resulting in increased enzyme activity and enhanced aerobic glycolysis in LUAD cells. Additionally, we provide preliminary evidence of its potential clinical value by investigating the therapeutic effects of tRF-Val-CAC-024 antagomir-loaded lipid nanoparticles (LNPs) in cell-line-derived xenograft models. These results could enhance our understanding of the regulatory mechanisms of tRFs in LUAD and provide a potential therapeutic target.

FIND-seq: high-throughput nucleic acid cytometry for rare single-cell transcriptomics

Rare cells have an important role in development and disease, and methods for isolating and studying cell subsets are therefore an essential part of biology research. Such methods traditionally rely on labeled antibodies targeted to cell surface proteins, but large public databases and sophisticated computational approaches increasingly define cell subsets on the basis of genomic, epigenomic and transcriptomic sequencing data. Methods for isolating cells on the basis of nucleic acid sequences powerfully complement these approaches by providing experimental access to cell subsets discovered in cell atlases, as well as those that cannot be otherwise isolated, including cells infected with pathogens, with specific DNA mutations or with unique transcriptional or splicing signatures. We recently developed a nucleic acid cytometry platform called 'focused interrogation of cells by nucleic acid detection and sequencing' (FIND-seq), capable of isolating rare cells on the basis of RNA or DNA markers, followed by bulk or single-cell transcriptomic analysis. This platform has previously been used to characterize the splicing-dependent activation of the transcription factor XBP1 in astrocytes and HIV persistence in memory CD4 T cells from people on long-term antiretroviral therapy. Here, we outline the molecular and microfluidic steps involved in performing FIND-seq, including protocol updates that allow detection and whole transcriptome sequencing of rare HIV-infected cells that harbor genetically intact virus genomes. FIND-seq requires knowledge of microfluidics, optics and molecular biology. We expect that FIND-seq, and this comprehensive protocol, will enable mechanistic studies of rare HIV+ cells, as well as other cell subsets that were previously difficult to recover and sequence.

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on drug response genes to predict prognosis and therapeutic response in ovarian cancer

Ovarian cancer represents a severe gynecological malignancy with a dire prognosis, underscoring the imperative need for dependable biomarkers that can accurately predict drug response and guide therapeutic choices. In this study, we harnessed online single-cell RNA sequencing (scRNAseq) and bulk RNA sequencing (RNAseq) datasets, applying the Scissor algorithm to identify cells responsive to paclitaxel. From these cells, we derived a gene signature, subsequently used to construct a prognostic model that demonstrated high sensitivity and specificity in predicting patient outcomes. Moreover, we conducted pathway and functional enrichment analyses to uncover potential molecular mechanisms driving the prognostic gene signature. This study illustrates the critical role of scRNAseq and bulk RNAseq in developing precise prognostic models for ovarian cancer, potentially transforming clinical decision-making.

An 18-gene signature of recurrence-associated endothelial cells predicts tumor progression and castration resistance in prostate cancer

BACKGROUND

The prognostic and therapeutic implications of endothelial cells (ECs) heterogeneity in prostate cancer (PCa) are poorly understood.

METHODS

We investigated associations of EC heterogeneity with PCa recurrence and castration resistance in 8 bulk transcriptomic and 4 single-cell RNA-seq cohorts. A recurrence-associated EC (RAEC) signature was constructed by comparing 11 machine learning algorithms through nested cross-validation. Functional relevances of RAEC-specific genes were also tested.

RESULTS

A subset of ECs was significantly associated with recurrence in primary PCa and named RAECs. RAECs were characteristic of tip and immature cells and were enriched in migration, angiogenesis, and collagen-related pathways. We then developed an 18-gene RAEC signature (RAECsig) representative of RAECs. Higher RAECsig scores independently predicted tumor recurrence and performed better or comparably compared to clinicopathological factors and commercial gene signatures in multiple PCa cohorts. Of the 18 RAECsig genes, FSCN1 was upregulated in ECs from PCa with higher Gleason scores; and the silencing of FSCN1, TMEME255B, or GABRD in ECs either attenuated tube formation or inhibited PCa cell proliferation. Finally, higher RAECsig scores predicted castration resistance in both primary and castration-resistant PCa.

CONCLUSION

This study establishes an endothelial signature that links a subset of ECs to prostate cancer recurrence and castration resistance.

RBMS1 reflects a distinct microenvironment and promotes tumor progression in ocular melanoma

Ocular melanoma, including uveal melanoma (UM) and conjunctival melanoma (CM), is the most common ocular cancer among adults with a high rate of recurrence and poor prognosis. Loss of epigenetic homeostasis disturbed gene expression patterns, resulting in oncogenesis. Herein, we comprehensively analyzed the DNA methylation, transcriptome profiles, and corresponding clinical information of UM patients through multiple machine-learning algorithms, finding that a methylation-driven gene RBMS1 was correlated with poor clinical outcomes of UM patients. RNA-seq and single-cell RNA-seq analyses revealed that RBMS1 reflected diverse tumor microenvironments, where high RBMS1 expression marked an immune active TME. Furthermore, we found that tumor cells were identified to have the higher communication probability in RBMS1+ state. The functional enrichment analysis revealed that RBMS1 was associated with pigment granule and melanosome, participating in cell proliferation as well as apoptotic signaling pathway. Biological experiments were performed and demonstrated that the silencing of RBMS1 inhibited ocular melanoma proliferation and promoted apoptosis. Our study highlighted that RBMS1 reflects a distinct microenvironment and promotes tumor progression in ocular melanoma, contributing to the therapeutic customization and clinical decision-making.

Integrative analysis identifies molecular features of fibroblast and the significance of fibrosis on neoadjuvant chemotherapy response in breast cancer

BACKGROUND

The molecular features of fibroblasts and the role of fibrosis in neoadjuvant chemotherapy (NAC) response and breast cancer (BRCA) prognosis remain unclear. Therefore, this study aimed to investigate the impact of interstitial fibrosis on the response and prognosis of patients with BRCA undergoing NAC treatment.

MATERIALS AND METHODS

The molecular characteristics of pathologic complete response (pCR) and non-pCR (npCR) in patients with BRCA were analyzed using multiomics analysis. A clinical cohort was collected to investigate the predictive value of fibrosis in patients with BRCA.

RESULTS

Fibrosis-related signaling pathways were significantly upregulated in patients with npCR. npCR may be associated with distinct and highly active fibroblast subtypes. Patients with high fibrosis had lower pCR rates. The fibrosis-dependent nomogram for pCR showed efficient predictive ability [training set: area under the curve [AUC]=0.871, validation set: AUC=0.792]. Patients with low fibrosis had a significantly better prognosis than those with high fibrosis, and those with a high fibrotic focus index had significantly shorter overall and recurrence-free survival. Therefore, fibrosis can be used to predict pCR. Our findings provide a basis for decision-making in the treatment of BRCA.

CONCLUSIONS

npCR is associated with a distinct and highly active fibroblast subtype. Furthermore, patients with high fibrosis have lower pCR rates and shorter long-term survival. Therefore, fibrosis can predict pCR. A nomogram that includes fibrosis can provide a basis for decision-making in the treatment of BRCA.

Multiomic analysis of monocyte-derived alveolar macrophages in idiopathic pulmonary fibrosis

BACKGROUND

Monocyte-derived alveolar macrophages (Mo_AMs) are increasingly recognised as potential pathogenic factors for idiopathic pulmonary fibrosis (IPF). While scRNAseq analysis has proven valuable in the transcriptome profiling of Mo_AMs, the integration analysis of multi-omics may provide additional dimensions of understanding of these cellular populations.

METHODS

We performed multi-omics analysis on 116 scRNAseq, 119 bulkseq and five scATACseq lung tissue samples from IPF. We built a large-scale IPF scRNAseq atlas and conducted the Monocle 2/3 as well as the Cellchat to explore the developmental path and intercellular communication on Mo_AMs. We also reported the difference in metabolisms, tissue repair and phagocytosis between Mo_AMs and tissue-resident alveolar macrophages (TRMs). To determine whether Mo_AMs affected pulmonary function, we projected clinical phenotypes (FVC%pred) from the bulkseq dataset onto the scRNAseq atlas. Finally, we used scATATCseq to uncover the upstream regulatory mechanisms and determine key drivers in Mo_AMs.

RESULTS

We identified three Mo_AMs clusters and the trajectory analysis further validated the origin of these clusters. Moreover, via the Cellchat analysis, the CXCL12/CXCR4 axis was found to be involved in the molecular basis of reciprocal interactions between Mo_AMs and fibroblasts through the activation of the ERK pathway in Mo_AMs. SPP1_RecMacs (RecMacs, recruited macrophages) were higher in the low-FVC group than in the high-FVC group. Specifically, compared with TRMs, the functions of lipid and energetic metabolism as well as tissue repair were higher in Mo_AMs than TRMs. But, TRMs may have higher level of phagocytosis than TRMs. SPIB (PU.1), JUNB, JUND, BACH2, FOSL2, and SMARCC1 showed stronger association with open chromatin of Mo_AMs than TRMs. Significant upregulated expression and deep chromatin accessibility of APOE were observed in both SPP1_RecMacs and TRMs.

CONCLUSION

Through trajectory analysis, it was confirmed that SPP1_RecMacs derived from Monocytes. Besides, Mo_AMs may influence FVC% pred and aggravate pulmonary fibrosis through the communication with fibroblasts. Furthermore, distinctive transcriptional regulators between Mo_AMs and TRMs implied that they may depend on different upstream regulatory mechanisms. Overall, this work provides a global overview of how Mo_AMs govern IPF and also helps determine better approaches and intervention therapies.

The association between vitamin D and the progression of diabetic nephropathy: insights into potential mechanisms

Aims

Vitamin D deficiency (VDD) is prevalent in the population, with inadequate intake, impaired absorption and metabolism as the main causative factors. VDD increases the risk of developing chronic diseases such as type 2 diabetes mellitus (T2DM) and diabetic nephropathy (DN), but the molecular mechanisms underlying this phenomenon are not known. The aim of this study was to investigate the association and potential mechanisms of vitamin D levels with the progression of DN by analyzing general clinical data and using bioinformatics methods.

Methods

The study included 567 diabetes mellitus type 2 (T2DM) patients from the Rocket Force Characteristic Medical Center as the case group and 221 healthy examinees as the normal control group. T2DM patients were categorized into T2DM, early diabetic nephropathy (EDN), and advanced diabetic nephropathy (ADN) based on the progression of diabetic nephropathy. The renal RNA-seq and scRNA-seq data of patients with DN were mined from public databases, and the differential expression of vitamin D-related genes in normal-EDN-ADN was analyzed by bioinformatics method, protein interaction network was constructed, immune infiltration was evaluated, single cell map was drawn, and potential mechanisms of VD and DN interaction were explored.

Results

Chi-square test showed that vitamin D level was significantly negatively correlated with DN progression (p < 0.001). Bioinformatics showed that the expression of vitamin D-related cytochrome P450 family genes was down-regulated, and TLR4 and other related inflammatory genes were abnormally up-regulated with the progression of DN. Vitamin D metabolism disturbance up-regulate "Nf-Kappa B signaling pathway," B cell receptor signaling pathway and other immune regulation and insulin resistance related pathways, and inhibit a variety of metabolic pathways. In addition, vitamin D metabolism disturbance are strongly associated with the development of diabetic cardiomyopathy and several neurological disease complications.

Conclusion

VDD or vitamin D metabolism disturbance is positively associated with the severity of renal injury. The mechanisms may involve abnormal regulation of the immune system by vitamin D metabolism disturbance, metabolic suppression, upregulation of insulin resistance and inflammatory signalling pathways.

scRank infers drug-responsive cell types from untreated scRNA-seq data using a target-perturbed gene regulatory network

Cells respond divergently to drugs due to the heterogeneity among cell populations. Thus, it is crucial to identify drug-responsive cell populations in order to accurately elucidate the mechanism of drug action, which is still a great challenge. Here, we address this problem with scRank, which employs a target-perturbed gene regulatory network to rank drug-responsive cell populations via in silico drug perturbations using untreated single-cell transcriptomic data. We benchmark scRank on simulated and real datasets, which shows the superior performance of scRank over existing methods. When applied to medulloblastoma and major depressive disorder datasets, scRank identifies drug-responsive cell types that are consistent with the literature. Moreover, scRank accurately uncovers the macrophage subpopulation responsive to tanshinone IIA and its potential targets in myocardial infarction, with experimental validation. In conclusion, scRank enables the inference of drug-responsive cell types using untreated single-cell data, thus providing insights into the cellular-level impacts of therapeutic interventions.

APOE2 protects against Aβ pathology by improving neuronal mitochondrial function through ERRα signaling

BACKGROUND

Alzheimer's disease (AD) is a progressive neurodegenerative disease and apolipoprotein E (APOE) genotypes (APOE2, APOE3, and APOE4) show different AD susceptibility. Previous studies indicated that individuals carrying the APOE2 allele reduce the risk of developing AD, which may be attributed to the potential neuroprotective role of APOE2. However, the mechanisms underlying the protective effects of APOE2 is still unclear.

METHODS

We analyzed single-nucleus RNA sequencing and bulk RNA sequencing data of APOE2 and APOE3 carriers from the Religious Orders Study and Memory and Aging Project (ROSMAP) cohort. We validated the findings in SH-SY5Y cells and AD model mice by evaluating mitochondrial functions and cognitive behaviors respectively.

RESULTS

The pathway analysis of six major cell types revealed a strong association between APOE2 and cellular stress and energy metabolism, particularly in excitatory and inhibitory neurons, which was found to be more pronounced in the presence of beta-amyloid (Aβ). Moreover, APOE2 overexpression alleviates Aβ1-42-induced mitochondrial dysfunction and reduces the generation of reactive oxygen species in SH-SY5Y cells. These protective effects may be due to ApoE2 interacting with estrogen-related receptor alpha (ERRα). ERRα overexpression by plasmids or activation by agonist was also found to show similar mitochondrial protective effects in Aβ1-42-stimulated SH-SY5Y cells. Additionally, ERRα agonist treatment improve the cognitive performance of Aβ injected mice in both Y maze and novel object recognition tests. ERRα agonist treatment increased PSD95 expression in the cortex of agonist-treated-AD mice.

CONCLUSIONS

APOE2 appears to enhance neural mitochondrial function via the activation of ERRα signaling, which may be the protective effect of APOE2 to treat AD.

Unveiling inflammatory and prehypertrophic cell populations as key contributors to knee cartilage degeneration in osteoarthritis using multi-omics data integration

OBJECTIVES

Single-cell and spatial transcriptomics analysis of human knee articular cartilage tissue to present a comprehensive transcriptome landscape and osteoarthritis (OA)-critical cell populations.

METHODS

Single-cell RNA sequencing and spatially resolved transcriptomic technology have been applied to characterise the cellular heterogeneity of human knee articular cartilage which were collected from 8 OA donors, and 3 non-OA control donors, and a total of 19 samples. The novel chondrocyte population and marker genes of interest were validated by immunohistochemistry staining, quantitative real-time PCR, etc. The OA-critical cell populations were validated through integrative analyses of publicly available bulk RNA sequencing data and large-scale genome-wide association studies.

RESULTS

We identified 33 cell population-specific marker genes that define 11 chondrocyte populations, including 9 known populations and 2 new populations, that is, pre-inflammatory chondrocyte population (preInfC) and inflammatory chondrocyte population (InfC). The novel findings that make this an important addition to the literature include: (1) the novel InfC activates the mediator MIF-CD74; (2) the prehypertrophic chondrocyte (preHTC) and hypertrophic chondrocyte (HTC) are potentially OA-critical cell populations; (3) most OA-associated differentially expressed genes reside in the articular surface and superficial zone; (4) the prefibrocartilage chondrocyte (preFC) population is a major contributor to the stratification of patients with OA, resulting in both an inflammatory-related subtype and a non-inflammatory-related subtype.

CONCLUSIONS

Our results highlight InfC, preHTC, preFC and HTC as potential cell populations to target for therapy. Also, we conclude that profiling of those cell populations in patients might be used to stratify patient populations for defining cohorts for clinical trials and precision medicine.

Establishment of a prognostic model for hypoxia-associated genes in OPSCC and revelation of intercellular crosstalk

Hypoxia exerts a profound influence on the tumor microenvironment and immune response, shaping treatment outcomes and prognosis. Utilizing consistency clustering, we discerned two hypoxia subtypes in OPSCC bulk sequencing data from GEO. Key modules within OPSCC were identified through weighted gene correlation network analysis (WGCNA). Core modules underwent CIBERSORT immune infiltration analysis and GSEA functional enrichment. Univariate Cox and LASSO analyses were employed to construct prognostic models for seven hypoxia-related genes. Further investigation into clinical characteristics, the immune microenvironment, and TIDE algorithm prediction for immunotherapy response was conducted in high- and low-risk groups. scRNA-seq data were visually represented through TSNE clustering, employing the scissors algorithm to map hypoxia phenotypes. Interactions among cellular subpopulations were explored using the Cellchat package, with additional assessments of metabolic and transcriptional activities. Integration with clinical data unveiled a prevalence of HPV-positive patients in the low hypoxia and low-risk groups. Immunohistochemical validation demonstrated low TDO2 expression in HPV-positive (P16-positive) patients. Our prediction suggested that HPV16 E7 promotes HIF-1α inhibition, leading to reduced glycolytic activity, ultimately contributing to better prognosis and treatment sensitivity. The scissors algorithm effectively segregated epithelial cells and fibroblasts into distinct clusters based on hypoxia characteristics. Cellular communication analysis illuminated significant crosstalk among hypoxia-associated epithelial, fibroblast, and endothelial cells, potentially fostering tumor proliferation and metastasis.

Integrative Single Cell Atlas Revealed Intratumoral Heterogeneity Generation from an Adaptive Epigenetic Cell State in Human Bladder Urothelial Carcinoma

Intratumor heterogeneity (ITH) of bladder cancer (BLCA) contributes to therapy resistance and immune evasion affecting clinical prognosis. The molecular and cellular mechanisms contributing to BLCA ITH generation remain elusive. It is found that a TM4SF1-positive cancer subpopulation (TPCS) can generate ITH in BLCA, evidenced by integrative single cell atlas analysis. Extensive profiling of the epigenome and transcriptome of all stages of BLCA revealed their evolutionary trajectories. Distinct ancestor cells gave rise to low-grade noninvasive and high-grade invasive BLCA. Epigenome reprograming led to transcriptional heterogeneity in BLCA. During early oncogenesis, epithelial-to-mesenchymal transition generated TPCS. TPCS has stem-cell-like properties and exhibited transcriptional plasticity, priming the development of transcriptionally heterogeneous descendent cell lineages. Moreover, TPCS prevalence in tumor is associated with advanced stage cancer and poor prognosis. The results of this study suggested that bladder cancer interacts with its environment by acquiring a stem cell-like epigenomic landscape, which might generate ITH without additional genetic diversification.

Single-cell transcriptome atlas reveals somatic cell embryogenic differentiation features during regeneration

Understanding somatic cell totipotency remains a challenge facing scientific inquiry today. Plants display remarkable cell totipotency expression, illustrated by single-cell differentiation during somatic embryogenesis (SE) for plant regeneration. Determining cell identity and exploring gene regulation in such complex heterogeneous somatic cell differentiation have been major challenges. Here, we performed high-throughput single-cell sequencing assays to define the precise cellular landscape and revealed the modulation mode of marker genes during embryogenic differentiation in cotton (Gossypium hirsutum L.) as the crop for biotechnology application. We demonstrated that nonembryogenic calli (NEC) and primary embryogenic calli (PEC) tissues were composed of heterogeneous cells that could be partitioned into four broad populations with six distinct cell clusters. Enriched cell clusters and cell states were identified in NEC and PEC samples, respectively. Moreover, a broad repertoire of new cluster-specific genes and associated expression modules were identified. The energy metabolism, signal transduction, environmental adaptation, membrane transport pathways, and a series of transcription factors were preferentially enriched in cell embryogenic totipotency expression. Notably, the SE-ASSOCIATED LIPID TRANSFER PROTEIN (SELTP) gene dose-dependently marked cell types with distinct embryogenic states and exhibited a parabolic curve pattern along the somatic cell embryogenic differentiation trajectory, suggesting that SELTP could serve as a favorable quantitative cellular marker for detecting embryogenic expression at the single-cell level. In addition, RNA velocity and Scissor analysis confirmed the pseudo-temporal model and validated the accuracy of the scRNA-seq data, respectively. This work provides valuable marker-genes resources and defines precise cellular taxonomy and trajectory atlases for somatic cell embryogenic differentiation in plant regeneration.

Bulk and single-cell RNA sequencing reveal the contribution of laminin γ2 -CD44 to the immune resistance in lymphocyte-infiltrated squamous lung cancer subtype

The high heterogeneity of lung squamous cell carcinomas (LUSC) and the complex tumor microenvironment lead to non-response to immunotherapy in many patients. Therefore, characterizing the heterogeneity of the tumor microenvironment in patients with LUSC and further exploring the immune features and molecular mechanisms that lead to immune resistance will help improve the efficacy of immunotherapy in such patients. Herein, we retrospectively analyzed the RNA sequencing (RNA-seq) data of 513 LUSC samples with other multiomics and single-cell RNA-seq data and validated key features using multiplex immunohistochemistry. We divided these samples into six subtypes (CS1-CS6) based on the RNA-seq data and found that CS3 activates the immune response with a high level of lymphocyte infiltration and gathers a large number of patients with advanced-stage disease but increases the expression of exhausted markers cytotoxic T-lymphocyte associated protein 4, lymphocyte-activation gene 3, and programmed death-1. The prediction of the response to immunotherapy showed that CS3 is potentially resistant to immune checkpoint blockade therapy, and multi-omic data analysis revealed that CS3 specifically expresses immunosuppression-related proteins B cell lymphoma 2, GRB2-associated binding protein, and dual-specificity phosphatase 4 and has a high mutation ratio of the driver gene ATP binding cassette subfamily A member 13. Furthermore, single-cell RNA-seq verified lymphocyte infiltration in the CS3 subtype and revealed a positive relationship between the expression of LAMC2-CD44 and immune resistance. LAMC2 and CD44 are epithelial-mesenchymal transition-associated genes that modulate tumor proliferation, and multicolor immunofluorescence validated the negative relationship between the expression of LAMC2-CD44 and immune infiltration. Thus, we identified a lymphocyte-infiltrated subtype (CS3) in patients with LUSC that exhibited resistance to immune checkpoint blockade therapy, and the co-hyperexpression of LAMC2-CD44 contributed to immune resistance, which could potentially improve immunological efficacy by targeting this molecule pair in combination with immunotherapy.

Single-Cell RNA Sequencing Identifies Crucial Genes Influencing the Polarization of Tumor-Associated Macrophages in Liver Cancer

Background

In the context of hepatocellular carcinoma (HCC), tumor-associated macrophages (TAMs) are pivotal for the immunosuppressive nature of the tumor microenvironment (TME). This investigation delves into the functional transformations of TAMs within the TME by leveraging single-cell transcriptomics to pinpoint critical genes influencing TAM subset polarization.

Methods

We procured single-cell and bulk transcriptomic data from the Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA), implementing quality assurance, dimensional reduction, clustering, and annotation on the single-cell sequencing data. To examine cellular interactions, CellChat was utilized, while single-cell regulatory network inference and clustering (SCENIC) was applied to deduce transcription factors (TFs) and their associated targets. Through gene enrichment, survival, and immune infiltration correlation analyses, we sought to pinpoint and validate influential genes. A TAM model under HCC conditions was then established to confirm the expression levels of these key genes.

Results

Our analysis encompassed 74,742 cells and 23,110 genes. Through postdimensional reduction and clustering, we identified seven distinct cell types and nine TAM subtypes. Analysis via CellChat highlighted a predominance of M2-phenotype-inclined TAM subsets within the tumor's core. SCENIC pinpointed the transcription factor PRDM1 and its target genes as pivotal in this region. Further analysis indicated these genes' involvement in macrophage polarization. Employing trajectory analysis, survival analysis, and immune infiltration correlation, we scrutinized and validated genes likely directing M2 polarization. Experimental validation confirmed PRDM1's heightened expression in TAMs conditioned by HCC.

Conclusions

Our findings suggest the PRDM1 gene is a key regulator of M2 macrophage polarization, contributing to the immunosuppressive TME in HCC.

PIPET: predicting relevant subpopulations in single-cell data using phenotypic information from bulk data

Single-cell RNA sequencing has revealed cellular heterogeneity in complex tissues, notably benefiting research on diseases such as cancer. However, the integration of single-cell data from small samples with extensive clinical features in bulk data remains underexplored. In this study, we introduce PIPET, an algorithmic method for predicting relevant subpopulations in single-cell data based on multivariate phenotypic information from bulk data. PIPET generates feature vectors for each phenotype from differentially expressed genes in bulk data and then identifies relevant cellular subpopulations by assessing the similarity between single-cell data and these vectors. Subsequently, phenotype-related cell states can be analyzed based on these subpopulations. In simulated datasets, PIPET showed robust performance in predicting multiclassification cellular subpopulations. Application of PIPET to lung adenocarcinoma single-cell RNA sequencing data revealed cellular subpopulations with poor survival and associations with TP53 mutations. Similarly, in breast cancer single-cell data, PIPET identified cellular subpopulations associated with the PAM50 clinical subtypes and triple-negative breast cancer subtypes. Overall, PIPET effectively identified relevant cellular subpopulations in single-cell data, guided by phenotypic information from bulk data. This approach comprehensively delineates the molecular characteristics of each cellular subpopulation, offering insights into disease-related subpopulations and guiding personalized treatment strategies.

Machine learning identifies the role of SMAD6 in the prognosis and drug susceptibility in bladder cancer

BACKGROUND

Bladder cancer (BCa) is among the most prevalent malignant tumors affecting the urinary system. Due to its highly recurrent nature, standard treatments such as surgery often fail to significantly improve patient prognosis. Our research aims to predict prognosis and identify precise therapeutic targets for novel treatment interventions.

METHODS

We collected and screened genes related to the TGF-β signaling pathway and performed unsupervised clustering analysis on TCGA-BLCA samples based on these genes. Our analysis revealed two novel subtypes of bladder cancer with completely different biological characteristics, including immune microenvironment, drug sensitivity, and more. Using machine learning classifiers, we identified SMAD6 as a hub gene contributing to these differences and further investigated the role of SMAD6 in bladder cancer in the single-cell transcriptome data. Additionally, we analyzed the relationship between SMAD6 and immune checkpoint genes. Finally, we performed a series of in vitro assays to verify the function of SMAD6 in bladder cancer cell lines.

RESULTS

We have revealed two novel subtypes of bladder cancer, among which C1 exhibits a worse prognosis, lower drug sensitivity, a more complex tumor microenvironment, and a 'colder' immune microenvironment compared to C2. We identified SMAD6 as a key gene responsible for the differences and further explored its impact on the molecular characteristics of bladder cancer. Through in vitro experiments, we found that SMAD6 promoted the prognosis of BCa patients by inhibiting the proliferation and migration of BCa cells.

CONCLUSION

Our study reveals two novel subtypes of BCa and identifies SMAD6 as a highly promising therapeutic target.

Tumor niche network-defined subtypes predict immunotherapy response of esophageal squamous cell cancer

Despite the promising outcomes of immune checkpoint inhibitors (ICIs), resistance to ICI presents a new challenge. Therefore, selecting patients for specific ICI applications is crucial for maximizing therapeutic efficacy. Herein, we curated 69 human esophageal squamous cell cancer (ESCC) patients' tumor microenvironment (TME) single-cell transcriptomic datasets to subtype ESCC. Integrative analyses of the cellular network and transcriptional signatures of T cells and myeloid cells define distinct ESCC subtypes characterized by T cell exhaustion, and interleukin (IL) and interferon (IFN) signaling. Furthermore, this approach classifies ESCC patients into ICI responders and non-responders, as validated by whole tumor transcriptomes and liquid biopsy-based single-cell transcriptomes of anti-PD-1 ICI responders and non-responders. Our study stratifies ESCC patients based on TME transcriptional network, providing novel insights into tumor niche remodeling and potentially predicting ICI responses in ESCC patients.

Mutant U2AF1-Induced Mis-Splicing of mRNA Translation Genes Confers Resistance to Chemotherapy in Acute Myeloid Leukemia

Patients with primary refractory acute myeloid leukemia (AML) have a dismal long-term prognosis. Elucidating the resistance mechanisms to induction chemotherapy could help identify strategies to improve AML patient outcomes. Herein, we retrospectively analyzed the multiomics data of more than 1,500 AML cases and found that patients with spliceosome mutations had a higher risk of developing refractory disease. RNA splicing analysis revealed that the mis-spliced genes in refractory patients converged on translation-associated pathways, promoted mainly by U2AF1 mutations. Integrative analyses of binding and splicing in AML cell lines substantiated that the splicing perturbations of mRNA translation genes originated from both the loss and gain of mutant U2AF1 binding. In particular, the U2AF1S34F and U2AF1Q157R mutants orchestrated the inclusion of exon 11 (encoding a premature termination codon) in the eukaryotic translation initiation factor 4A2 (EIF4A2). This aberrant inclusion led to reduced eIF4A2 protein expression via nonsense-mediated mRNA decay. Consequently, U2AF1 mutations caused a net decrease in global mRNA translation that induced the integrated stress response (ISR) in AML cells, which was confirmed by single-cell RNA sequencing. The induction of ISR enhanced the ability of AML cells to respond and adapt to stress, contributing to chemoresistance. A pharmacologic inhibitor of ISR, ISRIB, sensitized U2AF1 mutant cells to chemotherapy. These findings highlight a resistance mechanism by which U2AF1 mutations drive chemoresistance and provide a therapeutic approach for AML through targeting the ISR pathway.

SIGNIFICANCE

U2AF1 mutations induce the integrated stress response by disrupting splicing of mRNA translation genes that improves AML cell fitness to enable resistance to chemotherapy, which can be targeted to improve AML treatment.

Machine-learning developed an iron, copper, and sulfur-metabolism associated signature predicts lung adenocarcinoma prognosis and therapy response

BACKGROUND

Previous studies have largely neglected the role of sulfur metabolism in LUAD, and no study has combine iron, copper, and sulfur-metabolism associated genes together to create prognostic signatures.

METHODS

This study encompasses 1564 LUAD patients, 1249 NSCLC patients, and over 10,000 patients with various cancer types from diverse cohorts. We employed the R package ConsensusClusterPlus to separate patients into different ICSM (Iron, Copper, and Sulfur-Metabolism) subtypes. Various machine-learning methods were utilized to develop the ICSMI. Enrichment analyses were conducted using ClusterProfiler and GSVA, while IOBR quantified immune cell infiltration. GISTIC2.0 and maftools were utilized for CNV and SNV data analysis. The Oncopredict package predicted drug information based on GDSC1. TIDE algorithm and cohorts GSE91061 and IMvigor210 evaluated patient response to immunotherapy. Single-cell data was processed using the Seurat package, AUCell package calculated cells geneset activity scores, and the Scissor algorithm identified ICSMI-associated cells. In vitro experiments was conducted to explore the role of ICSMRGs in LUAD.

RESULTS

Unsupervised clustering identified two distinct ICSM subtypes of LUAD, each with unique clinical characteristics. The ICSMI, comprising 10 genes, was constructed using integrated machine-learning methods. Its prognostic power was validated in 10 independent datasets, revealing that LUAD patients with higher ICSMI levels had poorer prognoses. Furthermore, ICSMI demonstrated superior predictive abilities compared to 102 previously published signatures. A nomogram incorporating ICSMI and clinical features exhibited high predictive performance. ICSMI positively correlated with patients gene mutations, and integrated analysis of bulk and single-cell transcriptome data revealed its association with TME modulators. Cells representing the high-ICSMI phenotype exhibited more malignant features. LUAD patients with high ICSMI levels exhibited sensitivity to chemotherapy and targeted therapy but displayed resistance to immunotherapy. In a comprehensive analysis across various cancers, ICSMI retained significant prognostic value and emerged as a risk factor for the majority of cancer patients.

CONCLUSIONS

ICSMI provides critical prognostic insights for LUAD patients, offering valuable insights into the tumor microenvironment and predicting treatment responsiveness.

SCIPAC: quantitative estimation of cell-phenotype associations

Numerous algorithms have been proposed to identify cell types in single-cell RNA sequencing data, yet a fundamental problem remains: determining associations between cells and phenotypes such as cancer. We develop SCIPAC, the first algorithm that quantitatively estimates the association between each cell in single-cell data and a phenotype. SCIPAC also provides a p-value for each association and applies to data with virtually any type of phenotype. We demonstrate SCIPAC's accuracy in simulated data. On four real cancerous or noncancerous datasets, insights from SCIPAC help interpret the data and generate new hypotheses. SCIPAC requires minimum tuning and is computationally very fast.

HCCDB v2.0: Decompose Expression Variations by Single-cell RNA-seq and Spatial Transcriptomics in HCC

Large-scale transcriptomic data are crucial for understanding the molecular features of hepatocellular carcinoma (HCC). Integrated 15 transcriptomic datasets of HCC clinical samples, the first version of HCC database (HCCDB v1.0) was released in 2018. Through the meta-analysis of differentially expressed genes and prognosis-related genes across multiple datasets, it provides a systematic view of the altered biological processes and the inter-patient heterogeneities of HCC with high reproducibility and robustness. With four years having passed, the database now needs integration of recently published datasets. Furthermore, the latest single-cell and spatial transcriptomics have provided a great opportunity to decipher complex gene expression variations at the cellular level with spatial architecture. Here, we present HCCDB v2.0, an updated version that combines bulk, single-cell, and spatial transcriptomic data of HCC clinical samples. It dramatically expands the bulk sample size by adding 1656 new samples from 11 datasets to the existing 3917 samples, thereby enhancing the reliability of transcriptomic meta-analysis. A total of 182,832 cells and 69,352 spatial spots are added to the single-cell and spatial transcriptomics sections, respectively. A novel single-cell level and 2-dimension (sc-2D) metric is proposed as well to summarize cell type-specific and dysregulated gene expression patterns. Results are all graphically visualized in our online portal, allowing users to easily retrieve data through a user-friendly interface and navigate between different views. With extensive clinical phenotypes and transcriptomic data in the database, we show two applications for identifying prognosis-associated cells and tumor microenvironment. HCCDB v2.0 is available at http://lifeome.net/database/hccdb2.

Disruption of maternal vascular remodeling by a fetal endoretrovirus-derived gene in preeclampsia

BACKGROUND

Preeclampsia, one of the most lethal pregnancy-related diseases, is associated with the disruption of uterine spiral artery remodeling during placentation. However, the early molecular events leading to preeclampsia remain unknown.

RESULTS

By analyzing placentas from preeclampsia, non-preeclampsia, and twin pregnancies with selective intrauterine growth restriction, we show that the pathogenesis of preeclampsia is attributed to immature trophoblast and maldeveloped endothelial cells. Delayed epigenetic reprogramming during early extraembryonic tissue development leads to generation of excessive immature trophoblast cells. We find reduction of de novo DNA methylation in these trophoblast cells results in selective overexpression of maternally imprinted genes, including the endoretrovirus-derived gene PEG10 (paternally expressed gene 10). PEG10 forms virus-like particles, which are transferred from the trophoblast to the closely proximate endothelial cells. In normal pregnancy, only a low amount of PEG10 is transferred to maternal cells; however, in preeclampsia, excessive PEG10 disrupts maternal vascular development by inhibiting TGF-beta signaling.

CONCLUSIONS

Our study reveals the intricate epigenetic mechanisms that regulate trans-generational genetic conflict and ultimately ensure proper maternal-fetal interface formation.

Identification of a novel ADCC-related gene signature for predicting the prognosis and therapy response in lung adenocarcinoma

BACKGROUND

Previous studies have largely neglected the role of ADCC in LUAD, and no study has systematically compiled ADCC-associated genes to create prognostic signatures.

METHODS

In this study, 1564 LUAD patients, 2057 NSCLC patients, and more than 5000 patients with various cancer types from diverse cohorts were included. R package ConsensusClusterPlus was utilized to classify patients into different subtypes. A number of machine-learning algorithms were used to construct the ADCCRS. GSVA and ClusterProfiler were used for enrichment analyses, and IOBR was used to quantify immune cell infiltration level. GISTIC2.0 and maftools were used to analyze the CNV and SNV data. The Oncopredict package was used to predict drug information based on the GDSC1. Three immunotherapy cohorts were used to evaluate patient response to immunotherapy. The Seurat package was used to process single-cell data, the AUCell package was used to calculate cells' geneset activity scores, and the Scissor algorithm was used to identify ADCCRS-associated cells.

RESULTS

Through unsupervised clustering, two distinct subtypes of LUAD were identified, each exhibiting distinct clinical characteristics. The ADCCRS, consisted of 16 genes, was constructed by integrated machine-learning methods. The prognostic power of ADCCRS was validated in 28 independent datasets. Further, ADCCRS shows better predictive abilities than 102 previously published signatures in predicting LUAD patients' survival. A nomogram incorporating ADCCRS and clinical features was constructed, demonstrating high predictive performance. ADCCRS positively correlates with patients' gene mutation, and integrated analysis of bulk and single-cell transcriptome data revealed the association of ADCCRS with TME modulators. Cells representing high-ADCCRS phenotype exhibited more malignant features. LUAD patients with high ADCCRS levels exhibited sensitivity to chemotherapy and targeted therapy, while displaying resistance to immunotherapy. In pan-cancer analysis, ADCCRS still exhibited significant prognostic value and was found to be a risk factor for most cancer patients.

CONCLUSIONS

ADCCRS offers a critical prognostic insight for patients with LUAD, shedding light on the tumor microenvironment and forecasting treatment responsiveness.

Single-cell transcriptome analysis profiling lymphatic invasion-related TME in colorectal cancer

Lymphatic invasion (LI) is extremely aggressive and induces worse prognosis among patients with colorectal cancer (CRC). Thus, it is critical to characterize the cellular and molecular mechanisms underlying LI in order to establish novel and efficacious therapeutic targets that enhance the prognosis of CRC patients. RNA-seq data, clinical and survival information of colon adenocarcinoma (COAD) patients were obtained from the TCGA database. In addition, three scRNA-seq datasets of CRC patients were acquired from the GEO database. Data analyses were conducted with the R packages. We assessed the tumor microenvironment (TME) differences between LI+ and LI- based scRNA-seq data, LI+ cells exhibited augmented abundance of immunosuppression and invasive subset. Marked extracellular matrix network activation was also observed in LI+ cells within SPP1+ macrophages. We revealed that an immunosuppressive and pro-angiogenic TME strongly enhanced LI, as was evidenced by the CD4+ Tregs, CD8+ GZMK+, SPP1+ macrophages, e-myCAFs, and w-myCAFs subcluster infiltrations. Furthermore, we identified potential LI targets that influenced tumor development, metastasis, and immunotherapeutic response. Finally, a novel LIRS model was established based on the expression of 14 LI-related signatures, and in the two testing cohorts, LIRS was also proved to have accurate prognostic predictive ability. In this report, we provided a valuable resource and extensive insights into the LI of CRC. Our conclusions can potentially benefit the establishment of highly efficacious therapeutic targets as well as diagnostic biomarkers that improve patient outcomes.

Integrated analysis of single-cell RNA-seq and bulk RNA-seq reveals immune suppression subtypes and establishes a novel signature for determining the prognosis in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is the most common histological type of lung cancer with lower survival rates. Recent advancements in targeted therapies and immunotherapies targeting immune checkpoints have achieved remarkable success, there is still a large percentage of LUAD that lacks available therapeutic options. Due to tumor heterogeneity, the diagnosis and treatment of LUAD are challenging. Exploring the biology of LUAD and identifying new biomarker and therapeutic targets options are essential.

METHOD

We performed single-cell RNA sequencing (scRNA-seq) of 6 paired primary and adjacent LUAD tissues, and integrative omics analysis of the scRNA-seq, bulk RNA-seq and whole-exome sequencing data revealed molecular subtype characteristics. Our experimental results confirm that CDC25C gene can serve as a potential marker for poor prognosis in LUAD.

RESULTS

We investigated aberrant gene expression in diverse cell types in LUAD via the scRNA-seq data. Moreover, multi-omics clustering revealed four subgroups defined by transcriptional profile and molecular subtype 4 (MS4) with poor survival probability, and immune cell infiltration signatures revealed that MS4 tended to be the immunosuppressive subtype. Our study revealed that the CDC25C gene can be a distinct prognostic biomarker that indicates immune infiltration levels and response to immunotherapy in LUAD patients. Our experimental results concluded that CDC25C expression affects lung cancer cell invasion and migration, might play a key role in regulating Epithelial-Mesenchymal Transition (EMT) pathways.

CONCLUSIONS

Our multi-omics result revealed a comprehensive set of molecular attributes associated with prognosis-related genes in LUAD at the cellular and tissue level. Identification of a subtype of immunosuppressive TME and prognostic signature for LUAD. We identified the cell cycle regulation gene CDC25C affects lung cancer cell invasion and migration, which can be used as a potential biomarker for LUAD.

Integrating bulk and single-cell sequencing data to construct a Scissor+ dendritic cells prognostic model for predicting prognosis and immune responses in ESCC

Esophageal squamous cell carcinoma (ESCC) is characterized by molecular heterogeneity with various immune cell infiltration patterns, which have been associated with therapeutic sensitivity and resistance. In particular, dendritic cells (DCs) are recently discovered to be associated with prognosis and survival in cancer. However, how DCs differ among ESCC patients has not been fully comprehended. Recently, the advance of single-cell RNA sequencing (scRNA-seq) enables us to profile the cell types, states, and lineages in the heterogeneous ESCC tissues. Here, we dissect the ESCC tumor microenvironment at high resolution by integrating 192,078 single cells from 60 patients, including 4379 DCs. We then used Scissor, a method that identifies cell subpopulations from single-cell data that are associated bulk samples with genomic and clinical information, to stratify DCs into Scissorhi and Scissorlow subtypes. We applied the Scissorhi gene signature to stratify ESCC scRNAseq patient, and we found that PD-L1, TIGIT, PVR and IL6 ligand-receptor-mediated cell interactions existed mainly in Scissorhi patients. Finally, based on the Scissor results, we successfully developed a validated prognostic risk model for ESCC and further validated the reliability of the risk prediction model by recruiting 40 ESCC clinical patients. This information highlights the importance of these genes in assessing patient prognosis and may help in the development of targeted or personalized therapies for ESCC.

Integrated analysis of scRNA-seq and bulk RNA-seq identifies FBXO2 as a candidate biomarker associated with chemoresistance in HGSOC

Background

High-grade serous ovarian carcinoma (HGSOC) is the most prevalent and aggressive histological subtype of epithelial ovarian cancer. Around 80% of individuals will experience a recurrence within five years because of resistance to chemotherapy, despite initially responding well to platinum-based treatment. Biomarkers associated with chemoresistance are desperately needed in clinical practice.

Methods

We jointly analyzed the transcriptomic profiles of single-cell and bulk datasets of HGSOC to identify cell types associated with chemoresistance. Copy number variation (CNV) inference was performed to identify malignant cells. We subsequently analyzed the expression of candidate biomarkers and their relationship with patients' prognosis. The enrichment analysis and potential biological function of candidate biomarkers were explored. Then, we validated the candidate biomarker using in vitro experiments.

Results

We identified 8871 malignant epithelial cells in a single-cell RNA sequencing dataset, of which 861 cells were associated with chemoresistance. Among these malignant epithelial cells, FBXO2 (F-box protein 2) is highly expressed in cells related to chemoresistance. Moreover, FBXO2 expression was found to be higher in epithelial cells from chemoresistance samples compared to those from chemosensitivity samples in a separate single-cell RNA sequencing dataset. Patients exhibiting elevated levels of FBXO2 experienced poorer outcomes in terms of both overall survival (OS) and progression-free survival (PFS). FBXO2 could impact chemoresistance by influencing the PI3K-Akt signaling pathway, focal adhesion, and ECM-receptor interactions and regulating tumorigenesis. The 50% maximum inhibitory concentration (IC50) of cisplatin decreased in A2780 and SKOV3 ovarian carcinoma cell lines with silenced FBXO2 during an in vitro experiment.

Conclusions

We determined that FBXO2 is a potential biomarker linked to chemoresistance in HGSOC by combining single-cell RNA-seq and bulk RNA-seq dataset. Our results suggest that FBXO2 could serve as a valuable prognostic marker and potential target for drug development in HGSOC.

Single-cell analysis identified POSTN+ cells associated with the aggressive phenotype and risk of esophageal squamous cell carcinoma

Tumors are intricate and heterogeneous systems characterized by mosaic cancer cell populations with diverse expression profiles. Leveraging single-cell technologies, we employed the Scissor algorithm to delineate an epithelial subpopulation associated with the aggressive phenotype in esophageal squamous cell carcinoma (ESCC). This identified subpopulation exhibited elevated expression of genes involved in critical pathways, such as epithelial-mesenchymal transition and PI3K-Akt. Key signature genes within this subpopulation, namely CAV1, COL3A1, COL6A1, POSTN, and TAGLN, demonstrated significant upregulation concomitant with both tumorigenesis and tumor progression across independent single-cell datasets. Furthermore, we selected 1,450 expression quantitative trait loci of the top 62 signature genes of this cell subpopulation to investigate their potential in predicting ESCC risk. The results showed that the POSTN loci were predominantly associated with ESCC susceptibility. Through functional annotation and replication analyses, we identified that the rs1028728 in the POSTN promoter was significantly associated with increased ESCC risk in 7,049 ESCC cases and 8,063 controls (odds ratio = 1.29, 95% confidence interval: 1.18-1.42, p = 4.03 × 10-8). Subsequent biochemical experiments showed that the rs1028728[T] allele enhanced POSTN expression by affecting the binding of PRRX1 in the POSTN promoter. In summary, our meticulous single-cell analysis delineates an invasive epithelial subpopulation in ESCC, with POSTN emerging as an important marker for the aggressive phenotype. These findings offer more insights into potential strategies for the prevention and intervention of ESCC, enriching our understanding of this complex cancer landscape.

E-cadherin loss drives diffuse-type gastric tumorigenesis via EZH2-mediated reprogramming

Diffuse-type gastric adenocarcinoma (DGAC) is a deadly cancer often diagnosed late and resistant to treatment. While hereditary DGAC is linked to CDH1 mutations, the role of CDH1/E-cadherin inactivation in sporadic DGAC tumorigenesis remains elusive. We discovered CDH1 inactivation in a subset of DGAC patient tumors. Analyzing single-cell transcriptomes in malignant ascites, we identified two DGAC subtypes: DGAC1 (CDH1 loss) and DGAC2 (lacking immune response). DGAC1 displayed distinct molecular signatures, activated DGAC-related pathways, and an abundance of exhausted T cells in ascites. Genetically engineered murine gastric organoids showed that Cdh1 knock-out (KO), KrasG12D, Trp53 KO (EKP) accelerates tumorigenesis with immune evasion compared with KrasG12D, Trp53 KO (KP). We also identified EZH2 as a key mediator promoting CDH1 loss-associated DGAC tumorigenesis. These findings highlight DGAC's molecular diversity and potential for personalized treatment in CDH1-inactivated patients.

Applications of single‑cell omics and spatial transcriptomics technologies in gastric cancer (Review)

Gastric cancer (GC) is a prominent contributor to global cancer-related mortalities, and a deeper understanding of its molecular characteristics and tumor heterogeneity is required. Single-cell omics and spatial transcriptomics (ST) technologies have revolutionized cancer research by enabling the exploration of cellular heterogeneity and molecular landscapes at the single-cell level. In the present review, an overview of the advancements in single-cell omics and ST technologies and their applications in GC research is provided. Firstly, multiple single-cell omics and ST methods are discussed, highlighting their ability to offer unique insights into gene expression, genetic alterations, epigenomic modifications, protein expression patterns and cellular location in tissues. Furthermore, a summary is provided of key findings from previous research on single-cell omics and ST methods used in GC, which have provided valuable insights into genetic alterations, tumor diagnosis and prognosis, tumor microenvironment analysis, and treatment response. In summary, the application of single-cell omics and ST technologies has revealed the levels of cellular heterogeneity and the molecular characteristics of GC, and holds promise for improving diagnostics, personalized treatments and patient outcomes in GC.

Single-cell and bulk RNA sequencing highlights the role of M1-like infiltrating macrophages in antibody-mediated rejection after kidney transplantation

Background

Antibody-mediated rejection (ABMR) significantly affects transplanted kidney survival, yet the macrophage phenotype, ontogeny, and mechanisms in ABMR remain unclear.

Method

We analyzed post-transplant sequencing and clinical data from GEO and ArrayExpress. Using dimensionality reduction and clustering on scRNA-seq data, we identified macrophage subpopulations and compared their infiltration in ABMR and non-rejection cases. Cibersort quantified these subpopulations in bulk samples. Cellchat, SCENIC, monocle2, and monocle3 helped explore intercellular interactions, predict transcription factors, and simulate differentiation of cell subsets. The Scissor method linked macrophage subgroups with transplant prognosis. Furthermore, hdWGCNA, nichnet, and lasso regression identified key genes associated with core transcription factors in selected macrophages, validated by external datasets.

Results

Six macrophage subgroups were identified in five post-transplant kidney biopsies. M1-like infiltrating macrophages, prevalent in ABMR, correlated with pathological injury severity. MIF acted as a primary intercellular signal in these macrophages. STAT1 regulated monocyte-to-M1-like phenotype transformation, impacting transplant prognosis via the IFNγ pathway. The prognostic models built on the upstream and downstream genes of STAT1 effectively predicted transplant survival. The TLR4-STAT1-PARP9 axis may regulate the pro-inflammatory phenotype of M1-like infiltrating macrophages, identifying PARP9 as a potential target for mitigating ABMR inflammation.

Conclusion

Our study delineates the macrophage landscape in ABMR post-kidney transplantation, underscoring the detrimental impact of M1-like infiltrating macrophages on ABMR pathology and prognosis.

Identification of estrogen response-associated STRA6+ granulosa cells within high-grade serous ovarian carcinoma by single-cell sequencing

Background

High-grade serous ovarian carcinoma (HGSOC) is a pathologic subtype of ovarian cancer (OC) with a more lethal prognosis. Extensive heterogeneity results in HGSOC being more susceptible to treatment resistance and adverse treatment effects. Revealing the heterogeneity involved is crucial.

Methods

We downloaded the single-cell RNA-seq (scRNA) data from GEO database and performed a scRNA analysis for cell landscape of HGSOC by using the Seurat package. The highly expressed genes were uploaded into the DAVID and KEGG database for enrichment analysis, and the AUCell package was used to calculate cancer-associated hallmark score. The SCENIC analysis was used for key regulons, the estrogen response enrichment scores in TCGA-OV RNA-seq dataset were calculated by using the GSVA package. Besides, the expression of STRA6 and IRF1 and the cell invasion and migration in si-STRA6 OC cells were detected by using the quantitative reverse transcription (qRT)-PCR method and Transwell assay respectively.

Results

We successfully constructed a single-cell atlas of HGSOC and delineated the heterogeneity of epithelial cells therein. There were five epithelial cell subpopulations, GLDC + Epithelial cells, PEG3+ leydig cells, STRA6+ granulosa cells, POLE2+ Epithelial cells, and AURKA + Epithelial cells. STRA6+ granulosa cells have the potential to promote tumor growth as well as the highest estrogen response early activity through the biological pathways analysis of highly expressed genes and estrogen response score of ssGSEA. We found that IRF1 and STRA6 expression was remarkably upregulated in the OC cancer cell line HEY. Silencing of STRA6 markedly decreased the invasion and migration ability of the OC cancer cell line HEY.

Conclusion

There is extreme heterogeneity of epithelial cells in HGSOC, and STRA6+ granulosa cells may be able to promote cancer progression. Our findings are benefit to the heterogeneity identification of HGSOC and develop targeted therapy strategy for HGSOC patients.

The crucial prognostic signaling pathways of pancreatic ductal adenocarcinoma were identified by single-cell and bulk RNA sequencing data

Pancreatic ductal adenocarcinoma (PDAC) is a malignant tumor with poor prognosis and high mortality. Although a large number of studies have explored its potential prognostic markers using traditional RNA sequencing (RNA-Seq) data, they have not achieved good prediction effect. In order to explore the possible prognostic signaling pathways leading to the difference in prognosis, we identified differentially expressed genes from one scRNA-seq cohort and four GEO cohorts, respectively. Then Cox and Lasso regression analysis showed that 12 genes were independent prognostic factors for PDAC. AUC and calibration curve analysis showed that the prognostic model had good discrimination and calibration. Compared with the low-risk group, the high-risk group had a higher proportion of gene mutations than the low-risk group. Immune infiltration analysis revealed differences in macrophages and monocytes between the two groups. Prognosis related genes were mainly distributed in fibroblasts, macrophages and type 2 ducts. The results of cell communication analysis showed that there was a strong communication between cancer-associated fibroblasts (CAF) and type 2 ductal cells, and collagen formation was the main interaction pathway.

Engineering Tumor Stroma Morphogenesis Using Dynamic Cell-Matrix Spheroid Assembly

The tumor microenvironment consists of resident tumor cells organized within a compositionally diverse, three-dimensional (3D) extracellular matrix (ECM) network that cannot be replicated in vitro using bottom-up synthesis. We report a new self-assembly system to engineer ECM-rich 3D MatriSpheres wherein tumor cells actively organize and concentrate microgram quantities of decellularized ECM dispersions which modulate cell phenotype. 3D colorectal cancer (CRC) MatriSpheres were created using decellularized small intestine submucosa (SIS) as an orthotopic ECM source that had greater proteomic homology to CRC tumor ECM than traditional ECM formulations such as Matrigel. SIS ECM was rapidly concentrated from its environment and assembled into ECM-rich 3D stroma-like regions by mouse and human CRC cell lines within 4-5 days via a mechanism that was rheologically distinct from bulk hydrogel formation. Both ECM organization and transcriptional regulation by 3D ECM cues affected programs of malignancy, lipid metabolism, and immunoregulation that corresponded with an in vivo MC38 tumor cell subpopulation identified via single cell RNA sequencing. This 3D modeling approach stimulates tumor specific tissue morphogenesis that incorporates the complexities of both cancer cell and ECM compartments in a scalable, spontaneous assembly process that may further facilitate precision medicine.

A division-of-labor mode contributes to the cardioprotective potential of mesenchymal stem/stromal cells in heart failure post myocardial infarction

Background

Treatment of heart failure post myocardial infarction (post-MI HF) with mesenchymal stem/stromal cells (MSCs) holds great promise. Nevertheless, 2-dimensional (2D) GMP-grade MSCs from different labs and donor sources have different therapeutic efficacy and still in a low yield. Therefore, it is crucial to increase the production and find novel ways to assess the therapeutic efficacy of MSCs.

Materials and methods

hUC-MSCs were cultured in 3-dimensional (3D) expansion system for obtaining enough cells for clinical use, named as 3D MSCs. A post-MI HF mouse model was employed to conduct in vivo and in vitro experiments. Single-cell and bulk RNA-seq analyses were performed on 3D MSCs. A total of 125 combination algorithms were leveraged to screen for core ligand genes. Shinyapp and shinycell workflows were used for deploying web-server.

Result

3D GMP-grade MSCs can significantly and stably reduce the extent of post-MI HF. To understand the stable potential cardioprotective mechanism, scRNA-seq revealed the heterogeneity and division-of-labor mode of 3D MSCs at the cellular level. Specifically, scissor phenotypic analysis identified a reported wound-healing CD142+ MSCs subpopulation that is also associated with cardiac protection ability and CD142- MSCs that is in proliferative state, contributing to the cardioprotective function and self-renewal, respectively. Differential expression analysis was conducted on CD142+ MSCs and CD142- MSCs and the differentially expressed ligand-related model was achieved by employing 125 combination algorithms. The present study developed a machine learning predictive model based on 13 ligands. Further analysis using CellChat demonstrated that CD142+ MSCs have a stronger secretion capacity compared to CD142- MSCs and Flow cytometry sorting of the CD142+ MSCs and qRT-PCR validation confirmed the significant upregulation of these 13 ligand factors in CD142+ MSCs.

Conclusion

Clinical GMP-grade 3D MSCs could serve as a stable cardioprotective cell product. Using scissor analysis on scRNA-seq data, we have clarified the potential functional and proliferative subpopulation, which cooperatively contributed to self-renewal and functional maintenance for 3D MSCs, named as "division of labor" mode of MSCs. Moreover, a ligand model was robustly developed for predicting the secretory efficacy of MSCs. A user-friendly web-server and a predictive model were constructed and available (https://wangxc.shinyapps.io/3D_MSCs/).

Single-cell transcriptomic sequencing identifies subcutaneous patient-derived xenograft recapitulated medulloblastoma

BACKGROUND

Medulloblastoma (MB) is one of the most common malignant brain tumors that mainly affect children. Various approaches have been used to model MB to facilitate investigating tumorigenesis. This study aims to compare the recapitulation of MB between subcutaneous patient-derived xenograft (sPDX), intracranial patient-derived xenograft (iPDX), and genetically engineered mouse models (GEMM) at the single-cell level.

METHODS

We obtained primary human sonic hedgehog (SHH) and group 3 (G3) MB samples from six patients. For each patient specimen, we developed two sPDX and iPDX models, respectively. Three Patch+/- GEMM models were also included for sequencing. Single-cell RNA sequencing was performed to compare gene expression profiles, cellular composition, and functional pathway enrichment. Bulk RNA-seq deconvolution was performed to compare cellular composition across models and human samples.

RESULTS

Our results showed that the sPDX tumor model demonstrated the highest correlation to the overall transcriptomic profiles of primary human tumors at the single-cell level within the SHH and G3 subgroups, followed by the GEMM model and iPDX. The GEMM tumor model was able to recapitulate all subpopulations of tumor microenvironment (TME) cells that can be clustered in human SHH tumors, including a higher proportion of tumor-associated astrocytes and immune cells, and an additional cluster of vascular endothelia when compared to human SHH tumors.

CONCLUSIONS

This study was the first to compare experimental models for MB at the single-cell level, providing value insights into model selection for different research purposes. sPDX and iPDX are suitable for drug testing and personalized therapy screenings, whereas GEMM models are valuable for investigating the interaction between tumor and TME cells.

Discovering Root Causal Genes with High Throughput Perturbations

Root causal gene expression levels - or root causal genes for short - correspond to the initial changes to gene expression that generate patient symptoms as a downstream effect. Identifying root causal genes is critical towards developing treatments that modify disease near its onset, but no existing algorithms attempt to identify root causal genes from data. RNA-sequencing (RNA-seq) data introduces challenges such as measurement error, high dimensionality and non-linearity that compromise accurate estimation of root causal effects even with state-of-the-art approaches. We therefore instead leverage Perturb-seq, or high throughput perturbations with single cell RNA-seq readout, to learn the causal order between the genes. We then transfer the causal order to bulk RNA-seq and identify root causal genes specific to a given patient for the first time using a novel statistic. Experiments demonstrate large improvements in performance. Applications to macular degeneration and multiple sclerosis also reveal root causal genes that lie on known pathogenic pathways, delineate patient subgroups and implicate a newly defined omnigenic root causal model.

Characterizing neuroinflammation and identifying prenatal diagnostic markers for neural tube defects through integrated multi-omics analysis

BACKGROUND

Neural Tube Defects (NTDs) are congenital malformations of the central nervous system resulting from the incomplete closure of the neural tube during early embryonic development. Neuroinflammation refers to the inflammatory response in the nervous system, typically resulting from damage to neural tissue. Immune-related processes have been identified in NTDs, however, the detailed relationship and underlying mechanisms between neuroinflammation and NTDs remain largely unclear. In this study, we utilized integrated multi-omics analysis to explore the role of neuroinflammation in NTDs and identify potential prenatal diagnostic markers using a murine model.

METHODS

Nine public datasets from Gene Expression Omnibus (GEO) and ArrayExpress were mined using integrated multi-omics analysis to characterize the molecular landscape associated with neuroinflammation in NTDs. Special attention was given to the involvement of macrophages in neuroinflammation within amniotic fluid, as well as the dynamics of macrophage polarization and their interactions with neural cells at single-cell resolution. We also used qPCR assay to validate the key TFs and candidate prenatal diagnostic genes identified through the integrated analysis in a retinoic acid-induced NTDs mouse model.

RESULTS

Our analysis indicated that neuroinflammation is a critical pathological feature of NTDs, regulated both transcriptionally and epigenetically within central nervous system tissues. Key alterations in gene expression and pathways highlighted the crucial role of STATs molecules in the JAK-STAT signaling pathway in regulating NTDs-associated neuroinflammation. Furthermore, single-cell resolution analysis revealed significant polarization of macrophages and their interaction with neural cells in amniotic fluid, underscoring their central role in mediating neuroinflammation associated with NTDs. Finally, we identified a set of six potential prenatal diagnostic genes, including FABP7, CRMP1, SCG3, SLC16A10, RNASE6 and RNASE1, which were subsequently validated in a murine NTDs model, indicating their promise as prospective markers for prenatal diagnosis of NTDs.

CONCLUSIONS

Our study emphasizes the pivotal role of neuroinflammation in the progression of NTDs and underlines the potential of specific inflammatory and neural markers as novel prenatal diagnostic tools. These findings provide important clues for further understanding the underlying mechanisms between neuroinflammation and NTDs, and offer valuable insights for the future development of prenatal diagnostics.

Molecular classification reveals the sensitivity of lung adenocarcinoma to radiotherapy and immunotherapy: multi-omics clustering based on similarity network fusion

BACKGROUND

Due to individual differences in tumors and immune systems, the response rate to immunotherapy is low in lung adenocarcinoma (LUAD) patients. Combinations with other therapeutic strategies improve the efficacy of immunotherapy in LUAD patients. Although radioimmunotherapy has been demonstrated to effectively suppress tumors, the underlying mechanisms still need to be investigated.

METHODS

Total RNA from LUAD cells was sequenced before and after radiotherapy to identify differentially expressed radiation-associated genes. The similarity network fusion (SNF) algorithm was applied for molecular classification based on radiation-related genes, immune-related genes, methylation data, and somatic mutation data. The changes in gene expression, prognosis, immune cell infiltration, radiosensitivity, chemosensitivity, and sensitivity to immunotherapy were assessed for each subtype.

RESULTS

We used the SNF algorithm and multi-omics data to divide TCGA-LUAD patients into three subtypes. Patients with the CS3 subtype had the best prognosis, while those with the CS1 and CS2 subtypes had poorer prognoses. Among the strains tested, CS2 exhibited the most elevated immune cell infiltration and expression of immune checkpoint genes, while CS1 exhibited the least. Patients in the CS2 subgroup were more likely to respond to PD-1 immunotherapy. The CS2 patients were most sensitive to docetaxel and cisplatin, while the CS1 patients were most sensitive to paclitaxel. Experimental validation of signature genes in the CS2 subtype showed that inhibiting the expression of RHCG and TRPA1 could enhance the sensitivity of lung cancer cells to radiation.

CONCLUSIONS

In summary, this study identified a risk classifier based on multi-omics data that can guide treatment selection for LUAD patients.

Single-cell transcriptome analysis reveals the association between histone lactylation and cisplatin resistance in bladder cancer

Patients with bladder cancer (BCa) frequently acquires resistance to platinum-based chemotherapy, particularly cisplatin. This study centered on the mechanism of cisplatin resistance in BCa and highlighted the pivotal role of lactylation in driving this phenomenon. Utilizing single-cell RNA sequencing, we delineated the single-cell landscape of Bca, pinpointing a distinctive subset of BCa cells that exhibit marked resistance to cisplatin with association with glycolysis metabolism. Notably, we observed that H3 lysine 18 lactylation (H3K18la) plays a crucial role in activating the transcription of target genes by enriching in their promoter regions. Targeted inhibition of H3K18la effectively restored cisplatin sensitivity in these cisplatin-resistant epithelial cells. Furthermore, H3K18la-driven key transcription factors YBX1 and YY1 promote cisplatin resistance in BCa. These findings enhance our understanding of the mechanisms underlying cisplatin resistance, offering valuable insights for identifying novel intervention targets to overcome drug resistance in Bca.