Movatterモバイル変換

Title:	Extract, Analyze and Visualize Mutational Signatures for GenomicVariations
Version:	2.3.1
Description:	Genomic alterations including single nucleotide substitution, copy number alteration, etc. are the major force for cancer initialization and development. Due to the specificity of molecular lesions caused by genomic alterations, we can generate characteristic alteration spectra, called 'signature' (Wang, Shixiang, et al. (2021) <doi:10.1371/journal.pgen.1009557> & Alexandrov, Ludmil B., et al. (2020) <doi:10.1038/s41586-020-1943-3> & Steele Christopher D., et al. (2022) <doi:10.1038/s41586-022-04738-6>). This package helps users to extract, analyze and visualize signatures from genomic alteration records, thus providing new insight into cancer study.
License:	MIT + file LICENSE
URL:	https://github.com/ShixiangWang/sigminer,https://shixiangwang.github.io/sigminer/,https://shixiangwang.github.io/sigminer-book/
BugReports:	https://github.com/ShixiangWang/sigminer/issues
Depends:	R (≥ 3.5)
Imports:	cli (≥ 2.0.0), cowplot, data.table, dplyr, furrr (≥ 0.2.0),future, ggplot2 (≥ 3.3.0), ggpubr, maftools, magrittr,methods, NMF, purrr, Rcpp, rlang (≥ 0.1.2), stats, tidyr
Suggests:	Biobase, Biostrings, BSgenome, BSgenome.Hsapiens.UCSC.hg19,circlize, cluster, covr, digest, GenomicRanges, GenSA,ggalluvial, ggcorrplot, ggfittext, ggplotify, ggrepel, IRanges,knitr, lpSolve, markdown, matrixStats, nnls, parallel,patchwork, pheatmap, quadprog, R.utils, RColorBrewer,reticulate, rmarkdown, roxygen2, scales, synchronicity,testthat (≥ 3.0.0), tibble, UCSCXenaTools
LinkingTo:	Rcpp
VignetteBuilder:	knitr
Encoding:	UTF-8
LazyData:	true
RoxygenNote:	7.3.1
Config/testthat/edition:	3
NeedsCompilation:	yes
Packaged:	2024-05-11 04:07:34 UTC; wsx
Author:	Shixiang Wang [aut, cre], Ziyu Tao [aut], Huimin Li [aut], Tao Wu [aut], Xue-Song Liu [aut, ctb], Anand Mayakonda [ctb]
Maintainer:	Shixiang Wang <w_shixiang@163.com>
Repository:	CRAN
Date/Publication:	2024-05-11 08:50:02 UTC

sigminer: Extract, Analyze and Visualize Signatures for Genomic Variations

Description

• Author:Shixiang Wang (w_shixiang@163.com)

• Please go tohttps://shixiangwang.github.io/sigminer-doc/ for full vignette.

• Please go tohttps://shixiangwang.github.io/sigminer/reference/index.htmlfor organized documentation of functions and datasets.

• Result visualization forMAF is provide bymaftools package,please read itsvignette.

Author(s)

Maintainer: Shixiang Wangw_shixiang@163.com (ORCID)

Authors:

• Ziyu Taotaozy@shanghaitech.edu.cn (ORCID)

• Huimin Lilihm@shanghaitech.edu.cn (ORCID)

• Tao Wuwutao2@shanghaitech.edu.cn (ORCID)

• Xue-Song Liu (ORCID) [contributor]

Other contributors:

• Anand Mayakonda [contributor]

See Also

Useful links:

• https://github.com/ShixiangWang/sigminer

• https://shixiangwang.github.io/sigminer/

• https://shixiangwang.github.io/sigminer-book/

• Report bugs athttps://github.com/ShixiangWang/sigminer/issues

Pipe operator

Description

Seemagrittr::%>% for details.

Usage

lhs %>% rhs

Classification Table of Copy Number Features Devised by Wang et al. for Method 'W'

Description

Classification Table of Copy Number Features Devised by Wang et al. for Method 'W'

Format

Adata.table with "sigminer.features" class name

Source

Generate from code under data_raw/

Examples

data(CN.features)

Class CopyNumber

Description

S4 class for storing summarized absolute copy number profile.

Slots

data

data.table of absolute copy number calling.

summary.per.sample

data.table of copy number variation summary per sample.

genome_build

genome build version, should be one of 'hg19' or 'hg38'.

genome_measure

Set 'called' will use autosomo called segments size to compute total sizefor CNA burden calculation, this option is useful for WES and target sequencing.Set 'wg' will autosome size from genome build, this option is useful for WGS, SNP etc..

annotation

data.table of annotation for copy number segments.

dropoff.segs

data.table of copy number segments dropped from raw input.

Class MAF

Description

S4 class for storing summarized MAF. It is frommaftools package.

Details

More about MAF object please seemaftools.

Slots

data

data.table of MAF file containing all non-synonymous variants.

variants.per.sample

table containing variants per sample

variant.type.summary

table containing variant types per sample

variant.classification.summary

table containing variant classification per sample

gene.summary

table containing variant classification per gene

summary

table with basic MAF summary stats

maf.silent

subset of main MAF containing only silent variants

clinical.data

clinical data associated with each sample/Tumor_Sample_Barcode in MAF.

Add Horizontal Arrow with Text Label to a ggplot

Description

Add Horizontal Arrow with Text Label to a ggplot

Usage

add_h_arrow(  p,  x,  y,  label = "optimal number",  space = 0.01,  vjust = 0.3,  seg_len = 0.1,  arrow_len = unit(2, "mm"),  arrow_type = c("closed", "open"),  font_size = 5,  font_family = c("serif", "sans", "mono"),  font_face = c("plain", "bold", "italic"))

Arguments

Value

aggplot object.

Add Text Labels to a ggplot

Description

Add text labels to a ggplot object, such as the resultfromshow_sig_profile.

Usage

add_labels(  p,  x,  y,  y_end = NULL,  n_label = NULL,  labels = NULL,  revert_order = FALSE,  font_size = 5,  font_family = "serif",  font_face = c("plain", "bold", "italic"),  ...)

Arguments

Value

aggplot object.

Examples

# Load mutational signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Show signature profilep <- show_sig_profile(sig2, mode = "SBS")# Method 1p1 <- add_labels(p,  x = 0.75, y = 0.3, y_end = 0.9, n_label = 3,  labels = paste0("text", 1:3))p1# Method 2p2 <- add_labels(p,  x = c(0.15, 0.6, 0.75), y = c(0.3, 0.6, 0.9),  labels = paste0("text", 1:3))p2# Method 3sim <- get_sig_similarity(sig2)p3 <- add_labels(p,  x = c(0.15, 0.6, 0.75), y = c(0.25, 0.55, 0.8),  labels = sim, font_size = 2)p3

A Best Practice for Signature Extraction and Exposure (Activity) Attribution

Description

These functions are combined to provide a best practice for optimallyidentifying mutational signatures and attributing their activities (exposures)in tumor samples. They are listed in order to use.

• bp_extract_signatures() for extracting signatures.

• bp_show_survey() for showing measures change under differentsignature numbers to help user select optimal signature number.At default, an aggregated score (named score) is generated tosuggest the best solution.

• bp_show_survey2() for showing simplified signature number survey likeshow_sig_number_survey().

• bp_get_sig_obj() for get a (list of)Signature object which is commonused insigminer for analysis and visualization.

• bp_attribute_activity() for optimizing signature activities (exposures).NOTE: the activities from extraction step may be better!You can also usesig_extract to get optimal NMF result from multiple NMF runs.Besides, you can usesig_fit to quantify exposures based on signatures extractedfrombp_extract_signatures().

• bp_extract_signatures_iter() for extracting signature in a iteration way.

• bp_cluster_iter_list() for clustering (hclust with average linkage)iterated signatures to help collapsemultiple signatures into one. The result cluster can be visualized byplot() orfactoextra::fviz_dend().

• bp_get_clustered_sigs() for getting clustered (grouped) mean signatures from signature clusters.

• Extra:bp_get_stats() for obtaining stats for signatures and samples of a solution.These stats are aggregated (averaged) as the stats for a solution(specific signature number).

• Extra:bp_get_rank_score() for obtaining rank score for all signature numbers.

Usage

bp_extract_signatures(  nmf_matrix,  range = 2:5,  n_bootstrap = 20L,  n_nmf_run = 50,  RTOL = 0.001,  min_contribution = 0,  cores = min(4L, future::availableCores()),  cores_solution = min(cores, length(range)),  seed = 123456L,  handle_hyper_mutation = TRUE,  report_integer_exposure = FALSE,  only_core_stats = nrow(nmf_matrix) > 100,  cache_dir = file.path(tempdir(), "sigminer_bp"),  keep_cache = FALSE,  pynmf = FALSE,  use_conda = TRUE,  py_path = "/Users/wsx/anaconda3/bin/python")bp_extract_signatures_iter(  nmf_matrix,  range = 2:5,  sim_threshold = 0.95,  max_iter = 10L,  n_bootstrap = 20L,  n_nmf_run = 50,  RTOL = 0.001,  min_contribution = 0,  cores = min(4L, future::availableCores()),  cores_solution = min(cores, length(range)),  seed = 123456L,  handle_hyper_mutation = TRUE,  report_integer_exposure = FALSE,  only_core_stats = nrow(nmf_matrix) > 100,  cache_dir = file.path(tempdir(), "sigminer_bp"),  keep_cache = FALSE,  pynmf = FALSE,  use_conda = FALSE,  py_path = "/Users/wsx/anaconda3/bin/python")bp_cluster_iter_list(x, k = NULL, include_final_iteration = TRUE)bp_get_clustered_sigs(SigClusters, cluster_label)bp_get_sig_obj(obj, signum = NULL)bp_get_stats(obj)bp_get_rank_score(obj)bp_show_survey2(  obj,  x = "signature_number",  left_y = "silhouette",  right_y = "L2_error",  left_name = left_y,  right_name = right_y,  left_color = "black",  right_color = "red",  left_shape = 16,  right_shape = 18,  shape_size = 4,  highlight = NULL)bp_show_survey(  obj,  add_score = FALSE,  scales = c("free_y", "free"),  fixed_ratio = TRUE)bp_attribute_activity(  input,  sample_class = NULL,  nmf_matrix = NULL,  method = c("bt", "stepwise"),  bt_use_prop = FALSE,  return_class = c("matrix", "data.table"),  use_parallel = FALSE,  cache_dir = file.path(tempdir(), "sigminer_attribute_activity"),  keep_cache = FALSE)

Arguments

Details

The signature extraction approach is adopted from reference #1, #2, andthe whole best practice is adopted from the pipeline used by reference #3.I implement the whole procedure with R code based on the method descriptionof papers. The code is well organized, tested and documented so user willfind it pretty simple and useful. Besides, the structure of the results isvery clear to see and also visualize like other approaches provided bysigminer.

Value

It depends on the called function.

Measure Explanation in Survey Plot

The survey plot provides a pretty good way to facilitate the signature numberselection. Ascore measure is calculated as the weighted mean of selectedmeasures and visualized as the first sub-plot. The optimal number is highlightedwith red color dot and the best values for each measures are alsohighlighted with orange color dots. The detail of 6 measures shown in plot areexplained as below.

• score - an aggregated score based on rank scores from selected measures below.The higher, the better. When two signature numbers have the same score,the larger signature number is preferred (this is a rare situation, youhave to double check other measures).

• silhouette - the average silhouette width for signatures, also named as ASW in reference #2.The signature number with silhouette decreases sharply is preferred.

• distance - the average sample reconstructed cosine distance, the lower value is better.

• error - the average sample reconstructed error calculated with L2 formula(i.e. L2 error). This lower value is better. This measure represents asimilar concept likedistance above, they are all used to quantify how wellsample mutation profiles can be reconstructed from signatures, butdistancecares the whole mutation profile similarity whileerror here cares value difference.

• ⁠pos cor⁠ - the average positive signature exposure correlation coefficient.The lower value is better. This measure is constructed based on my understandingabout signatures: mutational signatures are typically treated as independentrecurrent patterns, so their activities are less correlated.

• similarity - the average similarity within in a signature cluster.Likesilhouette, the point decreases sharply is preferred.In the practice, results from multiple NMF runs are clusteredwith "clustering with match" algorithm proposed by reference #2. This valueindicates if the signature profiles extracted from different NMF runs are similar.

Author(s)

Shixiang Wangw_shixiang@163.com

References

Alexandrov, Ludmil B., et al. "Deciphering signatures of mutational processes operative in human cancer." Cell reports 3.1 (2013): 246-259.

Degasperi, Andrea, et al. "A practical framework and online tool for mutational signature analyses show intertissue variation and driver dependencies." Nature cancer 1.2 (2020): 249-263.

Alexandrov, Ludmil B., et al. “The repertoire of mutational signatures in human cancer.” Nature 578.7793 (2020): 94-101.

See Also

Seesig_estimate,sig_extract,sig_auto_extract,sigprofiler_extract for other approaches.

Examples

