Chapter 35 Grun mouse HSC (CEL-seq)

35.1 Introduction

This performs an analysis of the mouse haematopoietic stem cell (HSC) dataset generated with CEL-seq (Grun et al. 2016). Despite its name, this dataset actually contains both sorted HSCs and a population of micro-dissected bone marrow cells.

35.2 Data loading

After loading and annotation, we inspect the resulting SingleCellExperiment object:

## class: SingleCellExperiment 
## dim: 21817 1915 
## metadata(0):
## assays(1): counts
## rownames(21817): ENSMUSG00000109644 ENSMUSG00000007777 ...
##   ENSMUSG00000055670 ENSMUSG00000039068
## rowData names(3): GENEID SYMBOL SEQNAME
## colnames(1915): JC4_349_HSC_FE_S13_ JC4_350_HSC_FE_S13_ ...
##   JC48P6_1203_HSC_FE_S8_ JC48P6_1204_HSC_FE_S8_
## colData names(2): sample protocol
## reducedDimNames(0):
## altExpNames(0):

35.3 Quality control

For some reason, no mitochondrial transcripts are available, and we have no spike-in transcripts, so we only use the number of detected genes and the library size for quality control. We block on the protocol used for cell extraction, ignoring the micro-dissected cells when computing this threshold. This is based on our judgement that a majority of micro-dissected plates consist of a majority of low-quality cells, compromising the assumptions of outlier detection.

We examine the number of cells discarded for each reason.

##   low_lib_size low_n_features        discard 
##            465            482            488

We create some diagnostic plots for each metric (Figure 35.1). The library sizes are unusually low for many plates of micro-dissected cells; this may be attributable to damage induced by the extraction protocol compared to cell sorting.

Distribution of each QC metric across cells in the Grun HSC dataset. Each point represents a cell and is colored according to whether that cell was discarded.

Figure 35.1: Distribution of each QC metric across cells in the Grun HSC dataset. Each point represents a cell and is colored according to whether that cell was discarded.

35.4 Normalization

We examine some key metrics for the distribution of size factors, and compare it to the library sizes as a sanity check (Figure 35.2).

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   0.027   0.290   0.603   1.000   1.201  16.433
Relationship between the library size factors and the deconvolution size factors in the Grun HSC dataset.

Figure 35.2: Relationship between the library size factors and the deconvolution size factors in the Grun HSC dataset.

35.5 Variance modelling

We create a mean-variance trend based on the expectation that UMI counts have Poisson technical noise. We do not block on sample here as we want to preserve any difference between the micro-dissected cells and the sorted HSCs.

The lack of a typical “bump” shape in Figure 35.3 is caused by the low counts.

Per-gene variance as a function of the mean for the log-expression values in the Grun HSC dataset. Each point represents a gene (black) with the mean-variance trend (blue) fitted to the simulated Poisson-distributed noise.

Figure 35.3: Per-gene variance as a function of the mean for the log-expression values in the Grun HSC dataset. Each point represents a gene (black) with the mean-variance trend (blue) fitted to the simulated Poisson-distributed noise.

35.7 Clustering

## 
##   1   2   3   4   5   6   7   8   9  10  11  12 
## 259 148 221 103 177 108  48 122  98  63  62  18
Obligatory $t$-SNE plot of the Grun HSC dataset, where each point represents a cell and is colored according to the assigned cluster (left) or extraction protocol (right).

Figure 35.4: Obligatory \(t\)-SNE plot of the Grun HSC dataset, where each point represents a cell and is colored according to the assigned cluster (left) or extraction protocol (right).

35.8 Marker gene detection

To illustrate the manual annotation process, we examine the marker genes for one of the clusters. Upregulation of Camp, Lcn2, Ltf and lysozyme genes indicates that this cluster contains cells of neuronal origin.

Heatmap of the AUCs for the top marker genes in cluster 6 compared to all other clusters in the Grun HSC dataset.

Figure 35.5: Heatmap of the AUCs for the top marker genes in cluster 6 compared to all other clusters in the Grun HSC dataset.

Session Info

