Uncovering axes of variation among single-cell cancer specimens

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.

Mapping lineage-traced cells across time points with moslin

Simultaneous profiling of single-cell gene expression and lineage history holds enormous potential for studying cellular decision-making. Recent computational approaches combine both modalities into cellular trajectories; however, they cannot make use of all available lineage information in destructive time-series experiments. Here, we present moslin, a Gromov-Wasserstein-based model to couple cellular profiles across time points based on lineage and gene expression information. We validate our approach in simulations and demonstrate on Caenorhabditis elegans embryonic development how moslin predicts fate probabilities and putative decision driver genes. Finally, we use moslin to delineate lineage relationships among transiently activated fibroblast states during zebrafish heart regeneration.

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.

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.

Integrative, high-resolution analysis of single cells across experimental conditions with PARAFAC2

Effective tools for exploration and analysis are needed to extract insights from large-scale single-cell measurement data. However, current techniques for handling single-cell studies performed across experimental conditions (e.g., samples, perturbations, or patients) require restrictive assumptions, lack flexibility, or do not adequately deconvolute condition-to-condition variation from cell-to-cell variation. Here, we report that the tensor decomposition method PARAFAC2 (Pf2) enables the dimensionality reduction of single-cell data across conditions. We demonstrate these benefits across two distinct contexts of single-cell RNA-sequencing (scRNA-seq) experiments of peripheral immune cells: pharmacologic drug perturbations and systemic lupus erythematosus (SLE) patient samples. By isolating relevant gene modules across cells and conditions, Pf2 enables straightforward associations of gene variation patterns across specific patients or perturbations while connecting each coordinated change to certain cells without pre-defining cell types. The theoretical grounding of Pf2 suggests a unified framework for many modeling tasks associated with single-cell data. Thus, Pf2 provides an intuitive universal dimensionality reduction approach for multi-sample single-cell studies across diverse biological contexts.

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.

scPerturb: harmonized single-cell perturbation data

Analysis across a growing number of single-cell perturbation datasets is hampered by poor data interoperability. To facilitate development and benchmarking of computational methods, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform quality control pipelines and harmonize feature annotations. The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets. We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing, and demonstrate E-distance as a general distance measure between sets of single-cell expression profiles. We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations. The perturbation-response datasets and E-statistics computation software are publicly available at scperturb.org. This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.

QOT: Efficient Computation of Sample Level Distance Matrix from Single-Cell Omics Data through Quantized Optimal Transport

Single-cell technologies have emerged as a transformative technology enabling high-dimensional characterization of cell populations at an unprecedented scale. The data's innate complexity and voluminous nature 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 (OT) 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 enables 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 OT-based algorithms in terms of accuracy and robustness when obtaining a distance matrix at the sample level from high throughput single-cell measures. Moreover, the sample level distance matrix could be used in downstream analysis (i.e. uncover the trajectory of disease progression), highlighting its usage in biomedical informatics and data science.

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.

Trellis tree-based analysis reveals stromal regulation of patient-derived organoid drug responses

Patient-derived organoids (PDOs) can model personalized therapy responses; however, current screening technologies cannot reveal drug response mechanisms or how tumor microenvironment cells alter therapeutic performance. To address this, we developed a highly multiplexed mass cytometry platform to measure post-translational modification (PTM) signaling, DNA damage, cell-cycle activity, and apoptosis in >2,500 colorectal cancer (CRC) PDOs and cancer-associated fibroblasts (CAFs) in response to clinical therapies at single-cell resolution. To compare patient- and microenvironment-specific drug responses in thousands of single-cell datasets, we developed "Trellis"-a highly scalable, tree-based treatment effect analysis method. Trellis single-cell screening revealed that on-target cell-cycle blockage and DNA-damage drug effects are common, even in chemorefractory PDOs. However, drug-induced apoptosis is rarer, patient-specific, and aligns with cancer cell PTM signaling. We find that CAFs can regulate PDO plasticity-shifting proliferative colonic stem cells (proCSCs) to slow-cycling revival colonic stem cells (revCSCs) to protect cancer cells from chemotherapy.

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.