data("simulated_catalogs")# Here I reduce the values for n_bootstrap and n_nmf_run# for reducing the run time.# In practice, you should keep default or increase the values# for better estimation.## The input data here is simulated from 10 mutational signatures# e1 <- bp_extract_signatures(#   t(simulated_catalogs$set1),#   range = 8:12,#   n_bootstrap = 5,#   n_nmf_run = 10# )## To avoid computation in examples,# Here just load the result# (e1$signature and e1$exposure set to NA to reduce package size)load(system.file("extdata", "e1.RData", package = "sigminer"))# See the survey for different signature numbers# The suggested solution is marked as red dot# with highest integrated score.p1 <- bp_show_survey(e1)p1# You can also exclude plotting and highlighting the scorep2 <- bp_show_survey(e1, add_score = FALSE)p2# You can also plot a simplified versionp3 <- bp_show_survey2(e1, highlight = 10)p3# Obtain the suggested solution from extraction resultobj_suggested <- bp_get_sig_obj(e1, e1$suggested)obj_suggested# If you think the suggested signature number is not right# Just pick up the solution you wantobj_s8 <- bp_get_sig_obj(e1, 8)# Track the reconstructed profile similarityrec_sim <- get_sig_rec_similarity(obj_s8, t(simulated_catalogs$set1))rec_sim# After extraction, you can assign the signatures# to reference COSMIC signatures# More see ?get_sig_similaritysim <- get_sig_similarity(obj_suggested)# Visualize the match resultif (require(pheatmap)) {  pheatmap::pheatmap(sim$similarity)}# You already got the activities of signatures# in obj_suggested, however, you can still# try to optimize the result.# NOTE: the optimization step may not truly optimize the result!expo <- bp_attribute_activity(e1, return_class = "data.table")expo$abs_activity## Not run: # Iterative extraction:# This procedure will rerun extraction step# for those samples with reconstructed catalog similarity# lower than a threshold (default is 0.95)e2 <- bp_extract_signatures_iter(  t(simulated_catalogs$set1),  range = 9:11,  n_bootstrap = 5,  n_nmf_run = 5,  sim_threshold = 0.99)e2# When the procedure run multiple rounds# you can cluster the signatures from different rounds by# the following command# bp_cluster_iter_list(e2)## Extra utilitiesrank_score <- bp_get_rank_score(e1)rank_scorestats <- bp_get_stats(e2$iter1)# Get the mean reconstructed similarity1 - stats$stats_sample$cosine_distance_mean## End(Not run)

Location of Centromeres at Genome Build T2T

Description

Location of Centromeres at Genome Build T2T

Format

A data.frame

Source

from T2T study

Examples

data(centromeres.T2T)

Location of Centromeres at Genome Build hg19

Description

Location of Centromeres at Genome Build hg19

Format

A data.frame

Source

Generate from UCSC gold path

Examples

data(centromeres.hg19)

Location of Centromeres at Genome Build hg38

Description

Location of Centromeres at Genome Build hg38

Format

A data.frame

Source

Generate from Genome Reference Consortium

Examples

data(centromeres.hg38)

Location of Centromeres at Genome Build mm10

Description

Location of Centromeres at Genome Build mm10

Format

A data.frame

Source

Generate fromhttps://hgdownload.soe.ucsc.edu/goldenPath/mm10/database/gap.txt.gz

Examples

data(centromeres.mm10)

Location of Centromeres at Genome Build mm9

Description

Location of Centromeres at Genome Build mm9

Format

A data.frame

Source

Generate fromhttps://hgdownload.soe.ucsc.edu/goldenPath/mm9/database/with code:

for i in $(seq 1 19) X Y;dowget https://hgdownload.soe.ucsc.edu/goldenPath/mm9/database/chr${i}_gap.txt.gzdone

Examples

data(centromeres.mm9)

Chromosome Size of Genome Build T2T

Description

Chromosome Size of Genome Build T2T

Format

A data.frame

Source

from T2T study

Examples

data(chromsize.T2T)

Chromosome Size of Genome Build hg19

Description

Chromosome Size of Genome Build hg19

Format

A data.frame

Source

Generate from UCSC gold path

Examples

data(chromsize.hg19)

Chromosome Size of Genome Build hg38

Description

Chromosome Size of Genome Build hg38

Format

A data.frame

Source

Generate from UCSC gold path

Examples

data(chromsize.hg38)

Chromosome Size of Genome Build mm10

Description

Chromosome Size of Genome Build mm10

Format

A data.frame

Source

Generate from UCSC gold pathhttp://hgdownload.cse.ucsc.edu/goldenPath/mm10/bigZips/mm10.chrom.sizes

Examples

data(chromsize.mm10)

Chromosome Size of Genome Build mm9

Description

Chromosome Size of Genome Build mm9

Format

A data.frame

Source

Generate from UCSC gold pathhttp://hgdownload.cse.ucsc.edu/goldenPath/mm9/bigZips/mm9.chrom.sizes

Examples

data(chromsize.mm9)

Calculate Cosine Measures

Description

Calculate Cosine Measures

Usage

cosine(x, y)

Arguments

Value

a numeric value ormatrix.

Examples

x <- c(1, 1, 1, 0, 0, 0, 0, 0, 0, 0, 0, 0)y <- c(0, 0, 1, 1, 1, 1, 1, 0, 1, 0, 0, 0)z1 <- cosine(x, y)z1z2 <- cosine(matrix(x), matrix(y))z2

Location of Chromosome Cytobands at Genome Build T2T

Description

Location of Chromosome Cytobands at Genome Build T2T

Format

A data.frame

Source

from T2T study

Examples

data(cytobands.T2T)

Location of Chromosome Cytobands at Genome Build hg19

Description

Location of Chromosome Cytobands at Genome Build hg19

Format

A data.frame

Source

from UCSC

Examples

data(cytobands.hg19)

Location of Chromosome Cytobands at Genome Build hg38

Description

Location of Chromosome Cytobands at Genome Build hg38

Format

A data.frame

Source

from UCSC

Examples

data(cytobands.hg38)

Location of Chromosome Cytobands at Genome Build mm10

Description

Location of Chromosome Cytobands at Genome Build mm10

Format

A data.frame

Source

from UCSChttp://hgdownload.cse.ucsc.edu/goldenpath/mm10/database/cytoBand.txt.gz

Examples

data(cytobands.mm10)

Location of Chromosome Cytobands at Genome Build mm9

Description

Location of Chromosome Cytobands at Genome Build mm9

Format

A data.frame

Source

from UCSChttp://hgdownload.cse.ucsc.edu/goldenpath/mm9/database/cytoBand.txt.gz

Examples

data(cytobands.mm9)

Performs Strand Bias Enrichment Analysis for a Given Sample-by-Component Matrix

Description

Seesig_tally for examples.

Usage

enrich_component_strand_bias(mat)

Arguments

Value

adata.table sorted byp_value.

Get Aneuploidy Score from Copy Number Profile

Description

This implements a Cohen-Sharir method (see reference) like "Aneuploidy Score" computation.You can read the source code to see how it works. Basically, it followsthe logic of Cohen-Sharir method but with some difference in detail implementation.Their results should be counterpart, but with no data validation for now.Please raise an issue if you find problem/bugs in this function.

Usage

get_Aneuploidy_score(  data,  ploidy_df = NULL,  genome_build = "hg19",  rm_black_arms = FALSE)

Arguments

Value

Adata.frame

References

• Cohen-Sharir, Y., McFarland, J. M., Abdusamad, M., Marquis, C., Bernhard, S. V., Kazachkova, M., ... & Ben-David, U. (2021). Aneuploidy renders cancer cells vulnerable to mitotic checkpoint inhibition. Nature, 1-6.

• Logic reference:https://github.com/quevedor2/aneuploidy_score/.

• Taylor, Alison M., et al. "Genomic and functional approaches to understanding cancer aneuploidy." Cancer cell 33.4 (2018): 676-689.

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))df <- get_Aneuploidy_score(cn)dfdf2 <- get_Aneuploidy_score(cn@data)df2df3 <- get_Aneuploidy_score(cn@data,  ploidy_df = get_cn_ploidy(cn@data))df3

Get Adjust P Values from Group Comparison

Description

Settingaes(label=..p.adj..) inggpubr::compare_means() does notshow adjust p values. The returned result of this function can be combined withggpubr::stat_pvalue_manual() to fixthis problem.

Usage

get_adj_p(  data,  .col,  .grp = "Sample",  comparisons = NULL,  method = "wilcox.test",  p.adjust.method = "fdr",  p.digits = 3L,  ...)

Arguments

Details

More info seeggpubr::compare_means(),ggpubr::stat_compare_means() andstats::p.adjust().

Value

adata.frame containing comparison result

Source

https://github.com/kassambara/ggpubr/issues/143

Examples

library(ggpubr)# T-teststat.test <- compare_means(  len ~ dose,  data = ToothGrowth,  method = "t.test",  p.adjust.method = "fdr")stat.test# Create a simple box plotp <- ggboxplot(ToothGrowth, x = "dose", y = "len")p# Add p valuesmy_comparisons <- list(c("0.5", "1"), c("1", "2"), c("0.5", "2"))p + stat_compare_means(method = "t.test", comparisons = my_comparisons)# Try adding adjust p values# proposed by author of ggpubr# however it does not workp + stat_compare_means(aes(label = ..p.adj..), method = "t.test", comparisons = my_comparisons)# Solution:# calculate adjust p values and their location# then use stat_pvalue_manual() functionp_adj <- get_adj_p(ToothGrowth, .col = "len", .grp = "dose")p_adjp + stat_pvalue_manual(p_adj, label = "p.adj")# Show selected comparisons# Of note, p value is ajusted# for three comparisons, but only# two are showed in figurep_adj <- get_adj_p(ToothGrowth,  .col = "len", .grp = "dose",  comparisons = list(c("0.5", "1"), c("1", "2")))p + stat_pvalue_manual(p_adj, label = "p.adj")

Get Specified Bayesian NMF Result from Run

Description

Sometimes, we may want to use or inspect specified run result fromsig_auto_extract.This function is designed for this purpose.

Usage

get_bayesian_result(run_info)

Arguments

Value

alist.

Author(s)

Shixiang Wang

Examples

load(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))res <- sig_auto_extract(cn_tally_W$nmf_matrix, result_prefix = "Test_copynumber", nrun = 1)# All run info are stored in res$Raw$summary_run# Obtain result of run 1res_run1 <- get_bayesian_result(res$Raw$summary_run[1, ])

Get CNV Frequency Table

Description

Get CNV Frequency Table

Usage

get_cn_freq_table(  data,  genome_build = "hg19",  cutoff = 2L,  resolution_factor = 1L)

Arguments

Value

adata.table.

Get Ploidy from Absolute Copy Number Profile

Description

Get Ploidy from Absolute Copy Number Profile

Usage

get_cn_ploidy(data)

Arguments

Value

a value or adata.table

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))df <- get_cn_ploidy(cn)df

Get Genome Annotation

Description

Get Genome Annotation

Usage

