bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
bioRχiv Article Spatial tissue profling by imaging-free molecular tomography Halima Hannah Schede*,1, Christian G. Schneider*,1,2, Johanna Stergiadou3,4, Lars E. Borm3, Anurag Ranjak1, Tracy M. Yamawaki5,6, Fabrice P.A. David1, Peter Lönnerberg3, Gilles Laurent5, Maria Antonietta Tosches7, Simone Codeluppi3, Gioele La Manno1, 1Brain Mind Institute, School of Life Sciences, École Polytechnique Fédérale de Lausanne, CH-1015, Lausanne, Switzerland. 2NeuroCure Cluster of Excellence, Charité-Universitäts- medizin Berlin, corporate member of the Freie Universität Berlin and Humboldt-Universität zu Berlin, and Berlin Institute of Health, 10117 Berlin, Germany. 3Division of Molecular Neurobiology, Department of Medical Biochemistry and Biophysics, Karolinska Institute, S-17177, Stockholm, Sweden. 4Current address: Cartana AB, S-17165, Stockholm, Sweden. 5Max Planck Institute for Brain Research, 60438 Frankfurt am Main, Germany. 6Current Address: Amgen, Inc, South San Francisco, CA 94080 7Department of Biological Sciences, Columbia University, NY10027 New York, USA. *Tese authors contributed equally to this work. Genomics techniques are currently being adapted to provide spatially resolved omics profling. However, the adaptation of each new method typically requires the setup of specifc detection strategies or specialized in- strumentation. A generic approach to spatially resolve diferent types of high throughput data is missing. Here, we describe an imaging-free framework to localize high throughput readouts within a tissue by combining compressive sampling and image reconstruction. We implemented this framework to transform a low-input 5 RNA sequencing protocol into an imaging-free spatial transcriptomics technique (STRP-seq) and validated this method with a transcriptome profling of the murine brain. To verify the broad applicability of STRP-seq, we applied the technique on the brain of the Australian bearded dragon Pogona vitticeps. Our results reveal the molecular anatomy of the telencephalon of this lizard, providing evidence for a marked regionalization of the reptilian pallium and subpallium. Overall, the proposed framework constitutes a new approach that allows 10 upgrading in a generic fashion conventional genomic assays to spatially resolved techniques. Correspondence to G.L.M. ( [email protected] / @GioeleLaManno / gioelelamanno.com) Introduction captured on barcoded supports [7]–[17]. However, in-situ Spatially resolved molecular atlases constitute fundamental hybridization techniques require a preselection of target resources to the scientifc community [1]. In situ hybrid- genes, and sequencing methods sufer from resolution lim- 25 15 ization (ISH) compendia such as the Allen Brain Atlas or its and suboptimal capture efciency. Genepaint are extensively used as references to guide the An alternative approach for spatially resolving the distribu- mapping of newly discovered cell types [2]–[6]. Recently, tion of gene expression patterns relies on sampling a tissue the development of new technologies that can simultane- regularly on a grid by laser capture microdissection (LCM). ously map the expression of multiple genes on the same Afer LCM sampling, a next-generation sequencing (NGS) 30 20 tissue sample have accelerated the generation of new atlas- library can be prepared using a low-input protocol, and the es. Tese techniques are based on either quantitative in-si- gene expression pattern can be quantifed [18]. Overall, tu hybridization or on the sequencing of RNA molecules combining microdissection-based approaches with sin- a b c Secondary Slicing Angle α Strips Data for each gene Proj. α Proj. β Proj. γ Primary Sectioning
Angle α Probabilistic Model Strips Expectation Lysis/Capture P(Image | Data) RT-PCR Sections Library Prep Sequencing ... Iterations Angle β ≈ Strips ... ≈
Tissue Levels Voxel Genes Assumption Angle γ strip # Strips Parallel slice-projection α Final Maximum a posteriori Reconstruction estimation Figure 1 | Sampling and reconstruction approach to resolve the spatial localization of genomics data. a, Schematic representation of the tissue processing steps. The tissue is frst cut along the axis of interest to obtain consecutive 14 µm primary sections. The spatial pattern of the signal is assumed 35 constant between primary sections. Each section is then sliced at a diferent angle to generate tissue strips (secondary sections). Only three angles are shown for simplicity, but more are typically used. b, Library preparation and generation of raw data. Cells of tissue strips are lysed and mRNA is captured. Reverse transcription is followed by PCR, tagmentation and sequencing. For each primary section, the gene expression within each tissue strip can be concatenated as a parallel-slice projection vector. c, Reconstruction of gene expression patterns performed by the Tomographer algorithm. Gene expression count values for the strips are modeled as samples from a negative binomial distribution. The expected value is the sum of the intensity of the image to reconstruct. 40 The spatial localization of the gene of interest is obtained by maximum a posteriori estimation, considering probabilistic priors on the distribution of pixel intensities.
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 1 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license. gle-cell RNA sequencing has proved successful to delineate To this aim, we designed a user-friendly tissue sampling relationships between cells populating a niche and their po- scheme based on cryosectioning (Supplementary Fig. 1). 45 tency [19]. However, the sensitivity of LCM-based assays is Initially, the tissue is cut into consecutive thin slices, which limited by laser-induced degradation of nucleic-acids and we refer to as “primary sections”. Subsequently, primary the technical complexity of handling and preparing libraries sections are further sliced along an orthogonal plane at pre- 105 for multiple microscopic samples is considerable. defned orientations resulting in tissue strips (“secondary sections”) (Fig. 1a, see Online Methods). Here, we reason Tomo-seq is a more accessible and imaging-free spatial se- that primary sections are thin enough so that the diference 50 quencing method in which a tissue sample is cryosectioned in the signal distribution between them is negligible in re- into thin slices that are used to prepare NGS libraries [20], lation to the target resolution (i.e. we chose 14µm for the 110 [21]. Using this method, it is possible to generate a one-di- experiments below). Under this assumption, diferent sets of mensional transcriptomic pattern that allows the evaluation secondary sections obtained at diferent orientations can be of gene expression variation over a spatial axis. Te authors considered as a resampling of the same distribution. 55 also demonstrated the reconstruction of three-dimensional gene expression patterns by using the iterative proportional It is helpful to think of the set of measurements obtained ftting (IPF) algorithm to combine three orthogonal one-di- from one series of secondary sections as an analogue of a 115 mensional profles generated using multiple genetically parallel-beam sum projection in ray-based tomography equivalent samples [20]. However, IPF cannot solve this [28], [29]. Henceforth we refer to this set of measurements 60 reconstruction problem for complex anatomical confgu- for each particular angle as a “parallel-slice projection”. rations and lacks robustness to noise. Moreover, an entire However, in contrast to classical computer tomography, we specimen is required to obtain each of the 1D profles, mak- developed a probabilistic algorithm that requires only a few 120 ing 3D reconstructions impossible for many applications parallel-slice projections to reconstruct the spatial pattern where multiple genetically identical specimens are not avail- that generated the data. Tis is achieved by modeling the 65 able [22]–[24]. noise, setting priors on the degree of smoothness and spar- sity of the solution, and then fnding the maximum a poste- Even though the current atlassing eforts are focused on riori estimate (Fig. 1c; see Section 2 of the Online Methods). 125 mapping gene expression patterns, the possibility to create atlases probing other modalities will be benefcial for the We implemented this framework for the quantifcation of entire scientifc community [25], [26]. However, other NGS- gene expression and developed Spatial Transcriptomics by 70 based technologies, including those to measure chromatin Reoriented Projections and sequencing (STRP-seq), a meth- accessibility, histone modifcations, non-coding or nascent od that combines the sampling strategy presented above RNAs lack a spatial counterpart. Tus, the development of with a customized low-input RNA-seq protocol based on 130 a generalized approach to convert every existing NGS pro- STRT-seq chemistry (Fig. 1b, see Section 1 of Online Meth- tocol into a spatially resolved technique, without modifying ods) [30]. In this technique, we obtain a parallel-slice pro- 75 the detection strategy, would accelerate the transition to a jection for each gene by quantifying the reads that map to a spatial-omics paradigm. Motivated by this, we developed transcript in each of the secondary sections (Fig. 1b,c). a strategy to perform such an adaptation. Here, we present a framework that brings this transition to life by using a Tomographer’s maximum a posteriori estimates pro- 135 compressed sensing tissue sampling strategy based on mul- vide accurate reconstructions of complex spatial pat- 80 ti-angle-sectioning and an associated algorithm that enables terns. the reconstruction of complex spatial patterns (see Online To demonstrate the potential of the proposed framework for Methods). In this study, we demonstrate that our spatial-om- studying non-trivial spatial patterns, we carried out in-silico ics approach allows the reconstruction of spatial expression simulation experiments using the Adult Mouse Allen Brain 140 ISH Atlas profles as ground truth (see Online Methods Sec- patterns on a transcriptome-wide scale by benchmarking it 210 85 against the mouse Allen Institute ISH Mouse Brain atlas [1]. tion 3). We compared ground truth ISH profles to the re- Furthermore, we apply this framework to study the brain of constructions obtained from simulation experiments using the lizard non-model organism Pogona vitticeps. We reveal both our newly developed algorithm dubbed Tomographer quantitative aspects of the molecular anatomy of the pallium and the IPT-based reconstruction proposed for Tomo-seq 145 [20]. Since a 3D reconstruction problem is more strongly and comparative relationships between the reptilian and the 215 90 mammalian brain. underdetermined (i.e. signifcantly more unknowns than equations) than the one we propose to solve, we include a Results simple 2D adaptation of Tomo-seq that constitutes a fairer Te spatial distribution of a signal can be determined comparison (See Online Methods Section 3.2). 150
by compressed sampling and reconstruction. To assess the information that can be recovered by each of 220 Te objective of the proposed framework is to determine the sampling methods, we determined the reconstruction 95 the abundance of a target signal at a given tissue location accuracies of genes that are characterized by diferent spar- through the use of a compressed sensing approach [27]. Te sity and distribution patterns (Fig. 2a). We calculated the approach consists in measuring several locations of the tis- Pearson’s correlation coefcient, the total absolute diference 155
sue collectively so as to simplify tissue sampling while al- among pixels, and the relative error between the recon- 225 lowing the subsequent reconstruction of the signal at every structions and the ground truth. Evaluating these metrics 100 single point (Fig. 1). we concluded that our method outperforms Tomo-seq for
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 2 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
Iterative proportional fitting Tomographer Figure 2 | Reconstruction performance of Tomog- a Ground truth b 1.0 3D (sections) 2D (strips) 2D (strips) rapher on simulated data. a, Comparison of the nov- 160 el sampling and reconstruction algorithm (Tomogra- 0.6 pher) with existing methods: Iterative proportional
LAMP5 ftting 2D and 3D (Tomo-seq). Starting from gene ex- Correlation 0.2 pression patterns extracted from the Allen Brain Atlas, 1000 the simulated projections are computed and the im- 165
500 age is then reconstructed (noiseless case). b, Recon-
TAC1 struction accuracy for TAC1 evaluated using diferent 0
Different pixels metrics between the images: Pearson’s correlation co- 1.0 efcient, relative error and the number of pixels that
0.6 difer more than the range between the maximum 170 and minimum pixel intensity values are calculated PLP1 0.2 with respect to the ground truth. c, Resilience to low Relative error Signal intensity IPT3d IPT2D Tomog. signal to noise ratio of the Tomographer algorithm. min max Noisy data are simulated by projecting the expression c Expectation Reconstructions Parallel slice-projections and adding Poisson noise. d, Reconstruction accuracy 175 400 Simulated Data μ = 50 at diferent expression levels for a given gene. The rel- ative noise increases as the average expression level decays, causing a decrease in reconstruction accuracy 0 for lowly expressed genes. Pearson correlation coef- 20 μ = 5 cient (blue) and Peak signal-to-noise ratio (PSNR) (red) 180 shown for fve realizations with additive Poisson noise at various expression levels. The shaded area is the standard deviation over 10 realizations of the noisy 0 projections. e, Reconstruction accuracy across a sam- 4 μ = 0.5 ple of 100 genes of the Allen Brain Atlas with noise 185 simulated as in (c). Reconstructions were compared to the ground truth using Pearson’s correlation coef- Expression Level (molecule counts) fcient. Violin plots show the distribution of similarity 0 across all the genes. Reference violins show the simi- Angle 41 Angle 77 Angle 113 Angle 149 Angle 185 (146 strips) (130 strips) (133 strips) (142 strips) (128 strips) larity of each ground truth gene with the most similar 190 gene (green) and the gene median similarity (across Similarity to the ground truth d e 1.0 g Distance 0.3 0.78 1.27 1.75 the entire Allen Brain Atlas). f, Efect of the number 1.0 of projections (i.e. angles) on the reconstruction ac- 0.6 Ground truth 0.8 curacy. Reconstructions are scored using Pearson’s 0.6 0.2 Reconstriction correlation coefcient. Standard deviation is comput- 195 0.4 ed over diferent realizations of the noise for average Pearsons s
Pearson s -0.2 0.2 expression levels of 5 and 10. g, Assessment of two- 0.0 Tomographer Most similar ABA gene points discriminative power and resolution. The two- Reconstructed ABA gene with median 500 Two Gaussians (Avg.) similarity ±Std. deviation point efective resolution was defned as the distance 0.04 0.15 0.5 1.5 4.9 15.2 48 150 f One Gaussian (Ctrl.) 40 Similarity to the ground truth 400 ±Std. deviation at which the two-Gaussian model outperforms the 200 0.93 single-Gaussian model in explaining the reconstruct- 36 300 Resolution ed images, the Bayesian Information Criterion (BIC) 0.88 p < 0.05 was used as a metric. Resolution was determined to
PSNR 30 200
(BIC x1e3) (BIC be 1.15 times the strip thickness (n=16, p-value < 2 × 26 0.83 Mean expression level 100 10^-5). Standard deviation calculated among random Discrimination score Discrimination 205 Pearsons s 10 molecules realization of the location of the points, fxing the dis- 0.78 5 molecules 0 0.04 0.15 0.5 1.5 4.9 15.2 48 150 tance. Gene Expression Level 3 4 5 6 7 0.3 0.54 0.78 1.0 1.27 1.51 1.75 2.0 (log scale) Number of Angles Distances (relative to strip thickness)
the reconstruction of the spatial pattern of a given gene (Fig. squares approach (Supplementary Fig. 2a-d) 2b). To gauge the ability of the technique to generalize and resolve 210 We developed a probabilistic signal reconstruction strategy other spatial distributions, we simulated reconstructions for 230 taking into account that we expect read counts to be distrib- 100 randomly selected genes retrieved from the Allen Brain uted as a Negative Binomial (See Online Methods Section Atlas. Te resulting reconstructions were on average more 2; [31], [32]). In this way, we aim at bufering reconstruc- similar to the ground truth pattern of a gene than to the next tion errors that may arise because of the discrete nature of most similar gene available in the Allen Brain Atlas (Fig. 2e). 215 the data and the heteroscedasticity within each parallel-slice We then analyzed the efect of using diferent numbers of 235 span. Nevertheless, we expect the maximum achievable ac- secondary slicing angles on the reconstruction accuracy in curacy to remain dependent on average expression levels. To order to optimize the reconstruction while maintaining the test the extent of this dependence, we corrupted simulated feasibility of the protocol from an experimental viewpoint. parallel-slice projections drawing realizations from a Neg- We demonstrate that four cutting angles provide results that 220 ative Binomial distribution afer rescaling the average read are a fair compromise between the reconstruction quality 240 count to diferent levels (Fig. 2c-e and Supplementary Fig. and sample processing efort and cost (Fig. 2f). However, we 2). Analyzing the reconstructions, we found that greater ex- decided to use fve cutting angles in the implementation of pression levels correlate with the quality of the output and the protocol to increase resiliency to other technical sources that a minimum average expression level of 1 molecule per of error (e.g. cutting errors or lost strips). 225 strip is required to obtain reliable reconstructions. Notably, the reconstruction accuracy we obtained here was signif- Finally, we evaluated the two-point discriminative resolu- 245 icantly more robust to what we achieved by a naive least- tion of the technique in noisy conditions as a function of
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 3 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
the thickness of the strip generated by secondary slicing. By pixels using the Louvain community detection algorithm. 305 performing Monte Carlo simulations we demonstrate that Te clustering output identifed regions of neighboring pix- the approach is capable of discriminating two distinct points els that recapitulate gross anatomical structures defned by 365 250 provided that their distance is at least 1.15 times the second- the Allen Brain Reference Atlas (Fig. 3k). ary section width (Fig. 2g; 80 μm in this case; See Online Together, these analyses provide evidence that STRP-seq Methods section 3.6). data can be used to reveal the molecular organization of tis- 310 STRP-seq accurately recovers the spatial transcriptome sue heterogeneity. of the mouse brain. STRP-seq allows for the de novo reconstruction of the 255 To demonstrate the potential of STRP-seq to profle anatom- molecular anatomy of the lizard brain. ically complex tissues, we applied the technique to a coronal To demonstrate the applicability of STRP-seq for examining section of the mouse brain (Fig. 3). We cryosectioned the unconventional or rare samples we profled brain sections 315 brain to obtain fve primary and then 679 secondary sec- from Pogona vitticeps, a non-model organism lizard endem- tions, prepared the cDNA libraries and sequenced them ic to semi-arid regions of Australia (Fig. 4a,b; [33], [34]). 260 to a depth allowing a transcriptome-wide quantifcation We used STRP-seq to measure and reconstruct 8,183 anno- [8763 genes detected; an average of 40.3 counts per gene tated genes in the brain of P. vitticeps [35]. To validate our per secondary section]. We used this data to verify two key approach, we compared a set of images from our reconstruc- 320 assumptions for the reconstruction. First, we calculated the tions to ISH experiments and single-cell profling previously correlation of the gene expression over the primary sections reported, confrming that these genes were enriched in the 265 and confrmed their high similarity (median Pearson’s R = same locations. (Fig. 4c,d and Supplementary Fig. 4) [34]. 0.92, Supplementary Fig. 3a). Second, we evaluated the dis- tribution of the noise and found that the read counts closely According to classical neuroanatomical work, the lizard tel- matched a Poisson distribution at low levels of expression, encephalon includes the basal ganglia, derived from the ven- 325 but that they were overdispersed at higher levels (Supple- tral telencephalon or subpallium, a layered cortex and the 270 mentary Fig. 3b-c). non-layered dorsal-ventricular ridge (DVR), both derived from the dorsal telencephalon or pallium. While the ho- Next, we proceeded with the reconstructions (Fig. 3a-b). mology of reptilian and mammalian cortices is established, Based on the results of the previous simulations, we fltered the DVR, which is unique to birds and reptiles, remains an 330 out the genes that were too lowly expressed to allow for re- evolutionary dilemma. Recent work, based on single-cell constructions and retained the parallel-slice projections that RNA-seq but not spatial transcriptomics, suggested that an- 275 varied signifcantly from the trend of total read counts (3880 terior dorsal-ventricular ridge (aDVR) includes two distinct genes; see Online Methods Section 4.2). To assess whether regions, which have never been recognized in any previous we retrieved the correct spatial localization of RNA tran- work [36]. To validate this distinction and quantitatively 335 scripts, we compared 923 reconstructed genes to the ISH evaluate its relevance in a broader context, we studied the data from the Allen Brain Atlas using Pearson’s correlation molecular anatomy of the lizard telencephalon using STRP- 280 coefcient and found a median reconstruction accuracy of seq profles. R = 0.48 (Fig. 3c; See Online Methods). Notably, gene re- constructions with Pearson’s correlation coefcients as low First, we performed a multivariate analysis of STRP-seq data as 0.2 could already be recognized as similarly distributed examining single pixels. We began by calculating a pairwise 340 to the ground truth, indicating that anatomical localizations gene expression similarity matrix between diferent loca- 285 can be correctly determined even for low-scoring recon- tions of the tissue (Fig. 4e). Sorting the similarity matrix structions (Fig. 3d-e). In addition to this, across total recon- using the SPIN algorithm revealed sharply defned groups structions, 2.1% had parallel-slice projections that were in- of covarying pixels [37]. Visualizing the cell-to-cell similar- compatible to the Allen Brain Atlas ground truth, suggesting ity values in correspondence to the location in the brain, we 345 that there might be cases in which ISH and RNA-seq detect confrmed that the block-like covariance structure corre- 290 diferent signals, for example, capturing diferent splicing sponded to anatomical regions in the lizard brain. A more isoforms (Fig. 3g and Supplementary Fig. 3d). continuous variation was only observed at a fner level with- in those blocks (Fig. 4f). To further evaluate if STRP-seq profles are sufcient to de- lineate anatomical distinctions de novo, we computed pair- To explicitly delineate molecularly distinct regions, we per- 350 wise gene expression similarity scores between all the pixels. formed unbiased clustering of the pixels (Fig. 4g). We found 295 As predicted, pixels belonging to similar anatomical areas 15 clusters of contiguous pixels that explained most of the contained more similar gene expression levels. For example, variation present in the data (88.12 %) (Fig. 4g). Te 15 clus- the correlation between two points in the thalamus was 0.97 ters provided a nuanced description of the variability (cf. while it was only 0.20 between points in the thalamus and Fig 4e) and could be grouped to identify main regions of 355 the cortex (Fig. 3h). Additionally, displaying the similarity the lizard brain. Nonetheless, each cluster was defned by a 300 between a specifc point to the remaining points as an im- marked gene expression identity. Gene enrichment analysis age, we obtain a profle that demarcates diferent anatomical highlighted genes that were expressed specifcally in indi- regions (Fig. 3i,j). vidual areas such as SS2, FAT1 and NR4A2 (Fig. 4h,i and Supplementary Fig. 5). Other genes were enriched in mul- 360 Finally, to actively fnd those anatomical boundaries sug- tiple clusters, either marking a whole main area (DDX58, gested by the correlation structure of the data, we clustered BCL11B) or displaying a gradient-like pattern such as IL7R, Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 4 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
RTN4RL2 and BAALC (Fig. 4i). tween the subpallium and the other regions. Surprisingly, we 370 found that the dissimilarity between the long-believed ho- Representing the similarity relations between the 15 clusters mogenous medial aDVR and the lateral aDVR was as large 365 using hierarchical clustering revealed a molecular organiza- as the dissimilarity between the pallial cortex and the lateral tion of this part of the anterior telencephalon in four main aDVR, two functionally distinct regions. When comparing regions. Te regions corresponded to the pallial cortex, the the genes enriched in both regions of the aDVR, we found 375 subpallium, the medial aDVR and the lateral aDVR (Fig. the lateral part harboring genes encoding for membrane 4i). Te dendrogram suggested a marked dissimilarity be- Predicted Reference Reconstruction a Parallel slice-projection data Allen Brain Atlas (STRP-seq) 80 AMOTL1 Raw Data b AMOTL1
0 146 130 133 142 128 80 PDE10A PDE10A
0 146 130 133 142 128 80 NGEF NGEF Expression Level (molecule counts)
0 146 130 133 142 128 Angle 9˚ Angle 117˚ Angle 81˚ Angle 45˚ Angle 153˚
c 50 d 80th Percentile e 20th - 80th Percentile f 20th Percentile CPNE7 B3GALT2 40 AMPH
30
Frequency 20 0.788 0.491 0.273 10 ARC REPS2 IFNGR2 0 0.0 0.2 0.4 0.6 0.8 Pearson s Allen Brain g Predicted 25 IFNGR2 Data 0.815 0.472 -0.17
ITGAV RGS14 PNPO Molecule Count 0 Parallel slice-projections 0.813 0.577 0.237 Allen Brain Reconstruction Allen Brain Reconstruction Allen Brain Reconstruction
Thalamus-Thalamus Thalamus-Cortex h 0.7 0.7 i R = 0.97 R = 0.20
C C 1.0 A A Gene Expression
B Gene Expression B Pearson’s R Pixel C Pixel B 0.15 0 0 0 Pixel A Gene Expression 0.7 0 Pixel A Gene Expression 0.7 j Pixels pixel k location Location Allen Brain Reference Pixels 1 P earso Clustering results n ’ s R
0.4
Figure 3 | STRP-seq reconstructions of spatial expression profles for the mouse brain. a, Representative examples of parallel-slice projections at dif- ferent secondary sectioning angles for the genes AMOTL1, PDE10A and NGEF. Raw data points represented by green dots. Parallel-slice projections of the reconstruction as black lines, statistically those values represent the expectation on the counts measured by STRP-seq. b, Ground truth ISH data generated by Allen Brain Institute (left) and reconstructed spatial expression profles (right). c, Reconstruction accuracy of STRP-seq as evaluated against the Allen Brain 380 Atlas ISH data. Plot shows the distribution of Pearson’s correlation coefcients. The histogram is color-coded to indicate low accuracy ( <20th percentile), intermediate accuracy (20-80th percentile), and high accuracy (>80th percentile). d, Examples of high, e, intermediate and f, low accuracy reconstructions are shown side-to-side with the Allen Brain Atlas ISH data. g, Projection of a poorly reconstructed gene (IFNGR2). Ground truth projection depicted in yellow for comparison. h, Comparison of the expression profle of the medial thalamus, cortex and ventral thalamus. Markers (A, B, C) indicate diferent pixels of the image that are used to compute the scatter plots on the right. i, Brain section heat map representing the gene expression similarity (Pearson’s correlation 385 coefcient) of each point with respect to point A. Thalamic region outlined by central black contour. j, Pairwise correlation distance matrix showing the similarity of gene expression across the pixels of the section. Rows and columns are sorted by similarity using SPIN [37]. k, Top: Anatomical annotation by the Allen Brain Institute. Bottom: Unbiased determination of molecular anatomy obtained by pixels clustering. Louvain community detection algorithm was used to cluster the pixels on the basis of their gene expression, no neighborhood constraint was imposed.
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 5 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
Anatomical regions a Fishes Pogona vitticeps b Anterior Amphibians Subpallium Reptiles Tuatara Dorsomedial Cortex Lizards and Medial Cortex Snakes Turtles posterior Dorsal Cortex Crocodiles anterior Dorsal Cortex anterior Dorsoventricular Ridge Birds Lateral Cortex Mammals Posterior
c 600 PENK Parallel slice-projections d Reconstruction ISH 400
200
0 PENK 60 FIGN Reconstruction ISH 40 (molecule counts) Expression Level
20
0 Angle 9 Angle 117 Angle 81 Angle 45 Angle 153 FIGN (75 strips) (100 strips) (109 strips) (92 strips) (84 strips)
Correlation Pairwise similarity e f mapped on pixels g Clusters Main regions Pixels 13 Pallial Cortex 12 14 5 6 11 10 medial lateral 7 aDVR aDVR 15 9 8 2 4 3 1 Subpallum
Pixels h GRM4 IL7R
NT5E SS2
pixel location Pearson s R
colormap -0.1 1
0
G 6
i 2
1
3
5
5
C
2
1
L
A
9
B
4
D
C A
1 5
8
1
L 9
2 5
4
2 B 2 8
3
R Y
Genes C
1
3
4
2
P D
H
L
R
V N H 5
L S A
N 4
L
A X A 3 C
N 1
M E 1 H
P
1 N
B M
D B
R
B
I 4 K X C
L N T P t N X Y L
2 A
N 5 S
Q C R T
P
i N
G 7
Y
A R A
R 1 E R P A T D T D H
C F S v A G C N D C
F P I U A 1
Clusters L
N C R N K N E C M G R S Z P F B F C A N I K T A D C F F G P G B Subp1 1 Subp2 2 Subp3 3 Subp4 4 aDVRl1 5 806 aDVRl2 6 aDVRl3 7 aDVRl4 8 aDVRm1 9 503 aDVRm2 10 aDVRm3 11 424 DMC 12 aDC 13 pDC 14 LC 15
Distance Expression 0 Max
390 Figure 4 | Reconstruction of the molecular anatomy of the lizard brain de novo. a, Left: Phylogenetic tree displaying the relation of Lizards to other taxa. Reptilian family highlighted in blue and subsection to which Pogona vitticeps belongs in red. Right: Illustration of a Pogona vitticeps specimen and its geographical distribution. b, Left: Section of tissue in lizard brain that is investigated in our study. Right: A schematic of the anatomical reference annota- tion for the studied brain section reproduced from Tosches et al. [34]. c, Raw STRP-seq data for the genes PENK and FIGN in P. vitticeps brain. d, STRP-seq reconstruction of expression patterns (left) and validation by ISH (right). e, Multivariate analysis of the spatial transcriptome of P. vitticeps. The pairwise 395 correlation matrix computed between pixels contains information on the molecular anatomy as highlighted by a block-like structure. The matrix was sort- ed by the SPIN algorithm. f, Selected rows of the matrix are visualized at their anatomical locations the pixel against which the correlation is calculated is marked by the white stripe which runs through the matrix. g, Top: the pixels clustered by gene expression reveal the molecular anatomy of the P. vitticeps brain. Hierarchical clustering was performed and the tree was pruned at diferent levels to highlight main clusters and sub-clusters, represented by difer- ently colored regions and by gray lines respectively. Below: Visualization of a selection of genes enriched in distinct sub-clusters. h, Gene enrichment score 400 for a variety of genes with respect to individual sub-clusters. i, Left: Dendrogram displaying the gene expression distance between the clusters. Numbers indicate the distances between diferent branches of the tree. Right: Heatmap showing top genes enriched in individual clusters. Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 6 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
uptake and transport proteins such as ANXA9, TRARG1, of glutamate but not in neurotransmission, had a distinct KCNH5 and other solute carriers suggesting a possible func- expression pattern (Supplementary Fig. 6b). Furthermore, tional specialization (Fig. 4i and Supplementary Fig. 6a). In the genes for the monoamine associated neuropeptide 430 405 contrast, the diferent regions of the pallial cortex (medial, Neurotensin (NTS), the noradrenaline receptor ADRB3, dorsomedial, anterior dorsal, posterior dorsal and lateral ADRA2A, and the monoamine-degrading enzyme MAOA cortex), that have been extensively characterized as func- colocalized dorsally, suggesting that cortical neurons in this tionally distinct, displayed less marked dissimilarities. area receive monoaminergic inputs [40]. Genes expressed by cholinergic neurons were enriched in ventral regions, in 435 STRP-seq data allows for the analysis of regional iden- correspondence to the subpallium, including the co-trans- 410 tities in the lizard brain. mitter TAC3, TAC1, the nicotinic acetylcholine receptor Our clustering analysis of spatial data revealed an unex- CHRM4, and the enzyme CHAT, suggesting the presence pected molecular heterogeneity in the lizard anterior telen- of a population cholinergic cells in this region (Fig. 5a and cephalon. To better understand the molecular basis of this Supplementary Fig. 6b). 440 diference, we examined the expression pattern of genes 415 characteristic of diferent neuronal types (Fig. 5a and Sup- To interpret the heterogeneity in a border context, we used plementary Fig. 6). the gene expression data to support homology relationships with other vertebrate brain areas. First, we examined in the We found genes typically expressed by GABAergic neurons, lizard brain the expression of genes known to characterize such as the transcription factors TSHZ1, MEIS2, SIX3 and the avian and murine subpallium, a region of the vertebrate 445 the postsynaptic GABA receptor GABRG3 expressed ven- brain that has an essential role in motoric and sensory func- 420 trally [38]. Additionally, genes expressed by glutamatergic tions [41], [42]. We detected marker RNAs of the murine neurons, such as the transcription factors SATB2, NEU- subpallial striatum enriched in the lizard subpallium such ROD2 and LHX2 and the postsynaptic glutamate receptor as DRD1 and ISL1 (Fig. 5b and Supplementary Fig. 6e). In GRIN2B [39], were enriched in ventral and dorsomedial addition, we localized DLX5 and LHX8, which are known to 450 areas. GAD2 and ALDH5A1, which encode for enzymes characterize the entire subpallium and the palladial-preop- 425 involved in the biosynthesis of glutamate and GABA, were tic area respectively [41] (Fig. 5b). expressed consistently with the aforementioned markers. In contrast, GOT1, an enzyme involved in the metabolism Ten, we examined the expression of well-described mam-
a Gabaergic Glutamatergic Monoaminergic Cholinergic SIX3 SATB2 NTS TAC3
Transcription Factors/ Neuropeptides
GABRG3 GRIN2B ADRB3 CHRM4
Receptors
GAD2 ALDH5A1 MAOA CHAT
Enzymes
b c DRD1 LHX8 DLX5 HPCA GNG2 Subpallium medial aDVR d PTPRN2 JUN CPNE7 PLEKHG5 DUSP26 DMC/MC
Figure 5 | Identifcation of regional identities in the lizard brain. a, Localization of genes belonging to diferent ontology classes including neuropep- tides, transcription factors, neurotransmitter receptors and enzymes involved in neurotransmitter synthesis. Genes are arranged column wise by neuro- 455 transmitter type: GABAergic, glutamatergic, monoaminergic (monoamine-receiving neurons) and cholinergic genes. b, Visualization of genes expressed in the lizard brain known to be markers of the subpallium (DLX5), the striatal subdivision (DRD1) and the pallidal-preoptic subdivision (LHX8). c, Visual- ization of reptilian medial aDVR marker genes. d, Reptilian genes that are homologous of murine hippocampal markers are located in the dorsomedial cortex (DMC) and medial cortex (MC) (PTPRN2: pan-cornus ammonis, JUN: pan-dentate gyrus, CPNE7: granule cells in ventral dentate gyrus, PLEKHG5: granule cells in dorsal dentate gyrus, DUSP26: mossy fbers). 460
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 7 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
malian claustrum markers (RGS12, SYNPR, CNR1, HTR1B, We demonstrate how STRP-seq can strengthen our un- HTR1C, CACNA1B) in the lizard brain and found them to derstanding of evolutionary relationships between brain 520 be expressed in the medial aDVR, where they colocalized regions. We show this by supporting hypotheses on region with bona fde markers for this region (HPCA, GNG2, homology with molecular data and by presenting further ev- 580 465 ADARB2, RORB, CPNE4) (Fig. 5c and Supplementary Fig. idence for an anatomically divided aDVR [34], [36]. In par- 6c,d). Tis fnding confrms that the medial aDVR in lizards ticular, we suggest a starker molecular distinction than pre- is the claustrum homologue, as suggested by electrophysio- viously thought. In the medial aDVR, which is considered 525 logical characterization, axonal tracing and single-cell RNA- the reptilian claustrum homologue, we discovered a gene seq [36],[43]. expression signature related to mammalian telencephalic 585 inhibitory neurons, suggesting a diferential recruitment of 470 Lastly, to confrm the hypothesis that the reptilian pallial this set of genes to glutamatergic cells of the medial aDVR. cortex contains a region homologous to the hippocampus, In the lateral aDVR, we unexpectedly detected the expres- 530 we localized genes that are known to be expressed in dif- sion of genes related to non-cholinergic motor neurons, sug- ferent subregions of the murine hippocampus [34], [44]. 590 gesting that we capture mRNA of neurons projection from Expressed in the medial and dorsomedial cortex we found the brain stem to the lateral aDVR [46]. 475 markers of the cornus ammonis (CA) region 1-4, genes characterizing hippocampal mossy fbers and markers of Taken together, we provide a fundamental generic profling granule cells in the ventral and dorsal dentate gyrus (Fig. 5d scheme compatible with diferent types of genomics meas- 535 595 and Supplementary Fig. 6f). urements that can be set up at a low cost and without the need of specialized instrumentation. We envision our frame- Finally, we performed an independent comparison of the work will help bridge the gap between the various available 480 examined region with mammalian data, by both cross-ref- NGS techniques and their spatial counterpart. erencing the STRP-seq data with a recently proposed molec- 600 ular taxonomy of cell types in the mouse brain and scoring Acknowledgements 540 the correspondence between the expression pattern of a giv- We thank Sten Linnarsson (Karolinska Insitutet) for allow- en pixel and of the cell clades (See Online Methods; Supple- ing proof of principle tests in his laboratory, Noam Shental, 485 mentary Fig. 7; [45]). Importantly, the analysis confrmed Bar Shalem (Open University Israel) and Amit Zeisel (Tech- the localization of the dentate gyrus in the medial cortex nion) for stimulating discussions, Markus Schuelke (Char- 605 and the ventral positioning of the subpallium (Supplemen- ité) for enabling the participation of C.G.S. in this project. 545 tary Fig. 7). Interestingly, for the anterior dorsal-ventricular Tis work was supported by a grant from the Swiss National ridge, we found a set of genes compatible with telencephalic Science Foundation (CRSK-3_190495) to G.L.M. 490 inhibitory neurons and a non-cholinergic motor neuron sig- 610 nature in the medial and lateral aDVR respectively. Moreo- Author contributions ver, we observed an enrichment of astrocytic markers in the G.L.M. conceived the study design and supervised the pro- ventricular area, most likely marking reptilian ependymo- ject. G.L.M. H.H.S. and C.G.S. analyzed, annotated and in- 550 glial cells that surround the lateral ventricle (Supplementary terpreted the tomography data and wrote the manuscript. 495 Fig. 7; [34]). M.A.T. and T.J. performed P. vitticeps experiments, in situ 615 hybridization and contributed interpreting the results to- Discussion gether with G.L.. G.L.M., A.R., H.S designed and wrote the Here, we present a versatile framework that allows to trans- reconstruction algorithm. S.C. and L.E.B. designed the cry- 555 form existing low-input NGS techniques into methods ca- osectioning scheme. S.C. and J.S. performed the mouse ex- 620 pable of encoding spatial information through compressed periments and sectioning. P.L. ran the bioinformatics pipe- 500 sampling and a probabilistic image reconstruction algo- line. F.P.A.D. built the companion website. All the authors rithm [27]. We report an implementation of our framework critically reviewed the manuscript and approved the fnal for studying transcriptomic data, STRP-seq, and demon- version. 560 strate its application to profling rare samples such as the 625 brain of the Australian lizard P. vitticeps. To facilitate the ad- Code and data availability 505 aptation of the presented framework to applications beyond Source code and templates for customized 3D-printable cry- gene expression studies, we provide Tomographer, a versatile osectioning adaptor manifolds are available at https://github. sofware package to solve fast and parallelized image recon- com/lamanno-epf/tomographer. RNA-seq data is available 630 struction problems for genomic data. at the Gene Expression Omnibus repository (GEO; https:// 565 One limit of our approach is its relatively low resolution www.ncbi.nlm.nih.gov/geo) under accession GSE152989. Results of the reconstruction can be accessed at this paper 510 when compared to imaging based methods and its inability to reveal discontinuous, checkboard-like patterns [10]. For companion website https://strpseq-viewer.epf.ch/. 635 example, our method does not allow for discrimination of REFERENCES the signal at the single-cell level. Instead, it is designed to [1] E. S. Lein et al., “Genome-wide atlas of gene expression in the adult 570 capture broader tissue-level patterns of expression. For these mouse brain,” Nature, vol. 445, no. 7124, pp. 168–176, 2007. 515 reasons, we envision the method to be most useful as a way [2] S. L. Goldman, M. MacKay, E. Afshinnekoo, A. M. Melnick, S. Wu, to link single-cell resolved data to its spatial context, to pro- and C. E. Mason, “Te impact of heterogeneity on single-cell sequenc- 640 fle unique samples (e.g. biopsies), and to study the molecu- ing,” Front. Genet., vol. 10, no. MAR, pp. 1–8, 2019. lar anatomy of non-model organisms. [3] A. Haque, J. Engel, S. A. Teichmann, and T. Lönnberg, “A practical 575
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 8 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
guide to single-cell RNA-sequencing for biomedical research and clin- [26] E. Lundberg and G. H. H. Borner, “Spatial proteomics: a powerful 520 ical applications,” Genome Med., vol. 9, no. 1, pp. 1–12, 2017. discovery tool for cell biology,” Nat. Rev. Mol. Cell Biol., vol. 20, no. 5, [4] E. Lein, L. E. Borm, and S. Linnarsson, “Te promise of spatial pp. 285–302, May 2019. 645 transcriptomics for neuroscience in the era of molecular cell typing,” [27] D. L. Donoho, “Compressed sensing,” IEEE Trans. Inf. Teory, 580 Science (80-. )., vol. 358, no. 6359, pp. 64–69, 2017. vol. 52, no. 4, pp. 1289–1306, Apr. 2006. [5] S. Linnarsson and S. A. Teichmann, “Single-cell genomics: Coming [28] R. Cierniak and R. Cierniak, “Reconstruction from Parallel-beam of age,” Genome Biol., vol. 17, no. 1, pp. 16–18, 2016. Projections,” in X-Ray Computed Tomography in Biomedical Engi- 525 [6] B. Tasic et al., “Adult mouse cortical cell taxonomy revealed by sin- neering, Springer London, 2011, pp. 83–125. 650 gle cell transcriptomics,” Nat. Neurosci., vol. 19, no. 2, pp. 335–346, [29] G. T. Herman and A. Naparstek, “Fast Image Reconstruction 585 2016. Based on a Radon Inversion Formula Appropriate for Rapidly Col- [7] E. Berglund et al., “Spatial maps of prostate cancer transcriptomes lected Data Author ( s ): Gabor T . Herman and Abraham Naparstek reveal an unexplored landscape of heterogeneity,” Nat. Commun., vol. Source : SIAM Journal on Applied Mathematics , Vol . 33 , No . 3 ( Nov 9, no. 1, 2018. ., 1977 ), pp . 511-533 Publish,” SIAM J. Appl. Math., vol. 33, no. 3, pp. 655 530 [8] K. H. Chen, A. N. Boettiger, J. R. Moftt, S. Wang, and X. Zhuang, 511–533, 1977. 590 “Spatially resolved, highly multiplexed RNA profling in single cells,” [30] S. Islam et al., “Quantitative single-cell RNA-seq with unique mo- Science (80-. )., vol. 348, no. 6233, pp. 1360–1363, 2015. lecular identifers,” Nat. Methods, vol. 11, no. 2, pp. 163–166, 2014. [9] Q. Zhu, S. Shah, R. Dries, L. Cai, and G. C. Yuan, “Identifcation of [31] S. Anders and W. Huber, “Diferential expression analysis for se- spatially associated subpopulations by combining scRNAseq and se- quence count data,” Genome Biol., vol. 11, no. 10, p. R106, Oct. 2010. 660 quential fuorescence in situ hybridization data,” Nat. Biotechnol., vol. [32] D. Grün, L. Kester, and A. Van Oudenaarden, “Validation of noise 535 595 36, no. 12, pp. 1183–1190, 2018. models for single-cell transcriptomics,” Nat. Methods, vol. 11, no. 6, [10] S. Codeluppi et al., “Spatial organization of the somatosensory pp. 637–640, 2014. cortex revealed by osmFISH,” Nat. Methods, vol. 15, no. 11, pp. 932– [33] D. Hoops et al., A 3D MRI-based atlas of a lizard brain, vol. 526, 935, 2018. no. 16. 2018. 665 [11] C. H. L. Eng et al., “Transcriptome-scale super-resolved imaging [34] M. A. Tosches, T. M. Yamawaki, R. K. Naumann, A. A. Jacobi, 600 in tissues by RNA seqFISH+,” Nature, vol. 568, no. 7751, pp. 235–239, G. Tushev, and G. Laurent, “Evolution of pallium, hippocampus, and 540 2019. cortical cell types revealed by single-cell transcriptomics in reptiles,” [12] R. Ke et al., “In situ sequencing for RNA analysis in preserved Science (80-. )., vol. 360, no. 6391, pp. 881–888, 2018. tissue and cells,” Nat. Methods, vol. 10, no. 9, pp. 857–860, 2013. [35] A. Georges et al., “High-coverage sequencing and annotated as- 670 [13] S. G. Rodriques et al., “Slide-seq: A scalable technology for meas- sembly of the genome of the Australian dragon lizard Pogona vitti- 605 uring genome-wide expression at high spatial resolution,” Science (80- ceps,” Gigascience, vol. 4, no. 1, p. 45, Dec. 2015. . )., vol. 363, no. 6434, pp. 1463–1467, 2019. [36] H. Norimoto et al., “A claustrum in reptiles and its role in slow- 545 [14] F. Salmén et al., “Barcoded solid-phase RNA capture for Spatial wave sleep,” Nature, vol. 578, no. 7795, pp. 413–418, 2020. Transcriptomics profling in mammalian tissue sections,” Nat. Protoc., [37] D. Tsafrir, I. Tsafrir, O. Zuk, D. A. Notterman, and E. Domany, 675 vol. 13, no. 11, pp. 2501–2534, Nov. 2018. “Sorting points into neighborhoods ( SPIN ): data analysis and visual- 610 [15] S. Shah et al., “Dynamics and Spatial Genomics of the Nascent ization by ordering distance matrices,” vol. 21, no. 10, pp. 2301–2308, Transcriptome by Intron seqFISH,” Cell, vol. 174, no. 2, pp. 363-376. 2005. e16, 2018. [38] S.-M. Virolainen, K. Achim, P. Peltopuro, M. Salminen, and J. 550 [16] S. Vickovic et al., “High-defnition spatial transcriptomics for in Partanen, “Transcriptional regulatory mechanisms underlying the 680 situ tissue profling,” Nat. Methods, vol. 16, no. 10, pp. 987–990, 2019. GABAergic neuron fate in diferent diencephalic prosomeres,” Devel- 615 [17] G. Wang, J. R. Moftt, and X. Zhuang, “Multiplexed imaging of opment, vol. 139, no. 20, pp. 3795–3805, Oct. 2012. high-density libraries of RNAs with MERFISH and expansion micros- [39] S. Bell et al., “Disruption of GRIN2B Impairs Diferentiation in copy,” Sci. Rep., vol. 8, no. 1, pp. 1–13, 2018. Human Neurons,” Stem Cell Reports, vol. 11, no. 1, pp. 183–196, Jul. [18] S. Zechel, P. Zajac, P. Lönnerberg, C. F. Ibáñez, and S. Linnarsson, 2018. 685 555 “Topographical transcriptome mapping of the mouse medial gangli- [40] J. A. Garcia-Sevilla, T. Magnusson, A. Carlsson, J. Leban, and K. 620 onic eminence by spatially resolved RNA-seq,” Genome Biol., vol. 15, Folkers, “Neurotensin and its amide analogue [Gln4]-neurotensin: ef- no. 10, p. 486, 2014. fects on brain monoamine turnover,” Naunyn. Schmiedebergs. Arch. [19] A. E. Moor et al., “Spatial Reconstruction of Single Enterocytes Pharmacol., vol. 305, no. 3, pp. 213–218, Dec. 1978. Uncovers Broad Zonation along the Intestinal Villus Axis,” Cell, vol. [41] W. J. Kuenzel, L. Medina, A. Csillag, D. J. Perkel, and A. Reiner, 690 560 175, no. 4, pp. 1156-1167.e15, 2018. “Te avian subpallium: New insights into structural and functional 625 [20] J. P. Junker et al., “Genome-wide RNA Tomography in the Zebraf- subdivisions occupying the lateral subpallial wall and their embryo- ish Embryo,” Cell, vol. 159, no. 3, pp. 662–675, 2014. logical origins,” Brain Res., vol. 1424, pp. 67–101, Nov. 2011. [21] F. Kruse, J. P. Junker, A. van Oudenaarden, and J. Bakkers, “To- [42] J. Kaas, Evolution of Nervous Systems. 2016. mo-seq,” 2016, pp. 299–307. [43] Q. Wang et al., “Organization of the connections between claus- 695 [22] A. Ebbing et al., “Spatial Transcriptomics of C. elegans Males and trum and cortex in the mouse,” J. Comp. Neurol., vol. 525, no. 6, pp. 630 Hermaphrodites Identifes Sex-Specifc Diferences in Gene Expres- 1317–1346, 2017. 565 sion Patterns,” Dev. Cell, vol. 47, no. 6, pp. 801-813.e6, 2018. [44] G. F. Striedter, “Evolution of the hippocampus in reptiles and [23] G. P. A. Lacraz et al., “Tomo-Seq Identifes SOX9 as a Key Regu- birds,” J. Comp. Neurol., vol. 524, no. 3, pp. 496–517, Feb. 2016. lator of Cardiac Fibrosis during Ischemic Injury,” Circulation, vol. 136, [45] A. Zeisel et al., “Molecular Architecture of the Mouse Nervous 700 no. 15, pp. 1396–1409, 2017. System,” Cell, vol. 174, no. 4, pp. 999-1014.e22, 2018. 635 [24] C. C. Wu et al., “Spatially Resolved Genome-wide Transcriptional [46] H. J. ten Donkelaar and R. de Boer-van Huizen, “Brain stem af- Profling Identifes BMP Signaling as Essential Regulator of Zebraf- ferents to the anterior dorsal ventricular ridge in a lizard (Varanus ex- 570 ish Cardiomyocyte Regeneration,” Dev. Cell, vol. 36, no. 1, pp. 36–49, anthematicus),” Anat. Embryol. (Berl)., vol. 177, no. 5, pp. 465–475, 2016. 1988. 705 [25] I. Kaya et al., “Spatial Lipidomics Reveals Region and Long Chain 640 Base Specifc Accumulations of Monosialogangliosides in Amyloid Plaques in Familial Alzheimer’s Disease Mice (5xFAD) Brain,” ACS 575 Chem. Neurosci., vol. 11, no. 1, pp. 14–24, Jan. 2020.
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 9 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
a b
c cDNA Libraries from Strips
FU FU 100 200 50 100
0 0 Fragment Fragment 35 150 300 500 1000 7000 lenght (bp) 35 150 300 500 1000 7000 lenght (bp)
Supplementary Figure 1 | Generation of secondary sections and cDNA libraries. a, Representation of the device used for positioning primary sections to obtain secondary sections. The primary section is mounted on a slab of ice and positioned in a block according to a predefned angle using a 3D printed manifold. b, The slab is frozen at the predefned angle with the base (top). The base is mounted in a cryostat and sliced at a constant width to generate strips at the desired angle (bottom). c, Example of representative bioanalyzer profles of the cDNA library produced from two strips.
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 10 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
a Simulated noise Parallel slice-projections b Ground truth Expectation 400 μ = 50 LAMP5
0 NegBin MAP Least-squares 20 μ = 5 μ = 50
0
6 μ = 5 μ = 0.5 Expression Level (molecule counts)
μ = 0.5 0 Angle 1 Angle 2 Angle 3 Angle 4 Angle 5 (146 strips) (130 strips) (133 strips) (142 strips) (128 strips)
Simulated noise Ground truth c Expectation Parallel slice-projections d 600 μ = 50 PKP2
0 NegBin MAP Least-squares
40 μ = 5 μ = 50
0
6 μ = 5 μ = 0.5 Expression Level (molecule counts)
μ = 0.5 0 Angle 1 Angle 2 Angle 3 Angle 4 Angle 5 (146 strips) (130 strips) (133 strips) (142 strips) (128 strips)
e Reconstructions Ground Truth μ = 150 μ = 5.0 μ = 0.5 B3GALT2 PVALB PLP1 RGS14
Supplementary Figure 2 | Tomographer’s reconstructions at diferent signal-to-noise ratio conditions. a, Parallel-slice projections of simulated data 710 overlaid with expectation for LAMP5 for three degrees of noise. b, LAMP5 mathematical image reconstructions using Negative Binomial maximum a pos- teriori estimation and least-squares approach constrained to positive results. c, Parallel-slice projections of simulated data overlaid with expectation for PKP2. d, PKP2 mathematical image reconstructions using Negative Binomial Maximum A Posteriori Estimation and least-squares approach constrained to positive results. e, Ground truth and reconstructions using Negative Binomial maximum a posteriori estimation for B3GALT2, PVALB, PLP1 and RGS14.
Schede, Schneider et al. 2020 — Imaging-free tomography-based spatial molecular profling of tissues 11 bioRxiv preprint doi: https://doi.org/10.1101/2020.08.04.235655; this version posted August 4, 2020. The copyright holder for this preprint (which was not certified by peer review) is the author/funder, who has granted bioRxiv a license to display the preprint in perpetuity. It is made available under aCC-BY-ND 4.0 International license.
Supplementary Figure 3 | Verifcation of the assump- Expression Level - log2(Molecule counts) 715 a tions used in the Tomographer’s algorithm. a, Scatterplot 17 17 17 17 17 ° R=0.92 R=0.92 R=0.93 R=0.92 matrix comparing the total gene expression of the primary sections (sum of the counts sequenced across all the sec- ondary sections for each angle). Pearson’s R is indicated for Angle 9
0 0 0 0 0 Expression Level - log 720 each comparison. Histograms on the diagonal show the to- 0 1.0e3 0 17 0 17 0 17 number 0 17 tal molecule count distribution across genes. b, Evaluation of genes 17 17 17 R=0.92 17 R=0.92 R=0.92 of the over-dispersion of the count distribution for diferent Angle 117° genes across the parallel-slice projections. Scatterplots show the coefcient of variation as a function of the mean; the 0 0 0 0 relationship expected for a Poisson distribution is displayed 0 1.0e3 0 17 0 17 0 17 725 number as a reference. Note that in the scatter each gene appears of genes 17 17 17 R=0.92 R=0.92 multiple times: since genes are detected at diferent levels Angle 81° across the parallel-slice projection, we group points by the 2
(Molecule counts)(Molecule reconstruction expectation and, for each group, we compute 0 0 0 0 number 1.0e3 0 17 0 17 730 the coefcient of variation and the mean and plot a dot (See of genes 17 17 Online Methods). c, Molecule count distribution for high, R=0.92 medium and low gene expression. In red, the histogram of Angle 45° observed counts compared with, in gray, an expected Pois- son distribution with the same mean. d, Selection of genes 0 0 0 number 10e3 0 17 735 with low, medium and high correlation values between pro- of genes 17 vided data and Allen Brain ISH projection. Ground truth pro- jection (red) and predicted projection (blue) visualized. Angle 153°
0 0 number 1.0e3 of genes b 2 - - Poisson Distribution Gene 1
0