psupertime: supervised pseudotime analysis for time-series single-cell RNA-seq data

TSPAN4+ fibroblasts coordinate metastatic niche assembly through migrasome-driven metabolic reprogramming and stromal-immune crosstalk in pancreatic adenocarcinoma

Background

Pancreatic cancer (PC) is a highly aggressive pancreatic malignant tumor with poor prognosis due to its complex tumor microenvironment (TME) and metastatic potential. Fibroblasts play an important role in tumor progression and metastasis by remodeling the extracellular matrix (ECM) and influencing the immune response. This study explored migrasome-associated fibroblasts, especially TSPAN4 and ITGA5, as key regulators of pancreatic cancer metastasis, aiming to provide new ideas for therapeutic strategies targeting TME.

Methods

We employed single-cell RNA sequencing (scRNA-seq) and spatial transcriptomics to analyze pancreatic cancer tissues. Data from the GEO and TCGA databases were integrated and processed using batch correction techniques. Single-cell data were analyzed with Seurat and Monocle for dimensionality reduction and pseudotime trajectory analysis. Cell communication was assessed using CellCall and CellChat. Spatial transcriptomic analysis was conducted with RCTD and MISTy tools to investigate cell interactions within the TME. Additionally, gene enrichment, deconvolution, and prognostic analyses were performed, alongside experimental validation through siRNA transfection, qRT-PCR, and various functional assays to investigate the role of TSPAN4 in metastasis.

Results

Our results underscore the critical role of TSPAN4+ fibroblasts in pancreatic cancer. These fibroblasts were predominantly located at the tumor periphery and exhibited elevated migrasome gene expression, which was associated with enhanced ECM remodeling and immune suppression. Intercellular communication analysis revealed that TSPAN4+ fibroblasts engaged in extensive interactions with immune cells, such as macrophages and endothelial cells, facilitating metastasis and immune evasion. Moreover, the high expression of immune checkpoint markers indicated their involvement in modulating the immune response.

Conclusion

TSPAN4+ fibroblasts are key regulators of pancreatic cancer progression, contributing to metastasis, immune suppression, and ECM remodeling. Targeting these fibroblasts represents a promising therapeutic strategy to improve clinical outcomes and enhance the effectiveness of immunotherapies in pancreatic cancer.

Timestamp calibration for time-series single cell RNA-seq expression data

Timestamp automatic annotation (TAA) is a crucial procedure for analyzing time-series scRNA-seq data, as they unveil dynamic biological developments and cell regeneration processes. However, current TAA methods heavily rely on manual timestamps, often overlooking their reliability. This oversight can significantly degrade the performance of timestamp automatic annotation due to noisy timestamps. Nevertheless, the current approach for addressing this issue tends to select less critical cleaned samples for timestamp calibration. To tackle this challenge, we have developed a novel timestamp calibration model called ScPace for handling noisy labeled time-series scRNA-seq data. This approach incorporates a latent variable indicator within a base classifier instead of probability sampling to detect noisy samples effectively. To validate our proposed method, we conducted experiments on both simulated and real time-series scRNA-seq datasets. Cross validation experiments with different artificial mislabeling rates demonstrate that ScPace outperforms previous approaches. Furthermore, after calibrating the timestamps of the original time-series scRNA-seq data using our method, we performed supervised pseudotime analysis, revealing that ScPace enhances its performance significantly. These findings suggest that ScPace is an effective tool for timestamp calibration by enabling reclassification and deletion of detected noisy labeled samples while maintaining robustness across diverse ranges of time-series scRNA-seq datasets. The source code is available at https://github.com/OPUS-Lightphenexx/ScPace.

Advanced Carriers for Precise Delivery and Therapeutic Mechanisms of Traditional Chinese Medicines: Integrating Spatial Multi-Omics and Delivery Visualization

The complex composition of traditional Chinese medicines (TCMs) has posed challenges for in-depth study and global application, despite their abundance of bioactive compounds that make them valuable resources for disease treatment. To overcome these obstacles, it is essential to modernize TCMs by focusing on precise disease treatment. This involves elucidating the structure-activity relationships within their complex compositions, ensuring accurate in vivo delivery, and monitoring the delivery process. This review discusses the research progress of TCMs in precision disease treatment from three perspectives: spatial multi-omics technology for precision therapeutic activity, carrier systems for precise in vivo delivery, and medical imaging technology for visualizing the delivery process. The aim is to establish a novel research paradigm that advances the precision therapy of TCMs.

