Multicellular factor analysis of single-cell data for a tissue-centric understanding of disease

Multiple Sclerosis: Glial Cell Diversity in Time and Space

Multiple sclerosis (MS) is the most prevalent human inflammatory disease of the central nervous system with demyelination and glial scar formation as pathological hallmarks. Glial cells are key drivers of lesion progression in MS with roles in both tissue damage and repair depending on the surrounding microenvironment and the functional state of the individual glial subtype. In this review, we describe recent developments in the context of glial cell diversity in MS summarizing key findings with respect to pathological and maladaptive functions related to disease-associated glial subtypes. A particular focus is on the spatial and temporal dynamics of glial cells including subtypes of microglia, oligodendrocytes, and astrocytes. We contextualize recent high-dimensional findings suggesting that glial cells dynamically change with respect to epigenomic, transcriptomic, and metabolic features across the inflamed rim and during the progression of MS lesions. In summary, detailed knowledge of spatially restricted glial subtype functions is critical for a better understanding of MS pathology and its pathogenesis as well as the development of novel MS therapies targeting specific glial cell types.

Integrating Spatially-Resolved Transcriptomics Data Across Tissues and Individuals: Challenges and Opportunities

Advances in spatially-resolved transcriptomics (SRT) technologies have propelled the development of new computational analysis methods to unlock biological insights. The lowering cost of SRT data generation presents an unprecedented opportunity to create large-scale spatial atlases and enable population-level investigation, integrating SRT data across multiple tissues, individuals, species, or phenotypes. Here, unique challenges are described in the SRT data integration, where the analytic impact of varying spatial and biological resolutions is characterized and explored. A succinct review of spatially-aware integration methods and computational strategies is provided. Exciting opportunities to advance computational algorithms amenable to atlas-scale datasets along with standardized preprocessing methods, leading to improved sensitivity and reproducibility in the future are further highlighted.

snRNA-seq stratifies multiple sclerosis patients into distinct white matter glial responses

Poor understanding of the cellular and molecular basis of clinical and genetic heterogeneity in progressive multiple sclerosis (MS) has hindered the search for new effective therapies. To address this gap, we analyzed 632,000 single-nucleus RNA sequencing profiles from 156 brain tissue samples of MS and control donors to examine inter- and intra-donor heterogeneity. We found distinct cell type-specific gene expression changes between MS gray and white matter, highlighting clear pathology differences. MS lesion subtypes had different cellular compositions but surprisingly similar cell-type gene expression patterns both within and across patients, suggesting global changes. Most gene expression variability was instead explained by patient effects, allowing us to stratify patients and describe the different pathological processes occurring between patient subgroups. Future mapping of these brain molecular profiles with blood and/or CSF profiles from living MS patients will allow precision medicine approaches anchored in patient-specific pathological processes.

Unveiling adipose populations linked to metabolic health in obesity

Precision medicine is still not considered as a standard of care in obesity treatment, despite a large heterogeneity in the metabolic phenotype of individuals with obesity. One of the strongest factors influencing the variability in metabolic disease risk is adipose tissue (AT) dysfunction; however, there is little understanding of the link between distinct cell populations, cell-type-specific transcriptional programs, and disease severity. Here, we generated a comprehensive cellular map of subcutaneous and visceral AT of individuals with metabolically healthy and unhealthy obesity. By combining single-nucleus RNA-sequencing data with bulk transcriptomics and clinical parameters, we identified that mesothelial cells, adipocytes, and adipocyte-progenitor cells exhibit the strongest correlation with metabolic disease. Furthermore, we uncovered cell-specific transcriptional programs, such as the transitioning of mesothelial cells to a mesenchymal phenotype, that are involved in uncoupling obesity from metabolic disease. Together, these findings provide valuable insights by revealing biological drivers of clinical endpoints.

Cell type mapping reveals tissue niches and interactions in subcortical multiple sclerosis lesions

Multiple sclerosis (MS) is a chronic inflammatory disease of the central nervous system. Inflammation is gradually compartmentalized and restricted to specific tissue niches such as the lesion rim. However, the precise cell type composition of such niches, their interactions and changes between chronic active and inactive stages are incompletely understood. We used single-nucleus and spatial transcriptomics from subcortical MS and corresponding control tissues to map cell types and associated pathways to lesion and nonlesion areas. We identified niches such as perivascular spaces, the inflamed lesion rim or the lesion core that are associated with the glial scar and a cilia-forming astrocyte subtype. Focusing on the inflamed rim of chronic active lesions, we uncovered cell-cell communication events between myeloid, endothelial and glial cell types. Our results provide insight into the cellular composition, multicellular programs and intercellular communication in tissue niches along the conversion from a homeostatic to a dysfunctional state underlying lesion progression in MS.