LIDER: cell embedding based deep neural network classifier for supervised cell type identification

Background

Automatic cell type identification has been an urgent task for the rapid development of single-cell RNA-seq techniques. Generally, the current approach for cell type identification is to generate cell clusters by unsupervised clustering and later assign labels to each cell cluster with manual annotation.

Methods

Here, we introduce LIDER (celL embeddIng based Deep nEural netwoRk classifier), a deep supervised learning method that combines cell embedding and deep neural network classifier for automatic cell type identification. Based on a stacked denoising autoencoder with a tailored and reconstructed loss function, LIDER identifies cell embedding and predicts cell types with a deep neural network classifier. LIDER was developed upon a stacked denoising autoencoder to learn encoder-decoder structures for identifying cell embedding.

Results

LIDER accurately identifies cell types by using stacked denoising autoencoder. Benchmarking against state-of-the-art methods across eight types of single-cell data, LIDER achieves comparable or even superior enhancement performance. Moreover, LIDER suggests comparable robust to batch effects. Our results show a potential in deep supervised learning for automatic cell type identification of single-cell RNA-seq data. The LIDER codes are available at https://github.com/ShiMGLab/LIDER.

Multiscale geometric and topological analyses for characterizing and predicting immune responses from single cell data

Single cell genomics has revolutionized our ability to map immune heterogeneity and responses. With the influx of large-scale data sets from diverse modalities, the resolution achieved has supported the long-held notion that immune cells are naturally organized into hierarchical relationships, characterized at multiple levels. Such a multigranular structure corresponds to key geometric and topological features. Given that differences between an effective and ineffective immunological response may not be found at one level, there is vested interest in characterizing and predicting outcomes from such features. In this review, we highlight single cell methods and principles for learning geometric and topological properties of data at multiple scales, discussing their contributions to immunology. Ultimately, multiscale approaches go beyond classical clustering, revealing a more comprehensive picture of cellular heterogeneity.

Transcriptional space-time mapping identifies concerted immune and stromal cell patterns and gene programs in wound healing and cancer

Tissue repair responses in metazoans are highly coordinated by different cell types over space and time. However, comprehensive single-cell-based characterization covering this coordination is lacking. Here, we captured transcriptional states of single cells over space and time during skin wound closure, revealing choreographed gene-expression profiles. We identified shared space-time patterns of cellular and gene program enrichment, which we call multicellular "movements" spanning multiple cell types. We validated some of the discovered space-time movements using large-volume imaging of cleared wounds and demonstrated the value of this analysis to predict "sender" and "receiver" gene programs in macrophages and fibroblasts. Finally, we tested the hypothesis that tumors are like "wounds that never heal" and found conserved wound healing movements in mouse melanoma and colorectal tumor models, as well as human tumor samples, revealing fundamental multicellular units of tissue biology for integrative studies.

Unraveling non-genetic heterogeneity in cancer with dynamical models and computational tools

Individual cells within an otherwise genetically homogenous population constantly undergo fluctuations in their molecular state, giving rise to non-genetic heterogeneity. Such diversity is being increasingly implicated in cancer therapy resistance and metastasis. Identifying the origins of non-genetic heterogeneity is therefore crucial for making clinical breakthroughs. We discuss with examples how dynamical models and computational tools have provided critical multiscale insights into the nature and consequences of non-genetic heterogeneity in cancer. We demonstrate how mechanistic modeling has been pivotal in establishing key concepts underlying non-genetic diversity at various biological scales, from population dynamics to gene regulatory networks. We discuss advances in single-cell longitudinal profiling techniques to reveal patterns of non-genetic heterogeneity, highlighting the ongoing efforts and challenges in statistical frameworks to robustly interpret such multimodal datasets. Moving forward, we stress the need for data-driven statistical and mechanistically motivated dynamical frameworks to come together to develop predictive cancer models and inform therapeutic strategies.