Time-series single-cell transcriptomic profiling of luteal-phase endometrium uncovers dynamic characteristics and its dysregulation in recurrent implantation failures

Understanding human endometrial dynamics in the establishment of endometrial receptivity remains a challenge, which limits early diagnosis and treatment of endometrial-factor infertility. Here, we decode the endometrial dynamics of fertile women across the window of implantation and characterize the endometrial deficiency in women with recurrent implantation failure. A computational model capable of both temporal prediction and pattern discovery is used to analyze single-cell transcriptomic data from over 220,000 endometrial cells. The time-series atlas highlights a two-stage stromal decidualization process and a gradual transitional process of the luminal epithelial cells across the window of implantation. In addition, a time-varying gene set regulating epithelium receptivity is identified, based on which the recurrent implantation failure endometria are stratified into two classes of deficiencies. Further investigation uncovers a hyper-inflammatory microenvironment for the dysfunctional endometrial epithelial cells of recurrent implantation failure. The holistic characterization of the physiological and pathophysiological window of implantation and a computational tool trained on this temporal atlas provide a platform for future therapeutic developments.

A temporal cortex cell atlas highlights gene expression dynamics during human brain maturation

The human brain undergoes protracted postnatal maturation, guided by dynamic changes in gene expression. Most studies exploring these processes have used bulk tissue analyses, which mask cell-type-specific gene expression dynamics. Here, using single-nucleus RNA sequencing on temporal lobe tissue, including samples of African ancestry, we build a joint pediatric and adult atlas of 75 cell subtypes, which we verify with spatial transcriptomics. We explore the differences between pediatric and adult cell subtypes, revealing the genes and pathways that change during brain maturation. Our results highlight excitatory neuron subtypes, including the LTK and FREM subtypes, that show elevated expression of genes associated with cognition and synaptic plasticity in pediatric tissue. The resources we present here improve our understanding of the brain during its development and contribute to global efforts to build an inclusive brain cell map.

Constrained Pseudo-Time Ordering for Clinical Transcriptomics Data

Time series RNASeq studies can enable understanding of the dynamics of disease progression and treatment response in patients. They also provide information on biomarkers, activated and repressed pathways, and more. While useful, data from multiple patients is challenging to integrate due to the heterogeneity in treatment response among patients, and the small number of timepoints that are usually profiled. Due to the heterogeneity among patients, relying on the sampled time points to integrate data across individuals is challenging and does not lead to correct reconstruction of the response patterns. To address these challenges, we developed a new constrained based pseudo-time ordering method for analyzing transcriptomics data in clinical and response studies. Our method allows the assignment of samples to their correct placement on the response curve while respecting the individual patient order. We use polynomials to represent gene expression over the duration of the study and an EM algorithm to determine parameters and locations. Application to four treatment response datasets shows that our method improves on prior methods and leads to accurate orderings that provide new biological insight on the disease and response.

Quantifying cell-state densities in single-cell phenotypic landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive diverse biological processes. Here, we present Mellon, an algorithm for estimation of cell-state densities from high-dimensional representations of single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. We present evidence implicating enhancer priming and the activation of master regulators in emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during developmental processes. Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide insights into mechanisms guiding biological trajectories.

Integrating patients in time series clinical transcriptomics data

MOTIVATION

Analysis of time series transcriptomics data from clinical trials is challenging. Such studies usually profile very few time points from several individuals with varying response patterns and dynamics. Current methods for these datasets are mainly based on linear, global orderings using visit times which do not account for the varying response rates and subgroups within a patient cohort.

RESULTS

We developed a new method that utilizes multi-commodity flow algorithms for trajectory inference in large scale clinical studies. Recovered trajectories satisfy individual-based timing restrictions while integrating data from multiple patients. Testing the method on multiple drug datasets demonstrated an improved performance compared to prior approaches suggested for this task, while identifying novel disease subtypes that correspond to heterogeneous patient response patterns.

AVAILABILITY AND IMPLEMENTATION

