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

HMGB2 drives tumor progression and shapes the immunosuppressive microenvironment in hepatocellular carcinoma: insights from multi-omics analysis

Background

Hepatocellular carcinoma (HCC) poses a significant health burden globally, with high mortality rates despite various treatment options. Immunotherapy, particularly immune-checkpoint inhibitors (ICIs), has shown promise, but resistance and metastasis remain major challenges. Understanding the intricacies of the tumor microenvironment (TME) is imperative for optimizing HCC management strategies and enhancing patient prognosis.

Methods

This study employed a comprehensive approach integrating multi-omics approaches, including single-cell RNA sequencing (scRNA-seq), bulk RNA sequencing (Bulk RNA-seq), and validation in clinical samples using spatial transcriptomics (ST) and multiplex immunohistochemistry (mIHC). The analysis aimed to identify key factors influencing the immunosuppressive microenvironment associated with HCC metastasis and immunotherapy resistance.

Results

HMGB2 is significantly upregulated in HCCTrans, a transitional subgroup associated with aggressive metastasis. Furthermore, HMGB2 expression positively correlates with an immunosuppressive microenvironment, particularly evident in exhausted T cells. Notably, HMGB2 expression correlated positively with immunosuppressive markers and poor prognosis in HCC patients across multiple cohorts. ST combined with mIHC validated the spatial expression patterns of HMGB2 within the TME, providing additional evidence of its role in HCC progression and immune evasion.

Conclusion

HMGB2 emerges as a critical player of HCC progression, metastasis, and immunosuppression. Its elevated expression correlates with aggressive tumor behavior and poor patient outcomes, suggesting its potential as both a therapeutic target and a prognostic indicator in HCC management.

Microglia protect against age-associated brain pathologies

Microglia are brain-resident macrophages that contribute to central nervous system (CNS) development, maturation, and preservation. Here, we examine the consequences of permanent microglial deficiencies on brain aging using the Csf1rΔFIRE/ΔFIRE mouse model. In juvenile Csf1rΔFIRE/ΔFIRE mice, we show that microglia are dispensable for the transcriptomic maturation of other brain cell types. By contrast, with advancing age, pathologies accumulate in Csf1rΔFIRE/ΔFIRE brains, macroglia become increasingly dysregulated, and white matter integrity declines, mimicking many pathological features of human CSF1R-related leukoencephalopathy. The thalamus is particularly vulnerable to neuropathological changes in the absence of microglia, with atrophy, neuron loss, vascular alterations, macroglial dysregulation, and severe tissue calcification. We show that populating Csf1rΔFIRE/ΔFIRE brains with wild-type microglia protects against many of these pathological changes. Together with the accompanying study by Chadarevian and colleagues1, our results indicate that the lifelong absence of microglia results in an age-related neurodegenerative condition that can be counteracted via transplantation of healthy microglia.

The Impact of Low Protein Diet on the Molecular and Cellular Development of the Fetal Kidney

Background

Low nephron number has a direct impact on the development of hypertension and chronic kidney disease later in life. While intrauterine growth restriction caused by maternal low protein diet (LPD) is thought to be a significant cause of reduced nephron endowment in impoverished communities, its influence on the cellular and molecular processes which drive nephron formation are poorly understood.

Methods

We conducted a comprehensive characterization of the impact of LPD on kidney development using tomographic and confocal imaging to quantify changes in branching morphogenesis and the cellular and morphological features of nephrogenic niches across development. These analyses were paired with single-cell RNA sequencing to dissect the transcriptional changes that LPD imposes during renal development to affect nephron number.

Results

Single cell analysis at E14.5 and P0 revealed differences in the expression of genes and pathways involved in metabolism, cell cycle, epigenetic regulators and reciprocal inductive signals in most cell types analyzed, yielding imbalances and shifts in cellular energy production and cellular trajectories. In the nephron progenitor cells, LPD impeded cellular commitment and differentiation towards pre-tubular and renal vesicle structures. Confocal microscopy revealed a reduction in the number of pre-tubular aggregates and proliferation in nephron progenitor cells. We also found changes in branching morphogenesis, with a reduction in cell proliferation in the ureteric tips as well as reduced tip and tip parent lengths by optical projection tomography which causes patterning defects.

Conclusions

This unique profiling demonstrates how a fetal programming defect leads to low nephron endowment which is intricately linked to changes in both branching morphogenesis and the commitment of nephron progenitor cells. The commitment of progenitor cells is pivotal for nephron formation and is significantly influenced by nutritional factors, with a low protein diet driving alterations in this program which directly results in a reduced nephron endowment.

Significance Statement

While a mother's diet can negatively impact the number of nephrons in the kidneys of her offspring, the root cellular and molecular drivers of these deficits have not been rigorously explored. In this study we use advanced imaging and gene expression analysis in mouse models to define how a maternal low protein diet, analogous to that of impoverished communities, results in reduced nephron endowment. We find that low protein diet has pleiotropic effects on metabolism and the normal developmental programs of gene expression. These profoundly impact the process of branching morphogenesis necessary to establish niches for nephron generation and change cell behaviors which regulate how and when nephron progenitor cells commit to differentiation.

Elucidating the role of tumor-associated ALOX5+ mast cells with transformative function in cervical cancer progression via single-cell RNA sequencing

Background

Cervical cancer (CC) is the fourth most common malignancy among women globally and serves as the main cause of cancer-related deaths among women in developing countries. The early symptoms of CC are often not apparent, with diagnoses typically made at advanced stages, which lead to poor clinical prognoses. In recent years, numerous studies have shown that there is a close relationship between mast cells (MCs) and tumor development. However, research on the role MCs played in CC is still very limited at that time. Thus, the study conducted a single-cell multi-omics analysis on human CC cells, aiming to explore the mechanisms by which MCs interact with the tumor microenvironment in CC. The goal was to provide a scientific basis for the prevention, diagnosis, and treatment of CC, with the hope of improving patients' prognoses and quality of life.

Method

The present study acquired single-cell RNA sequencing data from ten CC tumor samples in the ArrayExpress database. Slingshot and AUCcell were utilized to infer and assess the differentiation trajectory and cell plasticity of MCs subpopulations. Differential expression analysis of MCs subpopulations in CC was performed, employing Gene Ontology, gene set enrichment analysis, and gene set variation analysis. CellChat software package was applied to predict cell communication between MCs subpopulations and CC cells. Cellular functional experiments validated the functionality of TNFRSF12A in HeLa and Caski cell lines. Additionally, a risk scoring model was constructed to evaluate the differences in clinical features, prognosis, immune infiltration, immune checkpoint, and functional enrichment across various risk scores. Copy number variation levels were computed using inference of copy number variations.

Result

The obtained 93,524 high-quality cells were classified into ten cell types, including T_NK cells, endothelial cells, fibroblasts, smooth muscle cells, epithelial cells, B cells, plasma cells, MCs, neutrophils, and myeloid cells. Furthermore, a total of 1,392 MCs were subdivided into seven subpopulations: C0 CTSG+ MCs, C1 CALR+ MCs, C2 ALOX5+ MCs, C3 ANXA2+ MCs, C4 MGP+ MCs, C5 IL32+ MCs, and C6 ADGRL4+ MCs. Notably, the C2 subpopulation showed close associations with tumor-related MCs, with Slingshot results indicating that C2 subpopulation resided at the intermediate-to-late stage of differentiation, potentially representing a crucial transition point in the benign-to-malignant transformation of CC. CNVscore and bulk analysis results further confirmed the transforming state of the C2 subpopulation. CellChat analysis revealed TNFRSF12A as a key receptor involved in the actions of C2 ALOX5+ MCs. Moreover, in vitro experiments indicated that downregulating the TNFRSF12A gene may partially inhibit the development of CC. Additionally, a prognosis model and immune infiltration analysis based on the marker genes of the C2 subpopulation provided valuable guidance for patient prognosis and clinical intervention strategies.

