Skip to main page content
U.S. flag

An official website of the United States government

Dot gov

The .gov means it’s official.
Federal government websites often end in .gov or .mil. Before sharing sensitive information, make sure you’re on a federal government site.

Https

The site is secure.
The https:// ensures that you are connecting to the official website and that any information you provide is encrypted and transmitted securely.

Access keys NCBI Homepage MyNCBI Homepage Main Content Main Navigation
. 2025 Jan 8;20(1):e0311296.
doi: 10.1371/journal.pone.0311296. eCollection 2025.

tomoseqr: A Bioconductor package for spatial reconstruction and visualization of 3D gene expression patterns based on RNA tomography

Ryosuke Matsuzawa  1   2 , Daichi Kawahara  3 , Makoto Kashima  3   4 , Hiromi Hirata  3 , Haruka Ozaki  2
Affiliations

tomoseqr: A Bioconductor package for spatial reconstruction and visualization of 3D gene expression patterns based on RNA tomography

Ryosuke Matsuzawa et al. PLoS One. .

Abstract

RNA tomography computationally reconstructs 3D spatial gene expression patterns genome-widely from 1D tomo-seq data, generated by RNA sequencing of cryosection samples along three orthogonal axes. We developed tomoseqr, an R package designed for RNA tomography analysis of tomo-seq data, to reconstruct and visualize 3D gene expression patterns through user-friendly graphical interfaces. We show the effectiveness of tomoseqr using simulated and real tomo-seq data, validating its utility for researchers. R package tomoseqr is available on Bioconductor (https://doi.org/doi:10.18129/B9.bioc.tomoseqr) and GitHub (https://github.com/bioinfo-tsukuba/tomoseqr).

PubMed Disclaimer

Conflict of interest statement

The authors have declared that no competing interests exist.

Figures

Fig 1
Fig 1. Overview of tomoseqr.
Based on 1D tomo-seq data along three mutually orthogonal axes and a mask data, tomoseqr performs RNA tomography using iterative proportional fitting and reconstructs a 3D expression pattern of each gene. A GUI called masker helps users to design and edit a mask data. Another GUI Image viewer visualizes the reconstructed gene expression patterns in 2D or 3D view.
Fig 2
Fig 2. User interface of masker.
Fig 3
Fig 3. User interface of Image viewer (2D view).
Fig 4
Fig 4. User interface of Image viewer (3D view).
Fig 5
Fig 5. Spatial expression patterns and reconstruction results for Gene1, Gene2, Gene3, Gene4, and Gene5.
Ground truth and reconstruction results are shown for (A) Gene1, (B) Gene2, (C) Gene3, (D) Gene4, and (E) Gene5. The dots are colored by gene expression levels, with white dots representing the sample shape (mask).
Fig 6
Fig 6. Results of reconstruction with zebrafish (Danio rerio, shield stage).
Reconstruction result for expression of (A) eve1 (ENSDARG00000056012), (B) bmp7a (ENSDARG00000018260), (C) ntla (ENSDARG00000009905), and (D) chd (ENSDARG00000006110). For each of the genes, 3D view, cross section perpendicular to the animal-vegetal axis (18 m × 50 sections), cross section perpendicular to the dorsal-ventral axis (18 m × 49 sections), and cross section perpendicular to the left-right axis (18 m × 56 sections) are shown from left to right. Parentheses after gene names indicate gene IDs. Colors represent the reconstructed expression levels.
Fig 7
Fig 7. Reconstruction results using planarian tomo-seq data.
Reconstruction result for expression of (A) piwiA (DjGI005146_001), (B) opsin (DjGI008464_001), (C) Djf-1 (DjGI006376_001), and (D) DjNp19 (DjGI017609_001). Parentheses after gene names indicate gene IDs. Colors represent the reconstructed expression levels.AP means anterior-posterior axis (80 μm × 78 sections), VD means ventral-dorsal axis,(20 μm × 28 sections) and LR means left-right axis (20 μm × 51 sections).
Fig 8
Fig 8. Systematic exploration of genes correlated with the reconstructed opsin expression pattern in planarian tomo-seq data.
(A) A histogram showing the Pearson correlation coefficients (PCCs) with the reconstructed opsin expression pattern for 18,768 expressed genes. The vertical lines indicate the top 4 genes with the highest PCCs. (B-E) Reconstruction results for the expression of (B) Arrb1 (DjGI008623_001), (C) Rnf13 (DjGI009330_001), (D) Cpne9 (DjGI006974_001), and (E) Tnnc1 (DjGI015938_001). The PCCs of these gene expression patterns with the opsin gene (Fig 7B) are shown. (F) A heatmap showing the expression of human homologs in human eye single-cell types using single-cell RNA sequencing data from the Human Protein Atlas. The gene expression levels, in normalized transcript per million (nTPM), were transformed into Z-scores.
Fig 9
Fig 9. Systematic exploration of genes with spatial autocorrelation in the reconstructed expression patterns in planarian tomo-seq data.
(A) A histogram showing the Moran’s I values of the reconstructed expression patterns for 214 genes encoding planarian homologs of mouse transcription factors. (B-F) Reconstruction results for the expression of (B) Creb3l1 (DjGI019398_001), (C) Kmt2d (DjGI001455_001), (D) Smad1 (DjGI002591_001), (E) Mzf1 (DjGI000932_002), and (F) Foxa2 (DjGI007454_001). The corresponding Moran’s I values are shown.
See this image and copyright information in PMC

References

    1. Junker J, Noël E, Guryev V, Peterson K, Shah G, Huisken J, et al.. Genome-wide RNA Tomography in the Zebrafish Embryo. Cell. 2014;159(3):662–675. doi: 10.1016/j.cell.2014.09.038 - DOI - PubMed
    1. Ebbing A, Ábel Vértesy, Betist MC, Spanjaard B, Junker JP, Berezikov E, et al.. Spatial Transcriptomics of C. elegans Males and Hermaphrodites Identifies Sex-Specific Differences in Gene Expression Patterns. Developmental cell. 2018;47(6):801–813.e6. doi: 10.1016/j.devcel.2018.10.016 - DOI - PubMed
    1. Mayeur H, Lanoizelet M, Quillien A, Menuet A, Michel L, Martin KJ, et al.. When Bigger Is Better: 3D RNA Profiling of the Developing Head in the Catshark Scyliorhinus canicula. Frontiers in cell and developmental biology. 2021;9:744982. doi: 10.3389/fcell.2021.744982 - DOI - PMC - PubMed
    1. Yvernogeau L, Klaus A, Maas J, Morin-Poulard I, Weijts B, Schulte-Merker S, et al.. Multispecies RNA tomography reveals regulators of hematopoietic stem cell birth in the embryonic aorta. Blood. 2020;136(7):831–844. doi: 10.1182/blood.2019004446 - DOI - PMC - PubMed
    1. Burkhard SB, Bakkers J. Spatially resolved RNA-sequencing of the embryonic heart identifies a role for Wnt/β-catenin signaling in autonomic control of heart rate. eLife. 2018;7. doi: 10.7554/eLife.31515 - DOI - PMC - PubMed

LinkOut - more resources