QOT: Quantized Optimal Transport for sample-level distance matrix in single-cell omics

Single-cell technologies have enabled the high-dimensional characterization of cell populations at an unprecedented scale. The innate complexity and increasing volume of data pose significant computational and analytical challenges, especially in comparative studies delineating cellular architectures across various biological conditions (i.e. generation of sample-level distance matrices). Optimal Transport is a mathematical tool that captures the intrinsic structure of data geometrically and has been applied to many bioinformatics tasks. In this paper, we propose QOT (Quantized Optimal Transport), a new method enabling efficient computation of sample-level distance matrix from large-scale single-cell omics data through a quantization step. We apply our algorithm to real-world single-cell genomics and pathomics datasets, aiming to extrapolate cell-level insights to inform sample-level categorizations. Our empirical study shows that QOT outperforms existing two OT-based algorithms in accuracy and robustness when obtaining a distance matrix from high throughput single-cell measures at the sample level. Moreover, the sample level distance matrix could be used in the downstream analysis (i.e. uncover the trajectory of disease progression), highlighting its usage in biomedical informatics and data science.

Uncovering disease-related multicellular pathway modules on large-scale single-cell transcriptomes with scPAFA

Pathway analysis is a crucial analytical phase in disease research on single-cell RNA sequencing (scRNA-seq) data, offering biological interpretations based on prior knowledge. However, currently available tools for generating cell-level pathway activity scores (PAS) exhibit computational inefficacy in large-scale scRNA-seq datasets. Additionally, disease-related pathways are often identified through cross-condition comparisons within specific cell types, overlooking potential patterns that involve multiple cell types. Here, we present single-cell pathway activity factor analysis (scPAFA), a Python library designed for large-scale single-cell datasets allowing rapid PAS computation and uncovering biologically interpretable disease-related multicellular pathway modules, which are low-dimensional representations of disease-related PAS alterations in multiple cell types. Application on colorectal cancer (CRC) datasets and large-scale lupus atlas over 1.2 million cells demonstrated that scPAFA can achieve over 40-fold reductions in the runtime of PAS computation and further identified reliable and interpretable multicellular pathway modules that capture the heterogeneity of CRC and transcriptional abnormalities in lupus patients, respectively. Overall, scPAFA presents a valuable addition to existing research tools in disease research, with the potential to reveal complex disease mechanisms and support biomarker discovery at the pathway level.

Visualizing scRNA-Seq data at population scale with GloScope

Increasingly, scRNA-Seq studies explore cell populations across different samples and the effect of sample heterogeneity on organism's phenotype. However, relatively few bioinformatic methods have been developed which adequately address the variation between samples for such population-level analyses. We propose a framework for representing the entire single-cell profile of a sample, which we call a GloScope representation. We implement GloScope on scRNA-Seq datasets from study designs ranging from 12 to over 300 samples and demonstrate how GloScope allows researchers to perform essential bioinformatic tasks at the sample-level, in particular visualization and quality control assessment.

Inferring pattern-driving intercellular flows from single-cell and spatial transcriptomics

From single-cell RNA-sequencing (scRNA-seq) and spatial transcriptomics (ST), one can extract high-dimensional gene expression patterns that can be described by intercellular communication networks or decoupled gene modules. These two descriptions of information flow are often assumed to occur independently. However, intercellular communication drives directed flows of information that are mediated by intracellular gene modules, in turn triggering outflows of other signals. Methodologies to describe such intercellular flows are lacking. We present FlowSig, a method that infers communication-driven intercellular flows from scRNA-seq or ST data using graphical causal modeling and conditional independence. We benchmark FlowSig using newly generated experimental cortical organoid data and synthetic data generated from mathematical modeling. We demonstrate FlowSig's utility by applying it to various studies, showing that FlowSig can capture stimulation-induced changes to paracrine signaling in pancreatic islets, demonstrate shifts in intercellular flows due to increasing COVID-19 severity and reconstruct morphogen-driven activator-inhibitor patterns in mouse embryogenesis.