Conclusions

We first identified the transformative tumor-associated MCs subpopulation C2 ALOX5+ MCs within CC, which was at a critical stage of tumor differentiation and impacted the progression of CC. In vitro experiments confirmed the inhibitory effect of knocking down the TNFRSF12A gene on the development of CC. The prognostic model constructed based on the C2 ALOX5+MCs subset demonstrated excellent predictive value. These findings offer a fresh perspective for clinical decision-making in CC.

Nucleotide depletion promotes cell fate transitions by inducing DNA replication stress

Control of cellular identity requires coordination of developmental programs with environmental factors such as nutrient availability, suggesting that perturbing metabolism can alter cell state. Here, we find that nucleotide depletion and DNA replication stress drive differentiation in human and murine normal and transformed hematopoietic systems, including patient-derived acute myeloid leukemia (AML) xenografts. These cell state transitions begin during S phase and are independent of ATR/ATM checkpoint signaling, double-stranded DNA break formation, and changes in cell cycle length. In systems where differentiation is blocked by oncogenic transcription factor expression, replication stress activates primed regulatory loci and induces lineage-appropriate maturation genes despite the persistence of progenitor programs. Altering the baseline cell state by manipulating transcription factor expression causes replication stress to induce genes specific for alternative lineages. The ability of replication stress to selectively activate primed maturation programs across different contexts suggests a general mechanism by which changes in metabolism can promote lineage-appropriate cell state transitions.

Respiratory infection with influenza A virus delays remyelination and alters oligodendrocyte metabolism

Peripheral viral infection disrupts oligodendrocyte (OL) homeostasis such that endogenous remyelination may be affected. Here, we demonstrate that influenza A virus infection perpetuated a demyelination- and disease-associated OL phenotype following cuprizone-induced demyelination that resulted in delayed OL maturation and remyelination in the prefrontal cortex. Furthermore, we assessed cellular metabolism ex vivo, and found that infection altered brain OL and microglia metabolism in a manner that opposed the metabolic profile induced by remyelination. Specifically, infection increased glycolytic capacity of OLs and microglia, an effect that was recapitulated by lipopolysaccharide (LPS) stimulation of mixed glia cultures. In contrast, mitochondrial dependence was increased in OLs during remyelination, which was similarly observed in OLs of myelinating P14 mice compared to adult and aged mice. Collectively, our data indicate that respiratory viral infection is capable of suppressing remyelination, and suggest that metabolic dysfunction of OLs is implicated in remyelination impairment.

Single-cell Transcriptional Landscape of Temporal Neutrophil Response to Burn Wound in Larval Zebrafish

Neutrophils accumulate early in tissue injury. However, the cellular and functional heterogeneity of neutrophils during homeostasis and in response to tissue damage remains unclear. In this study, we use larval zebrafish to understand neutrophil responses to thermal injury. Single-cell transcriptional mapping of myeloid cells during a 3-d time course in burn and control larvae revealed distinct neutrophil subsets and their cell-cell interactions with macrophages across time and conditions. The trajectory formed by three zebrafish neutrophil subsets resembles human neutrophil maturation, with varying transition patterns between conditions. Through ligand-receptor cell-cell interaction analysis, we found that neutrophils communicate more in burns in a pathway and temporal manner. Finally, we identified the correlation between zebrafish myeloid signatures and human burn severity, establishing GPR84+ neutrophils as a potential marker of early innate immune response in burns. This work builds a comparative single-cell transcriptomic framework to identify neutrophil markers of tissue damage using model organisms.

Associating transcription factors to single-cell trajectories with DREAMIT

Inferring gene regulatory networks from single-cell RNA-sequencing trajectories has been an active area of research yet methods are still needed to identify regulators governing cell transitions. We developed DREAMIT (Dynamic Regulation of Expression Across Modules in Inferred Trajectories) to annotate transcription-factor activity along single-cell trajectory branches, using ensembles of relations to target genes. Using a benchmark representing several different tissues, as well as external validation with ATAC-Seq and Perturb-Seq data on hematopoietic cells, the method was found to have higher tissue-specific sensitivity and specificity over competing approaches.

Spatially Resolved Niche and Tumor Microenvironmental Alterations in Gastric Cancer Peritoneal Metastases

BACKGROUND & AIMS

Peritoneal metastasis (PM) in gastric cancer (GC) is associated with poor prognosis and significant morbidity. We sought to understand the genomic, transcriptomic, and tumor microenvironment (TME) features that contribute to peritoneal organotropism in GC.

METHODS

We conducted a comprehensive multi-omic analysis of 548 samples from 326 patients, including primary tumors, matched normal tissues; peritoneal metastases, and adjacent-normal peritoneal tissues. We used whole exome sequencing, whole transcriptome sequencing, and digital spatial profiling to investigate molecular alterations, gene expression patterns, and TME characteristics associated with PM.

RESULTS

Our analysis identified specific genomic alterations in primary tumors, including mutations in ELF3, CDH1, and PIGR, and TME signatures, such as stromal infiltration and M2 macrophage enrichment, associated with increased risk of PM. We observed distinct transcriptional programs and immune compositions in GCPM compared with liver metastases, highlighting the importance of the TME in transcoelomic metastasis. We found differential expression of therapeutic targets between primary tumors and PM, with lower CLDN18.2 and FGFR2b expression in PM. We unravel the roles of the TME in niche reprogramming within the peritoneum, and provide evidence of pre-metastatic niche conditioning even in early GC without clinical PM. These findings were further validated using a humanized mouse model, which demonstrated niche remodeling in the peritoneum during transcoelomic metastasis.

CONCLUSION

Our study provides a comprehensive molecular characterization of GCPM and unveils key biological principles underlying transcoelomic metastasis. The identified predictive markers, therapeutic targets, and TME alterations offer potential avenues for targeted interventions and improved patient outcomes.

Inference of genomic lesions from single-cell RNA-seq in myeloma improves functional intraclonal and interclonal analysis

ABSTRACT

Smoldering multiple myeloma (SMM) is an asymptomatic plasma cell (PC) neoplasm that may evolve with variable frequency into multiple myeloma (MM). SMM is initiated by chromosomal translocations involving the immunoglobulin heavy-chain locus or by hyperdiploidy and evolves through acquisition of additional genetic lesions. In this scenario, we aimed at establishing a reliable analysis pipeline to infer genomic lesions from transcriptomic analysis, by combining single-cell RNA sequencing (scRNA-seq) with B-cell receptor sequencing and copy number abnormality (CNA) analysis to identify clonal PCs at the genetic level along their specific transcriptional landscape. We profiled 20 465 bone marrow PCs derived from 5 patients with SMM/MM and unbiasedly identified clonal and polyclonal PCs. Hyperdiploidy, t(11;14), and t(6;14) were identified at the scRNA level by analysis of chimeric reads. Subclone functional analysis was improved by combining transcriptome with CNA analysis. As examples, we illustrate the different functional properties of a light-chain escape subclone in SMM and of different B-cell and PC subclones in a patient affected by Wäldenstrom macroglobulinemia and SMM. Overall, our data provide a proof of principle for inference of clinically relevant genotypic data from scRNA-seq, which in turn will refine functional annotation of the clonal architecture of PC dyscrasias.

Building a learnable universal coordinate system for single-cell atlas with a joint-VAE model

