/
about.Rmd
129 lines (94 loc) · 7.5 KB
/
about.Rmd
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
---
title: "About"
author:
- Joseph Bergenstråhle, SciLifeLab, Royal Institute of Technology (KTH)
- Ludvig Larsson, SciLifeLab, Royal Institute of Technology (KTH)
date: ''
#output:
# html_document:
# theme: darkly
# css: style.css
# number_sections: false
# toc: yes
# toc_depth: 3
# toc_float:
# collapsed: false
# smooth_scroll: true
# pdf_document:
# number_sections: false
# toc: yes
# toc_depth: 3
---
<style type="text/css">
div.main-container {
background-color: #000000 !important;
max-width: 1400px;
margin-left: auto;
margin-right: auto;
}
</style>
<style>
#TOC {
background: url("https://www.spatialresearch.org/wp-content/uploads/2019/09/str-logo-spatial_research_3@2x.png");
background-size: contain;
padding-top: 100px !important;
background-repeat: no-repeat;
op: 5%;
opacity: 0.8;
width: 500px;
color: white;
border-color: #000000 !important;
}
</style>
<style> code, pre{
background-color: #000000 !important;
color: white !important;
}
</style>
<style>
body {
color: white
}
</style>
<style>
.list-group-item.active, .list-group-item.active:focus, .list-group-item.active:hover {
background-color: #375a7f;
}
</style>
```{r setup, include=FALSE}
knitr::opts_chunk$set(echo = TRUE, message = FALSE, warning = FALSE, autodep = FALSE)
```
# Background
STUtility R-package is an effort to create a user-friendly visualization and analysis tool for analysis of spatial transcriptomcis data. It's built around [Seurat](https://github.com/satijalab/seurat), which is a single-cell genomics toolkit.
### Spatial Transcriptomics (ST)
---
Spatial Transcriptomics is a method that allows visualization and quantitative analysis of the transcriptome in individual tissue sections by combining gene expression data and microscopy based image data. The invention was presented [in science](https://science.sciencemag.org/content/353/6294/78) in 2016. Prof Joakim Lundeberg (KTH Royal Institute of Technology) and Prof Jonas Frisén (Karolinska Institutet) received a key initial support from the Knut and Alice Wallenberg Foundation in 2012 to develop and use the Spatial Transcriptomics technology for analysis and discovery of transcriptional patterns in tissue, with a focus on the brain. The method has received increasing attention and is currently the basis of several national and international collaborations. The research is predominantly done at [Science for Life Laboratory](https://www.scilifelab.se/), Stockholm.
For details, see the publication. In short, the schematic below gives a brief overview of the concept.
An introductory animation is available on our website: http://www.spatialresearch.org/
![Schematic Spatial Transcriptomics](assets/st_method_1.png)
The array featured 1000 capture-spots, 100 µm in diameter and accordingly we refer this as the "1k" array in this tutorial and package parameters. The ST technology was futher developed, and the capture-spot number increased to 2000, hence the "2k" array.
### 10X Visium
In Dec 2018, 10X Genomics [aquired](https://www.10xgenomics.com/news/10x-genomics-acquires-spatial-transcriptomics/) Spatial Transcriptomics, and in Nov 2019 they [started shipping](https://investors.10xgenomics.com/news-releases/news-release-details/10x-genomics-begins-shipments-visium-spatial-gene-expression) the Visium array, which is a further development of the original ST array. This array features 5000 capture-spots 55µm in diameter.
# Notes about the tool
### Selecting spots - original ST arrays
The gene expression data consists of a count matrix with genes in rows and "capture-spots" in columns. Each spot represents a small area on an ST array from which the captured transcripts have been barcoded with a unique sequence. The unique barcode makes it possible to map the transcripts onto a spatial position on the tissue section and would be equivalent to a cell specific barcode in scRNA-seq data but can tag a mixture of transcripts from multiple cells. The spatial position of a spot is an (x, y) coordinate that defines the centroid of the spot area. These spatial coordinates are stored in the spot ids (column names) and allows us to visualize gene expression (and other spot features) in the array grid system. However, if you want to overlay a visualization on top the HE image you want to make sure that the spot coordinates are exact in relation to morphological features of the image. When the spots are printed onto the ST array surface, they will sometimes deviate from the (x, y) coordinates given by the spot ids and should therefore be adjusted. In addition to the spot adjustment, you will also need to label the spots that are located directly under the tissue. Spot adjustment and selection can be done automatically using our [ST spot detector](https://github.com/SpatialTranscriptomicsResearch/st_spot_detector) web tool which outputs a table of adjusted coordinates and labels for the spots under tissue.
### Selecting spots - 10X Visium arrays
10X Genomics has developed their own tool for visualization and spot selection called [SpaceRanger](www.10xgenomics.com) [fixa denna link]. In SpaceRanger you can .....
### Multiple samples
The STUtility tool was developed with the goal of multiple sample inputs. As with all biological data, using multiple samples add power to the analysis and is a necessity to enable comprehensive insight which otherwise suffers from stochastic uncertainty. Within this vignette, we display how you can input multiple samples, look for aggravating circumstances like batch effects and missing data, apply methods to correct such if present, get a holistic picture of your data as well as conduct more in depth analysis in various ways.
### Seurat workflow
We have extensively tried different methods and workflows for handling ST data. While all roads lead to Rome, as of the date of this writing we find the [Seurat approach](https://satijalab.org/seurat/) to be the most well suited for this type of data. Seurat is an R package designed for single-cell RNAseq data. Obviously, this deviates from the data that the ST technology currently produce, as the resolution on the array implies that each capture-spot consists of transcripts originating from multiple cells. Nevertheless, the characteristics of the ST data resembles that of scRNAseq to a large extent. Note that the STUtility package requires Seurat v3.0 or higher.
The data obtained from an ST experiment can treated like a scRNA-seq experiment and be processed and analyzed using the Seurat package. STUtility provides image processing and visualization functionallity on top of this framework.
### Naming conventions
For users familiar with the Seurat workflow, there are a number of Seruat plotting functions, e.g. `Seurat::FeaturePlot()`, those plotting functions all have a "ST version", which is called upon by adding "ST." prior to the original function name e.g. `STutility::ST.FeaturePlot()`.
The external STUtility functions are following a PascalCase convention.
<hr />
<p style="text-align: center;">A work by <a href="j.bergenstrahle@scilifelab.se">Joseph Bergenstråhle</a> and <a href="ludvig.larsson@scilifelab.se">Ludvig Larsson</a></p>
<!-- Add icon library -->
<link rel="stylesheet" href="https://cdnjs.cloudflare.com/ajax/libs/font-awesome/4.7.0/css/font-awesome.min.css">
<!-- Add font awesome icons -->
<p style="text-align: center;">
<a href="https://www.spatialresearch.org" class="fa fa-beer"></a>
</p>