Chapter 7 Lawlor human pancreas (SMARTer)

7.1 Introduction

This performs an analysis of the Lawlor et al. (2017) dataset, consisting of human pancreas cells from various donors.

7.3 Quality control

Distribution of each QC metric across cells from each donor of the Lawlor pancreas dataset. Each point represents a cell and is colored according to whether that cell was discarded.

Figure 7.1: Distribution of each QC metric across cells from each donor of the Lawlor pancreas dataset. Each point represents a cell and is colored according to whether that cell was discarded.

Percentage of mitochondrial reads in each cell in the 416B dataset compared to the total count. Each point represents a cell and is colored according to whether that cell was discarded.

Figure 7.2: Percentage of mitochondrial reads in each cell in the 416B dataset compared to the total count. Each point represents a cell and is colored according to whether that cell was discarded.

##              low_lib_size            low_n_features high_subsets_Mito_percent 
##                         9                         5                        25 
##                   discard 
##                        34

7.4 Normalization

##    Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
##   0.295   0.781   0.963   1.000   1.182   2.629
Relationship between the library size factors and the deconvolution size factors in the Lawlor pancreas dataset.

Figure 7.3: Relationship between the library size factors and the deconvolution size factors in the Lawlor pancreas dataset.

7.7 Clustering

##    
##     Acinar Alpha Beta Delta Ductal Gamma/PP None/Other Stellate
##   1      1     0    0    13      2       16          2        0
##   2      0     1   76     1      0        0          0        0
##   3      0   161    1     0      0        1          2        0
##   4      0     1    0     1      0        0          5       19
##   5      0     0  175     4      1        0          1        0
##   6     22     0    0     0      0        0          0        0
##   7      0    75    0     0      0        0          0        0
##   8      0     0    0     1     20        0          2        0
##    
##     ACCG268 ACCR015A ACEK420A ACEL337 ACHY057 ACIB065 ACIW009 ACJV399
##   1       8        2        2       4       4       4       9       1
##   2      14        3        2      33       3       2       4      17
##   3      36       23       14      13      14      14      21      30
##   4       7        1        0       1       0       4       9       4
##   5      34       10        4      39       7      23      24      40
##   6       0        2       13       0       0       0       5       2
##   7      32       12        0       5       6       7       4       9
##   8       1        1        2       1       2       1      12       3
Obligatory $t$-SNE plots of the Lawlor pancreas dataset. Each point represents a cell that is colored by cluster (left) or batch (right).

Figure 5.3: Obligatory \(t\)-SNE plots of the Lawlor pancreas dataset. Each point represents a cell that is colored by cluster (left) or batch (right).

Session Info

R version 4.1.1 (2021-08-10)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 20.04.3 LTS

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

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C              
 [3] LC_TIME=en_GB              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] stats4    stats     graphics  grDevices utils     datasets  methods  
[8] base     

other attached packages:
 [1] BiocSingular_1.10.0         scran_1.22.0               
 [3] scater_1.22.0               ggplot2_3.3.5              
 [5] scuttle_1.4.0               ensembldb_2.18.0           
 [7] AnnotationFilter_1.18.0     GenomicFeatures_1.46.0     
 [9] AnnotationDbi_1.56.0        AnnotationHub_3.2.0        
[11] BiocFileCache_2.2.0         dbplyr_2.1.1               
[13] scRNAseq_2.7.2              SingleCellExperiment_1.16.0
[15] SummarizedExperiment_1.24.0 Biobase_2.54.0             
[17] GenomicRanges_1.46.0        GenomeInfoDb_1.30.0        
[19] IRanges_2.28.0              S4Vectors_0.32.0           
[21] BiocGenerics_0.40.0         MatrixGenerics_1.6.0       
[23] matrixStats_0.61.0          BiocStyle_2.22.0           
[25] rebook_1.4.0               