A universal coordinate system that can ensemble the huge number of cells and capture their heterogeneities is of vital importance for constructing large-scale cell atlases as references for molecular and cellular studies. Studies have shown that cells exhibit multifaceted heterogeneities in their transcriptomic features at multiple resolutions. This nature of complexity makes it hard to design a fixed coordinate system through a combination of known features. It is desirable to build a learnable universal coordinate model that can capture major heterogeneities and serve as a controlled generative model for data augmentation. We developed UniCoord, a specially-tuned joint-VAE model to represent single-cell transcriptomic data in a lower-dimensional latent space with high interpretability. Each latent dimension can represent either discrete or continuous feature, and either supervised by prior knowledge or unsupervised. The latent dimensions can be easily reconfigured to generate pseudo transcriptomic profiles with desired properties. UniCoord can also be used as a pre-trained model to analyze new data with unseen cell types and thus can serve as a feasible framework for cell annotation and comparison. UniCoord provides a prototype for a learnable universal coordinate framework to enable better analysis and generation of cells with highly orchestrated functions and heterogeneities.

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

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

Deciphering tumour microenvironment and elucidating the origin of cancer cells in ovarian clear cell carcinoma

Ovarian clear cell carcinoma (CCC) has an East Asian preponderance. It is associated with endometriosis, a benign condition where endometrial (inner lining of the uterus) tissue is found outside the uterus and on the peritoneal surface, in the abdominal or pelvic space. CCC is relatively more resistant to conventional chemotherapy compared to other ovarian cancer subtypes and is associated with a poorer prognosis. In this study, we recruited and obtained tumour tissues from seven patients across the four stages of CCC. The tumour and the tumour microenvironment (TME) from 7 CCC patients spanning clinical stages 1-4 were transcriptionally profiled using high-resolution scRNA-seq to gain insight into CCC's biological mechanisms. Firstly, we built a scRNA-seq resource for the CCC tumour microenvironment (TME). Secondly, we identified the different cell type proportions and found high levels of immune infiltration in CCC. Thirdly, since CCC is associated with endometriosis, we compared CCC with two publicly available endometriosis scRNA-seq datasets. The CCC malignant cells showed similarities with glandular secretory and ciliated epithelial cells found in endometriosis. Finally, we determined the differences in cell-cell communication between various cell types present in CCC TME and endometriosis conditions to gain insights into the transformations in CCC.

Multi-omic analysis of Huntington's disease reveals a compensatory astrocyte state

The mechanisms underlying the selective regional vulnerability to neurodegeneration in Huntington's disease (HD) have not been fully defined. To explore the role of astrocytes in this phenomenon, we used single-nucleus and bulk RNAseq, lipidomics, HTT gene CAG repeat-length measurements, and multiplexed immunofluorescence on HD and control post-mortem brains. We identified genes that correlated with CAG repeat length, which were enriched in astrocyte genes, and lipidomic signatures that implicated poly-unsaturated fatty acids in sensitizing neurons to cell death. Because astrocytes play essential roles in lipid metabolism, we explored the heterogeneity of astrocytic states in both protoplasmic and fibrous-like (CD44+) astrocytes. Significantly, one protoplasmic astrocyte state showed high levels of metallothioneins and was correlated with the selective vulnerability of distinct striatal neuronal populations. When modeled in vitro, this state improved the viability of HD-patient-derived spiny projection neurons. Our findings uncover key roles of astrocytic states in protecting against neurodegeneration in HD.

Pregnancy Reduces Il33+ Hybrid Progenitor Accumulation in the Aged Mammary Gland

Aging increases breast cancer risk while an early first pregnancy reduces a woman's life-long risk. Several studies have explored the effect of either aging or pregnancy on mammary epithelial cells (MECs), but the combined effect of both remains unclear. Here, we interrogate the functional and transcriptomic changes at single cell resolution in the mammary gland of aged nulliparous and parous mice to discover that pregnancy normalizes age-related imbalances in lineage composition, while also inducing a differentiated cell state. Importantly, we uncover a minority population of Il33-expressing hybrid MECs with high cellular potency that accumulate in aged nulliparous mice but is significantly reduced in aged parous mice. Functionally, IL33 treatment of basal, but not luminal, epithelial cells from young mice phenocopies aged nulliparous MECs and promotes formation of organoids with Trp53 knockdown. Collectively, our study demonstrates that pregnancy blocks the age-associated loss of lineage integrity in the basal layer through a decrease in Il33+ hybrid MECs, potentially contributing to pregnancy-induced breast cancer protection.

Deciphering deep-sea chemosynthetic symbiosis by single-nucleus RNA-sequencing

Bathymodioline mussels dominate deep-sea methane seep and hydrothermal vent habitats and obtain nutrients and energy primarily through chemosynthetic endosymbiotic bacteria in the bacteriocytes of their gill. However, the molecular mechanisms that orchestrate mussel host-symbiont interactions remain unclear. Here, we constructed a comprehensive cell atlas of the gill in the mussel Gigantidas platifrons from the South China Sea methane seeps (1100 m depth) using single-nucleus RNA-sequencing (snRNA-seq) and whole-mount in situ hybridisation. We identified 13 types of cells, including three previously unknown ones, and uncovered unknown tissue heterogeneity. Every cell type has a designated function in supporting the gill's structure and function, creating an optimal environment for chemosynthesis, and effectively acquiring nutrients from the endosymbiotic bacteria. Analysis of snRNA-seq of in situ transplanted mussels clearly showed the shifts in cell state in response to environmental oscillations. Our findings provide insight into the principles of host-symbiont interaction and the bivalves' environmental adaption mechanisms.

Recruited atypical Ly6G+ macrophages license alveolar regeneration after lung injury

The lung is constantly exposed to airborne pathogens and particles that can cause alveolar damage. Hence, appropriate repair responses are essential for gas exchange and life. Here, we deciphered the spatiotemporal trajectory and function of an atypical population of macrophages after lung injury. Post-influenza A virus (IAV) infection, short-lived monocyte-derived Ly6G-expressing macrophages (Ly6G+ Macs) were recruited to the alveoli of lung perilesional areas. Ly6G+ Macs engulfed immune cells, exhibited a high metabolic potential, and clustered with alveolar type 2 epithelial cells (AT2s) in zones of active epithelial regeneration. Ly6G+ Macs were partially dependent on granulocyte-macrophage colony-stimulating factor and interleukin-4 receptor signaling and were essential for AT2-dependent alveolar regeneration. Similar macrophages were recruited in other models of injury and in the airspaces of lungs from patients with suspected pneumonia. This study identifies perilesional alveolar Ly6G+ Macs as a spatially restricted, short-lived macrophage subset promoting epithelial regeneration postinjury, thus representing an attractive therapeutic target for treating lung damage.

Navigating internal senses: A road map for the vagal interoceptive system

A road map for the vagal interoceptive system.

Revealing the crucial roles of suppressive immune microenvironment in cardiac myxoma progression

Cardiac myxoma is a commonly encountered tumor within the heart that has the potential to be life-threatening. However, the cellular composition of this condition is still not well understood. To fill this gap, we analyzed 75,641 cells from cardiac myxoma tissues based on single-cell sequencing. We defined a population of myxoma cells, which exhibited a resemblance to fibroblasts, yet they were distinguished by an increased expression of phosphodiesterases and genes associated with cell proliferation, differentiation, and adhesion. The clinical relevance of the cell populations indicated a higher proportion of myxoma cells and M2-like macrophage infiltration, along with their enhanced spatial interaction, were found to significantly contribute to the occurrence of embolism. The immune cells surrounding the myxoma exhibit inhibitory characteristics, with impaired function of T cells characterized by the expression of GZMK and TOX, along with a substantial infiltration of tumor-promoting macrophages expressed growth factors such as PDGFC. Furthermore, in vitro co-culture experiments showed that macrophages promoted the growth of myxoma cells significantly. In summary, this study presents a comprehensive single-cell atlas of cardiac myxoma, highlighting the heterogeneity of myxoma cells and their collaborative impact on immune cells. These findings shed light on the complex pathobiology of cardiac myxoma and present potential targets for intervention.