The source code and instructions to download the data have been deposited on GitHub at https://github.com/euxhenh/Truffle.

Single cell lineage tracing reveals clonal dynamics of anti-EGFR therapy resistance in triple negative breast cancer

BACKGROUND

Most primary Triple Negative Breast Cancers (TNBCs) show amplification of the Epidermal Growth Factor Receptor (EGFR) gene, leading to increased protein expression. However, unlike other EGFR-driven cancers, targeting this receptor in TNBC yields inconsistent therapeutic responses.

METHODS

To elucidate the underlying mechanisms of this variability, we employ cellular barcoding and single-cell transcriptomics to reconstruct the subclonal dynamics of EGFR-amplified TNBC cells in response to afatinib, a tyrosine kinase inhibitor (TKI) that irreversibly inhibits EGFR.

RESULTS

Integrated lineage tracing analysis revealed a rare pre-existing subpopulation of cells with distinct biological signature, including elevated expression levels of Insulin-Like Growth Factor Binding Protein 2 (IGFBP2). We show that IGFBP2 overexpression is sufficient to render TNBC cells tolerant to afatinib treatment by activating the compensatory insulin-like growth factor I receptor (IGF1-R) signalling pathway. Finally, based on reconstructed mechanisms of resistance, we employ deep learning techniques to predict the afatinib sensitivity of TNBC cells.

CONCLUSIONS

Our strategy proved effective in reconstructing the complex signalling network driving EGFR-targeted therapy resistance, offering new insights for the development of individualized treatment strategies in TNBC.

Single-cell time series analysis reveals the dynamics of HSPC response to inflammation

Hematopoietic stem and progenitor cells (HSPCs) are known to respond to acute inflammation; however, little is understood about the dynamics and heterogeneity of these stress responses in HSPCs. Here, we performed single-cell sequencing during the sensing, response, and recovery phases of the inflammatory response of HSPCs to treatment (a total of 10,046 cells from four time points spanning the first 72 h of response) with the pro-inflammatory cytokine IFNα to investigate the HSPCs' dynamic changes during acute inflammation. We developed the essential novel computational approaches to process and analyze the resulting single-cell time series dataset. This includes an unbiased cell type annotation and abundance analysis post inflammation, tools for identification of global and cell type-specific responding genes, and a semi-supervised linear regression approach for response pseudotime reconstruction. We discovered a variety of different gene responses of the HSPCs to the treatment. Interestingly, we were able to associate a global reduced myeloid differentiation program and a locally enhanced pyroptosis activity with reduced myeloid progenitor and differentiated cells after IFNα treatment. Altogether, the single-cell time series analyses have allowed us to unbiasedly study the heterogeneous and dynamic impact of IFNα on the HSPCs.

Mesenchymal ovarian cancer cells promote CD8+ T cell exhaustion through the LGALS3-LAG3 axis

Cancer cells often metastasize by undergoing an epithelial-mesenchymal transition (EMT). Although abundance of CD8+ T-cells in the tumor microenvironment correlates with improved survival, mesenchymal cancer cells acquire greater resistance to antitumor immunity in some cancers. We hypothesized the EMT modulates the immune response to ovarian cancer. Here we show that cancer cells from infiltrated/inflamed tumors possess more mesenchymal cells, than excluded and desert tumors. We also noted high expression of LGALS3 is associated with EMT in vivo, a finding validated with in vitro EMT models. Dissecting the cellular communications among populations in the tumor revealed that mesenchymal cancer cells in infiltrated tumors communicate through LGALS3 to LAG3 receptor expressed by CD8+ T cells. We found CD8+ T cells express high levels of LAG3, a marker of T cell exhaustion. The results indicate that EMT in ovarian cancer cells promotes interactions between cancer cells and T cells through the LGALS3 - LAG3 axis, which could increase T cell exhaustion in infiltrated tumors, dampening antitumor immunity.

Quantifying Cell-State Densities in Single-Cell Phenotypic Landscapes using Mellon

