159 |

This is the analysis of the Human squamous cell carcinoma (SCC) ST 160 | data analysis with IRIS.

161 |

Due to the storage issue of github, the required input dataset is 162 | stored here: on the google 163 | drive. Please see more details on the required data input from the 164 | main 165 | tutorial.

166 |
167 |

load the list of spatial transcriptomics data across multiple 168 | slices

169 |
library(IRIS)
170 | ### read in the list of spatial transcriptomics dataset across multiple tissue slices
171 | sp_input = readRDS("./countList_SCC_ST.RDS")
172 | #### extract the count data list and location list
173 | spatial_countMat_list = sp_input$spatial_countMat_list
174 | spatial_location_list = sp_input$spatial_location_list
175 |
176 |
177 |

load the reference scRNA-seq data

178 |
#### read in the example reference scRNA-seq data, 
179 | sc_input = readRDS("./scRef_input_SCC.RDS")
180 | #### extract the count data, metadata, column name indicating cell type annotation, and column name indicating samples/subjects.
181 | sc_count = sc_input$sc_count
182 | sc_meta = sc_input$sc_meta
183 | ct.varname = sc_input$ct.varname
184 | sample.varname = sc_input$sample.varname
185 |
186 |
187 |

integrative and reference-informed spatial domain detection with 188 | IRIS

189 |
##### create IRIS object
190 | IRIS_object <- createIRISObject(
191 |   spatial_countMat_list = spatial_countMat_list,
192 |   spatial_location_list = spatial_location_list,
193 |   sc_count = sc_count,
194 |   sc_meta = sc_meta,
195 |   ct.varname = ct.varname,
196 |   sample.varname = sample.varname,
197 |   minCountGene = 100,
198 |   minCountSpot =5) 
199 | #### preselect the number of spatial domains 
200 | numCluster = 12
201 | #### perform the integrative and refernce-informed domain detection
202 | IRIS_object <- IRIS_spatial(IRIS_object,numCluster = numCluster)
203 |

The results of detected domains are stored in 204 | IRIS_object@spatialDomain.

205 |
print(IRIS_object@spatialDomain[1:4,])
206 |        Slice spotName  x  y IRIS_domain
207 | 10x26 Slice1    10x26 10 26           6
208 | 10x28 Slice1    10x28 10 28           6
209 | 10x30 Slice1    10x30 10 30           6
210 | 10x32 Slice1    10x32 10 32           3
211 |
212 |