SnapShot: Spatial transcriptomics.

Spatially resolved transcriptomics methodologies using RNA sequencing principles have and will continue to contribute to decode the molecular landscape of tissues. Linking quantitative sequencing data with tissue morphology empowers profiling of cellular morphology and transcription over time and space in health and disease. To view this SnapShot, open or download the PDF.

Spatial Transcriptomics and In Situ Sequencing to Study Alzheimer's Disease.

Although complex inflammatory-like alterations are observed around the amyloid plaques of Alzheimer's disease (AD), little is known about the molecular changes and cellular interactions that characterize this response. We investigate here, in an AD mouse model, the transcriptional changes occurring in tissue domains in a 100-µm diameter around amyloid plaques using spatial transcriptomics. We demonstrate early alterations in a gene co-expression network enriched for myelin and oligodendrocyte genes (OLIGs), whereas a multicellular gene co-expression network of plaque-induced genes (PIGs) involving the complement system, oxidative stress, lysosomes, and inflammation is prominent in the later phase of the disease. We confirm the majority of the observed alterations at the cellular level using in situ sequencing on mouse and human brain sections. Genome-wide spatial transcriptomics analysis provides an unprecedented approach to untangle the dysregulated cellular network in the vicinity of pathogenic hallmarks of AD and other brain diseases.

Multimodal Analysis of Composition and Spatial Architecture in Human Squamous Cell Carcinoma.

To define the cellular composition and architecture of cutaneous squamous cell carcinoma (cSCC), we combined single-cell RNA sequencing with spatial transcriptomics and multiplexed ion beam imaging from a series of human cSCCs and matched normal skin. cSCC exhibited four tumor subpopulations, three recapitulating normal epidermal states, and a tumor-specific keratinocyte (TSK) population unique to cancer, which localized to a fibrovascular niche. Integration of single-cell and spatial data mapped ligand-receptor networks to specific cell types, revealing TSK cells as a hub for intercellular communication. Multiple features of potential immunosuppression were observed, including T regulatory cell (Treg) co-localization with CD8 T cells in compartmentalized tumor stroma. Finally, single-cell characterization of human tumor xenografts and in vivo CRISPR screens identified essential roles for specific tumor subpopulation-enriched gene networks in tumorigenesis. These data define cSCC tumor and stromal cell subpopulations, the spatial niches where they interact, and the communicating gene networks that they engage in cancer.

A Spatiotemporal Organ-Wide Gene Expression and Cell Atlas of the Developing Human Heart.

The process of cardiac morphogenesis in humans is incompletely understood. Its full characterization requires a deep exploration of the organ-wide orchestration of gene expression with a single-cell spatial resolution. Here, we present a molecular approach that reveals the comprehensive transcriptional landscape of cell types populating the embryonic heart at three developmental stages and that maps cell-type-specific gene expression to specific anatomical domains. Spatial transcriptomics identified unique gene profiles that correspond to distinct anatomical regions in each developmental stage. Human embryonic cardiac cell types identified by single-cell RNA sequencing confirmed and enriched the spatial annotation of embryonic cardiac gene expression. In situ sequencing was then used to refine these results and create a spatial subcellular map for the three developmental phases. Finally, we generated a publicly available web resource of the human developing heart to facilitate future studies on human cardiogenesis.

Reducing Pericyte-Derived Scarring Promotes Recovery after Spinal Cord Injury.

CNS injury often severs axons. Scar tissue that forms locally at the lesion site is thought to block axonal regeneration, resulting in permanent functional deficits. We report that inhibiting the generation of progeny by a subclass of pericytes led to decreased fibrosis and extracellular matrix deposition after spinal cord injury in mice. Regeneration of raphespinal and corticospinal tract axons was enhanced and sensorimotor function recovery improved following spinal cord injury in animals with attenuated pericyte-derived scarring. Using optogenetic stimulation, we demonstrate that regenerated corticospinal tract axons integrated into the local spinal cord circuitry below the lesion site. The number of regenerated axons correlated with improved sensorimotor function recovery. In conclusion, attenuation of pericyte-derived fibrosis represents a promising therapeutic approach to facilitate recovery following CNS injury.

B cell expansion hinders the stroma-epithelium regenerative cross talk during mucosal healing.

