Slingshot: cell lineage and pseudotime inference for single-cell transcriptomics

The pancreatic tumor microenvironment of treatment-naïve patients causes a functional shift in γδ T cells, impairing their anti-tumoral defense

Pancreatic ductal adenocarcinoma (PDAC) presents a unique challenge for researchers due to its late diagnosis caused by vague symptoms and lack of early detection markers. Additionally, PDAC is characterized by an immunosuppressive microenvironment (TME), making it a difficult tumor to treat. While γδ T cells have shown potential for anti-tumor activity, conflicting studies exist regarding their effectiveness in pancreatic cancer. This study aims to explore the hypothesis that the PDAC TME hinders the anti-tumor capabilities of γδ T cells through blockade of cytotoxic functions. For this reason, we chose to enroll PDAC treatment-naive patients to avoid the possibility of therapy modifying the TME. By flow cytometry, our research findings indicate that the presence of γδ T cells among CD45+ cells in tumor tissue is lower compared to CD66+ cells, but higher than in blood. Circulating Vδ1 T cells exhibit a terminal effector memory phenotype (TEMRA) more than Vδ2 T cells. Interestingly, Vδ1 and Vδ2 T cells appear to be more prevalent at different stages of tumor development. In our in vitro culture using conditioned medium derived from Patient-derived organoids ;(PDOs), we observed a shift in expression markers in γδ T cells of healthy individuals toward an activation and exhaustion phenotype, as confirmed by scRNA-seq analysis extracted from a public database. A deeper understanding of γδ T cells in PDAC could be valuable for developing novel therapies aimed at mitigating the impact of the pancreatic tumor microenvironment on this cell population.

Plasmodium falciparum infection induces T cell tolerance that is associated with decreased disease severity upon re-infection

Immunity to severe malaria is acquired quickly, operates independently of pathogen load, and represents a highly effective form of disease tolerance. The mechanism that underpins tolerance remains unknown. We used a human rechallenge model of falciparum malaria in which healthy adult volunteers were infected three times over a 12 mo period to track the development of disease tolerance in real-time. We found that parasitemia triggered a hardwired innate immune response that led to systemic inflammation, pyrexia, and hallmark symptoms of clinical malaria across the first three infections of life. In contrast, a single infection was sufficient to reprogram T cell activation and reduce the number and diversity of effector cells upon rechallenge. Crucially, this did not silence stem-like memory cells but instead prevented the generation of cytotoxic effectors associated with autoinflammatory disease. Tolerized hosts were thus able to prevent collateral tissue damage in the absence of antiparasite immunity.

Integration of Single-Cell and Spatial Transcriptomic Data Reveals Spatial Architecture and Potential Biomarkers in Alzheimer's Disease

Alzheimer's disease (AD) is a complex neurodegenerative disorder characterized by the gradual loss of neurons and the accumulation of amyloid plaques and neurofibrillary tangles. Despite advancements in the understanding of AD's pathophysiology, the cellular organization and interactions in the prefrontal cortex (PFC) remain elusive. Eight single-cell RNA sequencing (scRNA-seq) datasets from both normal controls and individuals with AD were harmonized. Stringent preprocessing protocols were implemented to uphold dataset integrity. Unsupervised clustering and annotation revealed 22 distinct cell clusters corresponding to 19 unique cell types. The spatial architecture of the PFC region was constructed using the CARD tool. Further analyses encompassed trajectory examination of Oligodendrocyte subtypes, evaluation of regulon activity scores, and spot clustering within white matter regions (WM). Differential expression analysis and functional enrichment assays unveiled molecular signatures linked to AD progression and were validated using microarray data sourced from neurodegenerative disorder patients. Our investigation employs scRNA-seq and spatial transcriptomics to uncover the cellular atlas and spatial architecture of the human PFC in AD. Moreover, our results indicate that Oligodendrocytes are more prevalent in AD patients, showcasing diverse subtypes and spatial organization within WM regions. Each subtype appears to be associated with distinct biological processes and transcriptional regulators, shedding light on their involvement in AD pathology. Notably, the Oligodendrocyte_C6 subtype is linked to neurological damage in AD patients, characterized by heightened expression of genes involved in cell-cell connections, cell membrane stability, and myelination. Additionally, 12 target genes regulated by NFIA were identified, which are upregulated in AD patients and associated with disease progression. Elevated PLXDC2 expression in peripheral blood was also identified, suggesting its potential as a non-invasive biomarker for early AD detection. Our study provides novel insights into the role of Oligodendrocytes in AD and highlights the potential of PLXDC2 as a blood biomarker for non-invasive diagnosis and monitoring of AD patients.

Circulating exosomes in pediatric obstructive sleep apnea with or without neurocognitive deficits and their effects on a 3D-blood-brain barrier spheroid model

Obstructive sleep apnea (OSA) in children is linked to cognitive impairments, potentially due to blood-brain barrier (BBB) dysfunction. Exosomes, small vesicles released by most cells, reflect cellular changes. This study examined the effects of exosomes from children with OSA, with or without cognitive deficits, on neurovascular unit (NVU) models. Twenty-six children were categorized into three groups: healthy controls (Cont, n = 6), OSA without cognitive deficits (OSA-NG, n = 10), and OSA with neurocognitive deficits (OSA-POS, n = 10). Plasma exosomes were characterized and applied to human 3D NVU spheroids for 24 h. Barrier integrity, permeability, and angiogenesis were assessed using trans-endothelial electrical resistance (TEER), tight junction integrity, and tube formation assays. Single-nucleus RNA sequencing (snRNA-seq) and bioinformatics, including CellChat analysis, identified intercellular signaling pathways. Results showed that exosomes from OSA-POS children disrupted TEER, increased permeability, and impaired ZO1 staining in spheroids, compared to the other groups. Both OSA-POS and OSA-NG exosomes increased permeability in NVU cells in monolayer and microfluidic BBB models. snRNA-seq analysis further revealed distinct cell clusters and pathways associated with the different groups. This 3D NVU spheroid model provides a robust platform to study BBB properties and the role of exosomes in OSA. These findings suggest that integrating snRNA-seq with exosome studies can uncover mechanisms underlying neurocognitive dysfunction in pediatric OSA, potentially leading to personalized therapeutic approaches.

Multi-modal analysis of human hepatic stellate cells identifies novel therapeutic targets for metabolic dysfunction-associated steatotic liver disease

BACKGROUND & AIMS

Metabolic dysfunction-associated steatotic liver disease ranges from metabolic dysfunction-associated steatotic liver (MASL) to metabolic dysfunction-associated steatohepatitis (MASH) with fibrosis. Transdifferentiation of hepatic stellate cells (HSCs) into fibrogenic myofibroblasts plays a critical role in the pathogenesis of MASH liver fibrosis. We compared transcriptome and chromatin accessibility of human HSCs from NORMAL, MASL, and MASH livers at single-cell resolution. We aimed to identify genes that are upregulated in activated HSCs and to determine which of these genes are key in the pathogenesis of MASH fibrosis.

METHODS

Eighteen human livers were profiled using single-nucleus (sn)RNA-seq and snATAC-seq. High priority targets were identified, then tested in 2D human HSC cultures, 3D human liver spheroids, and HSC-specific gene knockout mice.

RESULTS

MASH-enriched activated HSC subclusters are the major source of extracellular matrix proteins. We identified a set of concurrently upregulated and more accessible core genes (GAS7, SPON1, SERPINE1, LTBP2, KLF9, EFEMP1) that drive activation of HSCs. Expression of these genes was regulated via crosstalk between lineage-specific (JUNB/AP1), cluster-specific (RUNX1/2) and signal-specific (FOXA1/2) transcription factors. The pathological relevance of the selected targets, such as SERPINE1 (PAI-1), was demonstrated using dsiRNA-based HSC-specific gene knockdown or pharmacological inhibition of PAI-1 in 3D human MASH liver spheroids, and HSC-specific Serpine1 knockout mice.

CONCLUSION

This study identified novel gene targets and regulatory mechanisms underlying activation of fibrogenic HSCs in MASH, and demonstrated that genetic or pharmacological inhibition of select genes suppressed liver fibrosis.

IMPACT AND IMPLICATIONS

