Compressed sensing for highly efficient imaging transcriptomics

Interpretable deep learning for deconvolutional analysis of neural signals

The widespread adoption of deep learning to build models that capture the dynamics of neural populations is typically based on "black-box" approaches that lack an interpretable link between neural activity and network parameters. Here, we propose to apply algorithm unrolling, a method for interpretable deep learning, to design the architecture of sparse deconvolutional neural networks and obtain a direct interpretation of network weights in relation to stimulus-driven single-neuron activity through a generative model. We characterize our method, referred to as deconvolutional unrolled neural learning (DUNL), and show its versatility by applying it to deconvolve single-trial local signals across multiple brain areas and recording modalities. To exemplify use cases of our decomposition method, we uncover multiplexed salience and reward prediction error signals from midbrain dopamine neurons in an unbiased manner, perform simultaneous event detection and characterization in somatosensory thalamus recordings, and characterize the heterogeneity of neural responses in the piriform cortex and in the striatum during unstructured, naturalistic experiments. Our work leverages the advances in interpretable deep learning to gain a mechanistic understanding of neural activity.

Spatial Single-Cell Mapping of Transcriptional Differences Across Genetic Backgrounds in Mouse Brains

Genetic variation can alter brain structure and, consequently, function. Comparative statistical analysis of mouse brains across genetic backgrounds requires spatial, single-cell, atlas-scale data, in replicates-a challenge for current technologies. We introduce Atlas-scale Transcriptome Localization using Aggregate Signatures (ATLAS), a scalable tissue mapping method. ATLAS learns transcriptional signatures from scRNAseq data, encodes them in situ with tens of thousands of oligonucleotide probes, and decodes them to infer cell types and imputed transcriptomes. We validated ATLAS by comparing its cell type inferences with direct MERFISH measurements of marker genes and quantitative comparisons to four other technologies. Using ATLAS, we mapped the central brains of five male and five female C57BL/6J (B6) mice and five male BTBR T+ tf/J (BTBR) mice, an idiopathic model of autism, collectively profiling over 40 million cells across over 400 coronal sections. Our analysis revealed over 40 significant differences in cell type distributions and identified 16 regional composition changes across male-female and B6-BTBR comparisons. ATLAS thus enables systematic comparative studies, facilitating organ-level structure-function analysis of disease models.

Multiplexed, image-based pooled screens in primary cells and tissues with PerturbView

Optical pooled screening (OPS) is a scalable method for linking image-based phenotypes with cellular perturbations. However, it has thus far been restricted to relatively low-plex phenotypic readouts in cancer cell lines in culture due to limitations associated with in situ sequencing of perturbation barcodes. Here, we develop PerturbView, an OPS technology that leverages in vitro transcription to amplify barcodes before in situ sequencing, enabling screens with highly multiplexed phenotypic readouts across diverse systems, including primary cells and tissues. We demonstrate PerturbView in induced pluripotent stem cell-derived neurons, primary immune cells and tumor tissue sections from animal models. In a screen of immune signaling pathways in primary bone marrow-derived macrophages, PerturbView uncovered both known and novel regulators of NF-κB signaling. Furthermore, we combine PerturbView with spatial transcriptomics in tissue sections from a mouse xenograft model, paving the way to in situ screens with rich optical and transcriptomic phenotypes. PerturbView broadens the scope of OPS to a wide range of models and applications.

Toward a foundation model of causal cell and tissue biology with a Perturbation Cell and Tissue Atlas

Comprehensively charting the biologically causal circuits that govern the phenotypic space of human cells has often been viewed as an insurmountable challenge. However, in the last decade, a suite of interleaved experimental and computational technologies has arisen that is making this fundamental goal increasingly tractable. Pooled CRISPR-based perturbation screens with high-content molecular and/or image-based readouts are now enabling researchers to probe, map, and decipher genetically causal circuits at increasing scale. This scale is now eminently suitable for the deployment of artificial intelligence and machine learning (AI/ML) to both direct further experiments and to predict or generate information that was not-and sometimes cannot-be gathered experimentally. By combining and iterating those through experiments that are designed for inference, we now envision a Perturbation Cell Atlas as a generative causal foundation model to unify human cell biology.

Scalable genetic screening for regulatory circuits using compressed Perturb-seq

Pooled CRISPR screens with single-cell RNA sequencing readout (Perturb-seq) have emerged as a key technique in functional genomics, but they are limited in scale by cost and combinatorial complexity. In this study, we modified the design of Perturb-seq by incorporating algorithms applied to random, low-dimensional observations. Compressed Perturb-seq measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq with an order of magnitude cost reduction and greater power to learn genetic interactions. We identified known and novel regulators of immune responses and uncovered evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing genome-wide association studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

Highly sensitive spatial transcriptomics using FISHnCHIPs of multiple co-expressed genes