Cellular proliferation biases clonal lineage tracing and trajectory inference

MOTIVATION

Lineage tracing and trajectory inference from single-cell RNA-sequencing data hold tremendous potential for uncovering the genetic programs driving development and disease. Single cell datasets are thought to provide an unbiased view on the diverse cellular architecture of tissues. Sampling bias, however, can skew single cell datasets away from the cellular composition they are meant to represent.

RESULTS

We demonstrate a novel form of sampling bias, caused by a statistical phenomenon related to repeated sampling from a growing, heterogeneous population. Relative growth rates of cells influence the probability that they will be sampled in clones observed across multiple time points. We support our probabilistic derivations with a simulation study and an analysis of a real time-course of T-cell development. We find that this bias can impact fate probability predictions, and we explore how to develop trajectory inference methods which are robust to this bias.

AVAILABILITY AND IMPLEMENTATION

Source code for the simulated datasets and to create the figures in this manuscript is freely available in python at https://github.com/rbonhamcarter/simulate-clones. A python implementation of the extension of the LineageOT method is freely available at https://github.com/rbonhamcarter/LineageOT/tree/multi-time-clones.

A human pluripotent stem cell-based somitogenesis model using microfluidics

Emerging human pluripotent stem cell (hPSC)-based embryo models are useful for studying human embryogenesis. Particularly, there are hPSC-based somitogenesis models using free-floating culture that recapitulate somite formation. Somitogenesis in vivo involves intricately orchestrated biochemical and biomechanical events. However, none of the current somitogenesis models controls biochemical gradients or biomechanical signals in the culture, limiting their applicability to untangle complex biochemical-biomechanical interactions that drive somitogenesis. Herein, we develop a human somitogenesis model by confining hPSC-derived presomitic mesoderm (PSM) tissues in microfabricated trenches. Exogenous microfluidic morphogen gradients imposed on the PSM tissues cause axial patterning and trigger spontaneous rostral-to-caudal somite formation. A mechanical theory is developed to explain the size dependency between somites and the PSM. The microfluidic somitogenesis model is further exploited to reveal regulatory roles of cellular and tissue biomechanics in somite formation. This study presents a useful microengineered, hPSC-based model for understanding the biochemical and biomechanical events that guide somite formation.

A Cell Cycle-Aware Network for Data Integration and Label Transferring of Single-Cell RNA-Seq and ATAC-Seq

In recent years, the integration of single-cell multi-omics data has provided a more comprehensive understanding of cell functions and internal regulatory mechanisms from a non-single omics perspective, but it still suffers many challenges, such as omics-variance, sparsity, cell heterogeneity, and confounding factors. As it is known, the cell cycle is regarded as a confounder when analyzing other factors in single-cell RNA-seq data, but it is not clear how it will work on the integrated single-cell multi-omics data. Here, a cell cycle-aware network (CCAN) is developed to remove cell cycle effects from the integrated single-cell multi-omics data while keeping the cell type-specific variations. This is the first computational model to study the cell-cycle effects in the integration of single-cell multi-omics data. Validations on several benchmark datasets show the outstanding performance of CCAN in a variety of downstream analyses and applications, including removing cell cycle effects and batch effects of scRNA-seq datasets from different protocols, integrating paired and unpaired scRNA-seq and scATAC-seq data, accurately transferring cell type labels from scRNA-seq to scATAC-seq data, and characterizing the differentiation process from hematopoietic stem cells to different lineages in the integration of differentiation data.

Unraveling the connection between Hashimoto's Thyroiditis and non-alcoholic fatty liver disease: exploring the role of CD4+central memory T cells through integrated genetic approaches

PURPOSE

Exploring the connection between Hashimoto's thyroiditis (HT) and non-alcoholic fatty liver disease (NAFLD) through integrated genetic approaches.

METHODS

We utilized integrated genetic approaches, such as single-cell RNA sequencing (scRNA-seq) data analysis, Mendelian Randomization (MR), colocalization analysis, cell communication, and metabolic analyses, to investigate potential correlations between HT and NAFLD.

RESULTS

Through the integrated analysis of scRNA-seq data from individuals with HT, NAFLD, and healthy controls, we observed an upregulation in the proportion of CD4+central memory (CD4+CM) T cells among T cells in both diseases. A total of 63 differentially expressed genes (DEGs) were identified in the CD4+CM cells after the differential analysis. By using MR, 8 DEGs (MAGI3, CSGALNACT1, CAMK4, GRIP1, TRAT1, IL7R, ERN1, and MB21D2) were identified to have a causal relationship with HT, and 4 DEGs (MAGI3, RCAN3, DOCK10, and SAMD12) had a causal relationship with NAFLD. MAGI3 was found to be causally linked to both HT and NAFLD. Therefore, MAGI3 was designated as the marker gene. Reverse MR and Steiger filtering showed no evidence of reverse causality. Colocalization analyses further indicated close links between MAGI3 and HT as well as NAFLD. Finally, based on the expression levels of MAGI3, we stratified CD4+CM cells into two subsets: MAGI3+CD4+CM cells and MAGI3-CD4+CM cells. Functional analyses revealed significant differences between the two subsets, potentially related to the progression of the two diseases.

CONCLUSION

This study delves into the potential connections between HT and NAFLD through integrated genetic methods. Our research reveals an elevated proportion of CD4+CM cells within T cells in both HT and NAFLD. Through MR and colocalization analysis, we identify specific genes causally linked to HT and NAFLD, such as MAGI3. Ultimately, based on MAGI3 expression levels, we categorize CD4+CM cells into MAGI3+CD4+CM cells and MAGI3-CD4+CM cells, uncovering significant differences between them through functional analyses.

Immunosuppression and phenotypic plasticity in an atlas of human hepatocholangiocarcinoma

Background

Hepatocholangiocarcinoma (H-ChC) has the clinicopathological features of both hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) and is a more aggressive subtype of primary hepatic carcinoma than HCC or iCCA.

Methods

We sequenced 91,112 single-cell transcriptomes from 16 human samples to elucidate the molecular mechanisms underlying the coexistence of HCC and iCCA components in H-ChC.

Results

We observed two molecular subtypes of H-ChC at the whole-transcriptome level (CHP and CIP), where a metabolically active tumour cell subpopulation enriched in CHP was characterized by a cellular pre-differentiation property. To define the heterogeneity of tumours and their associated microenvironments, we observe greater tumour diversity in H-ChC than HCC and iCCA. H-ChC exhibits weaker immune cell infiltration and greater CD8+ exhausted T cell (Tex) dysfunction than HCC and iCCA. Then we defined two broad cell states of 6,852 CD8+ Tex cells: GZMK+ CD8+ Tex cells and terminal CD8+ Tex cells. GZMK+ CD8+ Tex cells exhibited higher infiltration of after treatment in H-ChC, the effector scores and expression of the immune checkpoints of them greatly increased after immunotherapy, which indicated that H-ChC might be more sensitive than HCC or iCCA to immunotherapy.

Conclusions

In this paper, H-ChC was explored, hoping to contribute to the study of mixed tumours in other cancers.

Tumour vasculature at single-cell resolution