Herein, we present the results of a multi-modal sequencing analysis of human hepatic stellate cells (HSCs) from NORMAL, MASL (metabolic dysfunction-associated steatotic liver), and metabolic dysfunction-associated steatohepatitis (MASH) livers. We identified additional subclusters that were not detected by previous studies and characterized the mechanism by which HSCs are activated in MASH livers, including the transcriptional machinery that induces the transdifferentiation of HSCs into myofibroblasts. For the first time, we described the pathogenic role of activated HSC-derived PAI-1 (a product of the SERPINE1 gene) in the development of MASH liver fibrosis. Targeting the RUNX1/2-SERPINE1 axis could be a novel strategy for the treatment of liver fibrosis in patients.

Mutation S139N on Zika virus prM protein shifts immune response from Asian to contemporary strain

Zika virus (ZIKV) has been associated with neurological diseases like microcephaly and Guillain-Barré syndrome. The S139N single mutation on the prM protein of the FSS13025 Asian strain increases the mortality rate in mice. Therefore, it is a valuable tool for studying the impact of immune responses on neural damage. Here, we used single-cell sequencing technology to systematically assess the immune response induced by three ZIKV strains: Asian ancestral strain FSS13025/2010, FSS13025 strain with S139N mutation (FSS13025-S139N), and contemporary strain GZ01/2016. By infecting 1-day-old mice, we observed that the immune spectrum elicited by FSS13025-S139N mutant resembled that induced by the contemporary strain. The FSS13025-S139N strain induces the proliferation of inflammatory microglial cells earlier than the FSS13025 strain, similar to GZ01. A specific cell cluster, Microglia_Ccr7, was induced by the S139N mutant strain and GZ01 strain, which suppresses T cell activation through the PDCD1LG2-PDCD1 signaling pathway. Furthermore, the proliferation of CD8+ T cells was weakened and prolonged in S139N strain-infected samples. Finally, we found that the S139N mutant strain causes more apoptosis of neurons compared to the FSS13025 strain. These results indicate that the S139N mutation plays an important role in the immune response pattern of ZIKV and prolongs the duration of neuroinflammation.

Distinct adipose progenitor cells emerging with age drive active adipogenesis

Starting at middle age, adults often suffer from visceral adiposity and associated adverse metabolic disorders. Lineage tracing in mice revealed that adipose progenitor cells (APCs) in visceral fat undergo extensive adipogenesis during middle age. Thus, despite the low turnover rate of adipocytes in young adults, adipogenesis is unlocked during middle age. Transplantations quantitatively showed that APCs in middle-aged mice exhibited high adipogenic capacity cell-autonomously. Single-cell RNA sequencing identified a distinct APC population, the committed preadipocyte, age-enriched (CP-A), emerging at this age. CP-As demonstrated elevated proliferation and adipogenesis activity. Pharmacological and genetic manipulations indicated that leukemia inhibitory factor receptor signaling was indispensable for CP-A adipogenesis and visceral fat expansion. These findings uncover a fundamental mechanism of age-dependent adipose remodeling, offering critical insights into age-related metabolic diseases.

α-Ketoglutarate promotes trophectoderm induction and maturation from naive human embryonic stem cells

Development and lineage choice are driven by interconnected transcriptional, epigenetic and metabolic changes. Specific metabolites, such as α-ketoglutarate (αKG), function as signalling molecules affecting the activity of chromatin-modifying enzymes. However, how metabolism coordinates cell-state changes, especially in human pre-implantation development, remains unclear. Here we uncover that inducing naive human embryonic stem cells towards the trophectoderm lineage results in considerable metabolic rewiring, characterized by αKG accumulation. Elevated αKG levels potentiate the capacity of naive embryonic stem cells to specify towards the trophectoderm lineage. Moreover, increased αKG levels promote blastoid polarization and trophectoderm maturation. αKG supplementation does not affect global histone methylation levels; rather, it decreases acetyl-CoA availability, reduces histone acetyltransferase activity and weakens the pluripotency network. We propose that metabolism functions as a positive feedback loop aiding in trophectoderm fate induction and maturation, highlighting that global metabolic rewiring can promote specificity in cell fate decisions through intricate regulation of signalling and chromatin.

STAT3 signaling enhances tissue expansion during postimplantation mouse development

Signal transducer and activator of transcription (STAT)3 signaling has been studied extensively using mouse embryonic stem cells. Zygotic deletion of Stat3 enables embryo implantation, but thereafter, mutants resemble non-affected littermates from the previous day until around mid-gestation. This probably results from the loss of serine-phosphorylated STAT3, the predominant form in early postimplantation embryonic tissues associated with mitochondrial activity. Bulk RNA sequencing of isolated mouse epiblasts confirmed developmental delay transcriptionally. Single-cell RNA sequencing revealed the exclusion of derivatives of Stat3 null embryonic stem cells exclusively from the erythroid lineage of mid-gestation chimeras. We show that Stat3 null embryonic stem cells can differentiate into erythroid and hematopoietic lineages in vitro but become outcompeted when mixed with wild-type cells. Our results implicate a role for STAT3 in the temporal control of embryonic progression, particularly in tissues requiring rapid cell division, such as postimplantation epiblast and hematopoietic lineages. Interestingly, mutations in STAT3 are associated with short stature in humans.

Integrated analysis of molecular atlases unveils modules driving developmental cell subtype specification in the human cortex

Human brain development requires generating diverse cell types, a process explored by single-cell transcriptomics. Through parallel meta-analyses of the human cortex in development (seven datasets) and adulthood (16 datasets), we generated over 500 gene co-expression networks that can describe mechanisms of cortical development, centering on peak stages of neurogenesis. These meta-modules show dynamic cell subtype specificities throughout cortical development, with several developmental meta-modules displaying spatiotemporal expression patterns that allude to potential roles in cell fate specification. We validated the expression of these modules in primary human cortical tissues. These include meta-module 20, a module elevated in FEZF2+ deep layer neurons that includes TSHZ3, a transcription factor associated with neurodevelopmental disorders. Human cortical chimeroid experiments validated that both FEZF2 and TSHZ3 are required to drive module 20 activity and deep layer neuron specification but through distinct modalities. These studies demonstrate how meta-atlases can engender further mechanistic analyses of cortical fate specification.

Single-cell analysis of the early Drosophila salivary gland reveals that morphogenetic control involves both the induction and exclusion of gene expression programs

How tissue shape and therefore function is encoded by the genome remains in many cases unresolved. The tubes of the salivary glands in the Drosophila embryo start from simple epithelial placodes, specified through the homeotic factors Scr/Hth/Exd. Previous work indicated that early morphogenetic changes are prepatterned by transcriptional changes, but an exhaustive transcriptional blueprint driving physical changes was lacking. We performed single-cell-RNAseq-analysis of FACS-isolated early placodal cells, making up less than 0.4% of cells within the embryo. Differential expression analysis in comparison to epidermal cells analyzed in parallel generated a repertoire of genes highly upregulated within placodal cells prior to morphogenetic changes. Furthermore, clustering and pseudotime analysis of single-cell-sequencing data identified dynamic expression changes along the morphogenetic timeline. Our dataset provides a comprehensive resource for future studies of a simple but highly conserved morphogenetic process of tube morphogenesis. Unexpectedly, we identified a subset of genes that, although initially expressed in the very early placode, then became selectively excluded from the placode but not the surrounding epidermis, including hth, grainyhead and tollo/toll-8. We show that maintaining tollo expression severely compromised the tube morphogenesis. We propose tollo is switched off to not interfere with key Tolls/LRRs that are expressed and function in the tube morphogenesis.

Granulosa cell transcription is similarly impacted by superovulation and aging and predicts early embryonic trajectories

In vitro fertilization efficiency is limited in part because a fraction of retrieved oocytes fails to fertilize. Accurately evaluating their quality could significantly improve in vitro fertilization efficiency, which would require better understanding how their maturation may be disrupted. Here, we quantitatively investigate the interplay between superovulation and aging in mouse oocytes and their paired granulosa cells using a newly adapted experimental methodology. We test the hypothesis that superovulation disrupts oocyte maturation, revealing the key intercellular communication pathways dysregulated at the transcriptional level by forced hormonal stimulation. We further demonstrate that granulosa cell transcriptional markers can prospectively predict an associated oocyte's early developmental potential. By using naturally ovulated old mice as a non-stimulated reference, we show that aging and superovulation dysregulate similar genes and interact with each other. By comparing mice and human transcriptional responses of granulosa cells, we find that age-related dysregulation of hormonal responses and cell cycle pathways are shared, though substantial divergence exists in other pathways.

Peripheral nervous system microglia-like cells regulate neuronal soma size throughout evolution