Targeting the Lysosomal Degradation of Rab22a-NeoF1 Fusion Protein for Osteosarcoma Lung Metastasis

Rab22a-NeoF fusion protein has recently been reported as a promising target for osteosarcoma lung metastasis. However, how this fusion protein is regulated in cells remains unknown. Here, using multiple screenings, it is reported that Rab22a-NeoF1 fusion protein is degraded by an E3 ligase STUB1 via the autophagy receptor NDP52-mediated lysosome pathway, which is facilitated by PINK1 kinase. Mechanistically, STUB1 catalyzes the K63-linked ubiquitin chains on lysine112 of Rab22a-NeoF1, which is responsible for the binding of Rab22a-NeoF1 to NDP52, resulting in lysosomal degradation of Rab22a-NeoF1. PINK1 is able to phosphorylate Rab22a-NeoF1 at serine120, which promotes ubiquitination and degradation of Rab22a-NeoF1. Consistently, by upregulating PINK1, Sorafenib and Regorafenib can inhibit osteosarcoma lung metastasis induced by Rab22a-NeoF1. These findings reveal that the lysosomal degradation of Rab22a-NeoF1 fusion protein is targetable for osteosarcoma lung metastasis, proposing that Sorafenib and Regorafenib may benefit cancer patients who are positive for the RAB22A-NeoF1 fusion gene.

Next-Generation Morphometry for pathomics-data mining in histopathology

Pathology diagnostics relies on the assessment of morphology by trained experts, which remains subjective and qualitative. Here we developed a framework for large-scale histomorphometry (FLASH) performing deep learning-based semantic segmentation and subsequent large-scale extraction of interpretable, quantitative, morphometric features in non-tumour kidney histology. We use two internal and three external, multi-centre cohorts to analyse over 1000 kidney biopsies and nephrectomies. By associating morphometric features with clinical parameters, we confirm previous concepts and reveal unexpected relations. We show that the extracted features are independent predictors of long-term clinical outcomes in IgA-nephropathy. We introduce single-structure morphometric analysis by applying techniques from single-cell transcriptomics, identifying distinct glomerular populations and morphometric phenotypes along a trajectory of disease progression. Our study provides a concept for Next-generation Morphometry (NGM), enabling comprehensive quantitative pathology data mining, i.e., pathomics.

Control of cell state transitions

Understanding cell state transitions and purposefully controlling them is a longstanding challenge in biology. Here we present cell state transition assessment and regulation (cSTAR), an approach for mapping cell states, modelling transitions between them and predicting targeted interventions to convert cell fate decisions. cSTAR uses omics data as input, classifies cell states, and develops a workflow that transforms the input data into mechanistic models that identify a core signalling network, which controls cell fate transitions by influencing whole-cell networks. By integrating signalling and phenotypic data, cSTAR models how cells manoeuvre in Waddington's landscape1 and make decisions about which cell fate to adopt. Notably, cSTAR devises interventions to control the movement of cells in Waddington's landscape. Testing cSTAR in a cellular model of differentiation and proliferation shows a high correlation between quantitative predictions and experimental data. Applying cSTAR to different types of perturbation and omics datasets, including single-cell data, demonstrates its flexibility and scalability and provides new biological insights. The ability of cSTAR to identify targeted perturbations that interconvert cell fates will enable designer approaches for manipulating cellular development pathways and mechanistically underpinned therapeutic interventions.

SCANCell reveals diverse inter-cluster interaction patterns in systemic lupus erythematosus across the disease spectrum

MOTIVATION

High-dimensional mass cytometry (CyTOF), which provides both cellular signatures and inter-cluster interactions like the antagonism between immune activation and suppression, and the pro-inflammatory synergy, sheds light on the cellular and molecular basis of disease pathogenesis. However, revealing the aberrance of inter-cluster communication networks in CyTOF datasets remains a significant challenge.

