Single-cell Transcriptome Analysis Research Articles

Prediction of clinical prognosis and drug sensitivity in hepatocellular carcinoma through the combination of multiple cell death pathways.

Hepatocellular carcinoma (HCC) is the sixth most common malignant tumor, highlighting a significant need for reliable predictive models to assess clinical prognosis, disease progression, and drug sensitivity. Recent studies have highlighted the critical role of various programmed cell death pathways, including apoptosis, necroptosis, pyroptosis, ferroptosis, cuproptosis, entotic cell death, NETotic cell death, parthanatos, lysosome-dependent cell death, autophagy-dependent cell death, alkaliptosis, oxeiptosis, and disulfidptosis, in tumor development. Therefore, by investigating these pathways, we aimed to develop a predictive model for HCC prognosis and drug sensitivity. We analyzed transcriptome, single-cell transcriptome, genomic, and clinical information using data from the TCGA-LIHC, GSE14520, GSE45436, and GSE166635 datasets. Machine learning algorithms were used to establish a cell death index (CDI) with seven gene signatures, which was validated across three independent datasets, showing that high CDI correlates with poorer prognosis. Unsupervised clustering revealed three molecular subtypes of HCC with distinct biological processes. Furthermore, a nomogram integrating CDI and clinical information demonstrated good predictive performance. CDI was associated with immune checkpoint genes and tumor microenvironment components using single-cell transcriptome analysis. Drug sensitivity analysis indicated that patients with high CDI may be resistant to oxaliplatin and cisplatin but sensitive to axitinib and sorafenib. In summary, our model offers a precise prediction of clinical outcomes and drug sensitivity for patients with HCC, providing valuable insights for personalized treatment strategies.

Cutaneous T cell lymphoma atlas reveals malignant TH2 cells supported by a B cell-rich tumor microenvironment

AbstractCutaneous T cell lymphoma (CTCL) is a potentially fatal clonal malignancy of T cells primarily affecting the skin. The most common form of CTCL, mycosis fungoides, can be difficult to diagnose, resulting in treatment delay. We performed single-cell and spatial transcriptomics analysis of skin from patients with mycosis fungoides-type CTCL and an integrated comparative analysis with human skin cell atlas datasets from healthy and inflamed skin. We revealed the co-optation of T helper 2 (TH2) cell-immune gene programs by malignant CTCL cells and modeling of the tumor microenvironment to support their survival. We identified MHC-II+ fibroblasts and dendritic cells that can maintain TH2 cell-like tumor cells. CTCL tumor cells are spatially associated with B cells, forming tertiary lymphoid structure-like aggregates. Finally, we validated the enrichment of B cells in CTCL and its association with disease progression across three independent patient cohorts. Our findings provide diagnostic aids, potential biomarkers for disease staging and therapeutic strategies for CTCL.

Abstract PR017: Anti-metastatic immunotherapy discovery using ex vivo lung tissue cultures and high-throughput single-cell chemical transcriptomics

