Unravelling the progression of the zebrafish primary body axis with reconstructed spatiotemporal transcriptomics

Yang Dong,Tao Cheng,Xiang Liu,Xin-xin Fu,Yang Hu,Xian-Fa Yang,Lingen Yang,Hao-Ran Li,Zhi-Wen Bian,NAIHE JING,Jie Liao,Xiaohui Fan,Peng-Fei Xu
DOI: https://doi.org/10.1101/2024.07.01.601472
2024-07-03
Abstract:Elucidating the spatiotemporal dynamics of gene expression is essential for understanding complex physiological and pathological processes. Traditional technologies like in situ hybridization (ISH) and immunostaining have been restricted to analyzing expression patterns of a limited number of genes. Spatial transcriptomics (ST) has emerged as a robust alternative, enabling the investigation of spatial patterns of thousands of genes simultaneously. However, current ST methods are hindered by low read depths and limited gene detection capabilities. Here, we introduce Palette, a pipeline that infers detailed spatial gene expression patterns from bulk RNA-seq data, utilizing existing ST data as only reference. This method identifies more precise expression patterns by smoothing, imputing and adjusting gene expressions. We applied Palette to construct the Danio rerio SpatioTemporal Expression Profiles (DreSTEP) by integrating 53-slice serial bulk RNA-seq data from three developmental stages with existing ST references and 3D zebrafish embryo images. DreSTEP provides a comprehensive cartographic resource for examining gene expression and spatial cell-cell interactions within zebrafish embryos. Utilizing machine learning-based screening, we identified key morphogens and transcription factors (TFs) essential for anteroposterior (AP) axis development and characterized their dynamic distribution throughout embryogenesis. In addition, among these TFs, Hox family genes were found to be pivotal in AP axis refinement. Their expression was closely correlated with cellular AP identities, and hoxb genes may act as central regulators in this process.
Developmental Biology
What problem does this paper attempt to address?
This paper aims to solve the problem of how to infer detailed gene spatial expression patterns from bulk RNA - sequencing data. Specifically, the researchers developed a pipeline named Palette, which uses existing spatial transcriptomics (ST) data as a reference to infer spatial gene expression patterns in bulk RNA - sequencing data. This method focuses particularly on improving spatial specificity and continuity, as well as the ability to effectively supplement in the case of a low number of detected genes. Through this method, the researchers constructed a three - dimensional spatio - temporal expression profile map (DreSTEP) of zebrafish embryos, which can more accurately reflect gene expression patterns and their spatial cell - cell interactions during embryonic development, especially during the development of the anterior - posterior axis (AP axis). In addition, the study also identified morphogens and transcription factors that are crucial for AP - axis refinement, especially the core regulatory role of Hox family genes in this process.