get_genome_annotation(  data_type = c("chr_size", "centro_loc", "cytobands", "transcript", "gene"),  chrs = paste0("chr", c(1:22, "X", "Y")),  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"))

Arguments

Value

adata.frame containing annotation data

Examples

df1 <- get_genome_annotation()df1df2 <- get_genome_annotation(genome_build = "hg38")df2df3 <- get_genome_annotation(data_type = "centro_loc")df3df4 <- get_genome_annotation(data_type = "centro_loc", genome_build = "hg38")df4df5 <- get_genome_annotation(data_type = "cytobands")df5df6 <- get_genome_annotation(data_type = "cytobands", genome_build = "hg38")df6

Get Comparison Result between Signature Groups

Description

Compare genotypes/phenotypes based on signature groups (samples are assigned toseveral groups). For categoricaltype, calculate fisher p value (usingstats::fisher.test) and count table.In larger than 2 by 2 tables, compute p-values by Monte Carlo simulation.For continuous type, calculate anova p value (usingstats::aov),summary table and Tukey Honest significant difference (usingstats::TukeyHSD).The result of this function can be plotted byshow_group_comparison().

Usage

get_group_comparison(  data,  col_group,  cols_to_compare,  type = "ca",  NAs = NA,  verbose = FALSE)

Arguments

Value

alist contains data, summary, p value etc..

Author(s)

Shixiang Wangw_shixiang@163.com

Examples

load(system.file("extdata", "toy_copynumber_signature_by_W.RData",  package = "sigminer", mustWork = TRUE))# Assign samples to clustersgroups <- get_groups(sig, method = "k-means")set.seed(1234)groups$prob <- rnorm(10)groups$new_group <- sample(c("1", "2", "3", "4", NA), size = nrow(groups), replace = TRUE)# Compare groups (filter NAs for categorical coloumns)groups.cmp <- get_group_comparison(groups[, -1],  col_group = "group",  cols_to_compare = c("prob", "new_group"),  type = c("co", "ca"), verbose = TRUE)# Compare groups (Set NAs of categorical columns to 'Rest')groups.cmp2 <- get_group_comparison(groups[, -1],  col_group = "group",  cols_to_compare = c("prob", "new_group"),  type = c("co", "ca"), NAs = "Rest", verbose = TRUE)

Get Sample Groups from Signature Decomposition Information

Description

One of key results from signature analysis is to cluster samples into differentgroups. This function takesSignature object as inputand return the membership in each cluster.

Usage

get_groups(  Signature,  method = c("consensus", "k-means", "exposure", "samples"),  n_cluster = NULL,  match_consensus = TRUE)

Arguments

Details

Users may find there are bigger differences between using method 'samples' and 'exposure' butthey use a similar idear to find dominant signature, here goes the reason:

Method 'samples' using data directly from NMF decomposition, this means the two matrixW (basis matrix or signature matrix) andH (coefficient matrix or exposure matrix) arethe results of NMF. For method 'exposure', it uses the signature exposure loading matrix.In this situation, each signture represents a number of mutations (alterations)about implementation please see source code ofsig_extract() function.

Value

adata.table object

See Also

NMF::predict(),show_groups.

Examples

# Load copy number prepare objectload(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))# Extract copy number signatureslibrary(NMF)sig <- sig_extract(cn_tally_W$nmf_matrix, 2,  nrun = 10)# Methods 'consensus' and 'samples' are from NMF::predict()g1 <- get_groups(sig, method = "consensus", match_consensus = TRUE)g1g2 <- get_groups(sig, method = "samples")g2# Use k-means clusteringg3 <- get_groups(sig, method = "k-means")g3

Get Overlap Size between Interval x and y

Description

Get Overlap Size between Interval x and y

Usage

get_intersect_size(x.start, x.end, y.start, y.end)

Arguments

Value

a numeric vector.

Examples

o1 <- get_intersect_size(1, 5, 3, 20)o1o2 <- get_intersect_size(3, 20, 1, 10)o2o3 <- get_intersect_size(c(1, 2, 1), c(10, 4, 6), c(4, 2, 5), c(10, 3, 22))o3

Get proportions of pLOH score from Allele Specific Copy Number Profile

Description

pLOH score represents the genome that displayed LOH.

Usage

get_pLOH_score(data, rm_chrs = c("chrX", "chrY"), genome_build = "hg19")

Arguments

Value

Adata.frame

References

Steele, Christopher D., et al. "Signatures of copy number alterations in human cancer." bioRxiv (2021).

Examples

# Load toy dataset of absolute copynumber profileload(system.file("extdata", "toy_segTab.RData",  package = "sigminer", mustWork = TRUE))set.seed(1234)segTabs$minor_cn <- sample(c(0, 1), size = nrow(segTabs), replace = TRUE)cn <- read_copynumber(segTabs,  seg_cols = c("chromosome", "start", "end", "segVal"),  genome_measure = "wg", complement = TRUE, add_loh = TRUE)df <- get_pLOH_score(cn)dfdf2 <- get_pLOH_score(cn@data)df2

Get Shannon Diversity Index for Signatures

Description

H = - \sum_{i=1}^n{p_i ln(p_i)}

wheren is the numberof signatures identified in the signature with exposure >cutoff,andpi is the normalized exposure of the ith signature withexposure >cutoff. Exposures of signatures were normalized tosum to1.

Usage

get_shannon_diversity_index(rel_expo, cutoff = 0.001)

Arguments

Value

adata.frame

References

Steele, Christopher D., et al. "Undifferentiated sarcomas develop through distinct evolutionary pathways." Cancer Cell 35.3 (2019): 441-456.

Examples

# Load mutational signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Get signature exposurerel_expo <- get_sig_exposure(sig2, type = "relative")rel_expodiversity_index <- get_shannon_diversity_index(rel_expo)diversity_index

Obtain Signature Index for Cancer Types

Description

Obtain Signature Index for Cancer Types

Usage

get_sig_cancer_type_index(  sig_type = c("legacy", "SBS", "DBS", "ID"),  seq_type = c("WGS", "WES"),  source = c("PCAWG", "TCGA", "nonPCAWG"),  keyword = NULL)

Arguments

Value

alist.

Examples

l1 <- get_sig_cancer_type_index()l2 <- get_sig_cancer_type_index(sig_type = "SBS")l3 <- get_sig_cancer_type_index(sig_type = "DBS", source = "PCAWG", seq_type = "WGS")l4 <- get_sig_cancer_type_index(sig_type = "ID")l5 <- get_sig_cancer_type_index(keyword = "breast")l1l2l3l4l5

Get Curated Reference Signature Database

Description

Reference mutational signatures and their aetiologies,mainly obtained from COSMIC database(SigProfiler results) and cleaned before saving intosigminer package. You can obtain:

• COSMIC legacy SBS signatures.

• COSMIC v3 SBS signatures.

• COSMIC v3 DBS signatures.

• COSMIC v3 ID (indel) signatures.

• SBS and RS (rearrangement) signatures from Nik lab 2020 Nature Cancer paper.

• RS signatures from BRCA560 and USARC cohorts.

• Copy number signatures from USARC cohort and TCGA.

• Copy number signatures from Liu lab 2023. It supports both PCAWG and TCGA cohort.

Usage

get_sig_db(sig_db = "legacy")

Arguments

Value

alist.

References

• Steele, Christopher D., et al. "Signatures of copy number alterations in human cancer." Nature 606.7916 (2022): 984-991.

• Alexandrov, Ludmil B., et al. "The repertoire of mutational signatures in human cancer." Nature 578.7793 (2020): 94-101.

• Steele, Christopher D., et al. "Undifferentiated sarcomas develop through distinct evolutionary pathways." Cancer Cell 35.3 (2019): 441-456.

• Ziyu Tao, et al. "The repertoire of copy number alteration signatures in human cancer." Briefings in Bioinformatics (2023): bbad053.

See Also

get_sig_similarity,sig_fit andshow_cosmic_sig_profile.

Examples

s1 <- get_sig_db()s2 <- get_sig_db("SBS")s3 <- get_sig_db("DBS")s4 <- get_sig_db("DBS_mm10")s5 <- get_sig_db("SBS_Nik_lab")s6 <- get_sig_db("ID")s7 <- get_sig_db("RS_BRCA560")s8 <- get_sig_db("RS_USARC")s9 <- get_sig_db("RS_Nik_lab")s10 <- get_sig_db("CNS_USARC")s11 <- get_sig_db("CNS_TCGA")s12 <- get_sig_db("CNS_TCGA176")s13 <- get_sig_db("CNS_PCAWG176")s1s2s3s4s5s6s7s8s9s10s11s12s13

Get Signature Exposure from 'Signature' Object

Description

The expected number of mutations (or copy number segment records) with each signature wasdetermined after a scaling transformation V ~ WH = W'H' where W' = WU' and H' = UH.The scaling matrix U is a KxK diagnal matrix (K is signature number, U' is the inverse of U)with the element corresponding to the L1-norm of column vectors of W(ie. the sum of the elements of the vector). As a result, the k-th row vector of the finalmatrix H' represents the absolute exposure (activity) of the k-th process across samples(e.g., for SBS, the estimated (or expected) number of mutations generated by the k-th process).Of note, for copy number signatures, only components of feature CN was used for calculating H'.

Usage

get_sig_exposure(  Signature,  type = c("absolute", "relative"),  rel_threshold = 0.01)

Arguments

Value

adata.table

Author(s)

Shixiang Wangw_shixiang@163.com

References

Kim, Jaegil, et al. "Somatic ERCC2 mutations are associated with a distinct genomic signature in urothelial tumors."Nature genetics 48.6 (2016): 600.

Examples

# Load mutational signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Get signature exposureexpo1 <- get_sig_exposure(sig2)expo1expo2 <- get_sig_exposure(sig2, type = "relative")expo2

Calculate Association between Signature Exposures and Other Features

Description

Association of signature exposures with other features will be performed using one of two procedures:for a continuous association variable (including ordinal variable), correaltion is performed;for a binary association variable, samples will be divided into two groups and Mann-Whitney U-testis performed to test for differences in signature exposure medians between the two groups.Seeget_tidy_association for cleaning association result.

Usage

get_sig_feature_association(  data,  cols_to_sigs,  cols_to_features,  type = "ca",  method_co = c("spearman", "pearson", "kendall"),  method_ca = stats::wilcox.test,  min_n = 0.01,  verbose = FALSE,  ...)

Arguments

Value

alist. For 'co' features, 'measure' means correlation coefficient.For 'ca' features, 'measure' means difference in means of signature exposure.

See Also

get_tidy_association

Get Reconstructed Profile Cosine Similarity, RSS, etc.

Description

Seebp_extract_signatures for examples.

Usage

get_sig_rec_similarity(Signature, nmf_matrix)

Arguments

Value

adata.table.

Calculate Similarity between Identified Signatures and Reference Signatures

Description

The reference signatures can be either aSignature object specified byRef argumentor known COSMIC signatures specified bysig_db argument.Two COSMIC databases are used for comparisons - "legacy" which includes 30 signaures,and "SBS" - which includes updated/refined 65 signatures. This function is modifiedfromcompareSignatures() inmaftools package.NOTE: all reference signatures are generated from gold standard tool:SigProfiler.

Usage

get_sig_similarity(  Signature,  Ref = NULL,  sig_db = c("SBS", "legacy", "DBS", "ID", "TSB", "SBS_Nik_lab", "RS_Nik_lab",    "RS_BRCA560", "RS_USARC", "CNS_USARC", "CNS_TCGA", "CNS_TCGA176", "CNS_PCAWG176",    "SBS_hg19", "SBS_hg38", "SBS_mm9", "SBS_mm10", "DBS_hg19", "DBS_hg38", "DBS_mm9",    "DBS_mm10", "SBS_Nik_lab_Organ", "RS_Nik_lab_Organ", "latest_SBS_GRCh37",    "latest_DBS_GRCh37", "latest_ID_GRCh37", "latest_SBS_GRCh38", "latest_DBS_GRCh38",    "latest_SBS_mm9", "latest_DBS_mm9", "latest_SBS_mm10", "latest_DBS_mm10",    "latest_SBS_rn6", "latest_DBS_rn6", "latest_CN_GRCh37",         "latest_RNA-SBS_GRCh37", "latest_SV_GRCh38"),  db_type = c("", "human-exome", "human-genome"),  method = "cosine",  normalize = c("row", "feature"),  feature_setting = sigminer::CN.features,  set_order = TRUE,  pattern_to_rm = NULL,  verbose = TRUE)

Arguments

Value

alist containing smilarities, aetiologies if available, best match and RSS.

Author(s)

Shixiang Wangw_shixiang@163.com

References

Alexandrov, Ludmil B., et al. "The repertoire of mutational signatures in human cancer." Nature 578.7793 (2020): 94-101.

Degasperi, Andrea, et al. "A practical framework and online tool for mutational signature analyses show intertissue variation and driver dependencies." Nature cancer 1.2 (2020): 249-263.

Steele, Christopher D., et al. "Undifferentiated sarcomas develop through distinct evolutionary pathways." Cancer Cell 35.3 (2019): 441-456.

Nik-Zainal, Serena, et al. "Landscape of somatic mutations in 560 breast cancer whole-genome sequences." Nature 534.7605 (2016): 47-54.

Steele, Christopher D., et al. "Signatures of copy number alterations in human cancer." Nature 606.7916 (2022): 984-991.

Examples

# Load mutational signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))s1 <- get_sig_similarity(sig2, Ref = sig2)s1s2 <- get_sig_similarity(sig2)s2s3 <- get_sig_similarity(sig2, sig_db = "SBS")s3# Set order for result similarity matrixs4 <- get_sig_similarity(sig2, sig_db = "SBS", set_order = TRUE)s4## Remove some components## in similarity calculations5 <- get_sig_similarity(sig2,  Ref = sig2,  pattern_to_rm = c("T[T>G]C", "T[T>G]G", "T[T>G]T"))s5## Same to DBS and ID signaturesx1 <- get_sig_db("DBS_hg19")x2 <- get_sig_db("DBS_hg38")s6 <- get_sig_similarity(x1$db, x2$db)s6

Get Tidy Signature Association Results

Description

Get Tidy Signature Association Results

Usage

get_tidy_association(cor_res, p_adjust = FALSE, method = "fdr")

Arguments

Value

adata.frame

See Also

get_sig_feature_association

General Group Enrichment Analysis

Description

This function takes adata.frame as input, compares proportion of positivecases or mean measure in one subgroup and the remaining samples.

Usage

group_enrichment(  df,  grp_vars = NULL,  enrich_vars = NULL,  cross = TRUE,  co_method = c("t.test", "wilcox.test"),  ref_group = NA)

Arguments

Value

adata.table with following columns:

• grp_var: group variable name.

• enrich_var: enrich variable (variable to be compared) name.

• grp1: the first group name, should be a member ingrp_var column.

• grp2: the remaining samples, marked as 'Rest'.

• grp1_size: sample size forgrp1.

• grp1_pos_measure: for binary variable, it stores the proportion ofpositive cases ingrp1; for continuous variable, it stores mean value.

• grp2_size: sample size forgrp2.

• grp2_pos_measure: same asgrp1_pos_measure but forgrp2.

• measure_observed: for binary variable, it stores odds ratio;for continuous variable, it stores scaled mean ratio.

• measure_tested: only for binary variable, it storesestimated odds ratio and its 95% CI fromfisher.test().

• p_value: for binary variable, it stores p value fromfisher.test();for continuous variable, it stores value fromwilcox.test() ort.test().

• type: one of "binary" and "continuous".

• method: one of "fish.test", "wilcox.test" and "t.test".

See Also

show_group_enrichment

Examples

set.seed(1234)df <- dplyr::tibble(  g1 = factor(abs(round(rnorm(99, 0, 1)))),  g2 = rep(LETTERS[1:4], c(50, 40, 8, 1)),  e1 = sample(c("P", "N"), 99, replace = TRUE),  e2 = rnorm(99))print(str(df))print(head(df))# Compare g1:e1, g1:e2, g2:e1 and g2:e2x1 <- group_enrichment(df, grp_vars = c("g1", "g2"), enrich_vars = c("e1", "e2"))x1# Only compare g1:e1, g2:e2x2 <- group_enrichment(df,  grp_vars = c("g1", "g2"),  enrich_vars = c("e1", "e2"),  co_method = "wilcox.test",  cross = FALSE)x2# Visualizationp1 <- show_group_enrichment(x1, fill_by_p_value = TRUE)p1p2 <- show_group_enrichment(x1, fill_by_p_value = FALSE)p2p3 <- show_group_enrichment(x1, return_list = TRUE)p3