RESULTS

Here, we developed Sample Classification and direct Association Network among Cell clusters (SCANCell) that quantifies the direct association (DA) network of cell clusters. SCANCell was applied to profile inter-cluster interaction patterns of a well-recruited systemic lupus erythematosus (SLE) cohort, including 8 healthy controls, 10 active SLE patients (APs) and 8 remission SLE patients (RPs). SCANCell identified decreased inter-cluster interactions of CD8+ T cells in APs compared with RPs, and enhanced DA of CD8+ T cells after stimulation with immunostimulatory cytokine interleukin-2 in vitro. These discoveries prove that SCANCell can uncover pathology- and drug stimulation-associated inter-cluster interactions, which potentially benefits understanding of pathogenesis and novel therapeutic strategies.

AVAILABILITY AND IMPLEMENTATION

The main processing scripts of SCNACell are available at https://github.com/Lxc417/SCANCell. Other codes for the following data statistics are available from the corresponding author upon request.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Alignment of single-cell trajectories by tuMap enables high-resolution quantitative comparison of cancer samples

Single-cell technologies allow characterization of cancer samples as continuous developmental trajectories. Yet, the obtained temporal resolution cannot be leveraged for a comparative analysis due to the large phenotypic heterogeneity existing between patients. Here, we present the tuMap algorithm that exploits high-dimensional single-cell data of cancer samples exhibiting an underlying developmental structure to align them with the healthy development, yielding the tuMap pseudotime axis that allows their systematic, meaningful comparison. We applied tuMap on single-cell mass cytometry data of acute lymphoblastic and myeloid leukemia to reveal associations between the tuMap pseudotime axis and clinics that outperform cellular assignment into developmental populations. Application of the tuMap algorithm on single-cell RNA sequencing data further identified gene signatures of stem cells residing at the very-early parts of the cancer trajectories. The quantitative framework provided by tuMap allows generation of metrics for cancer patients evaluation.

Can Systems Biology Advance Clinical Precision Oncology?

Precision oncology is perceived as a way forward to treat individual cancer patients. However, knowing particular cancer mutations is not enough for optimal therapeutic treatment, because cancer genotype-phenotype relationships are nonlinear and dynamic. Systems biology studies the biological processes at the systems' level, using an array of techniques, ranging from statistical methods to network reconstruction and analysis, to mathematical modeling. Its goal is to reconstruct the complex and often counterintuitive dynamic behavior of biological systems and quantitatively predict their responses to environmental perturbations. In this paper, we review the impact of systems biology on precision oncology. We show examples of how the analysis of signal transduction networks allows to dissect resistance to targeted therapies and inform the choice of combinations of targeted drugs based on tumor molecular alterations. Patient-specific biomarkers based on dynamical models of signaling networks can have a greater prognostic value than conventional biomarkers. These examples support systems biology models as valuable tools to advance clinical and translational oncological research.

Machine learning for perturbational single-cell omics