High-dimensional, spatially resolved analysis of intact tissue samples promises to transform biomedical research and diagnostics, but existing spatial omics technologies are costly and labor-intensive. We present Fluorescence In Situ Hybridization of Cellular HeterogeneIty and gene expression Programs (FISHnCHIPs) for highly sensitive in situ profiling of cell types and gene expression programs. FISHnCHIPs achieves this by simultaneously imaging ~2-35 co-expressed genes (clustered into modules) that are spatially co-localized in tissues, resulting in similar spatial information as single-gene Fluorescence In Situ Hybridization (FISH), but with ~2-20-fold higher sensitivity. Using FISHnCHIPs, we image up to 53 modules from the mouse kidney and mouse brain, and demonstrate high-speed, large field-of-view profiling of a whole tissue section. FISHnCHIPs also reveals spatially restricted localizations of cancer-associated fibroblasts in a human colorectal cancer biopsy. Overall, FISHnCHIPs enables fast, robust, and scalable cell typing of tissues with normal physiology or undergoing pathogenesis.

scGIST: gene panel design for spatial transcriptomics with prioritized gene sets

A critical challenge of single-cell spatial transcriptomics (sc-ST) technologies is their panel size. Being based on fluorescence in situ hybridization, they are typically limited to panels of about a thousand genes. This constrains researchers to build panels from only the marker genes of different cell types and forgo other genes of interest, e.g., genes encoding ligand-receptor complexes or those in specific pathways. We propose scGIST, a constrained feature selection tool that designs sc-ST panels prioritizing user-specified genes without compromising cell type detection accuracy. We demonstrate scGIST's efficacy in diverse use cases, highlighting it as a valuable addition to sc-ST's algorithmic toolbox.

Multiplex, Quantitative, High-Resolution Imaging of Protein:Protein Complexes via Hybridization Chain Reaction

Signal amplification based on the mechanism of hybridization chain reaction (HCR) facilitates spatial exploration of gene regulatory networks by enabling multiplex, quantitative, high-resolution imaging of RNA and protein targets. Here, we extend these capabilities to the imaging of protein:protein complexes, using proximity-dependent cooperative probes to conditionally generate a single amplified signal if and only if two target proteins are colocalized within the sample. HCR probes and amplifiers combine to provide automatic background suppression throughout the protocol, ensuring that even if reagents bind nonspecifically in the sample, they will not generate amplified background. We demonstrate protein:protein imaging with a high signal-to-background ratio in human cells, mouse proT cells, and highly autofluorescent formalin-fixed paraffin-embedded (FFPE) human breast tissue sections. Further, we demonstrate multiplex imaging of three different protein:protein complexes simultaneously and validate that HCR enables accurate and precise relative quantitation of protein:protein complexes with subcellular resolution in an anatomical context. Moreover, we establish a unified framework for simultaneous multiplex, quantitative, high-resolution imaging of RNA, protein, and protein:protein targets, with one-step, isothermal, enzyme-free HCR signal amplification performed for all target classes simultaneously.

Spatial profiling technologies illuminate the tumor microenvironment

The tumor microenvironment (TME) is composed of many different cellular and acellular components that together drive tumor growth, invasion, metastasis, and response to therapies. Increasing realization of the significance of the TME in cancer biology has shifted cancer research from a cancer-centric model to one that considers the TME as a whole. Recent technological advancements in spatial profiling methodologies provide a systematic view and illuminate the physical localization of the components of the TME. In this review, we provide an overview of major spatial profiling technologies. We present the types of information that can be extracted from these data and describe their applications, findings and challenges in cancer research. Finally, we provide a future perspective of how spatial profiling could be integrated into cancer research to improve patient diagnosis, prognosis, stratification to treatment and development of novel therapeutics.

Embracing enzyme promiscuity with activity-based compressed biosensing

The development of protease-activatable drugs and diagnostics requires identifying substrates specific to individual proteases. However, this process becomes increasingly difficult as the number of target proteases increases because most substrates are promiscuously cleaved by multiple proteases. We introduce a method-substrate libraries for compressed sensing of enzymes (SLICE)-for selecting libraries of promiscuous substrates that classify protease mixtures (1) without deconvolution of compressed signals and (2) without highly specific substrates. SLICE ranks substrate libraries using a compression score (C), which quantifies substrate orthogonality and protease coverage. This metric is predictive of classification accuracy across 140 in silico (Pearson r = 0.71) and 55 in vitro libraries (r = 0.55). Using SLICE, we select a two-substrate library to classify 28 samples containing 11 enzymes in plasma (area under the receiver operating characteristic curve [AUROC] = 0.93). We envision that SLICE will enable the selection of libraries that capture information from hundreds of enzymes using fewer substrates for applications like activity-based sensors for imaging and diagnostics.