Microglia, essential in the central nervous system (CNS), were historically considered absent from the peripheral nervous system (PNS). Here, we show a PNS-resident macrophage population that shares transcriptomic and epigenetic profiles as well as an ontogenetic trajectory with CNS microglia. This population (termed PNS microglia-like cells) enwraps the neuronal soma inside the satellite glial cell envelope, preferentially associates with larger neurons during PNS development, and is required for neuronal functions by regulating soma enlargement and axon growth. A phylogenetic survey of 24 vertebrates revealed an early origin of PNS microglia-like cells, whose presence is correlated with neuronal soma size (and body size) rather than evolutionary distance. Consistent with their requirement for soma enlargement, PNS microglia-like cells are maintained in vertebrates with large peripheral neuronal soma but absent when neurons evolve to have smaller soma. Our study thus reveals a PNS counterpart of CNS microglia that regulates neuronal soma size during both evolution and ontogeny.

Resolving floral development dynamics using genome and single-cell temporal transcriptome of Dendrobium devonianum

Dendrobium devonianum, a species of the Orchidaceae family, is notable for its unique floral characteristics, which include two yellow spots and purple tips on its labellum, as well as fringed edges. However, the molecular mechanisms underlying flower pattern formation in D. devonianum remain poorly understood, hindering advancements in its breeding process. Here, a chromosome-scale genome of D. devonianum was presented for the first time, revealing two significant polyploidization events. Additionally, a high-resolution single-cell transcriptomic atlas was constructed, capturing 11 distinct cell clusters. Expression patterns of MADS-box genes were identified through temporal and spatial bulk RNA-Seq, revealing alignment with the ABCDE model of flower formation. Meanwhile, mass spectrometry imaging and scRNA analyses showed that the yellow spots were primarily associated with carotenoid biosynthesis gene expression, while the purple colour is predominantly linked to anthocyanin biosynthesis gene expression. These genes were mainly expressed in the epidermis and vascular cells. Developmental trajectory analyses of epidermal cells further uncovered a gene regulatory network and several transcription factors likely responsible for fringes formation along the labellum margin. This study provides valuable insights into the molecular mechanisms driving floral colour differentiation and structural traits in D. devonianum, contributing to a deeper understanding of orchid evolution, diversification and breeding.

Integrating bulk, single-cell, and spatial transcriptomics to identify and functionally validate novel targets to enhance immunotherapy in NSCLC

Programmed cell deaths (PCDs) are crucial for tumor progression. By analyzing 18 PCDs, we generated a robust multigene signature, Combined Cell Death Index (CCDI), comprising necroptosis and autophagy genes for non-small cell lung cancer (NSCLC). The CCDI accurately stratified patients by survival prognosis and predicted immunotherapy responses. We validated CCDI and prioritized CCDI genes using five single-cell RNA sequencing and two spatial transcriptomics datasets. CCDI positively correlates with tumor malignancy, invasiveness, and immunotherapy resistance. Four necroptosis genes (PTGES3, MYO6, CCT6A, and CTSH) may affect cancer cell evolution. In vitro, CTSH overexpression or PTGES3 knockdown inhibited NSCLC cell proliferation and migration while inducing necroptosis with necrosome formation. Moreover, we observed diminished CTSH, heightened PTGES3, and low necroptosis activity in 12 pairs of NSCLC tumors and normal tissues. CTSH overexpression or PTGES3 knockdown induced necroptosis and improved anti-PD1 therapy efficiency in syngeneic cancer mouse models. These findings indicate necroptosis genes as potential therapeutic targets in cancer treatments.

Redefining the high variable genes by optimized LOESS regression with positive ratio

BACKGROUND

Single-cell RNA sequencing allows for the exploration of transcriptomic features at the individual cell level, but the high dimensionality and sparsity of the data pose substantial challenges for downstream analysis. Feature selection, therefore, is a critical step to reduce dimensionality and enhance interpretability.

RESULTS

We developed a robust feature selection algorithm that leverages optimized locally estimated scatterplot smoothing regression (LOESS) to precisely capture the relationship between gene average expression level and positive ratio while minimizing overfitting. Our evaluations showed that our algorithm consistently outperforms eight leading feature selection methods across three benchmark criteria and helps improve downstream analysis, thus offering a significant improvement in gene subset selection.

CONCLUSIONS

By preserving key biological information through feature selection, GLP provides informative features to enhance the accuracy and effectiveness of downstream analyses.

Single cell RNAseq to identify subpopulations of glial progenitors in iPSC-derived oligodendroglial lineage cultures

Cellular heterogeneity is a common issue in differentiation protocols of oligodendrocytes (OLs) from human induced pluripotent stem cells. Our previous work described a novel method to generate OLs and highlighted the presence of glial progenitors. Here, we unravel the glial heterogeneity and characterize the response of isolated subpopulations to differentiation. This study provides a novel tool for studying the dynamics of glial development in vitro and on a transcriptomic level.

Distinct CD8+ T cell dynamics associate with response to neoadjuvant cancer immunotherapies

We leverage a clinical trial (NCT04080804) that compared neoadjuvant anti-PD-1, anti-PD-1+CTLA-4, and anti-PD-1+LAG-3 therapies in head and neck squamous cell carcinoma patients. Combination therapies promote higher pathologic response rates versus monotherapy, and major pathologic response is associated with better survival. To address whether successful immune checkpoint inhibitor (ICI) regimens act through similar or distinct pathways, we robustly and longitudinally characterize transcriptional and proteomic dynamics of CD8+ tumor-infiltrating lymphocytes (TILs) in a clonal manner. Anti-PD-1+LAG-3 reprograms CD8+ TIL with type-I interferon response and exhaustion gene programs into effector memory and resident memory (TEM/TRM). In contrast, anti-PD-1+CTLA-4 activates and expands pre-existing TEM/TRM CD8+ TIL, but does not rejuvenate exhausted phenotypes into T effector cells. Anti-PD-1+LAG-3, but not anti-PD-1+CTLA-4, induces widespread TCR sharing among the different transcriptional states, as well as increased TCR diversity in responding patients. Our data suggest doublet regimen-specific transcriptional and clonal dynamics of tumor-reactive CD8+ T cells.

New horizons at the interface of artificial intelligence and translational cancer research

Artificial intelligence (AI) is increasingly being utilized in cancer research as a computational strategy for analyzing multiomics datasets. Advances in single-cell and spatial profiling technologies have contributed significantly to our understanding of tumor biology, and AI methodologies are now being applied to accelerate translational efforts, including target discovery, biomarker identification, patient stratification, and therapeutic response prediction. Despite these advancements, the integration of AI into clinical workflows remains limited, presenting both challenges and opportunities. This review discusses AI applications in multiomics analysis and translational oncology, emphasizing their role in advancing biological discoveries and informing clinical decision-making. Key areas of focus include cellular heterogeneity, tumor microenvironment interactions, and AI-aided diagnostics. Challenges such as reproducibility, interpretability of AI models, and clinical integration are explored, with attention to strategies for addressing these hurdles. Together, these developments underscore the potential of AI and multiomics to enhance precision oncology and contribute to advancements in cancer care.

Identification of osteoarthritis-associated chondrocyte subpopulations and key gene-regulating drugs based on multi-omics analysis

