Principles of Spatial Transcriptomics Analysis: A Practical Walk-Through in Kidney Tissue

Spatial transcriptomic technologies capture genome-wide readouts across biological tissue space. Moreover, recent advances in this technology, including Slide-seqV2, have achieved spatial transcriptomic data collection at a near-single cell resolution. To-date, a repertoire of computational tools has been developed to discern cell type classes given the transcriptomic profiles of tissue coordinates. Upon applying these tools, we can explore the spatial patterns of distinct cell types and characterize how genes are spatially expressed within different cell type contexts. The kidney is one organ whose function relies upon spatially defined structures consisting of distinct cellular makeup. Thus, the application of Slide-seqV2 to kidney tissue has enabled us to elucidate spatially characteristic cellular and genetic profiles at a scale that remains largely unexplored. Here, we review spatial transcriptomic technologies, as well as computational approaches for cell type mapping and spatial cell type and transcriptomic characterizations. We take kidney tissue as an example to demonstrate how the technologies are applied, while considering the nuances of this architecturally complex tissue.

A Model of Spatial Cell Development in Rat Hippocampus Based on Artificial Neural Network

Physiological studies have shown that the hippocampal structure of rats develops at different stages, in which the place cells continue to develop during the whole juvenile period of rats and mature after the juvenile period. As the main information source of place cells, grid cells should mature earlier than place cells. In order to make better use of the biological information exhibited by the rat brain hippocampus in the environment, we propose a position cognition model based on the spatial cell development mechanism of rat hippocampus. The model uses a recurrent neural network with parametric bias (RNNPB) to simulate changes in the discharge characteristics during the development of a single stripe cell. The oscillatory interference mechanism is able to fuse the developing stripe waves, thus indirectly simulating the developmental process of the grid cells. The output of the grid cells is then used as the information input of the place cells, whose development process is simulated by BP neural network. After the place cells matured, the position matrix generated by the place cell group was used to realize the position cognition of rats in a given spatial region. The experimental results show that this model can simulate the development process of grid cells and place cells, and it can realize high precision positioning in the given space area. Moreover, the experimental effect of cognitive map construction using this model is basically consistent with the effect of RatSLAM, which verifies the validity and accuracy of the model.

A Novel Averaging Principle Provides Insights in the Impact of Intratumoral Heterogeneity on Tumor Progression

Typically stochastic differential equations (SDEs) involve an additive or multiplicative noise term. Here, we are interested in stochastic differential equations for which the white noise is nonlinearly integrated into the corresponding evolution term, typically termed as random ordinary differential equations (RODEs). The classical averaging methods fail to treat such RODEs. Therefore, we introduce a novel averaging method appropriate to be applied to a specific class of RODEs. To exemplify the importance of our method, we apply it to an important biomedical problem, in particular, we implement the method to the assessment of intratumoral heterogeneity impact on tumor dynamics. Precisely, we model gliomas according to a well-known Go or Grow (GoG) model, and tumor heterogeneity is modeled as a stochastic process. It has been shown that the corresponding deterministic GoG model exhibits an emerging Allee effect (bistability). In contrast, we analytically and computationally show that the introduction of white noise, as a model of intratumoral heterogeneity, leads to monostable tumor growth. This monostability behavior is also derived even when spatial cell diffusion is taken into account.

The Hippocampus uses cryptographic principles rather than plasticity to enable animals to secretly cache and retrieve their food

For animals, the ability to hide and retrieve valuable information, such as the location of food, can mean the difference between life and death. Here, we propose that to achieve this, their brain uses spatial cells similarly to how we utilize encryption for data security. Some animals are able to cache hundreds of thousands of food items annually by each individual and later retrieve most of what they themselves stashed. Rather than memorizing their cache locations as previously suggested, we propose that they use a single cryptographic-like mechanism during both caching and retrieval. The model we developed is based on hippocampal spatial cells, which respond to an animal's positional attention, such as when the animal enters a specific region (place-cells) or gazes at a particular location (spatial-view-cells). We know that the region that activates each spatial cell remains consistent across subsequent visits to the same area but not between areas. This remapping, combined with the uniqueness of cognitive maps, produces a persistent crypto-hash function for both food caching and retrieval. We also show that the model stores temporal information that helps animals in food caching order preference as previously observed. This behavior, which we refer to as crypto-taxis, might also explain consistent differences in decision-making when animals are faced with a large number of alternatives such as in foraging.

How environmental movement constraints shape the neural code for space

Study of the neural code for space in rodents has many insights to offer for how mammals, including humans, construct a mental representation of space. This code is centered on the hippocampal place cells, which are active in particular places in the environment. Place cells are informed by numerous other spatial cell types including grid cells, which provide a signal for distance and direction and are thought to help anchor the place cell signal. These neurons combine self-motion and environmental information to create and update their map-like representation. Study of their activity patterns in complex environments of varying structure has revealed that this "cognitive map" of space is not a fixed and rigid entity that permeates space, but rather is variably affected by the movement constraints of the environment. These findings are pointing toward a more flexible spatial code in which the map is adapted to the movement possibilities of the space. An as-yet-unanswered question is whether these different forms of representation have functional consequences, as suggested by an enactivist view of spatial cognition.