Group Enrichment Analysis with Subsets

Description

More details seegroup_enrichment().

Usage

group_enrichment2(  df,  subset_var,  grp_vars,  enrich_vars,  co_method = c("t.test", "wilcox.test"),  ref_group = NA)

Arguments

See Also

show_group_enrichment

Handle Hypermutant Samples

Description

This can be used for SNV/INDEL count matrix. For copy number analysis,please skip it.

Usage

handle_hyper_mutation(nmf_matrix)

Arguments

Value

amatrix.

References

Kim, Jaegil, et al. "Somatic ERCC2 mutations are associated with a distinct genomic signature in urothelial tumors."Nature genetics 48.6 (2016): 600.

Say Hello to Users

Description

Say Hello to Users

Usage

hello()

Examples

hello()

Output Signature Bootstrap Fitting Results

Description

Output Signature Bootstrap Fitting Results

Usage

output_bootstrap(x, result_dir, mut_type = "SBS", sig_db = mut_type)

Arguments

Value

Nothing.

Output Signature Fitting Results

Description

Output Signature Fitting Results

Usage

output_fit(x, result_dir, mut_type = "SBS", sig_db = mut_type)

Arguments

Value

Nothing.

Output Signature Results

Description

Output Signature Results

Usage

output_sig(sig, result_dir, mut_type = "SBS", sig_db = mut_type)

Arguments

Value

Nothing.

Output Tally Result in Barplots

Description

Output Tally Result in Barplots

Usage

output_tally(x, result_dir, mut_type = "SBS")

Arguments

Value

Nothing.

Read Absolute Copy Number Profile

Description

Readabsolute copy number profile for preparing CNV signatureanalysis. See detail part ofsig_tally() to see how to handle sex to get correctsummary.

Usage

read_copynumber(  input,  pattern = NULL,  ignore_case = FALSE,  seg_cols = c("Chromosome", "Start.bp", "End.bp", "modal_cn"),  samp_col = "sample",  add_loh = FALSE,  loh_min_len = 10000,  loh_min_frac = 0.05,  join_adj_seg = TRUE,  skip_annotation = FALSE,  use_all = add_loh,  min_segnum = 0L,  max_copynumber = 20L,  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  genome_measure = c("called", "wg"),  complement = FALSE,  ...)

Arguments

Value

aCopyNumber object.

Author(s)

Shixiang Wangw_shixiang@163.com

See Also

read_maf for reading mutation data toMAF object.

Examples

# Load toy dataset of absolute copynumber profileload(system.file("extdata", "toy_segTab.RData",  package = "sigminer", mustWork = TRUE))cn <- read_copynumber(segTabs,  seg_cols = c("chromosome", "start", "end", "segVal"),  genome_build = "hg19", complement = FALSE)cncn_subset <- subset(cn, sample == "TCGA-DF-A2KN-01A-11D-A17U-01")# Add LOHset.seed(1234)segTabs$minor_cn <- sample(c(0, 1), size = nrow(segTabs), replace = TRUE)cn <- read_copynumber(segTabs,  seg_cols = c("chromosome", "start", "end", "segVal"),  genome_measure = "wg", complement = TRUE, add_loh = TRUE)# Use tally method "S" (Steele et al.)tally_s <- sig_tally(cn, method = "S")tab_file <- system.file("extdata", "metastatic_tumor.segtab.txt",  package = "sigminer", mustWork = TRUE)cn2 <- read_copynumber(tab_file)cn2

Read Copy Number Data from ASCAT Result Files

Description

Note, the result is not aCopyNumber object, you need to generate itby yourself.

Usage

read_copynumber_ascat(x)

Arguments

Value

a tidylist.

Read Absolute Copy Number Profile from Sequenza Result Directory

Description

Read Absolute Copy Number Profile from Sequenza Result Directory

Usage

read_copynumber_seqz(target_dir, return_df = FALSE, ...)

Arguments

Value

adata.frame or aCopyNumber object.

Read MAF Files

Description

This function is a wrapper ofmaftools::read.maf.Useless options inmaftools::read.maf are dropped here.You can also usemaftools::read.maf to read the data.All reference alleles and mutation alleles should be recorded inpositive strand format.

Usage

read_maf(maf, verbose = TRUE)read_maf_minimal(dt)

Arguments

Functions

• read_maf_minimal(): Read Maf data.frame from a minimal maf-like data

See Also

read_copynumber for reading copy number data toCopyNumber object.

Examples

laml.maf <- system.file("extdata", "tcga_laml.maf.gz", package = "maftools", mustWork = TRUE)if (!require("R.utils")) {  message("Please install 'R.utils' package firstly")} else {  laml <- read_maf(maf = laml.maf)  laml  laml_mini <- laml@data[, list(    Tumor_Sample_Barcode, Chromosome,    Start_Position, End_Position,    Reference_Allele, Tumor_Seq_Allele2  )]  laml2 <- read_maf_minimal(laml_mini)  laml2}

Read Structural Variation Data as RS object

Description

Read Structural Variation Data as RS object

Usage

read_sv_as_rs(input)

Arguments

Value

alist

Examples

sv <- readRDS(system.file("extdata", "toy_sv.rds", package = "sigminer", mustWork = TRUE))rs <- read_sv_as_rs(sv)# svclass is optionalrs2 <- read_sv_as_rs(sv[, setdiff(colnames(sv), "svclass")])identical(rs, rs2)## Not run: tally_rs <- sig_tally(rs)## End(Not run)

Read VCF Files as MAF Object

Description

MAF file is more recommended. In this function, we will mimicthe MAF object from the keyc(1, 2, 4, 5, 7) columns of VCF file.

Usage

read_vcf(  vcfs,  samples = NULL,  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  keep_only_pass = FALSE,  verbose = TRUE)

Arguments

Value

aMAF.

See Also

read_maf,read_copynumber

Examples

vcfs <- list.files(system.file("extdata", package = "sigminer"), "*.vcf", full.names = TRUE)maf <- read_vcf(vcfs)maf <- read_vcf(vcfs, keep_only_pass = TRUE)

Read UCSC Xena Variant Format Data as MAF Object

Description

Read UCSC Xena Variant Format Data as MAF Object

Usage

read_xena_variants(path)

Arguments

Value

aMAF object.

Examples

if (requireNamespace("UCSCXenaTools")) {  library(UCSCXenaTools)  options(use_hiplot = TRUE)  example_file <- XenaGenerate(subset = XenaDatasets == "mc3/ACC_mc3.txt") %>%    XenaQuery() %>%    XenaDownload()  x <- read_xena_variants(example_file$destfiles)  x@data  y <- sig_tally(x)  y}

Report P Values from bootstrap Results

Description

See examples insig_fit_bootstrap.

Usage

report_bootstrap_p_value(x, thresholds = c(0.01, 0.05, 0.1))

Arguments

Value

a (list of)matrix

Same Size Clustering

Description

This is a wrapper for several implementation that classify samples intosame size clusters, the details please seethis blog.The source code is modified based on code from the blog.

Usage

same_size_clustering(  mat,  diss = FALSE,  clsize = NULL,  algo = c("nnit", "hcbottom", "kmvar"),  method = c("maxd", "random", "mind", "elki", "ward.D", "average", "complete", "single"))

Arguments

Value

a vector.

Examples

set.seed(1234L)x <- rbind(  matrix(rnorm(100, sd = 0.3), ncol = 2),  matrix(rnorm(100, mean = 1, sd = 0.3), ncol = 2))colnames(x) <- c("x", "y")y1 <- same_size_clustering(x, clsize = 10)y11 <- same_size_clustering(as.matrix(dist(x)), clsize = 10, diss = TRUE)y2 <- same_size_clustering(x, clsize = 10, algo = "hcbottom", method = "ward.D")y3 <- same_size_clustering(x, clsize = 10, algo = "kmvar")y33 <- same_size_clustering(as.matrix(dist(x)), clsize = 10, algo = "kmvar", diss = TRUE)

Score Copy Number Profile

Description

Returns quantification of copy number profile and events includingtandem duplication and Chromothripisis etc.Only copy number data from autosome is used here.Some of the quantification methods are rough,you use at your risk. You should do some extra work to check theresult scores.

Usage

scoring(object, TD_size_cutoff = c(1000, 1e+05, 2e+06), TD_cn_cutoff = Inf)

Arguments

Value

adata.table with following scores:

• cnaBurden: CNA burden representing the altered genomic fraction as previously reported.

• cnaLoad: CNA load representing the quantity of copy number alteration.

• MACN: mean altered copy number (MACN) reflecting the property of altered copy number segments,calculated as

  MACN = \frac{\sum_{i} CN_i}{N_{cnv}}

  whereCN_i is the copy number of altered segmenti,N_{cnv} isthe number of CNV.

• weightedMACN: same as MACN but weighted with segment length.

  MACN_{weighted} = \frac{\sum_{i} (CN_i \times L_{i})}{ \sum_{i} L_{i} }

  whereL_{i} is the length of altered copy number segmenti.

• Ploidy: ploidy, the formula is same asweightedMACN but using all copy number segments instead ofaltered copy number segments.

• TDP_pnas: tandem duplication phenotype score from⁠https://www.pnas.org/doi/10.1073/pnas.1520010113⁠,the thresholdk in reference is omitted.

  TDP = - \frac{\sum_{chr} |TD_{obs}-TD_{exp}|}{TD_{total}}

  whereTD_{total} is the number of TD,TD_{obs} andTD_exp are observed number of TD and expected number of TD for each chromosome.

• TDP: tandem duplication score used defined by our group work,TD represents segment with copy number greater than 2.

  TD = \frac{TD_{total}}{\sum_{chr} |TD_{obs}-TD_{exp}|+1}

• sTDP: TDP score for short TD.

• lTDP: TDP score for long TD.

• TDP_size : TDP region size (Mb).

• sTDP_size: sTDP region size (Mb).

• lTDP_size: lTDP region size(Mb).

• Chromoth_state: chromothripsis state score,according to referencedoi:10.1016/j.cell.2013.02.023,chromothripsis frequently leads to massive loss of segments onthe affected chromosome with segmental losses being interspersed with regions displayingnormal (disomic) copy-number (e.g., copy-number states oscillating betweencopy-number = 1 and copy-number = 2), form tens to hundreds of locally clustered DNA rearrangements.Most of methods use both SV and CNV to infer chromothripsis, here we roughly quantify it with

  \sum_{chr}{N_{OsCN}^2}

  whereN_{OsCN} is the number of oscillating copy number pattern "2-1-2" for each chromosome.

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))d <- scoring(cn)dd2 <- scoring(cn, TD_cn_cutoff = 4L)d2

Show Alteration Catalogue Profile

Description

Show Alteration Catalogue Profile

Usage

show_catalogue(  catalogue,  mode = c("SBS", "copynumber", "DBS", "ID", "RS"),  method = "Wang",  normalize = c("raw", "row", "feature"),  style = c("default", "cosmic"),  samples = NULL,  samples_name = NULL,  x_lab = "Components",  y_lab = "Counts",  ...)

Arguments

Value

aggplot object

Examples

data("simulated_catalogs")p <- show_catalogue(simulated_catalogs$set1, style = "cosmic")p

Show Copy Number Profile in Circos

Description

Another visualization method for copy number profile likeshow_cn_profile.

Usage

show_cn_circos(  data,  samples = NULL,  show_title = TRUE,  chrs = paste0("chr", 1:22),  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  col = NULL,  side = "inside",  ...)

Arguments

Value

a circos plot

Examples

load(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))show_cn_circos(cn, samples = 1)show_cn_circos(cn, samples = "TCGA-99-7458-01A-11D-2035-01")## Remove titleshow_cn_circos(cn, samples = 1, show_title = FALSE)## Subset chromosomesshow_cn_circos(cn, samples = 1, chrs = c("chr1", "chr2", "chr3"))## Arrange plotslayout(matrix(1:4, 2, 2))show_cn_circos(cn, samples = 4)layout(1) # reset layout

Show Copy Number Components

Description

Show classified components ("Wang" ("W") method) for copy number data.

Usage

show_cn_components(  parameters,  method = "Wang",  show_weights = TRUE,  log_y = FALSE,  return_plotlist = FALSE,  base_size = 12,  nrow = 2,  align = "hv",  ...)

Arguments

Value

aggplot object

Author(s)

Shixiang Wangw_shixiang@163.com

Show Copy Number Distribution either by Length or Chromosome

Description

Visually summarize copy number distribution either by copy number segment lengthor chromosome. Input is aCopyNumber object,genome_build option willread fromgenome_build slot of object.

Usage

show_cn_distribution(  data,  rm_normal = TRUE,  mode = c("ld", "cd"),  fill = FALSE,  scale_chr = TRUE,  base_size = 14)

Arguments

Value

aggplot object

Author(s)

Shixiang Wangw_shixiang@163.com

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))# Plot distributionp1 <- show_cn_distribution(cn)p1p2 <- show_cn_distribution(cn, mode = "cd")p2p3 <- show_cn_distribution(cn, mode = "cd", fill = TRUE)p3

