Abstract 4708: Application of spatially resolved transcriptomics to screen multiple tumor biospecimens using tissue microarrays

The tumor microenvironment is composed of highly heterogeneous structures and cell types that dynamically influence and communicate with each other. Although examination of singular biospecimens is sufficient for diagnostic purposes, it is inadequate and cost prohibitive when scaling for complex and overarching studies. Thus, high density multi-tumor tissue microarrays (TMAs) have been a practical and effective solution for high-throughput molecular analysis of tissues. Introduced more than a decade ago, TMAs have been instrumental in the recent study of tumor biology, the development of diagnostic tests, the establishment of quality control, and the investigation and identification of oncological biomarkers. Here, we demonstrate the pairing of the 10x Genomics Visium Cytassist Spatial Gene Expression Solution and Xenium In-Situ Platform on multi-tumor TMAs to screen for common biomarkers among a cohort of samples. Spatial transcriptomics technology has proven valuable in mapping the whole transcriptome with spatial context (Visium), whereas In Situ (Xenium) enables high-throughput cellular characterization at single-cell resolution. With the addition of our CytAssist platform, we expand on the pre-existing standard Visium solution by facilitating the retrieval of RNA transcriptomic information from tissues placed on standard or archival slides. On the other hand, the novel Xenium platform compliments whole transcriptome Visium data by unlocking the potential to assign transcripts to a particular cell with spatial context and subcellular resolution. The combination of spatial transcriptomics and targeted in situ data with FFPE TMAs promotes a high-throughput method to accelerate the uncovering of molecular signatures suitable to understanding the tumor microenvironment. We showcase the ability to spatially and comprehensively resolve individual oncogene and tumor suppressor genes associated with multiple tumors from a cohort of cancer patients and from multiple different tumor samples. In addition, these markers are mapped back to distinct morphological features within each tissue core, and use differential gene expression data to identify distinct cell types throughout the different patient tissues. By combining the throughput of TMA samples and depth of the Visium and Xenium platforms, the strategy enables greater insights into cell-type specifics while also expanding the spectrum of biospecimen types that can be analyzed.

Citation Format: Syrus Mohabbat, Hardeep Singh, Stephen R. Williams, Lauren M M. Gutgesell, David J. Sukovich, Govinda M. Kamath, Hanyoup Kim, Amanda Janesick, Robert Shelansky, Ghezal Beliakoff, Augusto M. Tentori, Albert Kim, Cedric R. Uytingco, Sarah Taylor. Application of spatially resolved transcriptomics to screen multiple tumor biospecimens using tissue microarrays. [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 4708.

Genome-scale screens identify factors regulating tumor cell responses to natural killer cells

To systematically define molecular features in human tumor cells that determine their degree of sensitivity to human allogeneic natural killer (NK) cells, we quantified the NK cell responsiveness of hundreds of molecularly annotated 'DNA-barcoded' solid tumor cell lines in multiplexed format and applied genome-scale CRISPR-based gene-editing screens in several solid tumor cell lines, to functionally interrogate which genes in tumor cells regulate the response to NK cells. In these orthogonal studies, NK cell-sensitive tumor cells tend to exhibit 'mesenchymal-like' transcriptional programs; high transcriptional signature for chromatin remodeling complexes; high levels of B7-H6 (NCR3LG1); and low levels of HLA-E/antigen presentation genes. Importantly, transcriptional signatures of NK cell-sensitive tumor cells correlate with immune checkpoint inhibitor (ICI) resistance in clinical samples. This study provides a comprehensive map of mechanisms regulating tumor cell responses to NK cells, with implications for future biomarker-driven applications of NK cell immunotherapies.

Spectral Jaccard Similarity: A New Approach to Estimating Pairwise Sequence Alignments

Pairwise sequence alignment is often a computational bottleneck in genomic analysis pipelines, particularly in the context of third-generation sequencing technologies. To speed up this process, the pairwise k-mer Jaccard similarity is sometimes used as a proxy for alignment size in order to filter pairs of reads, and min-hashes are employed to efficiently estimate these similarities. However, when the k-mer distribution of a dataset is significantly non-uniform (e.g., due to GC biases and repeats), Jaccard similarity is no longer a good proxy for alignment size. In this work, we introduce a min-hash-based approach for estimating alignment sizes called Spectral Jaccard Similarity, which naturally accounts for uneven k-mer distributions. The Spectral Jaccard Similarity is computed by performing a singular value decomposition on a min-hash collision matrix. We empirically show that this new metric provides significantly better estimates for alignment sizes, and we provide a computationally efficient estimator for these spectral similarity scores.

Valid Post-clustering Differential Analysis for Single-Cell RNA-Seq

Single-cell computational pipelines involve two critical steps: organizing cells (clustering) and identifying the markers driving this organization (differential expression analysis). State-of-the-art pipelines perform differential analysis after clustering on the same dataset. We observe that because clustering "forces" separation, reusing the same dataset generates artificially low p values and hence false discoveries. We introduce a valid post-clustering differential analysis framework, which corrects for this problem. We provide software at https://github.com/jessemzhang/tn_test.

HINGE: long-read assembly achieves optimal repeat resolution

Long-read sequencing technologies have the potential to produce gold-standard de novo genome assemblies, but fully exploiting error-prone reads to resolve repeats remains a challenge. Aggressive approaches to repeat resolution often produce misassemblies, and conservative approaches lead to unnecessary fragmentation. We present HINGE, an assembler that seeks to achieve optimal repeat resolution by distinguishing repeats that can be resolved given the data from those that cannot. This is accomplished by adding “hinges” to reads for constructing an overlap graph where only unresolvable repeats are merged. As a result, HINGE combines the error resilience of overlap-based assemblers with repeat-resolution capabilities of de Bruijn graph assemblers. HINGE was evaluated on the long-read bacterial data sets from the NCTC project. HINGE produces more finished assemblies than Miniasm and the manual pipeline of NCTC based on the HGAP assembler and Circlator. HINGE also allows us to identify 40 data sets where unresolvable repeats prevent the reliable construction of a unique finished assembly. In these cases, HINGE outputs a visually interpretable assembly graph that encodes all possible finished assemblies consistent with the reads, while other approaches such as the NCTC pipeline and FALCON either fragment the assembly or resolve the ambiguity arbitrarily.

Fast and accurate single-cell RNA-seq analysis by clustering of transcript-compatibility counts

Current approaches to single-cell transcriptomic analysis are computationally intensive and require assay-specific modeling, which limits their scope and generality. We propose a novel method that compares and clusters cells based on their transcript-compatibility read counts rather than on the transcript or gene quantifications used in standard analysis pipelines. In the reanalysis of two landmark yet disparate single-cell RNA-seq datasets, we show that our method is up to two orders of magnitude faster than previous approaches, provides accurate and in some cases improved results, and is directly applicable to data from a wide variety of assays.