The mechanism by which chondrocytes respond to mechanical stress in joints significantly affects the balance and function of cartilage. This study aims to characterize osteoarthritis-associated chondrocyte subpopulations and key gene targets for regulatory drugs. To begin, single-cell and transcriptome datasets were obtained from the Gene Expression Omnibus (GEO) database. Cell communication and pseudo-temporal analysis, as well as High-dimensional Weighted Gene Co-expression Network Analysis (hdWGCNA), were conducted on the single-cell data to identify key chondrocyte subtypes and module genes. Subsequently, Consensus Cluster Plus analysis was utilized to identify distinct disease subgroups within the osteoarthritis (OA) training dataset based on the key module genes. Furthermore, differential gene expression analysis and GO/KEGG pathway enrichment analysis were performed on the identified subgroups. To screen for hub genes associated with OA, a combination of 10 machine learning algorithms and 113 algorithm compositions was integrated. Additionally, the immune and pathway scores of the training dataset samples were evaluated using the ESTIMATE, MCP-counter, and ssGSEA algorithms to establish the relationship between the hub genes and immune and pathways. Following this, a network depicting the interaction between the hub genes and transcription factors was constructed based on the Network Analyst database. Moreover, the hub genes were subjected to drug prediction and molecular docking using the RNAactDrug database and AutoDockTools. Finally, real-time fluorescence quantitative PCR (RT-qPCR) was employed to detect the expression of hub genes in the plasma samples collected from osteoarthritis patients and healthy adults. In the OA sample, there is a significant increase in the proportion of prehypertrophic chondrocytes (preHTC), particularly in subgroups 6, 7, and 9. We defined these subgroups as OA_PreHTC subgroups. The OA_PreHTC subgroup exhibits a higher communication intensity with proliferative-related pathways such as ANGPTL and TGF-β. Furthermore, two OA disease subgroups were identified in the training set samples. This led to the identification of 411 differentially expressed genes (DEGs) related to osteoarthritis, 2485 DEGs among subgroups, as well as 238 intersecting genes and 5 hub genes (MMP13, FAM26F, CHI3L1, TAC1, and CKS2). RT-qPCR results indicate significant differences in the expression levels of five hub genes and their related TFs in the clinical blood samples of OA patients compared to the healthy control group (NC). Moreover, these five hub genes are positively associated with inflammatory pathways such as TNF-α, JAK-STAT3, and inflammatory response, while being negatively associated with proliferation pathways like WNT and KRAS. Additionally, the five hub genes are positively associated with neutrophils, activated CD4 T cell, gamma delta T cell, and regulatory T cell, while being negatively associated with CD56dim natural killer cell and Type 17T helper cell. Molecular docking results reveal that CAY10603, Tenulin, T0901317, and Nonactin exhibit high binding activity to CHI3L1, suggesting their potential as therapeutic drugs for OA. The OA_PreHTC subgroups plays a crucial role in the occurrence and development of osteoarthritis (OA). Five hub genes may exert their effects on OA through interactions with PreHTC cells, other chondrocytes, and immune cells, playing a role in inhibiting cell proliferation and stimulating inflammation, thus having high diagnostic value for OA. Additionally, CAY10603, Tenulin, T0901317, and Nonactin have potential therapeutic effects for OA patients.

Combining spatial transcriptomics and ECM imaging in 3D for mapping cellular interactions in the tumor microenvironment

Tumors are complex ecosystems composed of malignant and non-malignant cells embedded in a dynamic extracellular matrix (ECM). In the tumor microenvironment, molecular phenotypes are controlled by cell-cell and ECM interactions in 3D cellular neighborhoods (CNs). While their inhibition can impede tumor progression, routine molecular tumor profiling fails to capture cellular interactions. Single-cell spatial transcriptomics (ST) maps receptor-ligand interactions but usually remains limited to 2D tissue sections and lacks ECM readouts. Here, we integrate 3D ST with ECM imaging in serial sections from one clinical lung carcinoma to systematically quantify molecular states, cell-cell interactions, and ECM remodeling in CN. Our integrative analysis pinpointed known immune escape and tumor invasion mechanisms, revealing several druggable drivers of tumor progression in the patient under study. This proof-of-principle study highlights the potential of in-depth CN profiling in routine clinical samples to inform microenvironment-directed therapies. A record of this paper's transparent peer review process is included in the supplemental information.

Characterisation of the avascular mesenchyme during digit outgrowth

The avascular mesenchyme at the tip of the developing digit contributes to digit outgrowth and patterning, however, it has been poorly characterised. Using newly developed fate mapping approaches, tissue manipulation and single-cell mRNA sequencing data, we explore the transcriptional nature and developmental potential of this tissue. We find that the avascular mesenchyme is essential to normal segmental patterning of the digit and has a distinct transcriptional identity. In addition, we uncover an unexpected relationship between the unspecified tissue of the avascular mesenchyme and the committed phalanx forming region, which controls patterning, but not outgrowth of the digit. This multifaceted approach provides insights into the mechanics and genetic pathways that regulate digit outgrowth and patterning.

Anti-vascular endothelial growth factor treatment potentiates immune checkpoint blockade through a BAFF- and IL-12-dependent reprogramming of the TME

Anti-vascular endothelial growth factor (VEGF) treatment has shown clinical activity together with immune checkpoint blockade (ICB), but the exact mechanism is not known. We show that VEGF blockade in combination with anti-cytotoxic T-lymphocyte associated protein 4 (CTLA4) + anti-programmed death-ligand 1 (PD-L1) in cholangiocarcinoma (CCA) potentiated a multimodal mechanism dependent on B cell activating factor (BAFF), leading to a proinflammatory B cell response. It led to a BAFF- and interleukin (IL)-12-dependent expansion and rewiring of T regulatory cells (Tregs) toward an anti-tumor T helper-1 (Th-1)-like fragile state. We translated this approach to the clinic and observed immunological changes characterized by Treg cell expansion and rewiring toward fragile and unstable states. We explored the effect of VEGF receptor 2 (VEGFR2) signaling on Treg cell transcriptional programming and established a mouse model ablating VEGFR2 expression on Treg cells. This study reveals the immunological interplay resulting from targeting VEGF together with CTLA-4 and PD-L1 blockade.

Mechanisms of hematopoietic clonal dominance in VEXAS syndrome

Clonal dominance characterizes hematopoiesis during aging and increases susceptibility to blood cancers and common nonmalignant disorders. VEXAS syndrome is a recently discovered, adult-onset, autoinflammatory disease burdened by a high mortality rate and caused by dominant hematopoietic clones bearing somatic mutations in the UBA1 gene. However, pathogenic mechanisms driving clonal dominance are unknown. Moreover, the lack of disease models hampers the development of disease-modifying therapies. In the present study, we performed immunophenotype characterization of hematopoiesis and single-cell transcriptomics in a cohort of nine male patients with VEXAS syndrome, revealing pervasive inflammation across all lineages. Hematopoietic stem and progenitor cells (HSPCs) in patients are skewed toward myelopoiesis and acquire senescence-like programs. Humanized models of VEXAS syndrome, generated by inserting the causative mutation in healthy HSPCs through base editing, recapitulated proteostatic defects, cytological alterations and senescence signatures of patients' cells, as well as hematological and inflammatory disease hallmarks. Competitive transplantations of human UBA1-mutant and wild-type HSPCs showed that, although mutant cells are more resilient to the inflammatory milieu, probably through the acquisition of the senescence-like state, wild-type ones are progressively exhausted and overwhelmed by VEXAS clones, overall impairing functional hematopoiesis and leading to bone marrow failure. Our study unveils the mechanism of clonal dominance and provides models for preclinical studies and preliminary insights that could inform therapeutic strategies.

Sphingolipid metabolism orchestrates establishment of the hair follicle stem cell compartment

Sphingolipids serve as building blocks of membranes to ensure subcellular compartmentalization and facilitate intercellular communication. How cell type-specific lipid compositions are achieved and what is their functional significance in tissue morphogenesis and maintenance has remained unclear. Here, we identify a stem cell-specific role for ceramide synthase 4 (CerS4) in orchestrating fate decisions in skin epidermis. Deletion of CerS4 prevents the proper development of the adult hair follicle bulge stem cell (HFSC) compartment due to altered differentiation trajectories. Mechanistically, HFSC differentiation defects arise from an imbalance of key ceramides and their derivate sphingolipids, resulting in hyperactivation of noncanonical Wnt signaling. This impaired HFSC compartment establishment leads to disruption of hair follicle architecture and skin barrier function, ultimately triggering a T helper cell 2-dominated immune infiltration resembling human atopic dermatitis. This work uncovers a fundamental role for a cell state-specific sphingolipid profile in stem cell homeostasis and in maintaining an intact skin barrier.

Multiomics Analysis Reveals Neuroblastoma Molecular Signature Predicting Risk Stratification and Tumor Microenvironment Differences

Neuroblastoma (NB) remains associated with high mortality and low initial response rate, especially for high-risk patients, thus warranting exploration of molecular markers for precision risk classifiers. Through integrating multiomics profiling, we identified a range of hub genes involved in cell cycle and associated with dismal prognosis and malignant cells. Single-cell transcriptome sequencing revealed that a subset of malignant cells, subcluster 1, characterized by high proliferation and dedifferentiation, was strongly correlated with the hub gene signature and orchestrated an immunosuppressive tumor microenvironment (TME). Furthermore, we constructed a robust malignant subcluster 1 related signature (MSRS), which was an independent prognostic factor and superior to other clinical characteristics and published signatures. Besides, TME differences conferred remarkably distinct therapeutic responses between high and low MSRS groups. Notably, polo-like kinase-1 (PLK1) was one of the most crucial contributors to MSRS and remarkably correlated with malignant subcluster 1, and PLK1 inhibition was effective for NB treatment as demonstrated by in silico analysis and in vitro experiments. Overall, our study constructs a novel molecular model to further guide the clinical classification and individualized treatment of NB.