Show Copy Number Feature Distributions

Description

Show Copy Number Feature Distributions

Usage

show_cn_features(  features,  method = "Wang",  rm_outlier = FALSE,  ylab = NULL,  log_y = FALSE,  return_plotlist = FALSE,  base_size = 12,  nrow = 2,  align = "hv",  ...)

Arguments

Value

aggplot object

Show Copy Number Variation Frequency Profile with Circos

Description

Show Copy Number Variation Frequency Profile with Circos

Usage

show_cn_freq_circos(  data,  groups = NULL,  cutoff = 2L,  resolution_factor = 1L,  title = c("AMP", "DEL"),  chrs = paste0("chr", 1:22),  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  cols = NULL,  plot_ideogram = TRUE,  track_height = 0.5,  ideogram_height = 1,  ...)

Arguments

Value

Nothing.

Examples

load(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))show_cn_freq_circos(cn)ss <- unique(cn@data$sample)show_cn_freq_circos(cn, groups = list(a = ss[1:5], b = ss[6:10]), cols = c("red", "green"))

Show Summary Copy Number Profile for Sample Groups

Description

Show Summary Copy Number Profile for Sample Groups

Usage

show_cn_group_profile(  data,  groups = NULL,  fill_area = TRUE,  cols = NULL,  chrs = paste0("chr", c(1:22, "X")),  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  cutoff = 2L,  resolution_factor = 1L,  force_y_limit = TRUE,  highlight_genes = NULL,  repel = FALSE,  nrow = NULL,  ncol = NULL,  return_plotlist = FALSE)

Arguments

Value

a (list of)ggplot object.

Examples

load(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))p1 <- show_cn_group_profile(cn)p1ss <- unique(cn@data$sample)p2 <- show_cn_group_profile(cn, groups = list(a = ss[1:5], b = ss[6:10]))p2p3 <- show_cn_group_profile(cn,  groups = list(g1 = ss[1:5], g2 = ss[6:10]),  force_y_limit = c(-1, 1), nrow = 2)p3## Set custom cutoff for custom datadata <- cn@datadata$segVal <- data$segVal - 2Lp4 <- show_cn_group_profile(data,  groups = list(g1 = ss[1:5], g2 = ss[6:10]),  force_y_limit = c(-1, 1), nrow = 2,  cutoff = c(0, 0))p4## Add highlight genep5 <- show_cn_group_profile(cn, highlight_genes = c("TP53", "EGFR"))p5

Show Sample Copy Number Profile

Description

Sometimes it is very useful to check details about copy number profile for one or multiplesamples. This function is designed to do this job and can be further modified byggplot2related packages.

Usage

show_cn_profile(  data,  samples = NULL,  show_n = NULL,  show_title = FALSE,  show_labels = NULL,  chrs = paste0("chr", 1:22),  position = NULL,  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  ylim = NULL,  nrow = NULL,  ncol = NULL,  return_plotlist = FALSE)

Arguments

Value

aggplot object or alist

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))p <- show_cn_profile(cn, nrow = 2, ncol = 1)pp2 <- show_cn_profile(cn,  nrow = 2, ncol = 1,  position = "chr1:3218923-116319008")p2

A Simple and General Way for Association Analysis

Description

All variables must be continuous.The matrix will be returned as an element ofggplot object.This is basically a wrapper of R packageggcorrplot.

Usage

show_cor(  data,  x_vars = colnames(data),  y_vars = x_vars,  cor_method = "spearman",  vis_method = "square",  lab = TRUE,  test = TRUE,  hc_order = FALSE,  p_adj = NULL,  ...)

Arguments

Value

aggplot object

See Also

show_sig_feature_corrplot for specific and more powerfulassociation analysis and visualization.

Examples

data("mtcars")p1 <- show_cor(mtcars)p2 <- show_cor(mtcars,  x_vars = colnames(mtcars)[1:4],  y_vars = colnames(mtcars)[5:8])p3 <- show_cor(mtcars, vis_method = "circle", p_adj = "fdr")p1p1$corp2p3## Auto detect problem variablesmtcars$xx <- 0Lp4 <- show_cor(mtcars)p4

Show Signature Information in Web Browser

Description

Show Signature Information in Web Browser

Usage

show_cosmic(x = "home")

Arguments

Value

Nothing.

Examples

## Not run: show_cosmic()show_cosmic("legacy")show_cosmic("SBS")show_cosmic("DBS")show_cosmic("ID")show_cosmic("SBS1")show_cosmic("DBS2")show_cosmic("ID3")## End(Not run)

Plot Reference (Mainly COSMIC) Signature Profile

Description

Plot Reference (Mainly COSMIC) Signature Profile

Usage

show_cosmic_sig_profile(  sig_index = NULL,  show_index = TRUE,  sig_db = "legacy",  ...)

Arguments

Value

aggplot object

Author(s)

Shixiang Wangw_shixiang@163.com

Examples

show_cosmic_sig_profile()show_cosmic_sig_profile(sig_db = "SBS")show_cosmic_sig_profile(sig_index = 1:5)show_cosmic_sig_profile(sig_db = "SBS", sig_index = c("10a", "17a"))gg <- show_cosmic_sig_profile(sig_index = 1:5)gg$aetiology

Plot Group Comparison Result

Description

Using result data fromget_group_comparison, this function plotsgenotypes/phenotypes comparison between signature groups usingggplot2 package and returna list ofggplot object contains individual and combined plots. The combinedplot is easily saved to local usingcowplot::save_plot(). Of note, default fishertest p values are shown for categorical data and fdr values are shown forcontinuous data.

Usage

show_group_comparison(  group_comparison,  xlab = "group",  ylab_co = NA,  legend_title_ca = NA,  legend_position_ca = "bottom",  set_ca_sig_yaxis = FALSE,  set_ca_custom_xlab = FALSE,  show_pvalue = TRUE,  ca_p_threshold = 0.01,  method = "wilcox.test",  p.adjust.method = "fdr",  base_size = 12,  font_size_x = 12,  text_angle_x = 30,  text_hjust_x = 0.2,  ...)

Arguments

Value

alist ofggplot objects.

Author(s)

Shixiang Wangw_shixiang@163.com

Examples

load(system.file("extdata", "toy_copynumber_signature_by_W.RData",  package = "sigminer", mustWork = TRUE))# Assign samples to clustersgroups <- get_groups(sig, method = "k-means")set.seed(1234)groups$prob <- rnorm(10)groups$new_group <- sample(c("1", "2", "3", "4", NA), size = nrow(groups), replace = TRUE)# Compare groups (filter NAs for categorical coloumns)groups.cmp <- get_group_comparison(groups[, -1],  col_group = "group",  cols_to_compare = c("prob", "new_group"),  type = c("co", "ca"), verbose = TRUE)# Compare groups (Set NAs of categorical columns to 'Rest')groups.cmp2 <- get_group_comparison(groups[, -1],  col_group = "group",  cols_to_compare = c("prob", "new_group"),  type = c("co", "ca"), NAs = "Rest", verbose = TRUE)show_group_comparison(groups.cmp)ggcomp <- show_group_comparison(groups.cmp2)ggcomp$co_combggcomp$ca_comb

Show Groupped Variable Distribution

Description

This is a general function, it can be used in any proper analysis.

Usage

show_group_distribution(  data,  gvar,  dvar,  fun = stats::median,  order_by_fun = FALSE,  alpha = 0.8,  g_label = "label",  g_angle = 0,  g_position = "top",  point_size = 1L,  segment_size = 1L,  segment_color = "red",  xlab = NULL,  ylab = NULL,  nrow = 1L,  background_color = c("#DCDCDC", "#F5F5F5"))

Arguments

Value

aggplot object.

Author(s)

Shixiang Wangw_shixiang@163.com

Examples

set.seed(1234)data <- data.frame(  yval = rnorm(120),  gr = c(rep("A", 50), rep("B", 40), rep("C", 30)))p <- show_group_distribution(data,  gvar = 2, dvar = 1,  g_label = "norm",  background_color = "grey")pp2 <- show_group_distribution(data,  gvar = "gr", dvar = "yval",  g_position = "bottom",  order_by_fun = TRUE,  alpha = 0.3)p2# Set custom group namesp3 <- show_group_distribution(data,  gvar = 2, dvar = 1,  g_label = c("A" = "X", "B" = "Y", "C" = "Z"))p3

Show Group Enrichment Result

Description

Seegroup_enrichment for examples.NOTE the box fill and the box text have different meanings.

Usage

show_group_enrichment(  df_enrich,  return_list = FALSE,  scales = "free",  add_text_annotation = TRUE,  fill_by_p_value = TRUE,  use_fdr = TRUE,  cut_p_value = FALSE,  cut_breaks = c(-Inf, -5, log10(0.05), -log10(0.05), 5, Inf),  cut_labels = c("↓ 1e-5", "↓ 0.05", "non-significant", "↑ 0.05", "↑ 1e-5"),  fill_scale = scale_fill_gradient2(low = "#08A76B", mid = "white", high = "red",    midpoint = ifelse(fill_by_p_value, 0, 1)),  cluster_row = FALSE,  cluster_col = FALSE,  ...)

Arguments

Value

a (list of)ggplot object.

Map Groups using Sankey

Description

This feature is designed for signature analysis. However, users can also useit in other similar situations.

Usage

show_group_mapping(  data,  col_to_flow,  cols_to_map,  include_sig = FALSE,  fill_na = FALSE,  title = NULL,  xlab = NULL,  ylab = NULL,  custom_theme = cowplot::theme_minimal_hgrid())

Arguments

Value

aggplot object

Examples

data <- dplyr::tibble(  Group1 = rep(LETTERS[1:5], each = 10),  Group2 = rep(LETTERS[6:15], each = 5),  zzzz = c(rep("xx", 20), rep("yy", 20), rep(NA, 10)))p1 <- show_group_mapping(data, col_to_flow = "Group1", cols_to_map = colnames(data)[-1])p1p2 <- show_group_mapping(data,  col_to_flow = "Group1", cols_to_map = colnames(data)[-1],  include_sig = TRUE)p2

Show Signature Contribution in Clusters

Description

See example section insig_fit() for an examples.

Usage

show_groups(grp_dt, ...)

Arguments

Value

nothing.

See Also

get_groups,sig_fit.

Show Signature Bootstrap Analysis Results

Description

See details for description.

Usage

show_sig_bootstrap_exposure(  bt_result,  sample = NULL,  signatures = NULL,  methods = "QP",  plot_fun = c("boxplot", "violin"),  agg_fun = c("mean", "median", "min", "max"),  highlight = "auto",  highlight_size = 4,  palette = "aaas",  title = NULL,  xlab = FALSE,  ylab = "Signature exposure",  width = 0.3,  dodge_width = 0.8,  outlier.shape = NA,  add = "jitter",  add.params = list(alpha = 0.3),  ...)show_sig_bootstrap_error(  bt_result,  sample = NULL,  methods = "QP",  plot_fun = c("boxplot", "violin"),  agg_fun = c("mean", "median"),  highlight = "auto",  highlight_size = 4,  palette = "aaas",  title = NULL,  xlab = FALSE,  ylab = "Reconstruction error (L2 norm)",  width = 0.3,  dodge_width = 0.8,  outlier.shape = NA,  add = "jitter",  add.params = list(alpha = 0.3),  legend = "none",  ...)show_sig_bootstrap_stability(  bt_result,  signatures = NULL,  measure = c("RMSE", "CV", "MAE", "AbsDiff"),  methods = "QP",  plot_fun = c("boxplot", "violin"),  palette = "aaas",  title = NULL,  xlab = FALSE,  ylab = "Signature instability",  width = 0.3,  outlier.shape = NA,  add = "jitter",  add.params = list(alpha = 0.3),  ...)

Arguments

Details

Functions:

• show_sig_bootstrap_exposure - this function plots exposures from bootstrap samples with both dotted boxplot.The optimal exposure (the exposure from original input) is shown as triangle point.Only one sample can be plotted.

• show_sig_bootstrap_error - this function plots decomposition errors from bootstrap samples with both dotted boxplot.The error from optimal solution (the decomposition error from original input) is shown as triangle point.Only one sample can be plotted.

• show_sig_bootstrap_stability - this function plots the signature exposure instability for specified signatures. Currently,the instability measure supports 3 types:

  • 'RMSE' for Mean Root Squared Error (default) of bootstrap exposures and original exposures for each sample.

  • 'CV' for Coefficient of Variation (CV) based on RMSE (i.e.RMSE / btExposure_mean).

  • 'MAE' for Mean Absolute Error of bootstrap exposures and original exposures for each sample.

  • 'AbsDiff' for Absolute Difference between mean bootstram exposure and original exposure.

Value

aggplot object

References

Huang X, Wojtowicz D, Przytycka TM. Detecting presence of mutational signatures in cancer with confidence. Bioinformatics. 2018;34(2):330–337. doi:10.1093/bioinformatics/btx604

See Also

sig_fit_bootstrap_batch,sig_fit,sig_fit_bootstrap

Examples