Cell biology is fundamentally limited in its ability to collect complete data on cellular phenotypes and the wide range of responses to perturbation. Areas such as computer vision and speech recognition have addressed this problem of characterizing unseen or unlabeled conditions with the combined advances of big data, deep learning, and computing resources in the past 5 years. Similarly, recent advances in machine learning approaches enabled by single-cell data start to address prediction tasks in perturbation response modeling. We first define objectives in learning perturbation response in single-cell omics; survey existing approaches, resources, and datasets (https://github.com/theislab/sc-pert); and discuss how a perturbation atlas can enable deep learning models to construct an informative perturbation latent space. We then examine future avenues toward more powerful and explainable modeling using deep neural networks, which enable the integration of disparate information sources and an understanding of heterogeneous, complex, and unseen systems.

Global absence and targeting of protective immune states in severe COVID-19

Although infection with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has pleiotropic and systemic effects in some individuals1-3, many others experience milder symptoms. Here, to gain a more comprehensive understanding of the distinction between severe and mild phenotypes in the pathology of coronavirus disease 2019 (COVID-19) and its origins, we performed a whole-blood-preserving single-cell analysis protocol to integrate contributions from all major immune cell types of the blood-including neutrophils, monocytes, platelets, lymphocytes and the contents of the serum. Patients with mild COVID-19 exhibit a coordinated pattern of expression of interferon-stimulated genes (ISGs)3 across every cell population, whereas these ISG-expressing cells are systemically absent in patients with severe disease. Paradoxically, individuals with severe COVID-19 produce very high titres of anti-SARS-CoV-2 antibodies and have a lower viral load compared to individuals with mild disease. Examination of the serum from patients with severe COVID-19 shows that these patients uniquely produce antibodies that functionally block the production of the ISG-expressing cells associated with mild disease, by activating conserved signalling circuits that dampen cellular responses to interferons. Overzealous antibody responses pit the immune system against itself in many patients with COVID-19, and perhaps also in individuals with other viral infections. Our findings reveal potential targets for immunotherapies in patients with severe COVID-19 to re-engage viral defence.

Discriminating the single-cell gene regulatory networks of human pancreatic islets: A novel deep learning application

Analysis of single-cell pancreatic data can play an important role in understanding various metabolic diseases and health conditions. Due to the sparsity and noise present in such single-cell gene expression data, inference of single-cell gene regulatory networks remains a challenge. Since recent studies have reported the reliable inference of single-cell gene regulatory networks (SCGRNs), the current study focused on discriminating the SCGRNs of T2D patients from those of healthy controls. By accurately distinguishing SCGRNs of healthy pancreas from those of T2D pancreas, it would be possible to annotate, organize, visualize, and identify common patterns of SCGRNs in metabolic diseases. Such annotated SCGRNs could play an important role in accelerating the process of building large data repositories. This study aimed to contribute to the development of a novel deep learning (DL) application. First, we generated a dataset consisting of 224 SCGRNs belonging to both T2D and healthy pancreas and made it freely available. Next, we chose seven DL architectures, including VGG16, VGG19, Xception, ResNet50, ResNet101, DenseNet121, and DenseNet169, trained each of them on the dataset, and checked their prediction based on a test set. Of note, we evaluated the DL architectures on a single NVIDIA GeForce RTX 2080Ti GPU. Experimental results on the whole dataset, using several performance measures, demonstrated the superiority of VGG19 DL model in the automatic classification of SCGRNs, derived from the single-cell pancreatic data.

Structure-Guided Computational Methods Predict Multiple Distinct Binding Modes for Pyrazoloquinolinones in GABAA Receptors

Pyrazoloquinolinones (PQs) are a versatile class of GABAA receptor ligands. It has been demonstrated that high functional selectivity for certain receptor subtypes can be obtained by specific substitution patterns, but so far, no clear SAR rules emerge from the studies. As is the case for many GABAA receptor targeting chemotypes, PQs can interact with distinct binding sites on a given receptor pentamer. In pentamers of αβγ composition, such as the most abundant α1β2γ2 subtype, many PQs are high affinity binders of the benzodiazepine binding site at the extracellular α+/γ2- interfaces. There they display a functionally near silent, flumazenil-like allosteric activity. More recently, interactions with extracellular α+/β- interfaces have been investigated, where strong positive modulation can be steered toward interesting subtype preferences. The most prominent examples are functionally α6-selective PQs. Similar to benzodiazepines, PQs also seem to interact with sites in the transmembrane domain, mainly the sites used by etomidate and barbiturates. This promiscuity leads to potential contributions from multiple sites to net modulation. Developing ligands that interact exclusively with the extracellular α+/β- interfaces would be desired. Correlating functional profiles with binding sites usage is hampered by scarce and heterogeneous experimental data, as shown in our meta-analysis of aggregated published data. In the absence of experimental structures, bound states can be predicted with pharmacophore matching methods and with computational docking. We thus performed pharmacophore matching studies for the unwanted sites, and computational docking for the extracellular α1,6+/β3- interfaces. The results suggest that PQs interact with their binding sites with diverse binding modes. As such, rational design of improved ligands needs to take a complex structure-activity landscape with branches between sub-series of derivatives into account. We present a workflow, which is suitable to identify and explore potential branching points on the structure-activity landscape of any small molecule chemotype.

Defining the Teratoma as a Model for Multi-lineage Human Development

We propose that the teratoma, a recognized standard for validating pluripotency in stem cells, could be a promising platform for studying human developmental processes. Performing single-cell RNA sequencing (RNA-seq) of 179,632 cells across 23 teratomas from 4 cell lines, we found that teratomas reproducibly contain approximately 20 cell types across all 3 germ layers, that inter-teratoma cell type heterogeneity is comparable with organoid systems, and teratoma gut and brain cell types correspond well to similar fetal cell types. Furthermore, cellular barcoding confirmed that injected stem cells robustly engraft and contribute to all lineages. Using pooled CRISPR-Cas9 knockout screens, we showed that teratomas can enable simultaneous assaying of the effects of genetic perturbations across all germ layers. Additionally, we demonstrated that teratomas can be sculpted molecularly via microRNA (miRNA)-regulated suicide gene expression to enrich for specific tissues. Taken together, teratomas are a promising platform for modeling multi-lineage development, pan-tissue functional genetic screening, and tissue engineering.

Global Absence and Targeting of Protective Immune States in Severe COVID-19

While SARS-CoV-2 infection has pleiotropic and systemic effects in some patients, many others experience milder symptoms. We sought a holistic understanding of the severe/mild distinction in COVID-19 pathology, and its origins. We performed a whole-blood preserving single-cell analysis protocol to integrate contributions from all major cell types including neutrophils, monocytes, platelets, lymphocytes and the contents of serum. Patients with mild COVID-19 disease display a coordinated pattern of interferon-stimulated gene (ISG) expression across every cell population and these cells are systemically absent in patients with severe disease. Severe COVID-19 patients also paradoxically produce very high anti-SARS-CoV-2 antibody titers and have lower viral load as compared to mild disease. Examination of the serum from severe patients demonstrates that they uniquely produce antibodies with multiple patterns of specificity against interferon-stimulated cells and that those antibodies functionally block the production of the mild disease-associated ISG-expressing cells. Overzealous and auto-directed antibody responses pit the immune system against itself in many COVID-19 patients and this defines targets for immunotherapies to allow immune systems to provide viral defense.

One Sentence Summary

In severe COVID-19 patients, the immune system fails to generate cells that define mild disease; antibodies in their serum actively prevents the successful production of those cells.

Context specificity of the EMT transcriptional response

Epithelial–mesenchymal plasticity contributes to many biological processes, including tumor progression. Various epithelial–mesenchymal transition (EMT) responses have been reported and no common, EMT-defining gene expression program has been identified. Here, we have performed a comparative analysis of the EMT response, leveraging highly multiplexed single-cell RNA sequencing (scRNA-seq) to measure expression profiles of 103,999 cells from 960 samples, comprising 12 EMT time course experiments and independent kinase inhibitor screens for each. We demonstrate that the EMT is vastly context specific, with an average of only 22% of response genes being shared between any two conditions, and over half of all response genes were restricted to 1–2 time course experiments. Further, kinase inhibitor screens revealed signaling dependencies and modularity of these responses. These findings suggest that the EMT is not simply a single, linear process, but is highly variable and modular, warranting quantitative frameworks for understanding nuances of the transition.

It is unclear if a common EMT expression program exists. Here, the authors perform multiplexed single-cell RNA sequencing across 12 EMT time courses and 16 kinase inhibitor screens, and find that EMT transcriptional responses are context specific and EMT is not a single, linear transition.