Coordinated, multicellular patterns of transcriptional variation that stratify patient cohorts are revealed by tensor decomposition

Tissue-level and organism-level biological processes often involve the coordinated action of multiple distinct cell types. The recent application of single-cell assays to many individuals should enable the study of how donor-level variation in one cell type is linked to that in other cell types. Here we introduce a computational approach called single-cell interpretable tensor decomposition (scITD) to identify common axes of interindividual variation by considering joint expression variation across multiple cell types. scITD combines expression matrices from each cell type into a higher-order matrix and factorizes the result using the Tucker tensor decomposition. Applying scITD to single-cell RNA-sequencing data on 115 persons with lupus and 83 persons with coronavirus disease 2019, we identify patterns of coordinated cellular activity linked to disease severity and specific phenotypes, such as lupus nephritis. scITD results also implicate specific signaling pathways likely mediating coordination between cell types. Overall, scITD offers a tool for understanding the covariation of cell states across individuals, which can yield insights into the complex processes that define and stratify disease.

LIANA+ provides an all-in-one framework for cell-cell communication inference

The growing availability of single-cell and spatially resolved transcriptomics has led to the development of many approaches to infer cell-cell communication, each capturing only a partial view of the complex landscape of intercellular signalling. Here we present LIANA+, a scalable framework built around a rich knowledge base to decode coordinated inter- and intracellular signalling events from single- and multi-condition datasets in both single-cell and spatially resolved data. By extending and unifying established methodologies, LIANA+ provides a comprehensive set of synergistic components to study cell-cell communication via diverse molecular mediators, including those measured in multi-omics data. LIANA+ is accessible at https://github.com/saezlab/liana-py with extensive vignettes ( https://liana-py.readthedocs.io/ ) and provides an all-in-one solution to intercellular communication inference.

Multimodal analysis unveils tumor microenvironment heterogeneity linked to immune activity and evasion

The cellular and molecular heterogeneity of tumors is a major obstacle to cancer immunotherapy. Here, we use a systems biology approach to derive a signature of the main sources of heterogeneity in the tumor microenvironment (TME) from lung cancer transcriptomics. We demonstrate that this signature, which we called iHet, is conserved in different cancers and associated with antitumor immunity. Through analysis of single-cell and spatial transcriptomics data, we trace back the cellular origin of the variability explaining the iHet signature. Finally, we demonstrate that iHet has predictive value for cancer immunotherapy, which can be further improved by disentangling three major determinants of anticancer immune responses: activity of immune cells, immune infiltration or exclusion, and cancer-cell foreignness. This work shows how transcriptomics data can be integrated to derive a holistic representation of the phenotypic heterogeneity of the TME and to predict its unfolding and fate during immunotherapy with immune checkpoint blockers.

Integrating spatially-resolved transcriptomics data across tissues and individuals: challenges and opportunities

Advances in spatially-resolved transcriptomics (SRT) technologies have propelled the development of new computational analysis methods to unlock biological insights. As the cost of generating these data decreases, these technologies provide an exciting opportunity to create large-scale atlases that integrate SRT data across multiple tissues, individuals, species, or phenotypes to perform population-level analyses. Here, we describe unique challenges of varying spatial resolutions in SRT data, as well as highlight the opportunities for standardized preprocessing methods along with computational algorithms amenable to atlas-scale datasets leading to improved sensitivity and reproducibility in the future.

The diversification of methods for studying cell-cell interactions and communication

No cell lives in a vacuum, and the molecular interactions between cells define most phenotypes. Transcriptomics provides rich information to infer cell-cell interactions and communication, thus accelerating the discovery of the roles of cells within their communities. Such research relies heavily on algorithms that infer which cells are interacting and the ligands and receptors involved. Specific pressures on different research niches are driving the evolution of next-generation computational tools, enabling new conceptual opportunities and technological advances. More sophisticated algorithms now account for the heterogeneity and spatial organization of cells, multiple ligand types and intracellular signalling events, and enable the use of larger and more complex datasets, including single-cell and spatial transcriptomics. Similarly, new high-throughput experimental methods are increasing the number and resolution of interactions that can be analysed simultaneously. Here, we explore recent progress in cell-cell interaction research and highlight the diversification of the next generation of tools, which have yielded a rich ecosystem of tools for different applications and are enabling invaluable discoveries.

Complementing Cell Taxonomies with a Multicellular Analysis of Tissues