Tumours can obtain nutrients and oxygen required to progress and metastasize through the blood supply1. Inducing angiogenesis involves the sprouting of established vessel beds and their maturation into an organized network2,3. Here we generate a comprehensive atlas of tumour vasculature at single-cell resolution, encompassing approximately 200,000 cells from 372 donors representing 31 cancer types. Trajectory inference suggested that tumour angiogenesis was initiated from venous endothelial cells and extended towards arterial endothelial cells. As neovascularization elongates (through angiogenic stages SI, SII and SIII), APLN+ tip cells at the SI stage (APLN+ TipSI) advanced to TipSIII cells with increased Notch signalling. Meanwhile, stalk cells, following tip cells, transitioned from high chemokine expression to elevated TEK (also known as Tie2) expression. Moreover, APLN+ TipSI cells not only were associated with disease progression and poor prognosis but also hold promise for predicting response to anti-VEGF therapy. Lymphatic endothelial cells demonstrated two distinct differentiation lineages: one responsible for lymphangiogenesis and the other involved in antigen presentation. In pericytes, endoplasmic reticulum stress was associated with the proangiogenic BASP1+ matrix-producing pericytes. Furthermore, intercellular communication analysis showed that neovascular endothelial cells could shape an immunosuppressive microenvironment conducive to angiogenesis. This study depicts the complexity of tumour vasculature and has potential clinical significance for anti-angiogenic therapy.

Mapping spatially resolved transcriptomes in human and mouse pulmonary fibrosis

Idiopathic pulmonary fibrosis (IPF) is a progressive lung disease with poor prognosis and limited treatment options. Efforts to identify effective treatments are thwarted by limited understanding of IPF pathogenesis and poor translatability of available preclinical models. Here we generated spatially resolved transcriptome maps of human IPF (n = 4) and bleomycin-induced mouse pulmonary fibrosis (n = 6) to address these limitations. We uncovered distinct fibrotic niches in the IPF lung, characterized by aberrant alveolar epithelial cells in a microenvironment dominated by transforming growth factor beta signaling alongside predicted regulators, such as TP53 and APOE. We also identified a clear divergence between the arrested alveolar regeneration in the IPF fibrotic niches and the active tissue repair in the acutely fibrotic mouse lung. Our study offers in-depth insights into the IPF transcriptional landscape and proposes alveolar regeneration as a promising therapeutic strategy for IPF.

Single-cell RNA-seq reveals the transcriptional program underlying tumor progression and metastasis in neuroblastoma

Neuroblastoma (NB) is one of the most common childhood malignancies. Sixty percent of patients present with widely disseminated clinical signs at diagnosis and exhibit poor outcomes. However, the molecular mechanisms triggering NB metastasis remain largely uncharacterized. In this study, we generated a transcriptomic atlas of 15 447 NB cells from eight NB samples, including paired samples of primary tumors and bone marrow metastases. We used time-resolved analysis to chart the evolutionary trajectory of NB cells from the primary tumor to the metastases in the same patient and identified a common 'starter' subpopulation that initiates tumor development and metastasis. The 'starter' population exhibited high expression levels of multiple cell cycle-related genes, indicating the important role of cell cycle upregulation in NB tumor progression. In addition, our evolutionary trajectory analysis demonstrated the involvement of partial epithelial-to-mesenchymal transition (p-EMT) along the metastatic route from the primary site to the bone marrow. Our study provides insights into the program driving NB metastasis and presents a signature of metastasis-initiating cells as an independent prognostic indicator and potential therapeutic target to inhibit the initiation of NB metastasis.

Single-cell transcriptomic analysis of B cells reveals new insights into atypical memory B cells in COVID-19

Here, we performed single-cell RNA sequencing of S1 and receptor binding domain protein-specific B cells from convalescent COVID-19 patients with different clinical manifestations. This study aimed to evaluate the role and developmental pathway of atypical memory B cells (MBCs) in response to severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection. The results revealed a proinflammatory signature across B cell subsets associated with disease severity, as evidenced by the upregulation of genes such as GADD45B, MAP3K8, and NFKBIA in critical and severe individuals. Furthermore, the analysis of atypical MBCs suggested a developmental pathway similar to that of conventional MBCs through germinal centers, as indicated by the expression of several genes involved in germinal center processes, including CXCR4, CXCR5, BCL2, and MYC. Additionally, the upregulation of genes characteristic of the immune response in COVID-19, such as ZFP36 and DUSP1, suggested that the differentiation and activation of atypical MBCs may be influenced by exposure to SARS-CoV-2 and that these genes may contribute to the immune response for COVID-19 recovery. Our study contributes to a better understanding of atypical MBCs in COVID-19 and the role of other B cell subsets across different clinical manifestations.

A group 3 medulloblastoma stem cell program is maintained by OTX2-mediated alternative splicing

OTX2 is a transcription factor and known driver in medulloblastoma (MB), where it is amplified in a subset of tumours and overexpressed in most cases of group 3 and group 4 MB. Here we demonstrate a noncanonical role for OTX2 in group 3 MB alternative splicing. OTX2 associates with the large assembly of splicing regulators complex through protein-protein interactions and regulates a stem cell splicing program. OTX2 can directly or indirectly bind RNA and this may be partially independent of its DNA regulatory functions. OTX2 controls a pro-tumorigenic splicing program that is mirrored in human cerebellar rhombic lip origins. Among the OTX2-regulated differentially spliced genes, PPHLN1 is expressed in the most primitive rhombic lip stem cells, and targeting PPHLN1 splicing reduces tumour growth and enhances survival in vivo. These findings identify OTX2-mediated alternative splicing as a major determinant of cell fate decisions that drive group 3 MB progression.

B-BIND: Biophysical Bayesian Inference for Neurodegenerative Dynamics

Throughout an organism's life, a multitude of complex and interdependent biological systems transition through biophysical processes that serve as indicators of the underlying biological states. Inferring these latent, unobserved states is a goal of modern biology and neuroscience. However, in many experimental setups, we can at best obtain discrete snapshots of the system at different times and for different individuals. This challenge is particularly relevant in the study of Alzheimer's Disease (AD) progression, where we observe the aggregation of pathology in brain donors, but the underlying disease state is unknown. This paper proposes a biophysically motivated Bayesian framework (B-BIND: Biophysical Bayesian Inference for Neurodegenerative Dynamics), where the disease state is modeled and continuously inferred from observed quantifications of multiple AD pathological proteins. Inspired by biophysical models, we describe pathological burden as an exponential process. The progression of AD is modeled by assigning a latent score, termed pseudotime, to each pathological state, creating a pseudotemporal order of donors based on their pathological burden. We study the theoretical properties of the model using linearization to reveal convergence and identifiability properties. We provide Markov chain Monte Carlo estimation algorithms, illustrating the effectiveness of our approach with multiple simulation studies across various data conditions. Applying this methodology to data from the Seattle Alzheimer's Disease Brain Cell Atlas, we infer the pseudotime ordering of donors. Finally, we analyze the information within each pathological feature to refine the model, focusing on the most informative pathologies. This framework lays the groundwork for continuous pseudotime modeling in the analysis of neurodegenerative diseases.

Oncostatin M-driven macrophage-fibroblast circuits as a drug target in autoimmune arthritis

BACKGROUND

Recent single-cell RNA sequencing (scRNA-seq) analysis revealed the functional heterogeneity and pathogenic cell subsets in immune cells, synovial fibroblasts and bone cells in rheumatoid arthritis (RA). JAK inhibitors which ameliorate joint inflammation and bone destruction in RA, suppress the activation of various types of cells in vitro. However, the key cellular and molecular mechanisms underlying the potent clinical effects of JAK inhibitors on RA remain to be determined. Our aim is to identify a therapeutic target for JAK inhibitors in vivo.

METHODS

We performed scRNA-seq analysis of the synovium of collagen-induced arthritis (CIA) mice treated with or without a JAK inhibitor, followed by a computational analysis to identify the drug target cells and signaling pathways. We utilized integrated human RA scRNA-seq datasets and genetically modified mice administered with the JAK inhibitor for the confirmation of our findings.

RESULTS