Transcriptional landscapes underlying Notch-induced lineage conversion and plasticity of mammary basal cells

The mammary epithelium derives from multipotent mammary stem cells (MaSCs) that engage into differentiation during embryonic development. However, adult MaSCs maintain the ability to reactivate multipotency in non-physiological contexts. We previously reported that Notch1 activation in committed basal cells triggers a basal-to-luminal cell fate switch in the mouse mammary gland. Here, we report conservation of this mechanism and found that in addition to the mammary gland, constitutive Notch1 signaling induces a basal-to-luminal cell fate switch in adult cells of the lacrimal gland, the salivary gland, and the prostate. Since the lineage transition is progressive in time, we performed single-cell transcriptomic analysis on index-sorted mammary cells at different stages of lineage conversion, generating a temporal map of changes in cell identity. Combining single-cell analyses with organoid assays, we demonstrate that cell proliferation is indispensable for this lineage conversion. We also reveal the individual transcriptional landscapes underlying the cellular plasticity switching of committed mammary cells in vivo with spatial and temporal resolution. Given the roles of Notch signaling in cancer, these results may help to better understand the mechanisms that drive cellular transformation.

Transcriptional conservation and evolutionary divergence of cell types across mammalian hypothalamus development

The hypothalamus, an "ancient" subcortical brain structure, maintains physiological homeostasis and controls native behaviors. The evolution of homeostatic regulation and behavioral control in mammals may rely on adaptable neuronal identity establishment but conserved neural patterning mechanisms during neurodevelopment. Here, we combined single-cell, single-nucleus, and spatial transcriptomic datasets to map the spatial patterning of diverse progenitor domains and reconstruct their neurogenic lineages in the developing human and mouse hypothalamus. While the regional organizers orchestrating neural patterning are conserved between primates and rodents, we identified a human-enriched neuronal subtype and found a substantial increase in neuromodulatory gene expression among human neurons. Furthermore, cross-species comparison demonstrated a potential redistribution of two neuroendocrine neuronal subtypes and a shift in inter-transmitter and transmitter-peptide coupling within hypothalamic dopamine neurons. Together, our study lays a critical foundation for understanding cellular development and evolution of the mammalian hypothalamus.

CXXC-finger protein 1 associates with FOXP3 to stabilize homeostasis and suppressive functions of regulatory T cells

FOXP3-expressing regulatory T (Treg) cells play a pivotal role in maintaining immune homeostasis and tolerance, with their activation being crucial for preventing various inflammatory responses. However, the mechanisms governing the epigenetic program in Treg cells during their dynamic activation remain unclear. In this study, we demonstrate that CXXC-finger protein 1 (CXXC1) interacts with the transcription factor FOXP3 and facilitates the regulation of target genes by modulating H3K4me3 deposition. Cxxc1 deletion in Treg cells leads to severe inflammatory disease and spontaneous T cell activation, with impaired immunosuppressive function. As a transcriptional regulator, CXXC1 promotes the expression of key Treg functional markers under steady-state conditions, which are essential for the maintenance of Treg cell homeostasis and their suppressive functions. Epigenetically, CXXC1 binds to the genomic regulatory regions of Treg program genes in mouse Treg cells, overlapping with FOXP3-binding sites. Given its critical role in Treg cell homeostasis, CXXC1 presents itself as a promising therapeutic target for autoimmune diseases.

Single-cell multi-omics profiling uncovers the immune heterogeneity in HIV-infected immunological non-responders

BACKGROUND

Immunological non-responders (INRs) are people living with HIV-1 who fail to achieve full immune reconstitution despite long-term effective antiretroviral therapy (ART). This incomplete recovery of CD4+ T cells increase the risk of opportunistic infections and non-AIDS-related morbidity and mortality. Understanding the mechanisms driving this immune dysfunction is critical for developing targeted therapies.

METHODS

We performed single-cell RNA sequencing (scRNA-seq) and single-cell VDJ sequencing (scVDJ-seq) on peripheral blood mononuclear cells (PBMCs) from INRs, immune responders (IRs), and healthy controls (HCs). We developed scGeneANOVA, a novel mixed model differential gene analysis tool, to detect differentially expressed genes and pathways. In addition, we developed the Viral Identification and Load Detection Analysis (VILDA) tool to quantify HIV-1 transcripts and investigate their relationship with interferon (IFN) pathway activation.

FINDINGS

Our analysis revealed that INRs exhibit a dysregulated IFN response, closely associated with CD4+ T cell exhaustion and immune recovery failure. The scGeneANOVA tool identified critical genes and pathways that were missed by traditional analysis methods, while VILDA showed higher levels of HIV-1 transcripts in INRs, which may drive the heightened IFN response. These findings support a potential contribution of IFN signalling in INR-related immune dysfunction.

INTERPRETATION

Our study provides new insights into the pathogenic mechanisms behind immune recovery failure in INRs, suggesting that IFN signalling might be involved in the development of CD4+ T cell exhaustion. The identification of key genes and pathways offers potential biomarkers and therapeutic targets for improving immune recovery in this vulnerable population.

FUNDING

This study was supported by the grants from Special Research Fund for the Central High-level Hospitals of Peking Union Medical College Hospital (Grant No. 2022-PUMCH-D-008), Chinese Academy of Medical Sciences (CAMS) Innovation Fund for Medical Sciences (Grant No. 2021-I2M-1-037), National Key Technologies R&D Program for the 13th Five-year Plan (Grant No. 2017ZX10202101-001). The funders played no role in the design, experiment conduction, data analysis and preparation of the manuscript of this work.

Integrative network analysis reveals novel moderators of Aβ-Tau interaction in Alzheimer's disease

BACKGROUND

Although interactions between amyloid-beta and tau proteins have been implicated in Alzheimer's disease (AD), the precise mechanisms by which these interactions contribute to disease progression are not yet fully understood. Moreover, despite the growing application of deep learning in various biomedical fields, its application in integrating networks to analyze disease mechanisms in AD research remains limited. In this study, we employed BIONIC, a deep learning-based network integration method, to integrate proteomics and protein-protein interaction data, with an aim to uncover factors that moderate the effects of the Aβ-tau interaction on mild cognitive impairment (MCI) and early-stage AD.

METHODS

Proteomic data from the ROSMAP cohort were integrated with protein-protein interaction (PPI) data using a Deep Learning-based model. Linear regression analysis was applied to histopathological and gene expression data, and mutual information was used to detect moderating factors. Statistical significance was determined using the Benjamini-Hochberg correction (p < 0.05).

RESULTS

Our results suggested that astrocytes and GPNMB + microglia moderate the Aβ-tau interaction. Based on linear regression with histopathological and gene expression data, GFAP and IBA1 levels and GPNMB gene expression positively contributed to the interaction of tau with Aβ in non-dementia cases, replicating the results of the network analysis.

CONCLUSIONS