Cell-state density characterizes the distribution of cells along phenotypic landscapes and is crucial for unraveling the mechanisms that drive cellular differentiation, regeneration, and disease. Here, we present Mellon, a novel computational algorithm for high-resolution estimation of cell-state densities from single-cell data. We demonstrate Mellon's efficacy by dissecting the density landscape of various differentiating systems, revealing a consistent pattern of high-density regions corresponding to major cell types intertwined with low-density, rare transitory states. Utilizing hematopoietic stem cell fate specification to B-cells as a case study, we present evidence implicating enhancer priming and the activation of master regulators in the emergence of these transitory states. Mellon offers the flexibility to perform temporal interpolation of time-series data, providing a detailed view of cell-state dynamics during the inherently continuous developmental processes. Scalable and adaptable, Mellon facilitates density estimation across various single-cell data modalities, scaling linearly with the number of cells. Our work underscores the importance of cell-state density in understanding the differentiation processes, and the potential of Mellon to provide new insights into the regulatory mechanisms guiding cellular fate decisions.

Single-cell RNA sequencing in skeletal muscle developmental biology

Skeletal muscle is the most extensive tissue in mammals, and they perform several functions; it is derived from paraxial mesodermal somites and undergoes hyperplasia and hypertrophy to form multinucleated, contractile, and functional muscle fibers. Skeletal muscle is a complex heterogeneous tissue composed of various cell types that establish communication strategies to exchange biological information; therefore, characterizing the cellular heterogeneity and transcriptional signatures of skeletal muscle is central to understanding its ontogeny's details. Studies of skeletal myogenesis have focused primarily on myogenic cells' proliferation, differentiation, migration, and fusion and ignored the intricate network of cells with specific biological functions. The rapid development of single-cell sequencing technology has recently enabled the exploration of skeletal muscle cell types and molecular events during development. This review summarizes the progress in single-cell RNA sequencing and its applications in skeletal myogenesis, which will provide insights into skeletal muscle pathophysiology.

Cell atlas of CCl 4-induced progressive liver fibrosis reveals stage-specific responses

Chronic liver injury leads to progressive liver fibrosis and ultimately cirrhosis, a major cause of morbidity and mortality worldwide. However, there are currently no effective anti-fibrotic therapies available, especially for late-stage patients, which is partly attributed to the major knowledge gap regarding liver cell heterogeneity and cell-specific responses in different fibrosis stages. To reveal the multicellular networks regulating mammalian liver fibrosis from mild to severe phenotypes, we generated a single-nucleus transcriptomic atlas encompassing 49 919 nuclei corresponding to all main liver cell types at different stages of murine carbon tetrachloride (CCl 4)-induced progressive liver fibrosis. Integrative analysis distinguished the sequential responses to injury of hepatocytes, hepatic stellate cells and endothelial cells. Moreover, we reconstructed cell-cell interactions and gene regulatory networks implicated in these processes. These integrative analyses uncovered previously overlooked aspects of hepatocyte proliferation exhaustion and disrupted pericentral metabolic functions, dysfunction for clearance by apoptosis of activated hepatic stellate cells, accumulation of pro-fibrotic signals, and the switch from an anti-angiogenic to a pro-angiogenic program during CCl 4-induced progressive liver fibrosis. Our dataset thus constitutes a useful resource for understanding the molecular basis of progressive liver fibrosis using a relevant animal model.

Single-cell gene set enrichment analysis and transfer learning for functional annotation of scRNA-seq data

Although an essential step, cell functional annotation often proves particularly challenging from single-cell transcriptional data. Several methods have been developed to accomplish this task. However, in most cases, these rely on techniques initially developed for bulk RNA sequencing or simply make use of marker genes identified from cell clustering followed by supervised annotation. To overcome these limitations and automatize the process, we have developed two novel methods, the single-cell gene set enrichment analysis (scGSEA) and the single-cell mapper (scMAP). scGSEA combines latent data representations and gene set enrichment scores to detect coordinated gene activity at single-cell resolution. scMAP uses transfer learning techniques to re-purpose and contextualize new cells into a reference cell atlas. Using both simulated and real datasets, we show that scGSEA effectively recapitulates recurrent patterns of pathways' activity shared by cells from different experimental conditions. At the same time, we show that scMAP can reliably map and contextualize new single-cell profiles on a breast cancer atlas we recently released. Both tools are provided in an effective and straightforward workflow providing a framework to determine cell function and significantly improve annotation and interpretation of scRNA-seq data.