scRNA-seq analysis revealed that oncostatin M (OSM) driven macrophage-fibroblast interaction is highly activated under arthritic conditions. OSM derived from macrophages, acts on OSM receptor (OSMR)-expressing synovial fibroblasts, activating both inflammatory and tissue-destructive subsets. Inflammatory synovial fibroblasts stimulate macrophages, mainly through IL-6, to exacerbate inflammation. Tissue-destructive synovial fibroblasts promote osteoclast differentiation by producing RANKL to accelerate bone destruction. scRNA-seq analysis also revealed that OSM-signaling in synovial fibroblasts is the main signaling pathway targeted by JAK inhibitors in vivo. Mice specifically lacking OSMR in synovial fibroblasts (Osmr∆Fibro) displayed ameliorated inflammation and joint destruction in arthritis. The JAK inhibitor was effective on the arthritis of the control mice while it had no effect on the arthritis of Osmr∆Fibro mice.

CONCLUSIONS

OSM functions as one of the key cytokines mediating pathogenic macrophage-fibroblast interaction. OSM-signaling in synovial fibroblasts is one of the main signaling pathways targeted by JAK inhibitors in vivo. The critical role of fibroblast-OSM signaling in autoimmune arthritis was shown by a combination of mice specifically deficient for OSMR in synovial fibroblasts and administration of the JAK inhibitor. Thus, the OSM-driven synovial macrophage-fibroblast circuit is proven to be a key driver of autoimmune arthritis, serving as a crucial drug target in vivo.

Metabolic reprogramming and heterogeneity during the decidualization process of endometrial stromal cells

The human endometrial decidualization is a transformative event in the pregnant uterus that involves the differentiation of stromal cells into decidual cells. While crucial to the establishment of a successful pregnancy, the metabolic characteristics of decidual cells in vivo remain largely unexplored. Here, we integrated the single-cell RNA sequencing (scRNA-seq) datasets on the endometrium of the menstrual cycle and the maternal-fetal interface in the first trimester to comprehensively decrypt the metabolic characteristics of stromal fibroblast cells. Our results revealed that the differentiation of stromal cells into decidual cells is accompanied by increased amino acid and sphingolipid metabolism. Furthermore, metabolic heterogeneity exists in decidual cells with differentiation maturity disparities. Decidual cells with high metabolism exhibit higher cellular activity and show a strong propensity for signaling. In addition, significant metabolic reprogramming in amino acids and lipids also occurs during the transition from non-pregnancy to pregnancy in the uteri of pigs, cattle, and mice. Our analysis provides comprehensive insights into the dynamic landscape of stromal fibroblast cell metabolism, contributing to our understanding of the metabolism at the molecular dynamics underlying the decidualization process in the human endometrium.

Mucosal Single-Cell Profiling of Crohn's-Like Disease of the Pouch Reveals Unique Pathogenesis and Therapeutic Targets

BACKGROUND & AIMS

The pathophysiology of Crohn's-like disease of the pouch (CDP) in patients with a history of ulcerative colitis (UC) is unknown. We examined mucosal cells from patients with and without CDP using single-cell analyses.

METHODS

Endoscopic samples were collected from pouch body and prepouch ileum (pouch/ileum) of 50 patients with an ileal pouch-anal anastomosis. Single-cell RNA sequencing was performed on pouch/ileal tissues of patients with normal pouch/ileum and CDP. Mass cytometry was performed on mucosal immune cells from patients with UC with normal pouch/ileum, CDP, pouchitis, and those with familial adenomatous polyposis after pouch formation. Findings were independently validated using immunohistochemistry.

RESULTS

The cell populations/states in the pouch body differed from those in the prepouch ileum, likely secondary to increased microbial burden. Compared with the familial adenomatous polyposis pouch, the UC pouch was enriched in colitogenic immune cells even without inflammation. CDP was characterized by increases in T helper 17 cells, inflammatory fibroblasts, inflammatory monocytes, TREM1+ monocytes, clonal expansion of effector T cells, and overexpression of T helper 17 cells-inducing cytokine genes such as IL23, IL1B, and IL6 by mononuclear phagocytes. Ligand-receptor analysis further revealed a stromal-mononuclear phagocytes-lymphocyte circuit in CDP. Integrated analysis showed that up-regulated immune mediators in CDP were similar to those in CD and pouchitis, but not UC. Additionally, CDP pouch/ileum exhibited heightened endoplasmic reticulum stress across all major cell compartments.

CONCLUSIONS

CDP likely represents a distinct entity of inflammatory bowel disease with heightened endoplasmic reticulum stress in both immune and nonimmune cells, which may become a novel diagnostic biomarker and therapeutic target for CDP.

Systematic analysis on the horse-shoe-like effect in PCA plots of scRNA-seq data

Motivation

In single-cell studies, principal component analysis (PCA) is widely used to reduce the dimensionality of dataset and visualize in 2D or 3D PC plots. Scientists often focus on different clusters within PC plot, overlooking the specific phenomenon, such as horse-shoe-like effect, that may reveal hidden knowledge about underlying biological dataset. This phenomenon remains largely unexplored in single-cell studies.

Results

In this study, we investigated into the horse-shoe-like effect in PC plots using simulated and real scRNA-seq datasets. We systematically explain horse-shoe-like phenomenon from various inter-related perspectives. Initially, we establish an intuitive understanding with the help of simulated datasets. Then, we generalized the acquired knowledge on real biological scRNA-seq data. Experimental results provide logical explanations and understanding for the appearance of horse-shoe-like effect in PC plots. Furthermore, we identify a potential problem with a well-known theory of 'distance saturation property' attributed to induce horse-shoe phenomenon. Finally, we analyse a mathematical model for horse-shoe effect that suggests trigonometric solutions to estimated eigenvectors. We observe significant resemblance after comparing the results of mathematical model with simulated and real scRNA-seq datasets.

Availability and implementation

The code for reproducing the results of this study is available at: https://github.com/najeebullahshah/PCA-Horse-Shoe.

Single-cell RNA sequencing reveals that MYBL2 in malignant epithelial cells is involved in the development and progression of ovarian cancer

Background

Ovarian carcinoma (OC) is a prevalent gynecological malignancy associated with high recurrence rates and mortality, often diagnosed at advanced stages. Despite advances in immunotherapy, immune exhaustion remains a significant challenge in achieving optimal tumor control. However, the exploration of intratumoral heterogeneity of malignant epithelial cells and the ovarian cancer tumor microenvironment is still limited, hindering our comprehensive understanding of the disease.

Materials and methods

Utilizing single-cell RNA sequencing (scRNA-seq), we comprehensively investigated the cellular composition across six ovarian cancer patients with omental metastasis. Our focus centered on analysis of the malignant epithelial cells. Employing CytoTRACE and slingshot pseudotime analyses, we identified critical subpopulations and explored associated transcription factors (TFs) influencing ovarian cancer progression. Furthermore, by integrating clinical factors from a large cohort of bulk RNA sequencing data, we have established a novel prognostic model to investigate the impact of the tumor immune microenvironment on ovarian cancer patients. Furthermore, we have investigated the condition of immunological exhaustion.

Results

Our study identified a distinct and highly proliferative subgroup of malignant epithelial cells, known as C2 TOP2A+ TCs. This subgroup primarily consisted of patients who hadn't received neoadjuvant chemotherapy. Ovarian cancer patients with elevated TOP2A expression exhibited heightened sensitivity to neoadjuvant chemotherapy (NACT). Moreover, the transcription factor MYBL2 in this subgroup played a critical role in ovarian cancer development. Additionally, we developed an independent prognostic indicator, the TOP2A TCs Risk Score (TTRS), which revealed a correlation between the High TTRS Group and unfavorable outcomes. Furthermore, immune infiltration and drug sensitivity analyses demonstrated increased responsiveness to Paclitaxel, Cisplatin, and Gemcitabine in the Low TTRS Group.

Conclusion