These findings suggest that GPNMB + microglia moderate the Aβ-tau interaction in early AD and therefore are a novel therapeutic target. To facilitate further research, we have made the integrated network available as a visualization tool for the scientific community (URL: https://igcore.cloud/GerOmics/AlzPPMap ).

Spatially Resolved Tumor Ecosystems and Cell States in Gastric Adenocarcinoma Progression and Evolution

SIGNIFICANCE

Integration of spatial transcriptomic (GeoMx Digital Spatial Profiler) and single-cell RNA sequencing data from multiple gastric cancers identifies spatially resolved expression-based intratumoral heterogeneity, associated with distinct immune microenvironments. We uncovered two separate evolutionary trajectories associated with specific molecular subtypes, clinical prognoses, stromal neighborhoods, and genetic drivers. Tumor-stroma interfaces emerged as a unique state of tumor ecology.

Vps34-orchestrated lipid signaling processes regulate the transitional heterogeneity and functional adaptation of effector regulatory T cells

Regulatory T cell (Treg) heterogeneity exists in lymphoid and non-lymphoid tissues, but we have limited understanding of context-dependent functions and spatiotemporal regulators of heterogenous Treg states, especially during perinatal life when immune tolerance is established. Here, we revealed that the class III PI3K Vps34 orchestrates effector Treg (eTreg) transitional heterogeneity during perinatal life. We found that loss of Vps34 reduced terminal eTreg accumulation in lymphoid tissues, associated with decreased Treg generation in non-lymphoid tissues and development of an early-onset autoimmune-like disease. After perinatal life, Vps34-deficient eTreg accumulation was further impaired due to reduced cell survival, highlighting temporal regulation of eTreg heterogeneity and maintenance by Vps34. Accordingly, inhibition of Vps34 in mature Tregs disrupted immune homeostasis but boosted anti-tumor immunity. Mechanistically, multiomics profiling approaches uncovered that Vps34-orchestrated transcriptional and epigenetic remodeling promotes terminal eTreg programming. Further, via genetic deletion of the Vps34-interacting proteins Atg14 or Uvrag in Tregs, we established that Atg14 but not Uvrag was required for the overall survival, but not terminal differentiation, of eTregs, suggesting that autophagy but not endocytosis partly contributed to Vps34-dependent effects. Accordingly, mice with Treg-specific loss of Atg14, but not Uvrag, had moderately disrupted immune homeostasis and reduced tumor growth, with Vps34- or Atg14-dependent gene signatures also being elevated in intratumoral Tregs from human cancer patients. Collectively, our study reveals distinct Vps34-orchestrated signaling events that regulate eTreg heterogeneity and functional adaptation and the pathophysiological consequences on autoimmunity versus anti-tumor immunity.

Origin and stepwise evolution of vertebrate lungs

Lungs are essential respiratory organs in terrestrial vertebrates, present in most bony fishes but absent in cartilaginous fishes, making them an ideal model for studying organ evolution. Here we analysed single-cell RNA sequencing data from adult and developing lungs across vertebrate species, revealing significant similarities in cell composition, developmental trajectories and gene expression patterns. Surprisingly, a large proportion of lung-related genes, coexpression patterns and many lung enhancers are present in cartilaginous fishes despite their lack of lungs, suggesting that a substantial genetic foundation for lung development existed in the last common ancestor of jawed vertebrates. In addition, the 1,040 enhancers that emerged since the last common ancestor of bony fishes probably contain lung-specific elements that led to the development of lungs. We further identified alveolar type 1 cells as a mammal-specific alveolar cell type, along with several mammal-specific genes, including ager and sfta2, that are highly expressed in lungs. Functional validation showed that deletion of sfta2 in mice leads to severe respiratory defects, highlighting its critical role in mammalian lung features. Our study provides comprehensive insights into the evolution of vertebrate lungs, demonstrating how both regulatory network modifications and the emergence of new genes have shaped lung development and specialization across species.

X-Linked Gene Dosage and SOX2 Act as Key Roadblocks for Human Germ Cell Specification in Klinefelter Syndrome

Klinefelter syndrome (KS), characterized by the presence of at least one extra X-chromosome, is a common cause of male infertility. However, the mechanism underlying the failure of germline specification is not well studied. Intriguingly, the differentiation efficiency of female human pluripotent stem cells (hPSCs) is often lower than that of male. This study investigates how X-linked gene dosage affects human primordial germ cell-like cells (hPGCLCs) specification in both healthy and diseased conditions. This work reveals that X-linked genes play a multifaceted role against the fate competency to hPGCLCs, with escape genes IGSF1 and CHRDL1 inhibiting the TGF-beta/Activin A and BMP pathways, respectively. Notably, this work identifies a previously unrecognized role of SOX2, upregulated by the escape gene USP9X, elucidating a species-specific function in the mammalian germline. The USP9X-SOX2 regulatory axis profoundly influenced cellular metabolism, mitochondrial morphology, and progenitor competence in hPGCLCs specification. Furthermore, the inability to downregulate SOX2 and upregulate SOX17 in response to BMP signaling impedes downstream gene activation due to motif binding competition. These findings shed novel insights into the human germline specification by elucidating the divergent roles of SOX2 versus SOX17 in mammals, influenced by X-linked gene dosage effects. These results offer potential applications for improving the induction efficiency of hPGCLCs, facilitating disease mechanistic studies.

Tumor Microenvironment On-A-Chip and Single-Cell Analysis Reveal Synergistic Stromal-Immune Crosstalk on Breast Cancer Progression

Solid tumors develop within a complex environment called the tumor microenvironment (TME), which is sculpted by the presence of other cells, such as cancer-associated fibroblasts (CAFs) and immune cells like macrophages (Mφs). Despite the presence of immune cells, tumor cells orchestrate a tumor-supportive environment through intricate interaction with the components of the TME. However, the specific mechanism by which this intercellular dialogue is regulated is not fully understood. To that end, the development of an organotypic 3D breast TME-on-a-chip (TMEC) model, integrated with single-cell RNA sequencing analysis, is reported to mechanistically evaluate the progression of triple-negative breast cancer (TNBC) cells in the presence of patient-derived CAFs and Mφs. Extensive functional assays, including invasion and morphometric characterization, reveal the synergistic influence of CAFs and Mφs on tumor cells. Furthermore, gene expression and pathway enrichment analyses identify the involvement of the KYNU gene, suggesting a potential immune evasion mechanism through the kynurenine pathway. Lastly, the pharmacological targeting of the identified pathway is investigated.

3D Bioprinted Coaxial Testis Model Using Human Induced Pluripotent Stem Cells:A Step Toward Bicompartmental Cytoarchitecture and Functionalization

Fertility preservation following pediatric cancer therapy programs has become a major avenue of infertility research. In vitro spermatogenesis (IVS) aims to generate sperm from banked prepubertal testicular tissues in a lab setting using specialized culture conditions. While successful using rodent tissues, progress with human tissues is limited by the scarcity of human prepubertal testicular tissues for research. This study posits that human induced pluripotent stem cells (hiPSCs) can model human prepubertal testicular tissue to facilitate the development of human IVS conditions. Testicular cells derived from hiPSCs are characterized for phenotype markers and profiled transcriptionally. HiPSC-derived testicular cells are bioprinted into core-shell constructs representative of testis cytoarchitecture and found to capture functional aspects of prepubertal testicular tissues within 7 days under xeno-free conditions. Moreover, hiPSC-derived Sertoli cells illustrate the capacity to mature under pubertal-like conditions. The utility of the model is tested by comparing 2 methods of supplementing retinoic acid (RA), the vitamin responsible for inducing spermatogenesis. The model reveals a significant gain in activity under microsphere-released RA compared to RA medium supplementation, indicating that the fragility of free RA in vitro may be a contributing factor to the molecular dysfunction observed in human IVS studies to date.

Intestinal secretory differentiation reflects niche-driven phenotypic and epigenetic plasticity of a common signal-responsive terminal cell

Enterocytes and four classic secretory cell types derive from intestinal epithelial stem cells. Based on morphology, location, and canonical markers, goblet and Paneth cells are considered distinct secretory types. Here, we report high overlap in their transcripts and sites of accessible chromatin, in marked contrast to those of their enteroendocrine or tuft cell siblings. Mouse and human goblet and Paneth cells express extraordinary fractions of few antimicrobial genes, which reflect specific responses to local niches. Wnt signaling retains some ATOH1+ secretory cells in crypt bottoms, where the absence of BMP signaling potently induces Paneth features. Cells that migrate away from crypt bottoms encounter BMPs and thereby acquire goblet properties. These phenotypes and underlying accessible cis-elements interconvert in post-mitotic cells. Thus, goblet and Paneth properties represent alternative phenotypic manifestations of a common signal-responsive terminal cell type. These findings reveal exquisite niche-dependent cell plasticity and cis-regulatory dynamics in likely response to antimicrobial needs.

Spatial transcriptomics reveals tryptophan metabolism restricting maturation of intratumoral tertiary lymphoid structures

Tertiary lymphoid structures (TLSs) are ectopic lymphoid aggregates found in numerous cancers, often linked to enhanced immunotherapy responses and better clinical outcomes. However, the factors driving TLS maturation are not fully understood. Using near single-cell spatial transcriptomic mapping, we comprehensively profile TLSs under various maturation stages and their microenvironment in hepatocellular carcinoma (HCC). Based on their developmental trajectories, we classify immature TLSs into two groups: conforming and deviating TLSs. Our findings indicate that conforming TLSs, similar to mature TLSs, possess a niche function for immunotherapy responses, while deviating TLSs do not. We discover that the tryptophan-enriched metabolic microenvironment shaped by malignant cells contributes to the deviation of TLS maturation. Inhibiting tryptophan metabolism promotes intratumoral TLS maturation and enhances tumor control, synergizing with anti-PD-1 treatments. Therefore, promoting TLS maturation represents a potential strategy to improve antitumor responses and immunotherapy outcomes.

How to view the female reproductive tract through single-cell looking glasses

Single-cell technologies have emerged as an unprecedented tool for biologists and clinicians, allowing them to assess organs and tissues at the level of individual cells. In the field of women's reproductive biology, single-cell studies have provided insights into the cellular and molecular processes that regulate reproductive and obstetrical functions in health and disease. The knowledge that these studies generate is helping clinicians to improve the understanding and diagnosis of infertility related issues or pregnancy complications and to find new avenues for their treatment. However, navigating the expansive landscape of this type of transcriptomic data analysis represents a pivotal challenge in current research. Single cell RNA sequencing involves isolating cells into droplets, reverse transcribing RNA to generate complementary DNA, with each droplet content uniquely labeled by a barcode. Upon sequencing the complementary DNAs, the barcodes enable the reassignment of sequencing reads to individual droplets, facilitating the reconstruction of the cellular landscape of the sample obtained from a tissue or organ and beyond. Researchers, equipped with the metaphorical 'single-cell glasses,' must adequately choose from a plethora of strategies to dissect and interpret cellular information. Sophisticated algorithms and the decision-making process are often underestimated, resulting in artefactual or cumbersome interpreted results. Computational biologists apply and innovate computational tools designed to process, model, and interpret expansive datasets. The ramifications of their work extend far beyond the realm of data processing; they give shape to the outcome of analyses, playing a pivotal role in drawing meaningful conclusions from the wealth of information garnered. In this review, we describe the wide variety of approaches and analytical steps available with enough detail to gain a concise picture of what a complete examination of a single-cell dataset would be. We commence with a discussion on key points in experimental design, highlighting crucial questions one should consider. Following this, we delve into the various preprocessing and quality control steps essential for any single-cell dataset. The subsequent section offers a detailed guide on constructing a single-cell atlas, exploring nuances such as differential characteristics in visualization and clustering techniques, as well as strategies for assigning identity to cell populations through gene marker annotations. Moving beyond the creation of an atlas, we explore methods for investigating pathological conditions. This involves conducting cell population comparison tests between conditions and analyzing specific cell-to-cell communications and cellular differentiation trajectories in both health and disease scenarios. This work aims to furnish a newcomer researcher and/or clinician with essential guidelines to embark on a single-cell adventure without succumbing to common pitfalls. By bridging the gap between theory and practice, it facilitates the translation of single-cell technologies into clinically relevant applications. Throughout the manuscript, practical examples of its usage in women's reproductive health studies are provided. Various sections delve into specific clinical scenarios, demonstrating how these guidelines can be instrumental in unraveling the molecular landscapes of diseases and physiological processes related to women's reproduction.

IFN-γ-Induced intestinal epithelial cell-type-specific programmed cell death: PANoptosis and its modulation in Crohn's disease

Background

Crohn's disease (CD) is a chronic inflammatory bowel disease (IBD) and is considered a Th1-mediated disease, supported by the over-expression of interferon-gamma (IFN-γ) in the intestinal lamina propria. IFN-γ has a pleiotropic effect on the intestinal epithelial cells (IECs), suggesting that IFN-γ-induced responses may differ between epithelial cell types.

Methods

We established human small intestinal organoids (enteroids) derived from non-IBD controls and CD patients. Using human enteroids, the major response of IECs induced by IFN-γ was evaluated, focusing on the IFN-γ-induced programmed cell death (PCD) pathway. Identified IFN-γ-induced responses were validated in surgically resected intestinal samples and publicly available single-cell RNA-sequencing datasets.

Results

IFN-γ stimulated programmed cell death (PCD) of IECs in both control and CD enteroids in a dose-dependent manner. Pyroptosis, apoptosis. and necroptosis were activated in enteroids, suggesting that PANoptosis was the main process of IFN-γ-induced PCD in IECs. The response to IFN-γ depends on the cell type of the IECs. IFN-γ induced depletion of enterocytes with upregulation of PANoptosis-associated genes, while leading to expansion of goblet cells without significant change in PANoptosis-associated gene expression. Individual PCD inhibitors were insufficient to block IFN-γ-induced cytotoxicity, whereas the selective JAK1 inhibitor (upadacitinib) effectively blocked IFN-γ-induced cytotoxicity and PANoptosis. Furthermore, PANoptosis was significantly activated in surgically resected tissues and in publicly available single-cell RNA-sequencing datasets of intestinal tissues from patients with CD.

Conclusion

IFN-γ induces PANoptosis in enterocytes, which can be treated with a selective JAK1 inhibitor in patients with CD.

Chromosome-level genome assembly and single-cell analysis unveil molecular mechanisms of arm regeneration in the ophiuroid Ophiura sarsii vadicola

BACKGROUND

Ophiuroids, belonging to Ophiuroidea in Echinodermata, possess remarkable regenerative capacities in their arms, relying on cellular recruitment and de-differentiation. However, limited high-quality genomic resources have hindered the investigation of the underlying molecular mechanisms of ophiuroid regeneration.

RESULTS

Here, we report a chromosome-level genome of Ophiura sarsii vadicola, 259.28 Mbp in length with a scaffold N50 length of 66.91 Mbp. We then perform bulk and single-cell RNA sequencing analysis to investigate gene expression and cellular dynamics during arm regeneration. We identify five distinct cellular clusters involved in the arm regeneration and infer the dynamic transformations from sensory stimulation to injury response, wound healing, and tissue regeneration. We find that progenitor cells derived from connective tissue cells differentiate into muscle, cartilage, endothelial, and epithelial cells. Pseudotime analysis indicates that muscle differentiation occurs early in the regeneration process.

CONCLUSIONS

Our genomic resource and single-cell atlas shed light on the mechanisms of organ regeneration in ophiuroids.

Reference-informed evaluation of batch correction for single-cell omics data with overcorrection awareness

Batch effect correction (BEC) is fundamental to integrate multiple single-cell RNA sequencing datasets, and its success is critical to empower in-depth interrogation for biological insights. However, no simple metric is available to evaluate BEC performance with sensitivity to data overcorrection, which erases true biological variations and leads to false biological discoveries. Here, we propose RBET, a reference-informed statistical framework for evaluating the success of BEC. Using extensive simulations and six real data examples including scRNA-seq and scATAC-seq datasets with different numbers of batches, batch effect sizes and numbers of cell types, we demonstrate that RBET evaluates the performance of BEC methods more fairly with biologically meaningful insights from data, while other methods may lead to false results. Moreover, RBET is computationally efficient, sensitive to overcorrection and robust to large batch effect sizes. Thus, RBET provides a robust guideline on selecting case-specific BEC method, and the concept of RBET is extendable to other modalities.

CytoSimplex: visualizing single-cell fates and transitions on a simplex

SUMMARY

Cells differentiate to their final fates along unique trajectories, often involving multi-potent progenitors that can produce multiple terminally differentiated cell types. Recent developments in single-cell transcriptomic and epigenomic measurement provide tremendous opportunities for mapping these trajectories. The visualization of single-cell data often relies on dimension reduction methods such as UMAP to simplify high-dimensional single-cell data down to an understandable 2D form. However, these dimension reduction methods are not constructed to allow direct interpretation of the reduced dimensions in terms of cell differentiation. To address these limitations, we developed a new approach that places each cell from a single-cell dataset within a simplex whose vertices correspond to terminally differentiated cell types. Our approach can quantify and visualize current cell fate commitment and future cell potential. We developed CytoSimplex, a standalone open-source package implemented in R and Python that provides simple and intuitive visualizations of cell differentiation in 2D ternary and 3D quaternary plots. We believe that CytoSimplex can help researchers gain a better understanding of cell type transitions in specific tissues and characterize developmental processes.

AVAILABILITY AND IMPLEMENTATION

The R version of CytoSimplex is available on Github at https://github.com/welch-lab/CytoSimplex. The Python version of CytoSimplex is available on Github at https://github.com/welch-lab/pyCytoSimplex.

Different fibroblast subtypes propel spatially defined ileal inflammation through TNFR1 signalling in murine ileitis

Crohn's disease (CD) is a persistent inflammatory disorder primarily affecting the terminal ileum. The TnfΔΑRE mice, which spontaneously develop CD-like ileitis due to TNF overexpression, represent a faithful model of the human disease. Here, via single-cell RNA sequencing in TnfΔΑRE mice, we show that murine TNF-dependent ileitis is characterized by cell expansion in tertiary lymphoid organs (TLO), T cell effector reprogramming, and accumulation of activated macrophages in the submucosal granulomas. Within the stromal cell compartment, fibroblast subsets (telocytes, trophocytes, PdgfraloCd81- cells) are less abundant while lymphatic endothelial cells (LEC) and fibroblastic reticular cells (FRC) show relative expansion compared to the wild type. All three fibroblast subsets show strong pro-inflammatory signature. TNFR1 loss or gain of function experiments in specific fibroblast subsets suggest that the TnfΔΑRE-induced ileitis is initiated in the lamina propria via TNF pathway activation in villus-associated fibroblasts (telocytes and PdgfraloCd81- cells), which are responsible for the organization of TLOs. Trophocytes drive disease progression in the submucosal layer, accompanied by the excessive formation of granulomas. These findings provide evidence for spatial regulation of inflammation by fibroblast subsets and underscore the pivotal role of fibroblasts in the inception and advancement of ileitis.

Cyborg organoids integrated with stretchable nanoelectronics can be functionally mapped during development

Organoids are in vitro miniaturized cellular models of organs that offer opportunities for studying organ development, disease mechanisms and drug screening. Understanding the complex processes governing organoid development and function requires methods suitable for the continuous, long-term monitoring of cell activities (for example, electrophysiological and mechanical activity) at single-cell resolution throughout the entire three-dimensional (3D) structure. Cyborg organoid technology addresses this need by seamlessly integrating stretchable mesh nanoelectronics with tissue-like properties, such as tissue-level flexibility, subcellular feature size and mesh-like networks, into 3D organoids through a 2D-to-3D tissue reconfiguration process during organogenesis. This approach enables longitudinal, tissue-wide, single-cell functional mapping, thereby overcoming the limitations of existing techniques including recording duration, spatial coverage, and the ability to maintain stable contact with the tissue during organoid development. This protocol describes the fabrication and characterization of stretchable mesh nanoelectronics, their electrical performance, their integration with organoids and the acquisition of long-term functional organoid activity requiring multimodal data analysis techniques. Cyborg organoid technology represents a transformative tool for investigating organoid development and function, with potential for improving in vitro disease models, drug screening and personalized medicine. The procedure is suitable for users within a multidisciplinary team with expertise in stem cell biology, tissue engineering, nanoelectronics fabrication, electrophysiology and data science.

TNF-α disrupts the malate-aspartate shuttle, driving metabolic rewiring in iPSC-derived enteric neural lineages from Parkinson's Disease patients

Gastrointestinal (GI) dysfunction emerges years before motor symptoms in Parkinson's disease (PD), implicating the enteric nervous system (ENS) in early disease progression. However, the mechanisms linking the PD hallmark protein, α-synuclein (α-syn), to ENS dysfunction - and whether these mechanisms are influenced by inflammation - remains elusive. Using iPSC-derived enteric neural lineages from patients with α-syn triplications, we reveal that TNF-α increases mitochondrial-α-syn interactions, disrupts the malate-aspartate shuttle, and forces a metabolic shift toward glutamine oxidation. These alterations drive mitochondrial dysfunction, characterizing metabolic impairment under cytokine stress. Interestingly, targeting glutamate metabolism with Chicago Sky Blue 6B restores mitochondrial function, reversing TNF-α-driven metabolic disruption. Our findings position the ENS as a central player in PD pathogenesis, establishing a direct link between cytokines, α-syn accumulation, metabolic stress and mitochondrial dysfunction. By uncovering a previously unrecognized metabolic vulnerability in the ENS, we highlight its potential as a therapeutic target for early PD intervention.

Not all controls are made equal: Definition of human kidney reference samples by single cell gene expression profiles

Identifying kidney disease mechanisms often requires comparing samples from disease states with healthy reference tissues. However, the effect of variations in sample procurement, storage and donor baseline characteristics of reference samples has thus far not been evaluated. Three distinct kidney reference sample types were evaluated for integrity and injury biomarkers and in their ability to define differentially expressed genes (DEGs) when compared to three different diabetic kidney disease (DKD) states. Unaffected parts of tumor nephrectomies (TN), pre-transplant living donor biopsies (LD), and percutaneous kidney research biopsies from healthy volunteers (HC) served as sources for reference tissue. Single cell gene expression profiles showed differences in the expression of injury or disease markers and the proportion of immune and proximal cell states. TN exhibited the highest expression of early stress response genes. A gene set associated with procurement effect in post-operative biopsies (LD and TN) was identified. An age-associated transcriptional signature was extracted from the reference data. Providing these tools to control for age and tissue procurement effects, immune-related pathways were found to be most enriched in DKD when compared to HC. Energy-related processes were enriched in DEGs from DKD versus LD. TN samples exhibited more underlying pathology than LD. The pathway analyses using the DEGs underscore the importance of accounting for appropriate confounding factors in differential expression analyses between disease and reference samples. Comparable controls are essential for appropriate molecular evaluation of pathologic tissues.

TRANSLATIONAL STATEMENT

Integrated single-cell data analysis of three reference sample types-needle biopsy from young healthy kidney tissue, pre-perfusion biopsy from transplant kidneys, and cancer-free tissue from tumor-nephrectomies-revealed distinct transcriptional profiles influenced by the biopsy procurement method and age. These differences impacted findings in diabetes-related kidney disease versus reference comparisons highlighting the need and providing the tools to account for these differences in interpreting analyses and identifying disease mechanisms.

A bHLH interaction code controls bipotential differentiation and self-renewal in the Drosophila gut

Multipotent adult stem cells balance self-renewal with differentiation into various cell types. How this balance is regulated at the transcriptional level is poorly understood. Here, we show that a network of basic helix-loop-helix (bHLH) transcription factors controls both stemness and bipotential differentiation in the Drosophila adult intestine. We find that homodimers of Daughterless (Da), a homolog of mammalian E proteins, maintain self-renewal of intestinal stem cells (ISCs), antagonizing the enteroendocrine fate promoted by heterodimers of Da and Scute (Sc; homolog of ASCL). The HLH factor Extramacrochaetae (Emc; homologous to Id proteins) promotes absorptive differentiation by titrating Da and Sc. Emc prevents the committed absorptive progenitor from dedifferentiating, underscoring the plasticity of these cells. Switching physical interaction partners in this way enables the active maintenance of stemness while priming stem cells for differentiation along two alternative fates. Such regulatory logic is likely operative in other bipotent stem cell systems.

Immunosuppressive JAG2+ tumor-associated neutrophils hamper PD-1 blockade response in ovarian cancer by mediating the differentiation of effector regulatory T cells

BACKGROUND

Tumor-associated neutrophils (TANs) play a critical role in modulating immune responses and exhibit significant heterogeneity. Our previous study demonstrated that jagged canonical Notch ligand 2 (JAG2)+ TANs were associated with an immunosuppressive microenvironment in high-grade serous ovarian cancer (HGSOC), but the underlying mechanism remains unclear. This study aimed to elucidate the role of JAG2+ TANs in tumor immunosuppressive microenvironment in HGSOC.

METHODS

HGSOC samples were collected, with 274 samples constituting two independent cohorts (training and validation cohorts) and an additional 30 samples utilized to establish patient-derived tumor organoids (PDTOs). We characterized the number and phenotype of JAG2+ TANs by multiplex immunohistochemistry, flow cytometry, and single-cell RNA sequencing (scRNA-seq). We investigated the biological functions of JAG2 in immune evasion using in vitro co-culture systems, flow cytometry, tumor-bearing mouse models, and PDTOs.

RESULTS

JAG2+ TANs expressed elevated levels of immunosuppressive molecules, including programmed cell death ligand 1 and CD14, and had independent prognostic value for the overall survival of patients with HGSOC. scRNA-seq analysis revealed that JAG2+ TANs exhibited a terminally mature phenotype. The infiltration of JAG2+ TANs was positively correlated with the abundance of effector regulatory T cells (eTregs). Interaction with JAG2+ TANs skewed CD4+ T cells towards an eTreg phenotype, a process that was suppressed by the Notch inhibitor LY3039478 and induced by recombinant Jagged2. Furthermore, we demonstrated that JAG2+ TANs enhanced Notch signaling activation, ultimately promoting recombination signal binding protein for immunoglobulin kappa J region (RBPJ)-induced differentiation of naïve CD4+ T cells into eTregs. Clinically, JAG2+ TANs could serve as a biomarker for assessing immunotherapy resistance in various solid tumors. Pharmacological targeting of Notch signaling with LY3039478 or JAG2 neutralization antibodies enhanced the efficacy of programmed cell death protein 1 (PD-1) monoclonal antibodies (mAbs) in both xenograft and PDTO models.

CONCLUSIONS

The emergence of JAG2+ TANs is crucial for the differentiation of eTregs, which triggers immune evasion and resistance to anti-PD-1 therapy. Inhibiting Notch signaling with LY3039478 or JAG2 neutralization antibodies may overcome this anti-PD-1 resistance in HGSOC.