Abstract Introduction: Tumors systemically remodel the immune system during metastasis. Developing anti-metastatic immunotherapies that target these tumor-immune interactions to make the metastatic niche hostile colonization would represent a paradigm shift in cancer treatment. However, traditional high-throughput screening (HTS) platforms that use simple read-outs and contrived in vitro systems are ill-suited to identify such therapies. In contrast, HTS platforms that accurately model the metastatic niche and use high-dimensional read-outs such as multiplexed single-cell RNA-sequencing (scRNA-seq) can accurately profile nuanced perturbation responses in individual immune cell types and therefore have the potential to discover anti-metastatic immunotherapies. In this study, we describe the development of an ex vivo lung tissue culture HTS platform that preserves immune gene expression signatures observed in the in vivo metastatic niche and is amenable to single-cell chemical transcriptomic analysis using MULTI- seq (McGinnis et al., Nature Methods, 2019). Inspired by our recent description of myeloid cell TLR-NFκB inflammation during breast cancer lung metastasis (McGinnis et al., Cancer Cell, 2024), we then performed the first anti-metastatic immunotherapy drug screen and identified TLR-NFκB inhibitors that effectively operate on all desired myeloid cell types in the lung metastatic niche. Methods and Results: Building upon existing techniques for culturing precision-cut lung slices (Wu et al., Journal of Experimental Medicine, 2024) and patient- derived tumor fragments (Voabil et al., Nature Medicine, 2021), we established and optimized an ex vivo lung tissue culture system that enables 48-hour culturing of lung slices isolated from 4T1 tumor-bearing mice. scRNA-seq profiling revealed successful capture of all relevant immune cell types and retention of metastasis-associated gene expression programs (e.g., myeloid TLR-NFκB inflammation and neutrophil degranulation) following culture. We then scaled our platform to address how myeloid cells in the lung metastatic niche respond to perturbations targeting every component of the canonical TLR-NFκB signaling cascade. Cell-type-specific perturbation modeling and quantification of TLR-NFκB signaling inhibition identified compounds that optimally perturb all intended myeloid cell types (e.g., tissue-resident macrophages, neutrophils, and monocytes) and represent prioritized candidates for future in vivo validation studies to assess impact on metastatic disease progression. Conclusions: In this study, we describe the development and application of a novel HTS platform that couples ex vivo lung tissue cultures with single-cell chemical transcriptomics to identify anti-metastatic immunotherapy candidates. We demonstrate how our platform can identify TLR-NFκB signaling inhibitors that optimally operate in the lung metastatic niche, paving the way for future interrogation of additional metastasis-associated immune cell signaling pathways in different disease and tissue backgrounds. Citation Format: Chris McGinnis, Winnie Yao, Ansuman Satpathy. Anti-metastatic immunotherapy discovery using ex vivo lung tissue cultures and high-throughput single-cell chemical transcriptomics [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Tumor-body Interactions: The Roles of Micro- and Macroenvironment in Cancer; 2024 Nov 17-20; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(22_Suppl):Abstract nr PR017.

Abstract A002: Dissecting the brain metastatic microenvironment to uncover immune and molecular correlates of response to immunotherapy

Abstract Despite promising results with immune checkpoint blockade (ICB) therapy, brain metastases (BrM) remain a deadly complication for cancer patients. Understanding the highly specialized brain metastatic tumor microenvironment (BrTME) is crucial for identifying determinants of response, however, progress has been hindered by limited access to patient samples and scarcity of relevant preclinical models. Here, we developed two mouse melanoma BrM models by intracardiac injection of cells derived from our previously established M4-B2905 melanoma mouse cell line. We characterized their mutational landscape and performed single-cell phenotypic and transcriptomic analysis, demonstrating that these models recapitulate the cellular and molecular features of human melanoma BrMs. Comparative and interactome analysis of our models revealed key factors contributing to ICB response and resistance. We found that the responsive model (BR1) was characterized by tumor cells that polarize microglia toward reactive states via inflammatory programs and response to IFN signaling. These microglia express high levels of antigen presentation and immunostimulatory molecules, triggering T cell recruitment and activation and promoting ICB therapy response. In contrast, in the resistant model (BR3), the tumor cells express neurological molecular signatures and ligands that sustain microglia homeostasis. This results in poor T cell infiltration predominantly composed of naïve cells, evading immune reaction and promoting ICB resistance. While we showed that systemic ICB therapy induced BrTME changes in both models, differences between responder and resistant models were maintained, emphasizing the importance of tumor cell-intrinsic programs in dictating the BrTME. We validated the translational relevance of our findings and demonstrated that BR1 signatures positively correlate with T cell infiltration and are associated with improved patient outcomes, whereas BR3 signatures negatively correlate with T cells and are found in patients with worse prognosis. Here we address an important challenge in the field by providing clinically relevant BrM models that recapitulate both the cellular and molecular features of human disease and the variability of ICB response seen in patients. We further reveal mechanistic insights into BrM ICB response and identify potential therapeutic targets to modulate the BrTME. Citation Format: Amelie Daugherty-Lopes, Eva Perez-Guijarro, Vishaka Gopalan, Jessica Rappaport, Quanyi Chen, April Huang, Khiem C. Lam, Sung Chin, Jessica Ebersole, Emily Wu, Gabriel Needle, Isabella Church, George Kyriakopoulos, Shaojun Xie, Yongmei Zaho, Charli Gruen, Antonella Sassano, Romina E Araya, Andres Thorkelsson, Cari Smith, Maxwell P. Lee, Sridhar Hannenhalli, Chi-Ping Day, Glenn Merlino, Romina S. Goldszmid, Shaojun Xie. Dissecting the brain metastatic microenvironment to uncover immune and molecular correlates of response to immunotherapy [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Tumor-body Interactions: The Roles of Micro- and Macroenvironment in Cancer; 2024 Nov 17-20; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(22_Suppl):Abstract nr A002.

Abstract A018: Endothelial pyrimidine synthesis deficiency promotes tumor growth

Abstract The tumor microenvironment plays a crucial role in cancer progression and treatment response. Systemically administered metabolic cancer therapies target not only malignant but also stromal cells, including immune cells - the only stromal population previously addressed in this context. However, how such treatments impact other non-malignant cell types in tumors remains poorly understood. Here we show that inhibition of de novo pyrimidine synthesis in endothelial cells accelerates tumor growth and alters tumor immune repertoire. We found that whole-body de novo pyrimidine synthesis deficiency in mice, caused by the inducible whole-body ablation of DHODH, accelerates the growth rate of orthotopic lung tumors. Single-cell transcriptomic analysis of tumor-bearing lungs revealed that DHODH deficiency in the stroma impacts multiple cell populations, including immune and, surprisingly, endothelial cells. To explore the endothelial-specific effects, we generated a mouse model with inducible DHODH deficiency restricted to endothelial cells. The endothelium-specific model recapitulates the accelerated lung tumor growth observed in the whole-body DHODH deficiency model. Single-cell transcriptomics analysis of tumor-bearing lungs in the endothelium-specific model pointed to changes in the immune repertoire, particularly an enrichment of monocytes. We confirmed these results on the protein level using spectral flow cytometry, and we are currently in the process of uncovering the mechanism by which endothelial deficiency of pyrimidine synthesis affects the immune landscape of tumors. Inhibitors of pyrimidine de novo synthesis have been tested in clinical trials but failed due to their low efficiency. Our findings indicate that systemic treatment targeting pyrimidine synthesis may be hampered by its pro-tumorigenic effects in the endothelium, highlighting the unexpected role of endothelial metabolism in this context. Citation Format: Petra Hyrossova, Isidora Milisav, Silvia Novais, Mirko Milosevic, Jakub Rohlena, Katerina Rohlenova. Endothelial pyrimidine synthesis deficiency promotes tumor growth [abstract]. In: Proceedings of the AACR Special Conference in Cancer Research: Tumor-body Interactions: The Roles of Micro- and Macroenvironment in Cancer; 2024 Nov 17-20; Boston, MA. Philadelphia (PA): AACR; Cancer Res 2024;84(22_Suppl):Abstract nr A018.

Single-cell data revealed the regulatory mechanism of TNK cell heterogeneity in liver metastasis from gastric cancer.

The present work set out to classify cell subpopulations related to liver metastasis from gastric cancer (GC) and the mechanisms of their interactions with other immune cell subpopulations. GC is characterized by a high degree of heterogeneity and liver metastasis. Exploring the mechanism of liver metastasis of GC from the perspective of heterogeneity of the tumor microenvironment (TME) might help improve the efficacy of GC treatment. Based on the cellular subpopulation characteristics of GC with liver metastasis, the regulatory mechanisms contributing to GC progression were analyzed, with special focuses on the roles of signaling pathways, transcription factors (TFs) and ligand-receptor pairs. The GSE163558 dataset was downloaded from the Gene Expression Omnibus (GEO) database to collect single-cell transcriptomic data of GC patients and their metastasis groups for cell clustering and relevant analyses. Differentially expressed genes (DEGs) in the GC and GC liver metastasis groups were screened and subjected to Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis. SCENIC analysis was used to mine TFs that affected cellular subpopulations during liver metastasis from GC. The relative expression levels of TFs in GC were determined using qRT-PCR. Transwell and wound healing assays were utilized to verify the regulation of the TFs on the migration and invasion of GC cells. Interaction network between the cellular subpopulations was developed applying CellChat. Single-cell clustering was performed to group six major cell subpopulations, namely, Myeloid cells, B cells, Mast cells, Epithelial cells, Fibroblasts, and TNK cells, among which the number of TNK cells was significantly increased in the GC liver metastasis group. Differentially enriched pathways of TNK cells between GC and GC liver metastasis groups mainly included IL-17 and Pi3k-Akt signaling pathways. TNK cell subsets could be further categorized into CD8 T cells, Exhausted T cells, NK cells, NKT cells, and Treg cells, with the GC liver metastasis group showing significantly more CD8 T cells and NKT cells. FOS and JUNB were the TFs of TNK cell marker genes that contributed to liver metastasis from GC and the invasion and migration of GC cell lines. Significant differences in immune cell communication ligand-receptor pairs existed between the GC and GC liver metastasis groups. This study revealed the critical role of TNK cell subsets in GC with liver metastasis applying single-cell transcriptomics analysis. The findings provided an important theoretical basis for developing novel therapies to inhibit liver metastasis from GC.

Single-cell transcriptomic analysis reveals the developmental trajectory and transcriptional regulatory networks of quinoa salt bladders.

Salt bladders, specialized structures on the surface of quinoa leaves, secrete Na+ to mitigate the effects of the plant from abiotic stresses, particularly salt exposure. Understanding the development of these structures is crucial for elucidating quinoa's salt tolerance mechanisms. In this study, we employed transmission electron microscopy to detail cellular differentiation across the developmental stages of quinoa salt bladders. To further explore the developmental trajectory and underlying molecular mechanisms, we conducted single-cell RNA sequencing on quinoa protoplasts derived from young leaves. This allowed us to construct a cellular atlas, identifying 13 distinct cell clusters. Through pseudotime analysis, we mapped the developmental pathways of salt bladders and identified regulatory factors involved in cell fate decisions. GO and KEGG enrichment analyses, as well as experimental results, revealed the impacts of salt stress and the deprivation of sulfur and nitrogen on the development of quinoa salt bladders. Analysis of the transcription factor interaction network in pre-stalk cells (pre-SC), stalk cells (SC), and epidermal bladder cells (EBCs) indicated that TCP5, YAB5, NAC078, SCL8, GT-3B, and T1P17.40 play crucial roles in EBC development. Based on our findings, we developed an informative model elucidating salt bladder formation. This study provides a vital resource for mapping quinoa leaf cells and contributes to our understanding of its salt tolerance mechanisms.

Single-cell and bulk transcriptome analysis reveals tumor cell heterogeneity and underlying molecular program in uveal melanoma.

Uveal melanoma (UM) is a rare and deadly eye cancer with high metastatic potential. Despite the predominance of malignant cells within the tumor microenvironment, the heterogeneity and underlying molecular features remain to be fully characterized. We analyzed single-cell transcriptomic profiling of 37,660 malignant cells from 17 UM tumors. A consensusnon-negative factorization algorithm was used to decipher transcriptional programs underlying tumor cell-intrinsic heterogeneity. Tumor-infiltrated immune cells were computationally estimated from bulk transcriptomes and bioinformatics methods. A gene signature was derived for subtyping and prognostic stratification and validated in multicenter cohorts. ScRNA-seq analysis revealed the existence of diverse subpopulations and transcriptional variability among malignant cells within and between tumors. Furthermore, we observed that the heterogeneity of malignant cell states and compositions correlated with genomic, immunological, and clinical characteristics. By identifying gene expression programs associated with malignant cell heterogeneity at the single cell level, UM samples were classified into two distinct intra-tumoral subtypes (ITMHlo and ITMHhi) with different prognoses and immune microenvironments. Finally, a machine learning-derived 9-gene signature was developed to translate single-cell information into bulk tissue transcriptomes for patient stratification and was validated in multicenter cohorts. Our study provides insight into understanding the intra-tumoral heterogeneity of UM and its potential impact on patient stratification.

TMIC-83. ABERRANT EGFR ACTIVATION SUPPORTS ENHANCED TUMOR-MYELOID CROSSTALK AND IMPAIRS EFFECTOR T CELL INFILTRATION AND FUNCTION IN GLIOBLASTOMA

Abstract Glioblastoma (GBM) is known to display a highly suppressive tumor immune microenvironment (TIME) dominated by myeloid cells with low T cell infiltration. Aberrant epidermal growth factor receptor (EGFR) activation is observed in ~50% of GBM patients via amplification and the constitutively active EGFRvIII variant, however the impact of EGFR activation on the brain TIME and T cell-mediated anti-tumor immunity remains poorly understood. Single cell transcriptomic and receptor-ligand interaction analyses of paired newly diagnosed GBM tumor and intratumoral CD45+ immune cells revealed increased crosstalk between EGFR activated tumors and EGFR ligand-expressing tumor-associated myeloid cells. Therefore, a novel orthotopic murine EGFRvIII-expressing glioma model wherein the mEGFRvIII gene is under a tetracycline-off system (MADR-mEGFRvIII) was leveraged to identify how EGFR signaling impacts overall TIME composition and the phenotype of intratumoral and systemic T cells in vivo. Temporal flow cytometric immunophenotyping following tumor or mock (saline) implantation identified no significant shift in T cell infiltration within the brain or cervical lymph nodes (cLNs). In contrast, one week of doxycycline-mediated genetic EGFR ablation conferred a significant survival benefit and was associated with enhanced intratumoral T cell infiltration as well as diminished tumor-associated F4/80+ myeloid cells compared to vehicle treated animals. Further characterization revealed genetic EGFR ablation enhanced recruitment of CXCR3+, TNFα+, and TIGIT+ CD4 T cells to the tumor, while the prevalence of these populations contracted in CD4 and CD8 T cells within the draining cLNs. In vitro, pharmacologic EGFR inhibition enhanced antigen-specific Pmel-1 T cell killing of gp100-expressing MADR-mEGFRvIII tumor cells and was associated with increased IFNγ secretion by ELISA. Collectively, our work demonstrates that aberrant EGFR signaling promotes an immunosuppressive brain TIME with reduced T cell infiltration potentially via tumor-myeloid signaling. Targeting this oncogenic axis may therefore reprogram the TIME and directly enable enhanced trafficking of effector T cells to the tumor.

EPCO-19. MULTIOMIC PROFILING OF EPIGENETIC HETEROGENEITY IN GLIOBLASTOMA AT SINGLE-CELL RESOLUTION USING MULTIOMIC SINGLE-CELL CUT&TAG

Abstract INTRODUCTION Recent evidence suggests epigenetic heterogeneity may be the primary driver of cell phenotype and treatment resistance in glioblastoma (GBM). Epigenetic heterogeneity is defined by alterations of chromatin histone modifications with activating marks such as H3K27 acetylation and repressive marks such as H3K27 trimethylation. Combinations of these histone marks govern the activity of functional genomic elements such as enhancers or promoters. Single-cell (“sc-”) technologies have only recently been adapted to the study of epigenetic cell states and significant technical challenges have limited their widespread adoption, especially in complex tissues. Here, we report the development of an optimized method for performing tissue processing and simultaneous single-cell epigenetic and transcriptomic analyses in flash frozen GBM samples. METHODS Building on our prior work, we developed a new protocol, multiomic single-cell cleavage-under-targets and tagmentation, or “mscCUT&TAG”. We optimized nuclei isolation, cleanup and permeabilization from flash frozen GBM and reaction conditions for the scCUT&TAG portion of the assay. Critically, these changes maintained high-quality, intact nuclei in suspension while limiting well-described technical issues such as nuclei clumping and sample loss. Furthermore, we incorporated the ability to multiplex both patient samples and histone marks. RESULTS We applied 10X scRNA-seq, 10X Multiome (RNA+ATAC-seq), and mscCUT&TAG for five histone marks to the same GBM sample, enabling an integrated evaluation of the tumor’s transcriptome, chromatin accessibility and epigenetic regulation at single-cell resolution. mscCUT&TAG enabled the first assessment of single-cell genome segmentation in GBM. We integrated across modalities to identify specific tumor cell populations, cell-type specific coordinated changes in epigenetics and transcriptome, and define functional regulatory elements. CONCLUSIONS We developed a novel multiomic protocol, mscCUT&TAG, that allows for a comprehensive analysis of epigenetic heterogeneity in GBM. We anticipate this protocol will enhance our understanding of fundamental biology and potential therapeutic avenues in GBM.

RBIO-02. DEVELOPMENT OF FLASH-CAR RADIOIMMUNOTHERAPY FOR PEDIATRIC BRAIN TUMORS

Abstract Pediatric brain tumors including medulloblastoma and high-grade glioma are generally refractory to T cell-based immunotherapy, largely due to an immune-hostile microenvironment infiltrated extensively with immunosuppressive macrophages. Ultra-high dose rate (FLASH) radiotherapy holds promise for treating solid tumors, given the potential lower toxicity in normal tissues and possible favorable impact on tumor immunity. Using a genetically engineered mouse model of medulloblastoma and glioma, we show that FLASH radiation stimulates pro-inflammatory polarization in tumor macrophages. Single-cell transcriptome analysis shows that FLASH proton beam radiation skews macrophages towards proinflammatory phenotypes and increases T cell infiltration in the medulloblastoma model. Further bulk transcriptome analyses reveal that FLASH radiation reduces peroxisome proliferator-activated receptor (PPAR)-γ and arginase-1 expression and inhibits immunosuppressive macrophage polarization under stimulus-inducible conditions. Mechanistically, FLASH radiation abrogates lipid oxidase expression and oxidized low-density lipid (oxLDL) generation to reduce PPARγ activity, while standard radiation induces reactive oxygen species (ROS)-dependent PPARγ activation in macrophages. Notably, FLASH radiotherapy improves infiltration and activation of chimeric antigen receptor (CAR) T cells and sensitizes tumors to GD2 CAR T immunotherapy in the autochthonous medulloblastoma and syngeneic glioma models. These findings suggest that FLASH radiotherapy may reprogram macrophage lipid metabolism to reverse tumor immunosuppression and that combination FLASH-CAR radioimmunotherapy may offer exciting opportunities for treatment of pediatric brain tumor.

IMMU-22. TWIST1-MEDIATED ENDOTHELIAL PLASTICITY DRIVES GLIOBLASTOMA RESISTANCE TO CAR T IMMUNOTHERAPY

Abstract Glioblastoma (GBM) is the most common and most aggressive malignant primary brain tumor in adults. Among the most lethal human malignancies, GBM is highly resistant to cytotoxic treatments and molecularly targeted therapies, and also generally refractory to T cell–based immunotherapies, largely due to their immunologically cold nature, i.e., characterized by a paucity of tumor T cell infiltrates that result from the immunosuppressive microenvironment. Here, we identify a distinct mechanism driven by spatiotemporal interaction of leukocytes with tumor stromal cells, namely, mesenchymal-like population of endothelial cells (ECs) form an immunosuppressive vascular niche that regulates macrophage polarization and induces GBM resistance to CAR T immunotherapy. Our single-cell transcriptome analyses of human and mouse GBM tumors identify a subpopulation of ECs acquire robust immunosuppressive phonotypes to drive alternative macrophage polarization. These ECs are spatially localized proximately to immunosuppressive macrophages and selectively express Twist1. Furthermore, we reveal a Twist1/SATB1-mediated sequential transcriptional activation mechanism, through which Twist1+ tumor ECs produce osteopontin to locally promote immunosuppressive macrophage phenotypes. Genetic or pharmacological ablation of Twist1 reverses macrophage-mediated immunosuppression and enhances T cell infiltration and activation, leading to reduced GBM growth and extended mouse survival. Importantly, pharmacological inhibition and lipid nanoparticle-mediated targeted inactivation of Twist1 overcome GBM resistance to Egfrviii CAR T cell immunotherapy. Thus, these findings uncover a spatially restricted mechanism controlling vasculature-mediated tumor immunity and suggest that targeting endothelial Twist1 may offer exciting opportunities for GBM immunotherapy.

Single-cell transcriptomic and cross-species comparison analyses reveal distinct molecular changes of porcine testes during puberty

The pig is an important model for studying human diseases and is also a significant livestock species, yet its testicular development remains underexplored. Here, we employ single-cell RNA sequencing to characterize the transcriptomic landscapes across multiple developmental stages of Bama pig testes from fetal stage through infancy, puberty to adulthood, and made comparisons with those of humans and mice. We reveal an exceptionally early onset of porcine meiosis shortly after birth, and identify a distinct subtype of porcine spermatogonia resembling transcriptome state 0 spermatogonial stem cells identified in humans, which were previously thought to be primate specific. We also discover the persistent presence of proliferating progenitors for myoid cells in postnatal testes. The regulatory roles of Leydig cell steroidogenesis and estrogen synthesis in supporting cell lineages are also explored, including the potential impact of estrogen on Sertoli cell maturation and spermatogenesis. Overall, this study offers valuable insights into porcine testicular development, paving the way for future research in reproductive biology, advancements in agricultural breeding, and potential applications in translational medicine.

Fibroblast Heterogeneity in Hepatocellular Carcinoma and Identification of Prognostic Markers Based on Single-cell Transcriptome Analysis.

HCC is a malignant tumor with high morbidity and mortality. Fibroblasts play a key role in the tumor microenvironment (TME). However, the transcriptional regulatory mechanisms of fibroblasts remained unclear in HCC. The aim of this study was to explore the complex role of fibroblasts in hepatocellular carcinoma (HCC) and to reveal their transcriptional regulatory mechanisms. The goal of this study was to discover potential prognostic markers for HCC by analyzing the genetic variations and differentiation process of fibroblasts. Single-cell transcriptome data from the non-tumor liver site and primary tumor site of HCC were acquired from GSE149614, processed, and clustered using the Seurat pipeline. The inferCNV algorithm was applied to infer copy number variations (CNVs) in fibroblasts. Subsequently, the mechanism underlying the interaction between fibroblasts and other cells in the TME of HCC was analyzed using CellChat software. The trajectory of cellular differentiation of fibroblasts from normal state to malignant state was examined using Monocle 2. SCENIC analysis was performed to identify key transcription factors (TFs) in fibroblasts and assess their correlation with HCC prognosis. Finally, qRT-PCR and Transwell assays were carried out to analyze the mRNA expression and cell metastasis. We identified a total of nine different cell types (B cells, cycling cells, endothelial cells, epithelial cells, fibroblasts, hepatocytes, macrophages, plasma cells, and T cells) based on the single-cell transcriptomic data of HCC. Among them, fibroblasts were highly enriched at the primary tumor site, and their number increased with advanced stages. In addition, significant deletions were detected on chromosome 6p of fibroblasts, and genes in this region were remarkably enriched in pathways associated with antigen processing and presentation. Intercellular communication showed that epithelial cells regulated fibroblasts the most. The differentiation of fibroblasts was mainly accompanied by a transition from normal to malignant state. Importantly, CEBPD and FOSB, the TFs most associated with the putative timing of fibroblasts, were under-expressed in human hepatocytes and showed a significant correlation with HCC prognosis. Overexpressed CEBPD inhibited HCC cell migration and invasion. In conclusion, our study revealed that fibroblast recruitment and differentiation, as well as copy number loss at chromosome 6p, were associated with a higher degree of malignancy and immune dysfunction in HCC. The current discoveries provided new insights into the clinical treatment and diagnosis of HCC.

Unveiling the influence of tumor and immune signatures on immune checkpoint therapy in advanced lung cancer.

This study investigates the variability among patients with non-small cell lung cancer (NSCLC) in their responses to immune checkpoint inhibitors (ICIs). Recognizing that patients with advanced-stage NSCLC rarely qualify for surgical interventions, it becomes crucial to identify biomarkers that influence responses to ICI therapy. We conducted an analysis of single-cell transcriptomes from 33 lung cancer biopsy samples, with a particular focus on 14 core samples taken before the initiation of palliative ICI treatment. Our objective was to link tumor and immune cell profiles with patient responses to ICI. We discovered that ICI non-responders exhibited a higher presence of CD4+ regulatory T cells, resident memory T cells, and TH17 cells. This contrasts with the diverse activated CD8+ T cells found in responders. Furthermore, tumor cells in non-responders frequently showed heightened transcriptional activity in the NF-kB and STAT3 pathways, suggesting a potential inherent resistance to ICI therapy. Through the integration of immune cell profiles and tumor molecular signatures, we achieved an discriminative power (area under the curve [AUC]) exceeding 95% in identifying patient responses to ICI treatment. These results underscore the crucial importance of the interplay between tumor and immune microenvironment, including within metastatic sites, in affecting the effectiveness of ICIs in NSCLC.

Unveiling the influence of tumor and immune signatures on immune checkpoint therapy in advanced lung cancer

This study investigates the variability among patients with non-small cell lung cancer (NSCLC) in their responses to immune checkpoint inhibitors (ICIs). Recognizing that patients with advanced-stage NSCLC rarely qualify for surgical interventions, it becomes crucial to identify biomarkers that influence responses to ICI therapy. We conducted an analysis of single-cell transcriptomes from 33 lung cancer biopsy samples, with a particular focus on 14 core samples taken before the initiation of palliative ICI treatment. Our objective was to link tumor and immune cell profiles with patient responses to ICI. We discovered that ICI non-responders exhibited a higher presence of CD4+ regulatory T cells, resident memory T cells, and TH17 cells. This contrasts with the diverse activated CD8+ T cells found in responders. Furthermore, tumor cells in non-responders frequently showed heightened transcriptional activity in the NF-kB and STAT3 pathways, suggesting a potential inherent resistance to ICI therapy. Through the integration of immune cell profiles and tumor molecular signatures, we achieved an discriminative power (area under the curve [AUC]) exceeding 95% in identifying patient responses to ICI treatment. These results underscore the crucial importance of the interplay between tumor and immune microenvironment, including within metastatic sites, in affecting the effectiveness of ICIs in NSCLC.

The molecular chronology of mammary epithelial cell fate switching.

The adult mammary gland is maintained by lineage-restricted progenitor cells through pregnancy, lactation, involution, and menopause. Injury resolution, transplantation-associated mammary gland reconstitution, and tumorigenesis are unique exceptions, wherein mammary basal cells gain the ability to reprogram to a luminal state. Here, we leverage newly developed cell-identity reporter mouse strains, and time-resolved single-cell epigenetic and transcriptomic analyses to decipher the molecular programs underlying basal-to-luminal fate switching in vivo . We demonstrate that basal cells rapidly reprogram toward plastic cycling intermediates that appear to hijack molecular programs we find in bipotent fetal mammary stem cells and puberty-associatiated cap cells. Loss of basal-cell specifiers early in dedifferentiation coincides with activation of Notch and BMP, among others. Pharmacologic blockade of each pathway disrupts basal-to-luminal transdifferentiation. Our studies provide a comprehensive map and resource for understanding the coordinated molecular changes enabling terminally differentiated epithelial cells to transition between cell lineages and highlights the stunning rapidity by which epigenetic reprogramming can occur in response to disruption of tissue structure.

Molecular characterization of human HSPCs with different cell fates in vivo using single-cell transcriptome analysis and lentiviral barcoding technology.

Hematopoietic stem and progenitor cells (HSPCs) possess the potential to produce all types of blood cells throughout their lives. It is well recognized that HSPCs are heterogeneous, which is of great significance for their clinical applications and the treatment of diseases associated with HSPCs. This study presents a novel technology called Single-Cell transcriptome Analysis and Lentiviral Barcoding (SCALeBa) to investigate the molecular mechanisms underlying the heterogeneity of human HSPCs in vivo. The SCALeBa incorporates a transcribed barcoding library and algorithm to analyze the individual cell fates and their gene expression profiles simultaneously. Our findings using SCALeBa reveal that HSPCs subset with stronger stemness highly expressed MYL6B, ATP2A2, MYO19, MDN1, ING3, and so on. The high expression of COA3, RIF1, RAB14, and GOLGA4 may contribute to the pluripotent-lineage differentiation of HSPCs. Moreover, the roles of the representative genes revealed in this study regarding the stemness of HPSCs were confirmed with biological experiments. HSPCs expressing MRPL23 and RBM4 genes may contribute to differentiation bias into myeloid and lymphoid lineage, respectively. In addition, transcription factor (TF) characteristics of lymphoid and myeloid differentiation bias HSPCs subsets were identified and linked to previously identified genes. Furthermore, the stemness, pluripotency, and differentiation-bias genes identified with SCALeBa were verified in another independent HSPCs dataset. Finally, this study proposes using the SCALeBa-generated tracking trajectory to improve the accuracy of pseudo-time analysis results. In summary, our study provides valuable insights for understanding the heterogeneity of human HSPCs in vivo and introduces a novel technology, SCALeBa, which holds promise for broader applications. KEY POINTS: SCALeBa and its algorithm are developed to study the molecular mechanism underlying human HSPCs identity and function. The human HSPCs expressing MYL6B, MYO19, ATP2A2, MDN1, ING3, and PHF20 may have the capability for high stemness. The human HSPCs expressing COA3, RIF1, RAB14, and GOLGA4 may have the capability for pluripotent-lineage differentiation. The human HSPCs expressing MRPL23 and RBM4 genes may have the capability to differentiate into myeloid and lymphoid lineage respectively in vivo. The legitimacy of the identified genes with SCALeBa was validated using biological experiments and a public human HSPCs dataset. SCALeBa improves the accuracy of differentiation trajectories in monocle2-based pseudo-time analysis.

Hidden liver-joint axis: HBV infection causes rheumatoid arthritis via TRAFD1 with imbalance of HBV X protein and trans-ferulic acid

ABSTRACT Liver metabolites are involved in the progression of rheumatoid arthritis (RA), indicating a connection between the liver and joints. However, the impact and mechanism of Hepatitis B virus (HBV), a hepatotropic virus, on RA are still unclear. We investigated the correlation between HBV and RA using Mendelian randomization analysis. Single-cell transcriptome analysis was conducted to investigate changes in cell subtypes in synovial tissue of HBV-RA patients. Fibroblast-like synoviocytes (FLS) were used to create a cell model, and the transcriptome was examined to identify the key downstream molecules of FMT regulated by HBx. CIA model was constructed using HBV transgenic, HBx transgenic, and TRADF1 knockout mice to investigate the impact and mechanism of HBV on CIA. The results of our study revealed a significant positive correlation between HBV and RA. The functional studies identified a crucial role of fibroblast-myofibroblast transition (FMT) in the progression of RA. The results suggest that HBV-encoded HBx may promote FMT in RA by upregulating TRAFD1. Furthermore, trans-ferulic acid (TFA) was identified by screening for common metabolites in the liver, joints, and peripheral blood using the metabolome and WGCNA. Interestingly, we found that TFA ameliorated HBx-induced RA by suppressing TRAFD1 expression. Our study demonstrates that hidden liver-joint axis, an imbalance between TFA and HBx, plays a critical role in HBV-induced RA, which could be a potential strategy for preventing RA development.

From bioinformatics to anti-inflammation and immune regulation

Background: ACT001 is a potent anti-inflammatory small-molecule drug. However, the single cell and spatial molecular basis of pyroptosis and whether ACT001 exerts a therapeutic effect by preventing pyroptosis on acute lung injury (ALI) remains unclear. Methods: The bioinformatics approach was employed to identify single cell and spatial landscape of nucleotide-binding domains and leucine-rich repeat protein 3 (NLRP3)-dependent pyroptosis in lipopolysaccharide (LPS) and influenza virus-induced ALI. Molecular docking was performed to elucidate the relationship between ACT001 and NLRP3. LPS-induced ALI mice model was established. Histopathological analysis and bronchoalveolar lavage fluid collection were conducted to investigate the anti-inflammatory and protective effects. In vitro experiments were also performed on bone marrow–derived macrophages to explore the effect of ACT001 on the balance of mitochondrial fusion and fission protein. Results: Single cell transcriptomic and spatial transcriptomic analysis predicted that NLRP3-dependent pyroptosis significantly correlated with the development of ALI both in single cell and spatial distribution. Molecular docking provided a stable and reliable docking between ACT001 and NLRP3. ACT001 improved the 7-day survival of mice by approximately 50% over the loading dose of LPS-induced ALI. ACT001 (5 uM) attenuated the disruption of mitochondrial integrity and reactive oxygen species. Further, ACT001 reduced the overexpression of the mitochondrial fission protein DRP1 without affecting fusion protein Mitofusin2 levels. Moreover, ACT001 exerted a similar protective effect of suppressing pyroptosis as the DRP1-inhibitor Mdivi-1. Conclusions: Our study revealed that pyroptosis genes were highly expressed in single-cell and spatial mapping along the first week of ALI occurrence. ACT001 attenuates ALI by reducing the NLRP3-dependent pyroptosis and balancing mitochondrial fission and fusion.