Visualizing somatic alterations in spatial transcriptomics data of skin cancer

Tools to visualize genetic alterations within tissues remain underdeveloped despite the growth of spatial transcriptomic technologies, which measure gene expression in different regions of tissues. Since genetic alterations can be detected in RNA-sequencing data, we explored the feasibility of observing somatic alterations in spatial transcriptomics data. Extracting genetic information from spatial transcriptomic data would illuminate the spatial distribution of clones and allow for correlations with regional changes in gene expression to support genotype-phenotype studies. Recent work demonstrates that copy number alterations can be inferred from spatial transcriptomics data1. Here, we describe new software to further enhance the inference of copy number from spatial transcriptomics data. Moreover, we demonstrate that single nucleotide variants are also detectable in spatial transcriptomic data. We applied these approaches to map the location of point mutations, copy number alterations, and allelic imbalances in spatial transcriptomic data of two cutaneous squamous cell carcinomas. We show that both tumors are dominated by a single clone of cells, suggesting that their regional variations in gene expression2 are likely driven by non-genetic factors. Furthermore, we observe mutant cells in histologically normal tissue surrounding one tumor, which were not discernible upon histopathologic evaluation. Finally, we detected mono-allelic expression of immunoglobulin heavy chains in B-cells, revealing clonal populations of plasma cells surrounding one tumor. In summary, we put forward solutions to add the genetic dimension to spatial transcriptomic datasets, augmenting the potential of this new technology.