This research deepens our understanding of malignant epithelial cells in ovarian cancer and enhances our knowledge of the ovarian cancer immune microenvironment and immune exhaustion. We have revealed the heightened susceptibility of the C2 TOP2A+ TCs subgroup to neoadjuvant chemotherapy and emphasized the role of MYBL2 within the C2 subgroup in promoting the occurrence and progression of ovarian cancer. These insights provide valuable guidance for the management of ovarian cancer treatment.

Mapping the molecular landscape of Lotus japonicus nodule organogenesis through spatiotemporal transcriptomics

Legumes acquire nitrogen-fixing ability by forming root nodules. Transferring this capability to more crops could reduce our reliance on nitrogen fertilizers, thereby decreasing environmental pollution and agricultural production costs. Nodule organogenesis is complex, and a comprehensive transcriptomic atlas is crucial for understanding the underlying molecular events. Here, we utilized spatial transcriptomics to investigate the development of nodules in the model legume, Lotus japonicus. Our investigation has identified the developmental trajectories of two critical regions within the nodule: the infection zone and peripheral tissues. We reveal the underlying biological processes and provide gene sets to achieve symbiosis and material exchange, two essential aspects of nodulation. Among the candidate regulatory genes, we illustrate that LjNLP3, a transcription factor belonging to the NIN-LIKE PROTEIN family, orchestrates the transition of nodules from the differentiation to maturation. In summary, our research advances our understanding of nodule organogenesis and provides valuable data for developing symbiotic nitrogen-fixing crops.

Endothelial LATS2 is a suppressor of bone marrow fibrosis

Myelofibrosis and osteosclerosis are fibrotic diseases disrupting bone marrow function that occur in various leukemias but also in response to non-malignant alterations in hematopoietic cells. Here we show that endothelial cell-specific inactivation of the Lats2 gene, encoding Hippo kinase large tumor suppressor kinase 2, or overexpression of the downstream effector YAP1 induce myofibroblast formation and lead to extensive fibrosis and osteosclerosis, which impair bone marrow function and cause extramedullary hematopoiesis in the spleen. Mechanistically, loss of LATS2 induces endothelial-to-mesenchymal transition, resulting in increased expression of extracellular matrix and secreted signaling molecules. Changes in endothelial cells involve increased expression of serum response factor target genes, and, strikingly, major aspects of the LATS2 mutant phenotype are rescued by inactivation of the Srf gene. These findings identify the endothelium as a driver of bone marrow fibrosis, which improves understanding of myelofibrotic and osteosclerotic diseases, for which drug therapies are currently lacking.

Timing and location dictate monocyte fate and their transition to tumor-associated macrophages

Tumor-associated macrophages (TAMs) are a heterogeneous population of cells whose phenotypes and functions are shaped by factors that are incompletely understood. Herein, we asked when and where TAMs arise from blood monocytes and how they evolve during tumor development. We initiated pancreatic ductal adenocarcinoma (PDAC) in inducible monocyte fate-mapping mice and combined single-cell transcriptomics and high-dimensional flow cytometry to profile the monocyte-to-TAM transition. We revealed that monocytes differentiate first into a transient intermediate population of TAMs that generates two longer-lived lineages of terminally differentiated TAMs with distinct gene expression profiles, phenotypes, and intratumoral localization. Transcriptome datasets and tumor samples from patients with PDAC evidenced parallel TAM populations in humans and their prognostic associations. These insights will support the design of new therapeutic strategies targeting TAMs in PDAC.

Cell cycle expression heterogeneity predicts degree of differentiation

Methods that predict fate potential or degree of differentiation from transcriptomic data have identified rare progenitor populations and uncovered developmental regulatory mechanisms. However, some state-of-the-art methods are too computationally burdensome for emerging large-scale data and all methods make inaccurate predictions in certain biological systems. We developed a method in R (stemFinder) that predicts single cell differentiation time based on heterogeneity in cell cycle gene expression. Our method is computationally tractable and is as good as or superior to competitors. As part of our benchmarking, we implemented four different performance metrics to assist potential users in selecting the tool that is most apt for their application. Finally, we explore the relationship between differentiation time and cell fate potential by analyzing a lineage tracing dataset with clonally labelled hematopoietic cells, revealing that metrics of differentiation time are correlated with the number of downstream lineages.

A single-cell pan-cancer analysis to show the variability of tumor-infiltrating myeloid cells in immune checkpoint blockade

Myeloid cells are vital components of the immune system and have pivotal functions in orchestrating immune responses. Understanding their functions within the tumor microenvironment and their interactions with tumor-infiltrating lymphocytes presents formidable challenges across diverse cancer types, particularly with regards to cancer immunotherapies. Here, we explore tumor-infiltrating myeloid cells (TIMs) by conducting a pan-cancer analysis using single-cell transcriptomics across eight distinct cancer types, encompassing a total of 192 tumor samples from 129 patients. By examining gene expression patterns and transcriptional activities of TIMs in different cancer types, we discern notable alterations in abundance of TIMs and kinetic behaviors prior to and following immunotherapy. We also identify specific cell-cell interaction targets in immunotherapy; unique and shared regulatory profiles critical for treatment response; and TIMs associated with survival outcomes. Overall, our study illuminates the heterogeneity of TIMs and improves our understanding of tissue-specific and cancer-specific myeloid subsets within the context of tumor immunotherapies.

Human lymph node fibroblastic reticular cells maintain heterogeneous characteristics in culture

Fibroblastic reticular cells (FRCs) are mesenchymal stromal cells in human lymph nodes (LNs) playing a pivotal role in adaptive immunity. Several FRC subsets have been identified, yet it remains to be elucidated if their heterogeneity is maintained upon culture. Here, we established a protocol to preserve and culture FRCs from human LNs and characterized their phenotypic profile in fresh LN suspensions and upon culture using multispectral flow cytometry. We found nine FRC subsets in fresh human LNs, independent of donor, of which four persisted in culture throughout several passages. Interestingly, the historically FRC-defining marker podoplanin (PDPN) was not present on all FRC subsets. Therefore, we propose that CD45negCD31neg human FRCs are not restricted by PDPN expression, as we found CD90, BST1, and CD146/MCAM to be more widely expressed. Together, our data provide insight into FRC heterogeneity in human LNs, enabling further investigation into the function of individual FRC subsets.

Transcriptional dynamics orchestrating the development and integration of neurons born in the adult hippocampus

The adult hippocampus generates new granule cells (aGCs) with functional capabilities that convey unique forms of plasticity to the preexisting circuits. While early differentiation of adult radial glia-like cells (RGLs) has been studied extensively, the molecular mechanisms guiding the maturation of postmitotic neurons remain unknown. Here, we used a precise birthdating strategy to study aGC differentiation using single-nuclei RNA sequencing. Transcriptional profiling revealed a continuous trajectory from RGLs to mature aGCs, with multiple immature stages bearing increasing levels of effector genes supporting growth, excitability, and synaptogenesis. Analysis of differential gene expression, pseudo-time trajectory, and transcription factors (TFs) revealed critical transitions defining four cellular states: quiescent RGLs, proliferative progenitors, immature aGCs, and mature aGCs. Becoming mature aGCs involved a transcriptional switch that shuts down pathways promoting cell growth, such SoxC TFs, to activate programs that likely control neuronal homeostasis. aGCs overexpressing Sox4 or Sox11 remained immature. Our results unveil precise molecular mechanisms driving adult RGLs through the pathway of neuronal differentiation.

Single-cell decoding of drug induced transcriptomic reprogramming in triple negative breast cancers

BACKGROUND

The encoding of cell intrinsic drug resistance states in breast cancer reflects the contributions of genomic and non-genomic variations and requires accurate estimation of clonal fitness from co-measurement of transcriptomic and genomic data. Somatic copy number (CN) variation is the dominant mutational mechanism leading to transcriptional variation and notably contributes to platinum chemotherapy resistance cell states. Here, we deploy time series measurements of triple negative breast cancer (TNBC) single-cell transcriptomes, along with co-measured single-cell CN fitness, identifying genomic and transcriptomic mechanisms in drug-associated transcriptional cell states.

