geneBasis: an iterative approach for unsupervised selection of targeted gene panels from scRNA-seq

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.

Gene panel selection for targeted spatial transcriptomics

Targeted spatial transcriptomics hold particular promise in analyzing complex tissues. Most such methods, however, measure only a limited panel of transcripts, which need to be selected in advance to inform on the cell types or processes being studied. A limitation of existing gene selection methods is their reliance on scRNA-seq data, ignoring platform effects between technologies. Here we describe gpsFISH, a computational method performing gene selection through optimizing detection of known cell types. By modeling and adjusting for platform effects, gpsFISH outperforms other methods. Furthermore, gpsFISH can incorporate cell type hierarchies and custom gene preferences to accommodate diverse design requirements.

Multimodal spatiotemporal phenotyping of human retinal organoid development

Organoids generated from human pluripotent stem cells provide experimental systems to study development and disease, but quantitative measurements across different spatial scales and molecular modalities are lacking. In this study, we generated multiplexed protein maps over a retinal organoid time course and primary adult human retinal tissue. We developed a toolkit to visualize progenitor and neuron location, the spatial arrangements of extracellular and subcellular components and global patterning in each organoid and primary tissue. In addition, we generated a single-cell transcriptome and chromatin accessibility timecourse dataset and inferred a gene regulatory network underlying organoid development. We integrated genomic data with spatially segmented nuclei into a multimodal atlas to explore organoid patterning and retinal ganglion cell (RGC) spatial neighborhoods, highlighting pathways involved in RGC cell death and showing that mosaic genetic perturbations in retinal organoids provide insight into cell fate regulation.

FEED: a feature selection method based on gene expression decomposition for single cell clustering

Single-cell clustering is a critical step in biological downstream analysis. The clustering performance could be effectively improved by extracting cell-type-specific genes. The state-of-the-art feature selection methods usually calculate the importance of a single gene without considering the information contained in the gene expression distribution. Moreover, these methods ignore the intrinsic expression patterns of genes and heterogeneity within groups of different mean expression levels. In this work, we present a Feature sElection method based on gene Expression Decomposition (FEED) of scRNA-seq data, which selects informative genes to enhance clustering performance. First, the expression levels of genes are decomposed into multiple Gaussian components. Then, a novel gene correlation calculation method is proposed to measure the relationship between genes from the perspective of distribution. Finally, a permutation-based approach is proposed to determine the threshold of gene importance to obtain marker gene subsets. Compared with state-of-the-art feature selection methods, applying FEED on various scRNA-seq datasets including large datasets followed by different common clustering algorithms results in significant improvements in the accuracy of cell-type identification. The source codes for FEED are freely available at https://github.com/genemine/FEED.

Single-cell RNA reveals a tumorigenic microenvironment in the interface zone of human breast tumors

BACKGROUND

The interface zone, area around invasive carcinoma, can be thought of as the actual tissue of the tumor microenvironment with precedent alterations for tumor invasion. However, the heterogeneity and characteristics of the microenvironment in the interface area have not yet been thoroughly explored.

METHODS

For in vitro studies, single-cell RNA sequencing (scRNA-seq) was used to characterize the cells from the tumor zone, the normal zone and the interface zone with 5-mm-wide belts between the tumor invasion front and the normal zone. Through scRNA-seq data analysis, we compared the cell types and their transcriptional characteristics in the different zones. Pseudotime, cell-cell communication and pathway analysis were performed to characterize the zone-specific microenvironment. Cell proliferation, wound healing and clone formation experiments explored the function of differentially expressed gene BMPR1B, which were confirmed by tumor models in vivo.

RESULTS

After screening, 88,548 high-quality cells were obtained and identified. Regulatory T cells, M2 macrophages, angiogenesis-related mast cells, stem cells with weak DNA repair ability, endothelial cells with angiogenic activity, fibroblasts with collagen synthesis and epithelial cells with proliferative activity form a unique tumorigenic microenvironment in the interface zone. Cell-cell communication analysis revealed that there are special ligand-receptor pairs between different cell types in the interface zone, which protects endothelial cell apoptosis and promotes epithelial cell proliferation and migration, compared to the normal zone. Compared with the normal zone, the highly expressed BMPR1B gene promotes the tumorigenic ability of cancer cells in the interface zone.

CONCLUSIONS

Our work identified a unique tumorigenic microenvironment of the interface zone and allowed for deeper insights into the tumor microenvironment of breast cancer that will serve as a helpful resource for advancing breast cancer diagnosis and therapy.