R version 4.0.3 (2020-10-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 20.04.1 LTS

Matrix products: default
BLAS:   /home/biocbuild/bbs-3.12-bioc/R/lib/libRblas.so
LAPACK: /home/biocbuild/bbs-3.12-bioc/R/lib/libRlapack.so

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_US.UTF-8        LC_COLLATE=C              
 [5] LC_MONETARY=en_US.UTF-8    LC_MESSAGES=en_US.UTF-8   
 [7] LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               LC_TELEPHONE=C            
[11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C       

attached base packages:
[1] parallel  stats4    stats     graphics  grDevices utils     datasets 
[8] methods   base     

other attached packages:
 [1] pheatmap_1.0.12             scran_1.18.0               
 [3] scater_1.18.0               ggplot2_3.3.2              
 [5] scuttle_1.0.0               AnnotationHub_2.22.0       
 [7] BiocFileCache_1.14.0        dbplyr_1.4.4               
 [9] ensembldb_2.14.0            AnnotationFilter_1.14.0    
[11] GenomicFeatures_1.42.0      AnnotationDbi_1.52.0       
[13] scRNAseq_2.4.0              SingleCellExperiment_1.12.0
[15] SummarizedExperiment_1.20.0 Biobase_2.50.0             
[17] GenomicRanges_1.42.0        GenomeInfoDb_1.26.0        
[19] IRanges_2.24.0              S4Vectors_0.28.0           
[21] BiocGenerics_0.36.0         MatrixGenerics_1.2.0       
[23] matrixStats_0.57.0          BiocStyle_2.18.0           
[25] rebook_1.0.0               

loaded via a namespace (and not attached):
  [1] Rtsne_0.15                    ggbeeswarm_0.6.0             
  [3] colorspace_1.4-1              ellipsis_0.3.1               
  [5] bluster_1.0.0                 XVector_0.30.0               
  [7] BiocNeighbors_1.8.0           farver_2.0.3                 
  [9] bit64_4.0.5                   interactiveDisplayBase_1.28.0
 [11] xml2_1.3.2                    codetools_0.2-16             
 [13] sparseMatrixStats_1.2.0       knitr_1.30                   
 [15] Rsamtools_2.6.0               graph_1.68.0                 
 [17] shiny_1.5.0                   BiocManager_1.30.10          
 [19] compiler_4.0.3                httr_1.4.2                   
 [21] dqrng_0.2.1                   assertthat_0.2.1             
 [23] Matrix_1.2-18                 fastmap_1.0.1                
 [25] lazyeval_0.2.2                limma_3.46.0                 
 [27] later_1.1.0.1                 BiocSingular_1.6.0           
 [29] htmltools_0.5.0               prettyunits_1.1.1            
 [31] tools_4.0.3                   igraph_1.2.6                 
 [33] rsvd_1.0.3                    gtable_0.3.0                 
 [35] glue_1.4.2                    GenomeInfoDbData_1.2.4       
 [37] dplyr_1.0.2                   rappdirs_0.3.1               
 [39] Rcpp_1.0.5                    vctrs_0.3.4                  
 [41] Biostrings_2.58.0             ExperimentHub_1.16.0         
 [43] rtracklayer_1.50.0            DelayedMatrixStats_1.12.0    
 [45] xfun_0.19                     stringr_1.4.0                
 [47] ps_1.4.0                      beachmat_2.6.0               
 [49] mime_0.9                      lifecycle_0.2.0              
 [51] irlba_2.3.3                   statmod_1.4.35               
 [53] XML_3.99-0.5                  edgeR_3.32.0                 
 [55] zlibbioc_1.36.0               scales_1.1.1                 
 [57] hms_0.5.3                     promises_1.1.1               
 [59] ProtGenerics_1.22.0           RColorBrewer_1.1-2           
 [61] yaml_2.2.1                    curl_4.3                     
 [63] gridExtra_2.3                 memoise_1.1.0                
 [65] biomaRt_2.46.0                stringi_1.5.3                
 [67] RSQLite_2.2.1                 highr_0.8                    
 [69] BiocVersion_3.12.0            BiocParallel_1.24.0          
 [71] rlang_0.4.8                   pkgconfig_2.0.3              
 [73] bitops_1.0-6                  evaluate_0.14                
 [75] lattice_0.20-41               purrr_0.3.4                  
 [77] labeling_0.4.2                GenomicAlignments_1.26.0     
 [79] CodeDepends_0.6.5             cowplot_1.1.0                
 [81] bit_4.0.4                     processx_3.4.4               
 [83] tidyselect_1.1.0              magrittr_1.5                 
 [85] bookdown_0.21                 R6_2.5.0                     
 [87] generics_0.1.0                DelayedArray_0.16.0          
 [89] DBI_1.1.0                     pillar_1.4.6                 
 [91] withr_2.3.0                   RCurl_1.98-1.2               
 [93] tibble_3.0.4                  crayon_1.3.4                 
 [95] rmarkdown_2.5                 viridis_0.5.1                
 [97] progress_1.2.2                locfit_1.5-9.4               
 [99] grid_4.0.3                    blob_1.2.1                   
[101] callr_3.5.1                   digest_0.6.27                
[103] xtable_1.8-4                  httpuv_1.5.4                 
[105] openssl_1.4.3                 munsell_0.5.0                
[107] viridisLite_0.3.0             beeswarm_0.2.3               
[109] vipor_0.4.5                   askpass_1.1                  

Bibliography

Grun, D., M. J. Muraro, J. C. Boisset, K. Wiebrands, A. Lyubimova, G. Dharmadhikari, M. van den Born, et al. 2016. “De Novo Prediction of Stem Cell Identity using Single-Cell Transcriptome Data.” Cell Stem Cell 19 (2): 266–77.