RESULTS

We present scRNA-seq data (53,641 filtered cells) from serial passaging TNBC patient-derived xenograft (PDX) experiments spanning 2.5 years, matched with genomic single-cell CN data from the same samples. Our findings reveal distinct clonal responses within TNBC tumors exposed to platinum. Clones with high drug fitness undergo clonal sweeps and show subtle transcriptional reversion, while those with weak fitness exhibit dynamic transcription upon drug withdrawal. Pathway analysis highlights convergence on epithelial-mesenchymal transition and cytokine signaling, associated with resistance. Furthermore, pseudotime analysis demonstrates hysteresis in transcriptional reversion, indicating generation of new intermediate transcriptional states upon platinum exposure.

CONCLUSIONS

Within a polyclonal tumor, clones with strong genotype-associated fitness under platinum remained fixed, minimizing transcriptional reversion upon drug withdrawal. Conversely, clones with weaker fitness display non-genomic transcriptional plasticity. This suggests CN-associated and CN-independent transcriptional states could both contribute to platinum resistance. The dominance of genomic or non-genomic mechanisms within polyclonal tumors has implications for drug sensitivity, restoration, and re-treatment strategies.

Intergenic risk variant rs56258221 skews the fate of naive CD4+ T cells via miR4464-BACH2 interplay in primary sclerosing cholangitis

Primary sclerosing cholangitis (PSC) is an immune-mediated liver disease of unknown pathogenesis, with a high risk to develop cirrhosis and malignancies. Functional dysregulation of T cells and association with genetic polymorphisms in T cell-related genes were previously reported for PSC. Here, we genotyped a representative PSC cohort for several disease-associated risk loci and identified rs56258221 (BACH2/MIR4464) to correlate with not only the peripheral blood T cell immunophenotype but also the functional capacities of naive CD4+ T (CD4+ TN) cells in people with PSC. Mechanistically, rs56258221 leads to an increased expression of miR4464, in turn causing attenuated translation of BACH2, a major gatekeeper of T cell quiescence. Thereby, the fate of CD4+ TN is skewed toward polarization into pro-inflammatory subsets. Clinically, people with PSC carrying rs56258221 show signs of accelerated disease progression. The data presented here highlight the importance of assigning functional outcomes to disease-associated genetic polymorphisms as potential drivers of diseases.

Directly selecting cell-type marker genes for single-cell clustering analyses

In single-cell RNA sequencing (scRNA-seq) studies, cell types and their marker genes are often identified by clustering and differentially expressed gene (DEG) analysis. A common practice is to select genes using surrogate criteria such as variance and deviance, then cluster them using selected genes and detect markers by DEG analysis assuming known cell types. The surrogate criteria can miss important genes or select unimportant genes, while DEG analysis has the selection-bias problem. We present Festem, a statistical method for the direct selection of cell-type markers for downstream clustering. Festem distinguishes marker genes with heterogeneous distribution across cells that are cluster informative. Simulation and scRNA-seq applications demonstrate that Festem can sensitively select markers with high precision and enables the identification of cell types often missed by other methods. In a large intrahepatic cholangiocarcinoma dataset, we identify diverse CD8+ T cell types and potential prognostic marker genes.

Fatecode enables cell fate regulator prediction using classification-supervised autoencoder perturbation

Cell reprogramming, which guides the conversion between cell states, is a promising technology for tissue repair and regeneration, with the ultimate goal of accelerating recovery from diseases or injuries. To accomplish this, regulators must be identified and manipulated to control cell fate. We propose Fatecode, a computational method that predicts cell fate regulators based only on single-cell RNA sequencing (scRNA-seq) data. Fatecode learns a latent representation of the scRNA-seq data using a deep learning-based classification-supervised autoencoder and then performs in silico perturbation experiments on the latent representation to predict genes that, when perturbed, would alter the original cell type distribution to increase or decrease the population size of a cell type of interest. We assessed Fatecode's performance using simulations from a mechanistic gene-regulatory network model and scRNA-seq data mapping blood and brain development of different organisms. Our results suggest that Fatecode can detect known cell fate regulators from single-cell transcriptomics datasets.

Early human fetal lung atlas reveals the temporal dynamics of epithelial cell plasticity

Studying human fetal lungs can inform how developmental defects and disease states alter the function of the lungs. Here, we sequenced >150,000 single cells from 19 healthy human pseudoglandular fetal lung tissues ranging between gestational weeks 10-19. We capture dynamic developmental trajectories from progenitor cells that express abundant levels of the cystic fibrosis conductance transmembrane regulator (CFTR). These cells give rise to multiple specialized epithelial cell types. Combined with spatial transcriptomics, we show temporal regulation of key signalling pathways that may drive the temporal and spatial emergence of specialized epithelial cells including ciliated and pulmonary neuroendocrine cells. Finally, we show that human pluripotent stem cell-derived fetal lung models contain CFTR-expressing progenitor cells that capture similar lineage developmental trajectories as identified in the native tissue. Overall, this study provides a comprehensive single-cell atlas of the developing human lung, outlining the temporal and spatial complexities of cell lineage development and benchmarks fetal lung cultures from human pluripotent stem cell differentiations to similar developmental window.

Citalopram exposure of hESCs during neuronal differentiation identifies dysregulated genes involved in neurodevelopment and depression

Selective serotonin reuptake inhibitors (SSRIs), including citalopram, are widely used antidepressants during pregnancy. However, the effects of prenatal exposure to citalopram on neurodevelopment remain poorly understood. We aimed to investigate the impact of citalopram exposure on early neuronal differentiation of human embryonic stem cells using a multi-omics approach. Citalopram induced time- and dose-dependent effects on gene expression and DNA methylation of genes involved in neurodevelopmental processes or linked to depression, such as BDNF, GDF11, CCL2, STC1, DDIT4 and GAD2. Single-cell RNA-sequencing analysis revealed distinct clusters of stem cells, neuronal progenitors and neuroblasts, where exposure to citalopram subtly influenced progenitor subtypes. Pseudotemporal analysis showed enhanced neuronal differentiation. Our findings suggest that citalopram exposure during early neuronal differentiation influences gene expression patterns associated with neurodevelopment and depression, providing insights into its potential neurodevelopmental impact and highlighting the importance of further research to understand the long-term consequences of prenatal SSRI exposure.

Endogenous FGFs drive ERK-dependent cell fate patterning in 2D human gastruloids

The role of FGF is the least understood of the morphogens driving mammalian gastrulation. Here we investigated the function of FGF in a stem cell model for human gastrulation known as a 2D gastruloid. We found a ring of FGF-dependent ERK activity that closely follows the emergence of primitive streak (PS)-like cells but expands further inward. We showed that this ERK activity pattern is required for PS-like differentiation and that loss of PS-like cells upon FGF receptor inhibition can be rescued by directly activating ERK. We further demonstrated that the ERK-ring depends on localized activation of basally localized FGF receptors (FGFR) by endogenous FGF gradients. We confirm and extend previous studies in analyzing expression of FGF pathway components, showing the main receptor to be FGFR1 and the key ligands FGF2/4/17, similar to the human and monkey embryo but different from the mouse. In situ hybridization and scRNA-seq revealed that FGF4 and FGF17 expression colocalize with the PS marker TBXT but only FGF17 is maintained in nascent mesoderm and endoderm. FGF4 and FGF17 reduction both reduced ERK activity and differentiation to PS-like cells and their derivatives, indicating overlapping function. Thus, we have identified a previously unknown role for FGF-dependent ERK signaling in 2D gastruloids and possibly the human embryo, driven by a mechanism where FGF4 and FGF17 signal through basally localized FGFR1 to induce PS-like cells.