Therapeutic promotion of intestinal regeneration holds great promise, but defining the cellular mechanisms that influence tissue regeneration remains an unmet challenge. To gain insight into the process of mucosal healing, we longitudinally examined the immune cell composition during intestinal damage and regeneration. B cells were the dominant cell type in the healing colon, and single-cell RNA sequencing (scRNA-seq) revealed expansion of an IFN-induced B cell subset during experimental mucosal healing that predominantly located in damaged areas and associated with colitis severity. B cell depletion accelerated recovery upon injury, decreased epithelial ulceration, and enhanced gene expression programs associated with tissue remodeling. scRNA-seq from the epithelial and stromal compartments combined with spatial transcriptomics and multiplex immunostaining showed that B cells decreased interactions between stromal and epithelial cells during mucosal healing. Activated B cells disrupted the epithelial-stromal cross talk required for organoid survival. Thus, B cell expansion during injury impairs epithelial-stromal cell interactions required for mucosal healing, with implications for the treatment of IBD.

Spatially resolved clonal copy number alterations in benign and malignant tissue.

Defining the transition from benign to malignant tissue is fundamental to improving early diagnosis of cancer1. Here we use a systematic approach to study spatial genome integrity in situ and describe previously unidentified clonal relationships. We used spatially resolved transcriptomics2 to infer spatial copy number variations in >120,000 regions across multiple organs, in benign and malignant tissues. We demonstrate that genome-wide copy number variation reveals distinct clonal patterns within tumours and in nearby benign tissue using an organ-wide approach focused on the prostate. Our results suggest a model for how genomic instability arises in histologically benign tissue that may represent early events in cancer evolution. We highlight the power of capturing the molecular and spatial continuums in a tissue context and challenge the rationale for treatment paradigms, including focal therapy.

Toward the single-hour high-quality genome.

Today, resequencing of a human genome can be performed in approximately a week using a single instrument. Thanks to a steady logarithmic rate of increase in performance for DNA sequencing platforms over the past seven years, DNA sequencing is one of the fastest developing technology fields. As the process becomes faster, it opens up possibilities within health care, diagnostics, and entirely new fields of research. Immediate genetic characterization of contagious outbreaks has been exemplified, and with such applications for the direct benefit of human health, expectations of future sensitive, rapid, high-throughput, and cost-effective technologies are steadily growing. Simultaneously, some of the limitations of a rapidly growing field have become apparent, and questions regarding the quality of some of the data deposited into databases have been raised. A human genome sequenced in only an hour is likely to become a reality in the future, but its definition may not be as certain.

Spatial landmark detection and tissue registration with deep learning.

Spatial landmarks are crucial in describing histological features between samples or sites, tracking regions of interest in microscopy, and registering tissue samples within a common coordinate framework. Although other studies have explored unsupervised landmark detection, existing methods are not well-suited for histological image data as they often require a large number of images to converge, are unable to handle nonlinear deformations between tissue sections and are ineffective for z-stack alignment, other modalities beyond image data or multimodal data. We address these challenges by introducing effortless landmark detection, a new unsupervised landmark detection and registration method using neural-network-guided thin-plate splines. Our proposed method is evaluated on a diverse range of datasets including histology and spatially resolved transcriptomics, demonstrating superior performance in both accuracy and stability compared to existing approaches.

A roadmap for the Human Developmental Cell Atlas.

The Human Developmental Cell Atlas (HDCA) initiative, which is part of the Human Cell Atlas, aims to create a comprehensive reference map of cells during development. This will be critical to understanding normal organogenesis, the effect of mutations, environmental factors and infectious agents on human development, congenital and childhood disorders, and the cellular basis of ageing, cancer and regenerative medicine. Here we outline the HDCA initiative and the challenges of mapping and modelling human development using state-of-the-art technologies to create a reference atlas across gestation. Similar to the Human Genome Project, the HDCA will integrate the output from a growing community of scientists who are mapping human development into a unified atlas. We describe the early milestones that have been achieved and the use of human stem-cell-derived cultures, organoids and animal models to inform the HDCA, especially for prenatal tissues that are hard to acquire. Finally, we provide a roadmap towards a complete atlas of human development.

Expansion spatial transcriptomics.

Capture array-based spatial transcriptomics methods have been widely used to resolve gene expression in tissues; however, their spatial resolution is limited by the density of the array. Here we present expansion spatial transcriptomics to overcome this limitation by clearing and expanding tissue prior to capturing the entire polyadenylated transcriptome with an enhanced protocol. This approach enables us to achieve higher spatial resolution while retaining high library quality, which we demonstrate using mouse brain samples.