loaded via a namespace (and not attached):
  [1] igraph_1.2.7                  lazyeval_0.2.2               
  [3] BiocParallel_1.28.0           digest_0.6.28                
  [5] htmltools_0.5.2               viridis_0.6.2                
  [7] fansi_0.5.0                   magrittr_2.0.1               
  [9] memoise_2.0.0                 ScaledMatrix_1.2.0           
 [11] cluster_2.1.2                 limma_3.50.0                 
 [13] Biostrings_2.62.0             prettyunits_1.1.1            
 [15] colorspace_2.0-2              blob_1.2.2                   
 [17] rappdirs_0.3.3                ggrepel_0.9.1                
 [19] xfun_0.27                     dplyr_1.0.7                  
 [21] crayon_1.4.1                  RCurl_1.98-1.5               
 [23] jsonlite_1.7.2                graph_1.72.0                 
 [25] glue_1.4.2                    gtable_0.3.0                 
 [27] zlibbioc_1.40.0               XVector_0.34.0               
 [29] DelayedArray_0.20.0           scales_1.1.1                 
 [31] edgeR_3.36.0                  DBI_1.1.1                    
 [33] Rcpp_1.0.7                    viridisLite_0.4.0            
 [35] xtable_1.8-4                  progress_1.2.2               
 [37] dqrng_0.3.0                   bit_4.0.4                    
 [39] rsvd_1.0.5                    metapod_1.2.0                
 [41] httr_1.4.2                    dir.expiry_1.2.0             
 [43] ellipsis_0.3.2                pkgconfig_2.0.3              
 [45] XML_3.99-0.8                  farver_2.1.0                 
 [47] CodeDepends_0.6.5             sass_0.4.0                   
 [49] locfit_1.5-9.4                utf8_1.2.2                   
 [51] tidyselect_1.1.1              labeling_0.4.2               
 [53] rlang_0.4.12                  later_1.3.0                  
 [55] munsell_0.5.0                 BiocVersion_3.14.0           
 [57] tools_4.1.1                   cachem_1.0.6                 
 [59] generics_0.1.1                RSQLite_2.2.8                
 [61] ExperimentHub_2.2.0           evaluate_0.14                
 [63] stringr_1.4.0                 fastmap_1.1.0                
 [65] yaml_2.2.1                    knitr_1.36                   
 [67] bit64_4.0.5                   purrr_0.3.4                  
 [69] KEGGREST_1.34.0               sparseMatrixStats_1.6.0      
 [71] mime_0.12                     xml2_1.3.2                   
 [73] biomaRt_2.50.0                compiler_4.1.1               
 [75] beeswarm_0.4.0                filelock_1.0.2               
 [77] curl_4.3.2                    png_0.1-7                    
 [79] interactiveDisplayBase_1.32.0 statmod_1.4.36               
 [81] tibble_3.1.5                  bslib_0.3.1                  
 [83] stringi_1.7.5                 highr_0.9                    
 [85] bluster_1.4.0                 lattice_0.20-45              
 [87] ProtGenerics_1.26.0           Matrix_1.3-4                 
 [89] vctrs_0.3.8                   pillar_1.6.4                 
 [91] lifecycle_1.0.1               BiocManager_1.30.16          
 [93] jquerylib_0.1.4               BiocNeighbors_1.12.0         
 [95] cowplot_1.1.1                 bitops_1.0-7                 
 [97] irlba_2.3.3                   httpuv_1.6.3                 
 [99] rtracklayer_1.54.0            R6_2.5.1                     
[101] BiocIO_1.4.0                  bookdown_0.24                
[103] promises_1.2.0.1              gridExtra_2.3                
[105] vipor_0.4.5                   codetools_0.2-18             
[107] assertthat_0.2.1              rjson_0.2.20                 
[109] withr_2.4.2                   GenomicAlignments_1.30.0     
[111] Rsamtools_2.10.0              GenomeInfoDbData_1.2.7       
[113] parallel_4.1.1                hms_1.1.1                    
[115] grid_4.1.1                    beachmat_2.10.0              
[117] rmarkdown_2.11                DelayedMatrixStats_1.16.0    
[119] Rtsne_0.15                    shiny_1.7.1                  
[121] ggbeeswarm_0.6.0              restfulr_0.0.13              

References

Lawlor, N., J. George, M. Bolisetty, R. Kursawe, L. Sun, V. Sivakamasundari, I. Kycia, P. Robson, and M. L. Stitzel. 2017. “Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes.” Genome Res. 27 (2): 208–22.