Unraveling Spatially Variable Genes: A Statistical Perspective on Spatial Transcriptomics


[

The article was written by Guanao Yan, Ph.D. student of Statistics and Data Science at UCLA. Guanao is the first author of the Nature Communications review article [1].

Spatially resolved transcriptomics (SRT) is revolutionizing Genomics by enabling the high-throughput measurement of gene expression whereas preserving spatial context. Not like single-cell RNA sequencing (scRNA-seq), which captures transcriptomes with out spatial location data, SRT permits researchers to map gene expression to specific places inside a tissue, offering insights into tissue group, mobile interactions, and spatially coordinated gene exercise. The growing quantity and complexity of SRT knowledge necessitate the event of sturdy statistical and computational strategies, making this discipline extremely related to knowledge scientists, statisticians, and machine studying (ML) professionals. Strategies corresponding to spatial statistics, graph-based fashions, and deep studying have been utilized to extract significant organic insights from these knowledge.

A key step in SRT evaluation is the detection of spatially variable genes (SVGs)—genes whose expression varies non-randomly throughout spatial places. Figuring out SVGs is essential for characterizing tissue structure, useful gene modules, and mobile heterogeneity. Nonetheless, regardless of the speedy growth of computational strategies for SVG detection, these strategies differ broadly of their definitions and statistical frameworks, resulting in inconsistent outcomes and challenges in interpretation.

In our latest assessment revealed in Nature Communications [1], we systematically examined 34 peer-reviewed SVG detection strategies and launched a classification framework that clarifies the organic significance of various SVG varieties. This text gives an outline of our findings, specializing in the three main classes of SVGs and the statistical rules underlying their detection.

SVG detection strategies purpose to uncover genes whose spatial expression displays organic patterns somewhat than technical noise. Based mostly on our assessment of 34 peer-reviewed strategies, we categorize SVGs into three teams: Total SVGs, Cell-Kind-Particular SVGs, and Spatial-Area-Marker SVGs (Determine 2).

Picture created by the authors, tailored from [1]. Publication timeline of 34 SVG detection strategies. Colours symbolize three SVG classes: total SVGs (inexperienced), cell-type-specific SVGs (pink), and spatial-domain-marker SVGs (purple).

Strategies for detecting the three SVG classes serve totally different functions (Fig. 3). First, the detection of total SVGs screens informative genes for downstream analyses, together with the identification of spatial domains and useful gene modules. Second, detecting cell-type-specific SVGs goals to disclose spatial variation inside a cell sort and assist determine distinct cell subpopulations or states inside cell varieties. Third, spatial-domain-marker SVG detection is used to seek out marker genes to annotate and interpret spatial domains already detected. These markers assist perceive the molecular mechanisms underlying spatial domains and help in annotating tissue layers in different datasets.

Picture created by the authors, tailored from [1]. Conceptual visualization of three SVG classes: total SVGs, cell-type-specific SVGs, and spatial-domain-marker SVGs. The left column reveals a tissue slice with two cell varieties and three spatial domains. The fitting column reveals exemplar genes with colours representing the expression ranges proven for an total SVG, a cell-type-specific SVG, and a spatial-domain-marker SVG, respectively.

The connection among the many three SVG classes will depend on the detection strategies, notably the null and different hypotheses they make use of. If an total SVG detection technique makes use of the null speculation {that a} non-SVG’s expression is impartial of spatial location and the choice speculation that any deviation from this independence signifies an SVG, then its SVGs ought to theoretically embody each cell-type-specific SVGs and spatial-domain-marker SVGs. For instance, DESpace [2] is a technique that detects each total SVGs and spatial-domain-marker SVGs, and its detected total SVGs have to be marker genes for some spatial domains. This inclusion relationship holds true besides in excessive eventualities, corresponding to when a gene displays reverse cell-type-specific spatial patterns that successfully cancel one another out. Nonetheless, if an total SVG detection technique’s different speculation is outlined for a particular spatial expression sample, then its SVGs might not embody some cell-type-specific SVGs or spatial-domain-marker SVGs.

To know how SVGs are detected, we categorized the statistical approaches into three main kinds of speculation exams: 

  1. Dependence Check – Examines the dependence between a gene’s expression stage and the spatial location. 
  2. Regression Fastened-Impact Check – Examines whether or not some or all the fixed-effect covariates, for example, spatial location, contribute to the imply of the response variable, i.e., a gene’s expression. 
  3. Regression Random-Impact Check (Variance Element Check) – Examines whether or not the random-effect covariates, for example, spatial location, contribute to the variance of the response variable, i.e., a gene’s expression.

To additional clarify how these exams are used for SVG detection, we denote Y as gene’s expression stage and S because the spatial places. Dependence check is probably the most common speculation check for SVG detection. For a given gene, it decides whether or not the gene’s expression stage Y is impartial of the spatial location S, i.e., the null speculation is:

There are two kinds of regression exams: fixed-effect exams, the place the impact of the spatial location is assumed to be mounted, and random-effect exams, which assume the impact of the spatial location as random. To clarify these two kinds of exams, we use a linear combined mannequin for a given gene for instance:

the place the response variable ( Y_i ) is the gene’s expression stage at spot ( i ), ( x_i ) ( epsilon ) ( R^p ) signifies the fixed-effect covariates of spot ( i ), ( z_i ) ( epsilon ) ( R^q ) denotes the random-effect covariates of spot ( i ), and ( epsilon_i ) is the random measurement error at spot ( i ) with zero imply. Within the mannequin parameters, ( beta_0 ) is the (mounted) intercept, ( beta ) ( epsilon ) ( R^p ) signifies the mounted results, and ( gamma ) ( epsilon ) ( R^q ) denotes the random results with zero means and the covariance matrix:

On this linear combined mannequin, independence is assumed between random impact and random errors and amongst random errors.

Fastened-effect exams look at whether or not some or all the fixed-effect covariates ( x_i ) (depending on spatial places S) contribute to the imply of the response variable. If all fixed-effect covariates make no contribution, then:

The null speculation

implies

Random-effect exams look at whether or not the random-effect covariates ( z_i ) (depending on spatial places S) contribute to the variance of the response variable Var⁡Yi, specializing in the decomposition:

and testing if the contribution of the random-effect covariates is zero. The null speculation:

implies

Among the many 23 strategies that use frequentist speculation exams, dependence exams and random-effect regression exams have been primarily utilized to detect total SVGs, whereas fixed-effect regression exams have been used throughout all three SVG classes. Understanding these distinctions is vital to deciding on the suitable technique for particular analysis questions.

Bettering SVG detection strategies requires balancing detection energy, specificity, and scalability whereas addressing key challenges in spatial transcriptomics evaluation. Future developments ought to concentrate on adapting strategies to totally different SRT applied sciences and tissue varieties, in addition to extending assist for multi-sample SRT knowledge to reinforce organic insights. Moreover, strengthening statistical rigor and validation frameworks shall be essential for making certain the reliability of SVG detection. Benchmarking research additionally want refinement, with clearer analysis metrics and standardized datasets to offer strong technique comparisons.

References

[1] Yan, G., Hua, S.H. & Li, J.J. (2025). Categorization of 34 computational strategies to detect spatially variable genes from spatially resolved transcriptomics knowledge. Nature Communication, 16, 1141. https://doi.org/10.1038/s41467-025-56080-w

[2] Cai, P., Robinson, M. D., & Tiberi, S. (2024). DESpace: spatially variable gene detection by way of differential expression testing of spatial clusters. Bioinformatics, 40(2). https://doi.org/10.1093/bioinformatics/btae027

]