Visualizing Large-scale Copy Number Variation in Single-Cell RNA-Seq Expression Data
2019-09-16
Abstract
InferCNV is used to explore tumor single cell RNA-Seq data to identify evidence for large-scale chromosomal copy number variations, such as gains or deletions of entire chromosomes or large segments of chromosomes. This is done by exploring expression intensity of genes across positions of the genome in comparison to the average or a set of reference ‘normal’ cells. A heatmap is generated illustrating the relative expression intensities across each chromosome, and it becomes readily apparent as to which regions of the genome are over-abundant or less-abundant as compared to normal cells (or the average, if reference normal cells are not provided).Package
infercnv 1.0.4
Contents
1 Installation
1.1 Required dependencies
inferCNV uses the R packages ape, BiocGenerics, binhf, caTools, coda, coin, dplyr, doparallel, edgeR, fastcluster, fitdistrplus, foreach, futile.logger, future, gplots, ggplot2, HiddenMarkov, reshape, rjags, RColorBrewer, SingleCellExperiment, SummarizedExperiment and imports functions from the archived GMD.
1.2 Installing
if (!requireNamespace("BiocManager", quietly = TRUE))
install.packages("BiocManager")
BiocManager::install("infercnv")1.3 Optional extension
If you want to use the interactive heatmap visualization, please check the add-on packge R inferCNV_NGCHM after installing the packages tibble, tsvio and NGCHMR. To install optional packages, type the following in an R command window:
install.packages("tibble")
install.packages("devtools")
devtools::install_github("bmbroom/tsvio")
devtools::install_github("bmbroom/NGCHMR", ref="stable")
devtools::install_github("broadinstitute/inferCNV_NGCHM")And download the NGCHM java application by typing the following in a regular shell:
wget http://tcga.ngchm.net/NGCHM/ShaidyMapGen.jar2 Running InferCNV
2.1 Create the InferCNV Object
Reading in the raw counts matrix and meta data, populating the infercnv object
infercnv_obj = CreateInfercnvObject(
raw_counts_matrix="../inst/extdata/oligodendroglioma_expression_downsampled.counts.matrix.gz",
annotations_file="../inst/extdata/oligodendroglioma_annotations_downsampled.txt",
delim="\t",
gene_order_file="../inst/extdata/gencode_downsampled.EXAMPLE_ONLY_DONT_REUSE.txt",
ref_group_names=c("Microglia/Macrophage","Oligodendrocytes (non-malignant)"))## INFO [2019-09-16 23:20:27] Parsing matrix: ../inst/extdata/oligodendroglioma_expression_downsampled.counts.matrix.gz
## INFO [2019-09-16 23:20:29] Parsing gene order file: ../inst/extdata/gencode_downsampled.EXAMPLE_ONLY_DONT_REUSE.txt
## INFO [2019-09-16 23:20:29] Parsing cell annotations file: ../inst/extdata/oligodendroglioma_annotations_downsampled.txt
## INFO [2019-09-16 23:20:30] ::order_reduce:Start.
## INFO [2019-09-16 23:20:30] .order_reduce(): expr and order match.
## INFO [2019-09-16 23:20:30] ::process_data:order_reduce:Reduction from positional data, new dimensions (r,c) = 10338,184 Total=18322440.6799817 Min=0 Max=34215.
## INFO [2019-09-16 23:20:30] num genes removed taking into account provided gene ordering list: 399 = 3.8595473012188% removed.
## INFO [2019-09-16 23:20:30] validating infercnv_obj2.2 Running the full default analysis
out_dir = tempfile()
infercnv_obj_default = infercnv::run(
infercnv_obj,
cutoff=1, # cutoff=1 works well for Smart-seq2, and cutoff=0.1 works well for 10x Genomics
out_dir=out_dir,
cluster_by_groups=TRUE,
plot_steps=FALSE,
denoise=TRUE,
HMM=FALSE,
no_prelim_plot=TRUE,
png_res=60
)Basic ouput from running inferCNV.
3 Additional Information
3.1 Online Documentation
For additional explanations on files, usage, and a tutorial please visit the wiki.
3.2 TrinityCTAT
This tool is a part of the TrinityCTAT toolkit focused on leveraging the use of RNA-Seq to better understand cancer transcriptomes. To find out more please visit TrinityCTAT
3.3 Applications
This methodology was used in:
4 Session info
## R version 3.6.1 (2019-07-05)
## Platform: x86_64-pc-linux-gnu (64-bit)
## Running under: Ubuntu 18.04.3 LTS
##
## Matrix products: default
## BLAS: /home/biocbuild/bbs-3.9-bioc/R/lib/libRblas.so
## LAPACK: /home/biocbuild/bbs-3.9-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] stats graphics grDevices utils datasets methods base
##
## other attached packages:
## [1] infercnv_1.0.4 BiocStyle_2.12.0
##
## loaded via a namespace (and not attached):
## [1] nlme_3.1-141 bitops_1.0-6
## [3] matrixStats_0.55.0 doParallel_1.0.15
## [5] RColorBrewer_1.1-2 GenomeInfoDb_1.20.0
## [7] backports_1.1.4 tools_3.6.1
## [9] R6_2.4.0 KernSmooth_2.23-15
## [11] lazyeval_0.2.2 BiocGenerics_0.30.0
## [13] colorspace_1.4-1 npsurv_0.4-0
## [15] gridExtra_2.3 tidyselect_0.2.5
## [17] compiler_3.6.1 argparse_2.0.1
## [19] Biobase_2.44.0 formatR_1.7
## [21] DelayedArray_0.10.0 sandwich_2.5-1
## [23] bookdown_0.13 caTools_1.17.1.2
## [25] scales_1.0.0 mvtnorm_1.0-11
## [27] stringr_1.4.0 digest_0.6.20
## [29] rmarkdown_1.15 XVector_0.24.0
## [31] pkgconfig_2.0.2 htmltools_0.3.6
## [33] limma_3.40.6 rlang_0.4.0
## [35] zoo_1.8-6 jsonlite_1.6
## [37] BiocParallel_1.18.1 gtools_3.8.1
## [39] dplyr_0.8.3 RCurl_1.95-4.12
## [41] magrittr_1.5 modeltools_0.2-22
## [43] GenomeInfoDbData_1.2.1 futile.logger_1.4.3
## [45] Matrix_1.2-17 Rcpp_1.0.2
## [47] munsell_0.5.0 S4Vectors_0.22.1
## [49] ape_5.3 lifecycle_0.1.0
## [51] stringi_1.4.3 multcomp_1.4-10
## [53] yaml_2.2.0 edgeR_3.26.8
## [55] MASS_7.3-51.4 SummarizedExperiment_1.14.1
## [57] zlibbioc_1.30.0 plyr_1.8.4
## [59] gplots_3.0.1.1 grid_3.6.1
## [61] parallel_3.6.1 gdata_2.18.0
## [63] listenv_0.7.0 crayon_1.3.4
## [65] lattice_0.20-38 splines_3.6.1
## [67] findpython_1.0.5 locfit_1.5-9.1
## [69] zeallot_0.1.0 knitr_1.24
## [71] pillar_1.4.2 fastcluster_1.1.25
## [73] GenomicRanges_1.36.1 codetools_0.2-16
## [75] stats4_3.6.1 futile.options_1.0.1
## [77] glue_1.3.1 evaluate_0.14
## [79] lsei_1.2-0 lambda.r_1.2.3
## [81] BiocManager_1.30.4 vctrs_0.2.0
## [83] foreach_1.4.7 tidyr_1.0.0
## [85] gtable_0.3.0 purrr_0.3.2
## [87] reshape_0.8.8 future_1.14.0
## [89] assertthat_0.2.1 ggplot2_3.2.1
## [91] xfun_0.9 coin_1.3-1
## [93] libcoin_1.0-5 coda_0.19-3
## [95] rjags_4-9 survival_2.44-1.1
## [97] SingleCellExperiment_1.6.0 tibble_2.1.3
## [99] iterators_1.0.12 IRanges_2.18.2
## [101] globals_0.12.4 fitdistrplus_1.0-14
## [103] TH.data_1.0-10