The application of single-cell molecular profiling coupled with spatial technologies has enabled charting of cellular heterogeneity in reference tissues and in disease. This new wave of molecular data has highlighted the expected diversity of single-cell dynamics upon shared external queues and spatial organizations. However, little is known about the relationship between single-cell heterogeneity and the emergence and maintenance of robust multicellular processes in developed tissues and its role in (patho)physiology. Here, we present emerging computational modeling strategies that use increasingly available large-scale cross-condition single-cell and spatial datasets to study multicellular organization in tissues and complement cell taxonomies. This perspective should enable us to better understand how cells within tissues collectively process information and adapt synchronized responses in disease contexts and to bridge the gap between structural changes and functions in tissues.

Combining LIANA and Tensor-cell2cell to decipher cell-cell communication across multiple samples

In recent years, data-driven inference of cell-cell communication has helped reveal coordinated biological processes across cell types. Here, we integrate two tools, LIANA and Tensor-cell2cell, which, when combined, can deploy multiple existing methods and resources to enable the robust and flexible identification of cell-cell communication programs across multiple samples. In this work, we show how the integration of our tools facilitates the choice of method to infer cell-cell communication and subsequently perform an unsupervised deconvolution to obtain and summarize biological insights. We explain how to perform the analysis step by step in both Python and R and provide online tutorials with detailed instructions available at https://ccc-protocols.readthedocs.io/. This workflow typically takes ∼1.5 h to complete from installation to downstream visualizations on a graphics processing unit-enabled computer for a dataset of ∼63,000 cells, 10 cell types, and 12 samples.

Integrating single-cell multi-omics and prior biological knowledge for a functional characterization of the immune system

The immune system comprises diverse specialized cell types that cooperate to defend the host against a wide range of pathogenic threats. Recent advancements in single-cell and spatial multi-omics technologies provide rich information about the molecular state of immune cells. Here, we review how the integration of single-cell and spatial multi-omics data with prior knowledge-gathered from decades of detailed biochemical studies-allows us to obtain functional insights, focusing on gene regulatory processes and cell-cell interactions. We present diverse applications in immunology and critically assess underlying assumptions and limitations. Finally, we offer a perspective on the ongoing technological and algorithmic developments that promise to get us closer to a systemic mechanistic understanding of the immune system.

Detection of PatIent-Level distances from single cell genomics and pathomics data with Optimal Transport (PILOT)

Although clinical applications represent the next challenge in single-cell genomics and digital pathology, we still lack computational methods to analyze single-cell or pathomics data to find sample-level trajectories or clusters associated with diseases. This remains challenging as single-cell/pathomics data are multi-scale, i.e., a sample is represented by clusters of cells/structures, and samples cannot be easily compared with each other. Here we propose PatIent Level analysis with Optimal Transport (PILOT). PILOT uses optimal transport to compute the Wasserstein distance between two individual single-cell samples. This allows us to perform unsupervised analysis at the sample level and uncover trajectories or cellular clusters associated with disease progression. We evaluate PILOT and competing approaches in single-cell genomics or pathomics studies involving various human diseases with up to 600 samples/patients and millions of cells or tissue structures. Our results demonstrate that PILOT detects disease-associated samples from large and complex single-cell or pathomics data. Moreover, PILOT provides a statistical approach to find changes in cell populations, gene expression, and tissue structures related to the trajectories or clusters supporting interpretation of predictions.

Multicellular factor analysis of single-cell data for a tissue-centric understanding of disease

Biomedical single-cell atlases describe disease at the cellular level. However, analysis of this data commonly focuses on cell-type-centric pairwise cross-condition comparisons, disregarding the multicellular nature of disease processes. Here, we propose multicellular factor analysis for the unsupervised analysis of samples from cross-condition single-cell atlases and the identification of multicellular programs associated with disease. Our strategy, which repurposes group factor analysis as implemented in multi-omics factor analysis, incorporates the variation of patient samples across cell-types or other tissue-centric features, such as cell compositions or spatial relationships, and enables the joint analysis of multiple patient cohorts, facilitating the integration of atlases. We applied our framework to a collection of acute and chronic human heart failure atlases and described multicellular processes of cardiac remodeling, independent to cellular compositions and their local organization, that were conserved in independent spatial and bulk transcriptomics datasets. In sum, our framework serves as an exploratory tool for unsupervised analysis of cross-condition single-cell atlases and allows for the integration of the measurements of patient cohorts across distinct data modalities.