A radical switch in clonality reveals a stem cell niche in the epiphyseal growth plate.

Longitudinal bone growth in children is sustained by growth plates, narrow discs of cartilage that provide a continuous supply of chondrocytes for endochondral ossification1. However, it remains unknown how this supply is maintained throughout childhood growth. Chondroprogenitors in the resting zone are thought to be gradually consumed as they supply cells for longitudinal growth1,2, but this model has never been proved. Here, using clonal genetic tracing with multicolour reporters and functional perturbations, we demonstrate that longitudinal growth during the fetal and neonatal periods involves depletion of chondroprogenitors, whereas later in life, coinciding with the formation of the secondary ossification centre, chondroprogenitors acquire the capacity for self-renewal, resulting in the formation of large, stable monoclonal columns of chondrocytes. Simultaneously, chondroprogenitors begin to express stem cell markers and undergo symmetric cell division. Regulation of the pool of self-renewing progenitors involves the hedgehog and mammalian target of rapamycin complex 1 (mTORC1) signalling pathways. Our findings indicate that a stem cell niche develops postnatally in the epiphyseal growth plate, which provides a continuous supply of chondrocytes over a prolonged period.

The spatial landscape of gene expression isoforms in tissue sections.

In situ capturing technologies add tissue context to gene expression data, with the potential of providing a greater understanding of complex biological systems. However, splicing variants and full-length sequence heterogeneity cannot be characterized at spatial resolution with current transcriptome profiling methods. To that end, we introduce spatial isoform transcriptomics (SiT), an explorative method for characterizing spatial isoform variation and sequence heterogeneity using long-read sequencing. We show in mouse brain how SiT can be used to profile isoform expression and sequence heterogeneity in different areas of the tissue. SiT reveals regional isoform switching of Plp1 gene between different layers of the olfactory bulb, and the use of external single-cell data allows the nomination of cell types expressing each isoform. Furthermore, SiT identifies differential isoform usage for several major genes implicated in brain function (Snap25, Bin1, Gnas) that are independently validated by in situ sequencing. SiT also provides for the first time an in-depth A-to-I RNA editing map of the adult mouse brain. Data exploration can be performed through an online resource (https://www.isomics.eu), where isoform expression and RNA editing can be visualized in a spatial context.

Semla: a versatile toolkit for spatially resolved transcriptomics analysis and visualization.

SUMMARY: Spatially resolved transcriptomics technologies generate gene expression data with retained positional information from a tissue section, often accompanied by a corresponding histological image. Computational tools should make it effortless to incorporate spatial information into data analyses and present analysis results in their histological context. Here, we present semla, an R package for processing, analysis, and visualization of spatially resolved transcriptomics data generated by the Visium platform, that includes interactive web applications for data exploration and tissue annotation. AVAILABILITY AND IMPLEMENTATION: The R package semla is available on GitHub (https://github.com/ludvigla/semla), under the MIT License, and deposited on Zenodo (https://doi.org/10.5281/zenodo.8321645). Documentation and tutorials with detailed descriptions of usage can be found at https://ludvigla.github.io/semla/.

Spatial transcriptomic analysis of virtual prostate biopsy reveals confounding effect of tissue heterogeneity on genomic signatures.

Genetic signatures have added a molecular dimension to prognostics and therapeutic decision-making. However, tumour heterogeneity in prostate cancer and current sampling methods could confound accurate assessment. Based on previously published spatial transcriptomic data from multifocal prostate cancer, we created virtual biopsy models that mimic conventional biopsy placement and core size. We then analysed the gene expression of different prognostic signatures (OncotypeDx®, Decipher®, Prostadiag®) using a step-wise approach with increasing resolution from pseudo-bulk analysis of the whole biopsy, to differentiation by tissue subtype (benign, stroma, tumour), followed by distinct tumour grade and finally clonal resolution. The gene expression profile of virtual tumour biopsies revealed clear differences between grade groups and tumour clones, compared to a benign control, which were not reflected in bulk analyses. This suggests that bulk analyses of whole biopsies or tumour-only areas, as used in clinical practice, may provide an inaccurate assessment of gene profiles. The type of tissue, the grade of the tumour and the clonal composition all influence the gene expression in a biopsy. Clinical decision making based on biopsy genomics should be made with caution while we await more precise targeting and cost-effective spatial analyses.

The discovAIR project: a roadmap towards the Human Lung Cell Atlas.

The Human Cell Atlas (HCA) consortium aims to establish an atlas of all organs in the healthy human body at single-cell resolution to increase our understanding of basic biological processes that govern development, physiology and anatomy, and to accelerate diagnosis and treatment of disease. The Lung Biological Network of the HCA aims to generate the Human Lung Cell Atlas as a reference for the cellular repertoire, molecular cell states and phenotypes, and cell-cell interactions that characterise normal lung homeostasis in healthy lung tissue. Such a reference atlas of the healthy human lung will facilitate mapping the changes in the cellular landscape in disease. The discovAIR project is one of six pilot actions for the HCA funded by the European Commission in the context of the H2020 framework programme. discovAIR aims to establish the first draft of an integrated Human Lung Cell Atlas, combining single-cell transcriptional and epigenetic profiling with spatially resolving techniques on matched tissue samples, as well as including a number of chronic and infectious diseases of the lung. The integrated Human Lung Cell Atlas will be available as a resource for the wider respiratory community, including basic and translational scientists, clinical medicine, and the private sector, as well as for patients with lung disease and the interested lay public. We anticipate that the Human Lung Cell Atlas will be the founding stone for a more detailed understanding of the pathogenesis of lung diseases, guiding the design of novel diagnostics and preventive or curative interventions.

High-definition spatial transcriptomics for in situ tissue profiling.

Spatial and molecular characteristics determine tissue function, yet high-resolution methods to capture both concurrently are lacking. Here, we developed high-definition spatial transcriptomics, which captures RNA from histological tissue sections on a dense, spatially barcoded bead array. Each experiment recovers several hundred thousand transcript-coupled spatial barcodes at 2-µm resolution, as demonstrated in mouse brain and primary breast cancer. This opens the way to high-resolution spatial analysis of cells and tissues.

sepal: identifying transcript profiles with spatial patterns by diffusion-based modeling.

MOTIVATION: Collection of spatial signals in large numbers has become a routine task in multiple omics-fields, but parsing of these rich datasets still pose certain challenges. In whole or near-full transcriptome spatial techniques, spurious expression profiles are intermixed with those exhibiting an organized structure. To distinguish profiles with spatial patterns from the background noise, a metric that enables quantification of spatial structure is desirable. Current methods designed for similar purposes tend to be built around a framework of statistical hypothesis testing, hence we were compelled to explore a fundamentally different strategy. RESULTS: We propose an unexplored approach to analyze spatial transcriptomics data, simulating diffusion of individual transcripts to extract genes with spatial patterns. The method performed as expected when presented with synthetic data. When applied to real data, it identified genes with distinct spatial profiles, involved in key biological processes or characteristic for certain cell types. Compared to existing methods, ours seemed to be less informed by the genes' expression levels and showed better time performance when run with multiple cores. AVAILABILITYAND IMPLEMENTATION: Open-source Python package with a command line interface (CLI), freely available at https://github.com/almaan/sepal under an MIT licence. A mirror of the GitHub repository can be found at Zenodo, doi: 10.5281/zenodo.4573237. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Novel loss-of-function variant in DENND5A impedes melanosomal cargo transport and predisposes to familial cutaneous melanoma.

PURPOSE: More than half of the familial cutaneous melanomas have unknown genetic predisposition. This study aims at characterizing a novel melanoma susceptibility gene. METHODS: We performed exome and targeted sequencing in melanoma-prone families without any known melanoma susceptibility genes. We analyzed the expression of candidate gene DENND5A in melanoma samples in relation to pigmentation and UV signature. Functional studies were carried out using microscopic approaches and zebrafish model. RESULTS: We identified a novel DENND5A truncating variant that segregated with melanoma in a Swedish family and 2 additional rare DENND5A variants, 1 of which segregated with the disease in an American family. We found that DENND5A is significantly enriched in pigmented melanoma tissue. Our functional studies show that loss of DENND5A function leads to decrease in melanin content in vitro and pigmentation defects in vivo. Mechanistically, harboring the truncating variant or being suppressed leads to DENND5A losing its interaction with SNX1 and its ability to transport the SNX1-associated vesicles from melanosomes. Consequently, untethered SNX1-premelanosome protein and redundant tyrosinase are redirected to lysosomal degradation by default, causing decrease in melanin content. CONCLUSION: Our findings provide evidence of a physiological role of DENND5A in the skin context and link its variants to melanoma susceptibility.