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

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.