Compressed Perturb-seq: highly efficient screens for regulatory circuits using random composite perturbations

Pooled CRISPR screens with single-cell RNA-seq readout (Perturb-seq) have emerged as a key technique in functional genomics, but are limited in scale by cost and combinatorial complexity. Here, we reimagine Perturb-seq's design through the lens of algorithms applied to random, low-dimensional observations. We present compressed Perturb-seq, which measures multiple random perturbations per cell or multiple cells per droplet and computationally decompresses these measurements by leveraging the sparse structure of regulatory circuits. Applied to 598 genes in the immune response to bacterial lipopolysaccharide, compressed Perturb-seq achieves the same accuracy as conventional Perturb-seq at 4 to 20-fold reduced cost, with greater power to learn genetic interactions. We identify known and novel regulators of immune responses and uncover evolutionarily constrained genes with downstream targets enriched for immune disease heritability, including many missed by existing GWAS or trans-eQTL studies. Our framework enables new scales of interrogation for a foundational method in functional genomics.

Pooled genetic screens with image-based profiling

Spatial structure in biology, spanning molecular, organellular, cellular, tissue, and organismal scales, is encoded through a combination of genetic and epigenetic factors in individual cells. Microscopy remains the most direct approach to exploring the intricate spatial complexity defining biological systems and the structured dynamic responses of these systems to perturbations. Genetic screens with deep single-cell profiling via image features or gene expression programs have the capacity to show how biological systems work in detail by cataloging many cellular phenotypes with one experimental assay. Microscopy-based cellular profiling provides information complementary to next-generation sequencing (NGS) profiling and has only recently become compatible with large-scale genetic screens. Optical screening now offers the scale needed for systematic characterization and is poised for further scale-up. We discuss how these methodologies, together with emerging technologies for genetic perturbation and microscopy-based multiplexed molecular phenotyping, are powering new approaches to reveal genotype-phenotype relationships.

Computational multiplex panel reduction to maximize information retention in breast cancer tissue microarrays

Recent state-of-the-art multiplex imaging techniques have expanded the depth of information that can be captured within a single tissue sample by allowing for panels with dozens of markers. Despite this increase in capacity, space on the panel is still limited due to technical artifacts, tissue loss, and long imaging acquisition time. As such, selecting which markers to include on a panel is important, since removing important markers will result in a loss of biologically relevant information, but identifying redundant markers will provide a room for other markers. To address this, we propose computational approaches to determine the amount of shared information between markers and select an optimally reduced panel that captures maximum amount of information with the fewest markers. Here we examine several panel selection approaches and evaluate them based on their ability to reconstruct the full panel images and information within breast cancer tissue microarray datasets using cyclic immunofluorescence as a proof of concept. We show that all methods perform adequately and can re-capture cell types using only 18 of 25 markers (72% of the original panel size). The correlation-based selection methods achieved the best single-cell marker mean intensity predictions with a Spearman correlation of 0.90 with the reduced panel. Using the proposed methods shown here, it is possible for researchers to design more efficient multiplex imaging panels that maximize the amount of information retained with the limited number of markers with respect to certain evaluation metrics and architecture biases.

The emerging landscape of spatial profiling technologies

Improved scale, multiplexing and resolution are establishing spatial nucleic acid and protein profiling methods as a major pillar for cellular atlas building of complex samples, from tissues to full organisms. Emerging methods yield omics measurements at resolutions covering the nano- to microscale, enabling the charting of cellular heterogeneity, complex tissue architectures and dynamic changes during development and disease. We present an overview of the developing landscape of in situ spatial genome, transcriptome and proteome technologies, exemplify their impact on cell biology and translational research, and discuss current challenges for their community-wide adoption. Among many transformative applications, we envision that spatial methods will map entire organs and enable next-generation pathology.

CoSpar identifies early cell fate biases from single-cell transcriptomic and lineage information

A goal of single-cell genome-wide profiling is to reconstruct dynamic transitions during cell differentiation, disease onset and drug response. Single-cell assays have recently been integrated with lineage tracing, a set of methods that identify cells of common ancestry to establish bona fide dynamic relationships between cell states. These integrated methods have revealed unappreciated cell dynamics, but their analysis faces recurrent challenges arising from noisy, dispersed lineage data. In this study, we developed coherent, sparse optimization (CoSpar) as a robust computational approach to infer cell dynamics from single-cell transcriptomics integrated with lineage tracing. Built on assumptions of coherence and sparsity of transition maps, CoSpar is robust to severe downsampling and dispersion of lineage data, which enables simpler experimental designs and requires less calibration. In datasets representing hematopoiesis, reprogramming and directed differentiation, CoSpar identifies early fate biases not previously detected, predicting transcription factors and receptors implicated in fate choice. Documentation and detailed examples for common experimental designs are available at https://cospar.readthedocs.io/ .