The dawn of spatial omics

Spatial omics has been widely heralded as the new frontier in life sciences. This term encompasses a wide range of techniques that promise to transform many areas of biology and eventually revolutionize pathology by measuring physical tissue structure and molecular characteristics at the same time. Although the field came of age in the past 5 years, it still suffers from some growing pains: barriers to entry, robustness, unclear best practices for experimental design and analysis, and lack of standardization. In this Review, we present a systematic catalog of the different families of spatial omics technologies; highlight their principles, power, and limitations; and give some perspective and suggestions on the biggest challenges that lay ahead in this incredibly powerful-but still hard to navigate-landscape.

Predictive and robust gene selection for spatial transcriptomics

A prominent trend in single-cell transcriptomics is providing spatial context alongside a characterization of each cell's molecular state. This typically requires targeting an a priori selection of genes, often covering less than 1% of the genome, and a key question is how to optimally determine the small gene panel. We address this challenge by introducing a flexible deep learning framework, PERSIST, to identify informative gene targets for spatial transcriptomics studies by leveraging reference scRNA-seq data. Using datasets spanning different brain regions, species, and scRNA-seq technologies, we show that PERSIST reliably identifies panels that provide more accurate prediction of the genome-wide expression profile, thereby capturing more information with fewer genes. PERSIST can be adapted to specific biological goals, and we demonstrate that PERSIST's binarization of gene expression levels enables models trained on scRNA-seq data to generalize with to spatial transcriptomics data, despite the complex shift between these technologies.

Gene panel selection for targeted spatial transcriptomics

Targeted spatial transcriptomics hold particular promise in analysis of complex tissues. Most such methods, however, measure only a limited panel of transcripts, which need to be selected in advance to inform on the cell types or processes being studied. A limitation of existing gene selection methods is that they rely on scRNA-seq data, ignoring platform effects between technologies. Here we describe gpsFISH, a computational method to perform gene selection through optimizing detection of known cell types. By modeling and adjusting for platform effects, gpsFISH outperforms other methods. Furthermore, gpsFISH can incorporate cell type hierarchies and custom gene preferences to accommodate diverse design requirements.

Spatial RNA sequencing methods show high resolution of single cell in cancer metastasis and the formation of tumor microenvironment

Cancer metastasis often leads to death and therapeutic resistance. This process involves the participation of a variety of cell components, especially cellular and intercellular communications in the tumor microenvironment (TME). Using genetic sequencing technology to comprehensively characterize the tumor and TME is therefore key to understanding metastasis and therapeutic resistance. The use of spatial transcriptome sequencing enables the localization of gene expressions and cell activities in tissue sections. By examining the localization change as well as gene expression of these cells, it is possible to characterize the progress of tumor metastasis and TME formation. With improvements of this technology, spatial transcriptome sequencing technology has been extended from local regions to whole tissues, and from single sequencing technology to multimodal analysis combined with a variety of datasets. This has enabled the detection of every single cell in tissue slides, with high resolution, to provide more accurate predictive information for tumor treatments. In this review, we summarize the results of recent studies dealing with new multimodal methods and spatial transcriptome sequencing methods in tumors to illustrate recent developments in the imaging resolution of micro-tissues.

Concordance of MERFISH spatial transcriptomics with bulk and single-cell RNA sequencing

Spatial transcriptomics extends single-cell RNA sequencing (scRNA-seq) by providing spatial context for cell type identification and analysis. Imaging-based spatial technologies such as multiplexed error-robust fluorescence in situ hybridization (MERFISH) can achieve single-cell resolution, directly mapping single-cell identities to spatial positions. MERFISH produces a different data type than scRNA-seq, and a technical comparison between the two modalities is necessary to ascertain how to best integrate them. We performed MERFISH on the mouse liver and kidney and compared the resulting bulk and single-cell RNA statistics with those from the Tabula Muris Senis cell atlas and from two Visium datasets. MERFISH quantitatively reproduced the bulk RNA-seq and scRNA-seq results with improvements in overall dropout rates and sensitivity. Finally, we found that MERFISH independently resolved distinct cell types and spatial structure in both the liver and kidney. Computational integration with the Tabula Muris Senis atlas did not enhance these results. We conclude that MERFISH provides a quantitatively comparable method for single-cell gene expression and can identify cell types without the need for computational integration with scRNA-seq atlases.