if (require("BSgenome.Hsapiens.UCSC.hg19")) {  laml.maf <- system.file("extdata", "tcga_laml.maf.gz", package = "maftools")  laml <- read_maf(maf = laml.maf)  mt_tally <- sig_tally(    laml,    ref_genome = "BSgenome.Hsapiens.UCSC.hg19",    use_syn = TRUE  )  library(NMF)  mt_sig <- sig_extract(mt_tally$nmf_matrix,    n_sig = 3,    nrun = 2,    cores = 1  )  mat <- t(mt_tally$nmf_matrix)  mat <- mat[, colSums(mat) > 0]  bt_result <- sig_fit_bootstrap_batch(mat, sig = mt_sig, n = 10)  ## Parallel computation  ## bt_result = sig_fit_bootstrap_batch(mat, sig = mt_sig, n = 10, use_parallel = TRUE)  ## At default, mean bootstrap exposure for each sample has been calculated  p <- show_sig_bootstrap_exposure(bt_result, methods = c("QP"))  ## Show bootstrap exposure (optimal exposure is shown as triangle)  p1 <- show_sig_bootstrap_exposure(bt_result, methods = c("QP"), sample = "TCGA-AB-2802")  p1  p2 <- show_sig_bootstrap_exposure(bt_result,    methods = c("QP"),    sample = "TCGA-AB-3012",    signatures = c("Sig1", "Sig2")  )  p2  ## Show bootstrap error  ## Similar to exposure above  p <- show_sig_bootstrap_error(bt_result, methods = c("QP"))  p  p3 <- show_sig_bootstrap_error(bt_result, methods = c("QP"), sample = "TCGA-AB-2802")  p3  ## Show exposure (in)stability  p4 <- show_sig_bootstrap_stability(bt_result, methods = c("QP"))  p4  p5 <- show_sig_bootstrap_stability(bt_result, methods = c("QP"), measure = "MAE")  p5  p6 <- show_sig_bootstrap_stability(bt_result, methods = c("QP"), measure = "AbsDiff")  p6  p7 <- show_sig_bootstrap_stability(bt_result, methods = c("QP"), measure = "CV")  p7} else {  message("Please install package 'BSgenome.Hsapiens.UCSC.hg19' firstly!")}

Show Signature Consensus Map

Description

This function is a wrapper ofNMF::consensusmap().

Usage

show_sig_consensusmap(  sig,  main = "Consensus matrix",  tracks = c("consensus:", "silhouette:"),  lab_row = NA,  lab_col = NA,  ...)

Arguments

Value

nothing

Plot Signature Exposure

Description

Currently support copy number signatures and mutational signatures.

Usage

show_sig_exposure(  Signature,  sig_names = NULL,  groups = NULL,  grp_order = NULL,  grp_size = NULL,  samps = NULL,  cutoff = NULL,  style = c("default", "cosmic"),  palette = use_color_style(style),  base_size = 12,  font_scale = 1,  rm_space = FALSE,  rm_grid_line = TRUE,  rm_panel_border = FALSE,  hide_samps = TRUE,  legend_position = "top")

Arguments

Value

aggplot object

Author(s)

Shixiang Wang

Examples

# Load mutational signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Show signature exposurep1 <- show_sig_exposure(sig2)p1# Load copy number signatureload(system.file("extdata", "toy_copynumber_signature_by_W.RData",  package = "sigminer", mustWork = TRUE))# Show signature exposurep2 <- show_sig_exposure(sig)p2

Draw Corrplot for Signature Exposures and Other Features

Description

This function is for association visualization. Of note,the parametersp_val anddrop will affect the visualizationof association results under p value threshold.

Usage

show_sig_feature_corrplot(  tidy_cor,  feature_list,  sort_features = FALSE,  sig_orders = NULL,  drop = TRUE,  return_plotlist = FALSE,  p_val = 0.05,  xlab = "Signatures",  ylab = "Features",  co_gradient_colors = scale_color_gradient2(low = "blue", mid = "white", high = "red",    midpoint = 0),  ca_gradient_colors = co_gradient_colors,  plot_ratio = "auto",  breaks_count = NULL)

Arguments

Value

aggplot2 object

See Also

get_tidy_association andget_sig_feature_association

Examples

# The data is generated from Wang, Shixiang et al.load(system.file("extdata", "asso_data.RData",  package = "sigminer", mustWork = TRUE))p <- show_sig_feature_corrplot(            tidy_data.seqz.feature,            p_val = 0.05,            breaks_count = c(0L,200L, 400L, 600L, 800L, 1020L))p

Show Signature Fit Result

Description

Seesig_fit for examples.

Usage

show_sig_fit(  fit_result,  samples = NULL,  signatures = NULL,  plot_fun = c("boxplot", "violin", "scatter"),  palette = "aaas",  title = NULL,  xlab = FALSE,  ylab = "Signature exposure",  legend = "none",  width = 0.3,  outlier.shape = NA,  add = "jitter",  add.params = list(alpha = 0.3),  ...)

Arguments

Value

aggplot object.

See Also

sig_fit,show_sig_bootstrap_exposure,sig_fit_bootstrap,sig_fit_bootstrap_batch

Show Signature Profile

Description

Who don't like to show a barplot for signature profile? This is for it.

Usage

show_sig_profile(  Signature,  mode = c("SBS", "copynumber", "DBS", "ID", "RS"),  method = "Wang",  by_context = FALSE,  normalize = c("row", "column", "raw", "feature"),  y_tr = NULL,  filters = NULL,  feature_setting = sigminer::CN.features,  style = c("default", "cosmic"),  palette = use_color_style(style, ifelse(by_context, "SBS", mode), method),  set_gradient_color = FALSE,  free_space = "free_x",  rm_panel_border = style == "cosmic",  rm_grid_line = style == "cosmic",  rm_axis_text = FALSE,  bar_border_color = ifelse(style == "default", "grey50", "white"),  bar_width = 0.7,  paint_axis_text = TRUE,  x_label_angle = ifelse(mode == "copynumber" & !(startsWith(method, "T") | method ==    "X"), 60, 90),  x_label_vjust = ifelse(mode == "copynumber" & !(startsWith(method, "T") | method ==    "X"), 1, 0.5),  x_label_hjust = 1,  x_lab = "Components",  y_lab = "auto",  y_limits = NULL,  params = NULL,  show_cv = FALSE,  params_label_size = 3,  params_label_angle = 60,  y_expand = 1,  digits = 2,  base_size = 12,  font_scale = 1,  sig_names = NULL,  sig_orders = NULL,  check_sig_names = TRUE)

Arguments

Value

aggplot object

Author(s)

Shixiang Wang

See Also

show_sig_profile_loop,show_sig_profile_heatmap

Examples

# Load SBS signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Show signature profilep1 <- show_sig_profile(sig2, mode = "SBS")p1# Use 'y_tr' option to transform values in y axisp11 <- show_sig_profile(sig2, mode = "SBS", y_tr = function(x) x * 100)p11# Load copy number signature from method "W"load(system.file("extdata", "toy_copynumber_signature_by_W.RData",  package = "sigminer", mustWork = TRUE))# Show signature profilep2 <- show_sig_profile(sig,  style = "cosmic",  mode = "copynumber",  method = "W",  normalize = "feature")p2# Visualize rearrangement signaturess <- get_sig_db("RS_Nik_lab")ss <- s$db[, 1:3]colnames(ss) <- c("Sig1", "Sig2", "Sig3")p3 <- show_sig_profile(ss, mode = "RS", style = "cosmic")p3

Show Signature Profile with Heatmap

Description

This is a complementary function toshow_sig_profile(), it is used for visualizingsome big signatures, i.e. SBS-1536, not all signatures are supported. See details forcurrent supported signatures.

Usage

show_sig_profile_heatmap(  Signature,  mode = c("SBS", "DBS"),  normalize = c("row", "column", "raw"),  filters = NULL,  x_lab = NULL,  y_lab = NULL,  legend_name = "auto",  palette = "red",  x_label_angle = 90,  x_label_vjust = 1,  x_label_hjust = 0.5,  y_label_angle = 0,  y_label_vjust = 0.5,  y_label_hjust = 1,  flip_xy = FALSE,  sig_names = NULL,  sig_orders = NULL,  check_sig_names = TRUE)

Arguments

Details

Support:

• SBS-24

• SBS-96

• SBS-384

• SBS-1536

• SBS-6144

• DBS-78

• DBS-186

Value

aggplot object.

Examples

# Load SBS signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Show signature profilep1 <- show_sig_profile_heatmap(sig2, mode = "SBS")p1

Show Signature Profile with Loop Way

Description

Show Signature Profile with Loop Way

Usage

show_sig_profile_loop(  Signature,  sig_names = NULL,  ncol = 1,  nrow = NULL,  x_lab = "Components",  ...)

Arguments

Value

aggplot result fromcowplot::plot_grid().

See Also

show_sig_profile

Examples

load(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Show signature profilep1 <- show_sig_profile_loop(sig2, mode = "SBS")p1p2 <- show_sig_profile_loop(sig2, mode = "SBS", style = "cosmic", sig_names = c("A", "B", "C"))p2

Extract Signatures through the Automatic Relevance Determination Technique

Description

A bayesian variant of NMF algorithm to enable optimal inferences for thenumber of signatures through the automatic relevance determination technique.This functions delevers highly interpretable and sparse representations forboth signature profiles and attributions at a balance between data fitting andmodel complexity (this method may introduce more signatures than expected,especially for copy number signatures (thusI don't recommend you to use this featureto extract copy number signatures)). See detail part and references for more.

Usage

sig_auto_extract(  nmf_matrix = NULL,  result_prefix = "BayesNMF",  destdir = tempdir(),  method = c("L1W.L2H", "L1KL", "L2KL"),  strategy = c("stable", "optimal", "ms"),  ref_sigs = NULL,  K0 = 25,  nrun = 10,  niter = 2e+05,  tol = 1e-07,  cores = 1,  optimize = FALSE,  skip = FALSE,  recover = FALSE)

Arguments

Details

There are three methods available in this function: "L1W.L2H", "L1KL" and "L2KL".They use different priors for the bayesian variant of NMF algorithm(seemethod parameter) written by reference #1 and implemented inSignatureAnalyzer software(reference #2).

I copied source code for the three methods from Broad Institute and supplementaryfiles of reference #3, and wrote this higher function. It is more friendly for usersto extract, visualize and analyze signatures by combining with other powerful functionsinsigminer package. Besides, I implemented parallel computation to speed upthe calculation process and a similar input and output structure likesig_extract().

Value

alist withSignature class.

Author(s)

Shixiang Wang

References

Tan, Vincent YF, and Cédric Févotte. "Automatic relevance determination in nonnegative matrix factorization with the/spl beta/-divergence."IEEE Transactions on Pattern Analysis and Machine Intelligence 35.7 (2012): 1592-1605.

Kim, Jaegil, et al. "Somatic ERCC2 mutations are associated with a distinct genomic signature in urothelial tumors."Nature genetics 48.6 (2016): 600.

Alexandrov, Ludmil, et al. "The repertoire of mutational signatures in human cancer." BioRxiv (2018): 322859.

See Also

sig_tally for getting variation matrix,sig_extract for extracting signatures usingNMF package,sig_estimate forestimating signature number forsig_extract.

Examples

load(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))res <- sig_auto_extract(cn_tally_W$nmf_matrix, result_prefix = "Test_copynumber", nrun = 1)# At default, all run files are stored in tempdir()dir(tempdir(), pattern = "Test_copynumber")laml.maf <- system.file("extdata", "tcga_laml.maf.gz", package = "maftools")laml <- read_maf(maf = laml.maf)mt_tally <- sig_tally(  laml,  ref_genome = "BSgenome.Hsapiens.UCSC.hg19",  use_syn = TRUE)x <- sig_auto_extract(mt_tally$nmf_matrix,  strategy = "ms", nrun = 3, ref_sigs = "legacy")x

Convert Signatures between different Genomic Distribution of Components

Description

Converts signatures between two representations relative to different sets of mutational opportunities.Currently, only SBS signature is supported.

Usage

sig_convert(sig, from = "human-genome", to = "human-exome")

Arguments

Details

The default opportunity matrix for "human-genome" and "human-exome" comes from COSMICsignature database v2 and v3.

Value

amatrix.

References

convert_signatures function from sigfit package.

Examples

# Load SBS signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))# Exome-relative to Genome-relativesig_converted <- sig_convert(sig2,  from = "human-exome",  to = "human-genome")sig_convertedshow_sig_profile(sig2, style = "cosmic")show_sig_profile(sig_converted, style = "cosmic")

Estimate Signature Number

Description

UseNMF package to evaluate the optimal number of signatures.This is used along withsig_extract.Users shouldlibrary(NMF) firstly. If NMF objects are returned,the result can be further visualized by NMF plot methods likeNMF::consensusmap() andNMF::basismap().

sig_estimate() shows comprehensive rank survey generated byNMF package, sometimesit is hard to consider all measures.show_sig_number_survey() provides aone or two y-axis visualization method to help users determinethe optimal signature number (showing bothstability ("cophenetic") and error (RSS) at default).Users can also set custom measures to show.

show_sig_number_survey2() is modified fromNMF package tobetter help users to explore survey of signature number.

Usage

sig_estimate(  nmf_matrix,  range = 2:5,  nrun = 10,  use_random = FALSE,  method = "brunet",  seed = 123456,  cores = 1,  keep_nmfObj = FALSE,  save_plots = FALSE,  plot_basename = file.path(tempdir(), "nmf"),  what = "all",  verbose = FALSE)show_sig_number_survey(  object,  x = "rank",  left_y = "cophenetic",  right_y = "rss",  left_name = left_y,  right_name = toupper(right_y),  left_color = "black",  right_color = "red",  left_shape = 16,  right_shape = 18,  shape_size = 4,  highlight = NULL)show_sig_number_survey2(  x,  y = NULL,  what = c("all", "cophenetic", "rss", "residuals", "dispersion", "evar", "sparseness",    "sparseness.basis", "sparseness.coef", "silhouette", "silhouette.coef",    "silhouette.basis", "silhouette.consensus"),  na.rm = FALSE,  xlab = "Total signatures",  ylab = "",  main = "Signature number survey using NMF package")

Arguments

Details

The most common approach is to choose the smallest rank for which cophenetic correlation coefficientstarts decreasing (Used by this function). Another approach is to choose the rank for which the plotof the residual sum of squares (RSS) between the input matrix and its estimate shows an inflection point.More custom features please directly useNMF::nmfEstimateRank.

Value

• sig_estimate: alist contains information of NMF run and rank survey.

• show_sig_number_survey: aggplot object

• show_sig_number_survey2: aggplot object

Author(s)

Shixiang Wang

References

Gaujoux, Renaud, and Cathal Seoighe. "A flexible R package for nonnegative matrix factorization." BMC bioinformatics 11.1 (2010): 367.

See Also

sig_extract for extracting signatures usingNMF package,sig_auto_extract forextracting signatures using automatic relevance determination technique.

sig_estimate for estimating signature number forsig_extract,show_sig_number_survey2 for more visualization method.

Examples

load(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))library(NMF)cn_estimate <- sig_estimate(cn_tally_W$nmf_matrix,  cores = 1, nrun = 5,  verbose = TRUE)p <- show_sig_number_survey2(cn_estimate$survey)p# Show two measuresshow_sig_number_survey(cn_estimate)# Show one measurep1 <- show_sig_number_survey(cn_estimate, right_y = NULL)p1p2 <- add_h_arrow(p, x = 4.1, y = 0.953, label = "selected number")p2# Show data from a data.framep3 <- show_sig_number_survey(cn_estimate$survey)p3# Show other measureshead(cn_estimate$survey)p4 <- show_sig_number_survey(cn_estimate$survey,  right_y = "dispersion",  right_name = "dispersion")p4p5 <- show_sig_number_survey(cn_estimate$survey,  right_y = "evar",  right_name = "evar")p5

Extract Signatures through NMF

Description

Do NMF de-composition and then extract signatures.

Usage

sig_extract(  nmf_matrix,  n_sig,  nrun = 10,  cores = 1,  method = "brunet",  optimize = FALSE,  pynmf = FALSE,  use_conda = TRUE,  py_path = "/Users/wsx/anaconda3/bin/python",  seed = 123456,  ...)

Arguments

Value

alist withSignature class.

Author(s)

Shixiang Wang

References

Gaujoux, Renaud, and Cathal Seoighe. "A flexible R package for nonnegative matrix factorization." BMC bioinformatics 11.1 (2010): 367.

Mayakonda, Anand, et al. "Maftools: efficient and comprehensive analysis of somatic variants in cancer." Genome research 28.11 (2018): 1747-1756.

See Also

sig_tally for getting variation matrix,sig_estimate for estimating signature number forsig_extract,sig_auto_extract forextracting signatures using automatic relevance determination technique.

Examples

load(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))# Extract copy number signaturesres <- sig_extract(cn_tally_W$nmf_matrix, 2, nrun = 1)

Fit Signature Exposures with Linear Combination Decomposition

Description

The function performs a signatures decomposition of a given mutationalcatalogueV with known signaturesW by solving the minimization problem⁠min(||W*H - V||)⁠ where W and V are known.

Usage

sig_fit(  catalogue_matrix,  sig,  sig_index = NULL,  sig_db = c("legacy", "SBS", "DBS", "ID", "TSB", "SBS_Nik_lab", "RS_Nik_lab",    "RS_BRCA560", "RS_USARC", "CNS_USARC", "CNS_TCGA", "CNS_TCGA176", "CNS_PCAWG176",    "SBS_hg19", "SBS_hg38", "SBS_mm9", "SBS_mm10", "DBS_hg19", "DBS_hg38", "DBS_mm9",    "DBS_mm10", "SBS_Nik_lab_Organ", "RS_Nik_lab_Organ", "latest_SBS_GRCh37",    "latest_DBS_GRCh37", "latest_ID_GRCh37", "latest_SBS_GRCh38", "latest_DBS_GRCh38",    "latest_SBS_mm9", "latest_DBS_mm9", "latest_SBS_mm10", "latest_DBS_mm10",    "latest_SBS_rn6", "latest_DBS_rn6", "latest_CN_GRCh37",         "latest_RNA-SBS_GRCh37", "latest_SV_GRCh38"),  db_type = c("", "human-exome", "human-genome"),  show_index = TRUE,  method = c("QP", "NNLS", "SA"),  auto_reduce = FALSE,  type = c("absolute", "relative"),  return_class = c("matrix", "data.table"),  return_error = FALSE,  rel_threshold = 0,  mode = c("SBS", "DBS", "ID", "copynumber"),  true_catalog = NULL,  ...)

Arguments

Details

The method 'NNLS' solves the minimization problem with nonnegative least-squares constraints.The method 'QP' and 'SA' are modified from SignatureEstimation package.See references for details.Of note, when fitting exposures for copy number signatures, only components offeature CN is used.

Value

The exposure result either inmatrix ordata.table format.Ifreturn_error setTRUE, alist is returned.

References

Daniel Huebschmann, Zuguang Gu and Matthias Schlesner (2019). YAPSA: Yet Another Package for Signature Analysis. R package version 1.12.0.

Huang X, Wojtowicz D, Przytycka TM. Detecting presence of mutational signatures in cancer with confidence. Bioinformatics. 2018;34(2):330–337. doi:10.1093/bioinformatics/btx604

Kim, Jaegil, et al. "Somatic ERCC2 mutations are associated with a distinct genomic signature in urothelial tumors."Nature genetics 48.6 (2016): 600.

See Also

sig_extract,sig_auto_extract,sig_fit_bootstrap,sig_fit_bootstrap_batch

Examples

# For mutational signatures ----------------# SBS is used for illustration, similar# operations can be applied to DBS, INDEL, CN, RS, etc.# Load simulated datadata("simulated_catalogs")data = simulated_catalogs$set1data[1:5, 1:5]# Fitting with all COSMIC v2 reference signaturessig_fit(data, sig_index = "ALL")# Check ?sig_fit for sig_db options# e.g., use the COSMIC SBS v3sig_fit(data, sig_index = "ALL", sig_db = "SBS")# Fitting with specified signatures# opt 1. use selected reference signaturessig_fit(data, sig_index = c(1, 5, 9, 2, 13), sig_db = "SBS")# opt 2. use user specified signaturesref = get_sig_db()$dbref[1:5, 1:5]ref = ref[, 1:10]# The `sig` used here can be result object from `sig_extract`# or any reference matrix with similar structure (96-motif)v1 = sig_fit(data, sig = ref)v1# If possible, auto-reduce the reference signatures# for better fitting data from a samplev2 = sig_fit(data, sig = ref, auto_reduce = TRUE)v2all.equal(v1, v2)# Some samples reported signatures dropped# but its original activity values are 0s,# so the data remain same (0 -> 0)all.equal(v1[, 2], v2[, 2])# For COSMIC_10, 6.67638 -> 0v1[, 4]; v2[, 4]all.equal(v1[, 4], v2[, 4])# For general purpose -----------------------W <- matrix(c(1, 2, 3, 4, 5, 6), ncol = 2)colnames(W) <- c("sig1", "sig2")W <- apply(W, 2, function(x) x / sum(x))H <- matrix(c(2, 5, 3, 6, 1, 9, 1, 2), ncol = 4)colnames(H) <- paste0("samp", 1:4)V <- W %*% HVif (requireNamespace("quadprog", quietly = TRUE)) {  H_infer <- sig_fit(V, W, method = "QP")  H_infer  H  H_dt <- sig_fit(V, W, method = "QP", auto_reduce = TRUE, return_class = "data.table")  H_dt  ## Show results  show_sig_fit(H_infer)  show_sig_fit(H_dt)  ## Get clusters/groups  H_dt_rel <- sig_fit(V, W, return_class = "data.table", type = "relative")  z <- get_groups(H_dt_rel, method = "k-means")  show_groups(z)}# if (requireNamespace("GenSA", quietly = TRUE)) {#   H_infer <- sig_fit(V, W, method = "SA")#   H_infer#   H##   H_dt <- sig_fit(V, W, method = "SA", return_class = "data.table")#   H_dt##   ## Modify arguments to method#   sig_fit(V, W, method = "SA", maxit = 10, temperature = 100)##   ## Show results#   show_sig_fit(H_infer)#   show_sig_fit(H_dt)# }

Obtain Bootstrap Distribution of Signature Exposures of a Certain Tumor Sample

Description

This can be used to obtain the confidence of signature exposures or searchthe suboptimal decomposition solution.

Usage

sig_fit_bootstrap(  catalog,  sig,  n = 100L,  sig_index = NULL,  sig_db = "legacy",  db_type = c("", "human-exome", "human-genome"),  show_index = TRUE,  method = c("QP", "NNLS", "SA"),  auto_reduce = FALSE,  SA_not_bootstrap = FALSE,  type = c("absolute", "relative"),  rel_threshold = 0,  mode = c("SBS", "DBS", "ID", "copynumber"),  find_suboptimal = FALSE,  suboptimal_ref_error = NULL,  suboptimal_factor = 1.05,  ...)

Arguments

Value

alist

References

Huang X, Wojtowicz D, Przytycka TM. Detecting presence of mutational signatures in cancer with confidence. Bioinformatics. 2018;34(2):330–337. doi:10.1093/bioinformatics/btx604

See Also

report_bootstrap_p_value,sig_fit,sig_fit_bootstrap_batch

Examples

# This function is designed for processing# one sample, thus is not very useful in practice# please check `sig_fit_bootstrap_batch`# For general purpose -------------------W <- matrix(c(1, 2, 3, 4, 5, 6), ncol = 2)colnames(W) <- c("sig1", "sig2")W <- apply(W, 2, function(x) x / sum(x))H <- matrix(c(2, 5, 3, 6, 1, 9, 1, 2), ncol = 4)colnames(H) <- paste0("samp", 1:4)V <- W %*% HVif (requireNamespace("quadprog", quietly = TRUE)) {  H_bootstrap <- sig_fit_bootstrap(V[, 1], W, n = 10, type = "absolute")  ## Typically, you have to run many times to get close to the answer  boxplot(t(H_bootstrap$expo))  H[, 1]  ## Return P values  ## In practice, run times >= 100  ## is recommended  report_bootstrap_p_value(H_bootstrap)  ## For multiple samples  ## Input a list  report_bootstrap_p_value(list(samp1 = H_bootstrap, samp2 = H_bootstrap))  #   ## Find suboptimal decomposition  #   H_suboptimal <- sig_fit_bootstrap(V[, 1], W,  #     n = 10,  #     type = "absolute",  #     method = "SA",  #     find_suboptimal = TRUE  #   )}

Exposure Instability Analysis of Signature Exposures with Bootstrapping

Description

Readsig_fit_bootstrap for more option setting.

Usage

sig_fit_bootstrap_batch(  catalogue_matrix,  methods = c("QP"),  n = 100L,  min_count = 1L,  p_val_thresholds = c(0.05),  use_parallel = FALSE,  seed = 123456L,  job_id = NULL,  result_dir = tempdir(),  ...)

Arguments

Value

alist ofdata.table.

See Also

sig_fit,sig_fit_bootstrap

Examples

# For mutational signatures ----------------# SBS is used for illustration, similar# operations can be applied to DBS, INDEL, CN, RS, etc.# Load simulated datadata("simulated_catalogs")data = simulated_catalogs$set1data[1:5, 1:5]# Fitting with COSMIC reference signatures# Generally set n = 100rv = sig_fit_bootstrap_batch(data,  sig_index = c(1, 5, 9, 2, 13),   sig_db = "SBS", n = 10)rv# For general purpose --------------------W <- matrix(c(1, 2, 3, 4, 5, 6), ncol = 2)colnames(W) <- c("sig1", "sig2")W <- apply(W, 2, function(x) x / sum(x))H <- matrix(c(2, 5, 3, 6, 1, 9, 1, 2), ncol = 4)colnames(H) <- paste0("samp", 1:4)V <- W %*% HVif (requireNamespace("quadprog")) {  z10 <- sig_fit_bootstrap_batch(V, sig = W, n = 10)  z10}

Obtain or Modify Signature Information

Description

Obtain or Modify Signature Information

Usage

sig_names(sig)sig_modify_names(sig, new_names)sig_number(sig)sig_attrs(sig)sig_signature(sig, normalize = c("row", "column", "raw", "feature"))sig_exposure(sig, type = c("absolute", "relative"))

Arguments

Value

aSignature object or data.

Examples

## Operate signature namesload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))sig_names(sig2)cc <- sig_modify_names(sig2, new_names = c("Sig2", "Sig1", "Sig3"))sig_names(cc)# The older names are stored in tags.print(attr(cc, "tag"))## Get signature numbersig_number(sig2)## Get signature attributessig_number(sig2)## Get signature matrixz <- sig_signature(sig2)z <- sig_signature(sig2, normalize = "raw")## Get exposure matrix## Of note, this is different from get_sig_exposure()## it returns a matrix instead of data table.z <- sig_exposure(sig2) # it is same as sig$Exposurez <- sig_exposure(sig2, type = "relative") # it is same as sig2$Exposure.norm

Tally a Genomic Alteration Object

Description

Tally a variation object likeMAF,CopyNumber and return a matrix for NMF de-composition and more.This is a generic function,so it can be further extended to other mutation cases.Please read details about how to set sex for identifying copy number signatures.Please readhttps://osf.io/s93d5/ for the generation of SBS, DBS and ID (INDEL)components.

Usage

sig_tally(object, ...)## S3 method for class 'CopyNumber'sig_tally(  object,  method = "Wang",  ignore_chrs = NULL,  indices = NULL,  add_loh = FALSE,  feature_setting = sigminer::CN.features,  cores = 1,  keep_only_matrix = FALSE,  ...)## S3 method for class 'RS'sig_tally(object, keep_only_matrix = FALSE, ...)## S3 method for class 'MAF'sig_tally(  object,  mode = c("SBS", "DBS", "ID", "ALL"),  ref_genome = "BSgenome.Hsapiens.UCSC.hg19",  genome_build = NULL,  add_trans_bias = FALSE,  ignore_chrs = NULL,  use_syn = TRUE,  keep_only_matrix = FALSE,  ...)

Arguments

Details

For identifying copy number signatures, we have to derive copy numberfeatures firstly. Due to the difference of copy number values in sex chromosomesbetween male and female, we have to do an extra stepif we don't want toignore them.

I create two options to control this, the default values are shown asthe following, you can use the same way to set (per R session).

options(sigminer.sex = "female", sigminer.copynumber.max = NA_integer_)

• If your cohort are all females, you can totally ignore this.

• If your cohort are all males, setsigminer.sex to 'male' andsigminer.copynumber.max to a proper value (the best is consistentwithread_copynumber).

• If your cohort contains both males and females, setsigminer.sexas adata.frame with two columns "sample" and "sex". Andsetsigminer.copynumber.max to a proper value (the best is consistentwithread_copynumber).

Value

alist contains amatrix used for NMF de-composition.

Methods (by class)

• sig_tally(CopyNumber): Returns copy number features, components and component-by-sample matrix

• sig_tally(RS): Returns genome rearrangement sample-by-component matrix

• sig_tally(MAF): Returns SBS mutation sample-by-component matrix and APOBEC enrichment

Author(s)

Shixiang Wang

References

Wang, Shixiang, et al. "Copy number signature analyses in prostate cancer revealdistinct etiologies and clinical outcomes." medRxiv (2020).

Steele, Christopher D., et al. "Undifferentiated sarcomas develop throughdistinct evolutionary pathways." Cancer Cell 35.3 (2019): 441-456.

Mayakonda, Anand, et al. "Maftools: efficient and comprehensive analysis of somatic variants in cancer." Genome research 28.11 (2018): 1747-1756.

Roberts SA, Lawrence MS, Klimczak LJ, et al. An APOBEC Cytidine Deaminase Mutagenesis Pattern is Widespread in Human Cancers. Nature genetics. 2013;45(9):970-976. doi:10.1038/ng.2702.

Bergstrom EN, Huang MN, Mahto U, Barnes M, Stratton MR, Rozen SG, Alexandrov LB: SigProfilerMatrixGenerator: a tool for visualizing and exploring patterns of small mutational events. BMC Genomics 2019, 20:685 https://bmcgenomics.biomedcentral.com/articles/10.1186/s12864-019-6041-2

See Also

sig_estimate for estimating signature number forsig_extract,sig_auto_extract for extracting signatures using automatic relevance determination technique.

Examples

# Load copy number objectload(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))# Use method designed by Wang, Shixiang et al.cn_tally_W <- sig_tally(cn, method = "W")# Use method designed by Steele et al.# See example in read_copynumber# Prepare SBS signature analysislaml.maf <- system.file("extdata", "tcga_laml.maf.gz", package = "maftools")laml <- read_maf(maf = laml.maf)if (require("BSgenome.Hsapiens.UCSC.hg19")) {  mt_tally <- sig_tally(    laml,    ref_genome = "BSgenome.Hsapiens.UCSC.hg19",    use_syn = TRUE  )  mt_tally$nmf_matrix[1:5, 1:5]  ## Use strand bias categories  mt_tally <- sig_tally(    laml,    ref_genome = "BSgenome.Hsapiens.UCSC.hg19",    use_syn = TRUE, add_trans_bias = TRUE  )  ## Test it by enrichment analysis  enrich_component_strand_bias(mt_tally$nmf_matrix)  enrich_component_strand_bias(mt_tally$all_matrices$SBS_24)} else {  message("Please install package 'BSgenome.Hsapiens.UCSC.hg19' firstly!")}

An Unified Interface to Extract Signatures

Description

This function provides an unified interface to signature extractorimplemented insigminer. If you determine a specificapproach,please also read the documentation of corresponding extractor.See "Arguments" part.

Usage

sig_unify_extract(  nmf_matrix,  range = 2:5,  nrun = 10,  approach = c("bayes_nmf", "repeated_nmf", "bootstrap_nmf", "sigprofiler"),  cores = 1L,  ...)

Arguments

Value

Result dependent on theapproach setting.

See Also

sig_extract,sig_auto_extract,bp_extract_signatures,sigprofiler

Examples

load(system.file("extdata", "toy_copynumber_tally_W.RData",  package = "sigminer", mustWork = TRUE))# Extract signatures# It is same as sig_extract(cn_tally_W$nmf_matrix, 2, nrun = 1)res <- sig_unify_extract(cn_tally_W$nmf_matrix, 2,  nrun = 1,  approach = "repeated_nmf")# Auto-extract signatures based on bayesian NMFres2 <- sig_unify_extract(cn_tally_W$nmf_matrix,  nrun = 1,  approach = "bayes_nmf")

Extract Signatures with SigProfiler

Description

This function provides an interface to software SigProfiler.More please seehttps://github.com/AlexandrovLab/SigProfilerExtractor.Typically, a reference genome is not required because the input is a matrix (my understanding).If you are using refitting result by SigProfiler, please make sure you have input the matrix same order as examples athttps://github.com/AlexandrovLab/SigProfilerMatrixGenerator/tree/master/SigProfilerMatrixGenerator/references/matrix/BRCA_example. If not, usesigprofiler_reorder() firstly.

Usage

sigprofiler_extract(  nmf_matrix,  output,  output_matrix_only = FALSE,  range = 2:5,  nrun = 10L,  refit = FALSE,  refit_plot = FALSE,  is_exome = FALSE,  init_method = c("random", "nndsvd_min", "nndsvd", "nndsvda", "nndsvdar"),  cores = -1L,  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  use_conda = FALSE,  py_path = NULL,  sigprofiler_version = "1.1.3")sigprofiler_import(  output,  order_by_expo = FALSE,  type = c("suggest", "refit", "all"))sigprofiler_reorder(  nmf_matrix,  type = c("SBS96", "SBS6", "SBS12", "SBS192", "SBS1536", "SBS3072", "DBS78", "DBS312",    "DBS1248", "DBS4992"))

Arguments

Value

Forsigprofiler_extract(), returns nothing. Seeoutput directory.

Forsigprofiler_import(), alist containingSignature object.

A NMF matrix for input ofsigprofiler_extract().

Examples

if (FALSE) {  load(system.file("extdata", "toy_copynumber_tally_W.RData",    package = "sigminer", mustWork = TRUE  ))  reticulate::conda_list()  sigprofiler_extract(cn_tally_W$nmf_matrix, "~/test/test_sigminer",    use_conda = TRUE  )  sigprofiler_extract(cn_tally_W$nmf_matrix, "~/test/test_sigminer",    use_conda = FALSE, py_path = "/Users/wsx/anaconda3/bin/python"  )}data("simulated_catalogs")sigprofiler_reorder(t(simulated_catalogs$set1))

A List of Simulated SBS-96 Catalog Matrix

Description

Data fromdoi:10.1038/s43018-020-0027-5.5 simulated mutation catalogs are used by the paper but only 4 are available.The data are simulated from COSMIC mutational signatures 1, 2, 3, 5, 6, 8,12, 13, 17 and 18. Each sample is a linear combination of 5 randomly selectedsignatures with the addiction of Poisson noise. The number of mutation ineach sample is randomly selected between 1,000 and 50,000 mutations, in logscale so that a lower number of mutations is more likely to be selected.The proportion of each signature in each sample is also random.

Format

A list of matrix

Source

Generate from code under data_raw/

Examples

data(simulated_catalogs)

Simulation Analysis

Description

• simulate_signature() - Simulate signatures from signature pool.

• simulate_catalogue() - Simulate catalogs from signature/catalog pool.

• simulate_catalogue_matrix() - Simulate a bootstrapped catalog matrix.

Usage

simulate_signature(x, weights = NULL)simulate_catalogue(x, n, weights = NULL)simulate_catalogue_matrix(x)

Arguments

Value

amatrix.

Examples

# Generate a catalogset.seed(1234)catalog <- as.integer(table(sample(1:96, 1000, replace = TRUE)))names(catalog) <- paste0("comp", 1:96)# Generate a signaturesig <- catalog / sum(catalog)# Simulate catalogsx1 <- simulate_catalogue(catalog, 10) # 10 mutationsx1x2 <- simulate_catalogue(catalog, 100) # 100 mutationsx2x3 <- simulate_catalogue(catalog, 1000) # 1000 mutationsx3# Similar with a signaturex4 <- simulate_catalogue(sig, 10) # 10 mutationsx4# Load SBS signatureload(system.file("extdata", "toy_mutational_signature.RData",  package = "sigminer", mustWork = TRUE))s <- t(sig2$Signature.norm)# Generate a signature from multiple signatures/catalogss1 <- simulate_signature(s)s1s2 <- simulate_signature(s, weights = 1:3)s2# Generate a catalog from multiple signatures/catalogsc1 <- simulate_catalogue(s, 100, weights = 1:3)c1

Subsetting CopyNumber object

Description

Subsetdata slot ofCopyNumber object, un-selected rows will move todropoff.segs slot, annotation slot will update in the same way.

Usage

## S3 method for class 'CopyNumber'subset(x, subset = TRUE, ...)

Arguments

Value

aCopyNumber object

Author(s)

Shixiang Wang

Tidy eval helpers

Description

• sym() creates a symbol from a string andsyms() creates a list of symbols from acharacter vector.

• enquo() andenquos() delay the execution of one orseveral function arguments.enquo() returns a single quotedexpression, which is like a blueprint for the delayed computation.enquos() returns a list of such quoted expressions.

• expr() quotes a new expressionlocally. Itis mostly useful to build new expressions around argumentscaptured withenquo() orenquos():expr(mean(!!enquo(arg), na.rm = TRUE)).

• as_name() transforms a quoted variable nameinto a string. Supplying something else than a quoted variablename is an error.

  That's unlikeas_label() which also returnsa single string but supports any kind of R object as input,including quoted function calls and vectors. Its purpose is tosummarise that object into a single label. That label is oftensuitable as a default name.

  If you don't know what a quoted expression contains (for instanceexpressions captured withenquo() could be a variablename, a call to a function, or an unquoted constant), then useas_label(). If you know you have quoted a simple variablename, or would like to enforce this, useas_name().

To learn more about tidy eval and how to use these tools, visithttps://dplyr.tidyverse.org/articles/programming.html and theMetaprogrammingsection ofAdvanced R.

Merged Transcript Location at Genome Build T2T

Description

Merged Transcript Location at Genome Build T2T

Format

Adata.table

Source

from T2T study.

Examples

data(transcript.T2T)

Merged Transcript Location at Genome Build hg19

Description

Merged Transcript Location at Genome Build hg19

Format

Adata.table

Source

from GENCODE release v33.

Examples

data(transcript.hg19)

Merged Transcript Location at Genome Build hg38

Description

Merged Transcript Location at Genome Build hg38

Format

Adata.table

Source

from GENCODE release v33.

Examples

data(transcript.hg38)

Merged Transcript Location at Genome Build mm10

Description

Merged Transcript Location at Genome Build mm10

Format

Adata.table

Source

from GENCODE release M25.

Examples

data(transcript.mm10)

Merged Transcript Location at Genome Build mm9

Description

Merged Transcript Location at Genome Build mm9

Format

Adata.table

Source

from UCSChttp://hgdownload.cse.ucsc.edu/goldenPath/mm9/database/transcriptome.txt.gz

Examples

data(transcript.mm9)

Transform Copy Number Table

Description

Transform Copy Number Table

Usage

transform_seg_table(  data,  genome_build = c("hg19", "hg38", "T2T", "mm10", "mm9", "ce11"),  ref_type = c("cytoband", "gene"),  values_fill = NA,  values_fn = function(x, ...) {     round(mean(x, ...)) },  resolution_factor = 1L)

Arguments

Value

adata.table.

Examples

load(system.file("extdata", "toy_copynumber.RData",  package = "sigminer", mustWork = TRUE))# Compute the mean segVal in each cytobandx <- transform_seg_table(cn, resolution_factor = 1)x# Compute the mean segVal in each half-cytobandx2 <- transform_seg_table(cn, resolution_factor = 2)x2

Set Color Style for Plotting

Description

Set Color Style for Plotting

Usage

use_color_style(  style,  mode = c("SBS", "copynumber", "DBS", "ID", "RS"),  method = "Wang")

Arguments

Value

color values.

Examples

use_color_style("default")use_color_style("cosmic")