Next-generation sequencing: unraveling genetic mechanisms that shape cancer immunotherapy efficacy

Immunity is governed by fundamental genetic processes. These processes shape the nature of immune cells and set the rules that dictate the myriad complex cellular interactions that power immune systems. Everything from the generation of T cell receptors and antibodies, control of epitope presentation, and recognition of pathogens by the immunoediting of cancer cells is, in large part, made possible by core genetic mechanisms and the cellular machinery that they encode. In the last decade, next-generation sequencing has been used to dissect the complexities of cancer immunity with potent effect. Sequencing of exomes and genomes has begun to reveal how the immune system recognizes “foreign” entities and distinguishes self from non-self, especially in the setting of cancer. High-throughput analyses of transcriptomes have revealed deep insights into how the tumor microenvironment affects immunotherapy efficacy. In this Review, we discuss how high-throughput sequencing has added to our understanding of how immune systems interact with cancer cells and how cancer immunotherapies work.

Direct RNA targeted in situ sequencing for transcriptomic profiling in tissue

Highly multiplexed spatial mapping of transcripts within tissues allows for investigation of the transcriptomic and cellular diversity of mammalian organs previously unseen. Here we explore a direct RNA (dRNA) detection approach incorporating the use of padlock probes and rolling circle amplification in combination with hybridization-based in situ sequencing chemistry. We benchmark a High Sensitivity Library Preparation Kit from CARTANA that circumvents the reverse transcription needed for cDNA-based in situ sequencing (ISS) via direct RNA detection. We found a fivefold increase in transcript detection efficiency when compared to cDNA-based ISS and also validated its multiplexing capability by targeting a curated panel of 50 genes from previous publications on mouse brain sections, leading to additional data interpretation such as de novo cell clustering. With this increased efficiency, we also found to maintain specificity, multiplexing capabilities and ease of implementation. Overall, the dRNA chemistry shows significant improvements in target detection efficiency, closing the gap to other fluorescent in situ hybridization-based technologies and opens up possibilities to explore new biological questions previously not possible with cDNA-based ISS.

Enhancing the Throughput of FT Mass Spectrometry Imaging Using Joint Compressed Sensing and Subspace Modeling

Mass spectrometry imaging (MSI) allows for untargeted mapping of the chemical composition of tissues with attomole detection limits. MSI using Fourier transform (FT)-based mass spectrometers, such as FT-ion cyclotron resonance (FT-ICR), grants the ability to examine the chemical space with unmatched mass resolution and mass accuracy. However, direct imaging of large tissue samples using FT-ICR is slow. In this work, we present an approach that combines the subspace modeling of ICR temporal signals with compressed sensing to accelerate high-resolution FT-ICR MSI. A joint subspace and spatial sparsity constrained model computationally reconstructs high-resolution MSI data from the sparsely sampled transients with reduced duration, allowing a significant reduction in imaging time. Simulation studies and experimental implementation of the proposed method in investigation of brain tissues demonstrate a 10-fold enhancement in throughput of FT-ICR MSI, without the need for instrumental or hardware modifications.

Spatial components of molecular tissue biology

Methods for profiling RNA and protein expression in a spatially resolved manner are rapidly evolving, making it possible to comprehensively characterize cells and tissues in health and disease. To maximize the biological insights obtained using these techniques, it is critical to both clearly articulate the key biological questions in spatial analysis of tissues and develop the requisite computational tools to address them. Developers of analytical tools need to decide on the intrinsic molecular features of each cell that need to be considered, and how cell shape and morphological features are incorporated into the analysis. Also, optimal ways to compare different tissue samples at various length scales are still being sought. Grouping these biological problems and related computational algorithms into classes across length scales, thus characterizing common issues that need to be addressed, will facilitate further progress in spatial transcriptomics and proteomics.

Zero-preserving imputation of single-cell RNA-seq data

A key challenge in analyzing single cell RNA-sequencing data is the large number of false zeros, where genes actually expressed in a given cell are incorrectly measured as unexpressed. We present a method based on low-rank matrix approximation which imputes these values while preserving biologically non-expressed genes (true biological zeros) at zero expression levels. We provide theoretical justification for this denoising approach and demonstrate its advantages relative to other methods on simulated and biological datasets.

Evaluating microglial phenotypes using single-cell technologies

Recent single-cell technologies have enabled researchers to simultaneously assess the transcriptomes and other modalities of thousands of cells within their spatial context. Here, we have summarized available single-cell methods for dissociated tissues and tissue slides with respect to the specifics of microglial biology. We have focused on next-generation-based technologies. We review the potential of these single-cell sequencing methods and newer multiomics approaches to extend the understanding of microglia function beyond the status quo.