Cell segmentation-free inference of cell types from in situ transcriptomics data

Multiplexed fluorescence in situ hybridization techniques have enabled cell-type identification, linking transcriptional heterogeneity with spatial heterogeneity of cells. However, inaccurate cell segmentation reduces the efficacy of cell-type identification and tissue characterization. Here, we present a method called Spot-based Spatial cell-type Analysis by Multidimensional mRNA density estimation (SSAM), a robust cell segmentation-free computational framework for identifying cell-types and tissue domains in 2D and 3D. SSAM is applicable to a variety of in situ transcriptomics techniques and capable of integrating prior knowledge of cell types. We apply SSAM to three mouse brain tissue images: the somatosensory cortex imaged by osmFISH, the hypothalamic preoptic region by MERFISH, and the visual cortex by multiplexed smFISH. Here, we show that SSAM detects regions occupied by known cell types that were previously missed and discovers new cell types.

Spatial cell type mapping of the oligodendrocyte lineage in the mouse juvenile and adult CNS with in situ sequencing

Oligodendrocytes show transcriptional heterogeneity but the regional and functional implications of this heterogeneity are less clear. Here, we apply in situ sequencing (ISS) to simultaneously probe the expression of 124 marker genes of distinct oligodendrocyte populations, providing comprehensive maps of corpus callosum, cingulate, motor and somatosensory cortex in the brain, as well as gray (GM) and white matter (WM) regions in the spinal cord, at juvenile and adult stages. We systematically compare abundances of these populations and investigate the neighboring preference of distinct oligodendrocyte populations. As previously described, we observed that oligodendrocyte lineage progression is more advanced in the juvenile spinal cord compared to the brain. Additionally, myelination is ongoing in the adult corpus callosum while it is mostly completed in the cortex. Interestingly, we found a medial-to-lateral gradient of oligodendrocyte lineage progression in the juvenile cortex, which could be linked to arealization, as well as a deep-to-superficial gradient with mature oligodendrocytes preferentially accumulating in the deeper layers of the cortex. We observed differences in abundances and population dynamics over time between GM and WM regions in the brain and spinal cord, indicating regional differences within GM and WM. We also found that oligodendroglia populations' neighboring preferences are altered from the juvenile to the adult CNS. Thus, our ISS dataset reveals spatial heterogeneity of the oligodendrocyte lineage progression in the brain and spinal cord, which could be relevant to further investigate functional heterogeneity of oligodendroglia.

stPlus: a reference-based method for the accurate enhancement of spatial transcriptomics

Motivation: Single-cell RNA sequencing (scRNA-seq) techniques have revolutionized the investigation of transcriptomic landscape in individual cells. Recent advancements in spatial transcriptomic technologies further enable gene expression profiling and spatial organization mapping of cells simultaneously. Among the technologies, imaging-based methods can offer higher spatial resolutions, while they are limited by either the small number of genes imaged or the low gene detection sensitivity. Although several methods have been proposed for enhancing spatially resolved transcriptomics, inadequate accuracy of gene expression prediction and insufficient ability of cell-population identification still impede the applications of these methods. Results: We propose stPlus, a reference-based method that leverages information in scRNA-seq data to enhance spatial transcriptomics. Based on an auto-encoder with a carefully tailored loss function, stPlus performs joint embedding and predicts spatial gene expression via a weighted k-NN. stPlus outperforms baseline methods with higher gene-wise and cell-wise Spearman correlation coefficients. We also introduce a clustering-based approach to assess the enhancement performance systematically. Using the data enhanced by stPlus, cell populations can be better identified than using the measured data. The predicted expression of genes unique to scRNA-seq data can also well characterize spatial cell heterogeneity. Besides, stPlus is robust and scalable to datasets of diverse gene detection sensitivity levels, sample sizes, and number of spatially measured genes. We anticipate stPlus will facilitate the analysis of spatial transcriptomics. Availability: stPlus with detailed documents is freely accessible at http://health.tsinghua.edu.cn/software/stPlus/ and the source code is openly available on https://github.com/xy-chen16/stPlus.

Optimization and parallelization of the discrete ordinate method for radiation transport simulation in OpenFOAM: Hierarchical combination of shared and distributed memory approaches

This paper describes the reduction in memory and computational time for the simulation of complex radiation transport problems with the discrete ordinate method (DOM) model in the open-source computational fluid dynamics platform OpenFOAM. Finite volume models require storage of vector variables in each spatial cell; DOM introduces two additional discretizations, in direction and wavelength, making memory a limiting factor. Using specific classes for radiation sources data, changing the store of fluxes and other minor changes allowed a reduction of 75% in memory requirements. Besides, a hierarchical parallelization was developed, where each node of the standard parallelization uses several computing threads, allowing higher speed and scalability of the problem. This architecture, combined with optimization of some parts of the code, allowed a global speedup of x15. This relevant reduction in time and memory of radiation transport opens a new horizon of applications previously unaffordable.