Cancer genomes are altered by various mutational processes and, like palimpsests, bear the signatures of these successive processes. The Palimpsest R package provides a complete workflow for the characterisation and visualisation of mutational signatures, including their evolution along tumour development. The package includes a wide range of functions for extracting single base substitution (SBS), double base substitution (DBS) and indel mutational signatures as well as structural variant (SV) signatures. Palimpsest estimates the probability of each mutation being due to each signature, which allows the clonality of each alteration to be calculated, and the mechanism at the origin of each driver event to be predicted. In short, Palimpsest is an easy-to-use toolset for the reconstruction of the natural history of tumours using whole exome or whole genome sequencing data.
Install from the GitHub repository using devtools:
install.packages("devtools")
library(devtools)
devtools::install_github("FunGeST/Palimpsest")
To add indel mutation categories we use a python script that is embedded in the R function annotate_VCF(). For this to work the function must be run in a Unix environment (i.e. Mac or Linux) with python 2.7 installed. The other aspects of the annotate_VCF() function, and indeed all other functions, work on a Windows operating system. The indel aspect of this function also requires you to have a FASTA file compatable with the input VCF genome (including position and chromosome names) accessible in your local environment. If you only wish to work on SBS/DBS/SV signatures you can skip this step.
The R package bedr is required to perform structural variant signature analysis. The bedr API gives access to "BEDTools" and offers additional utilities for genomic region processing. To gain the functionality of bedr package you will need to have the BEDTools program installed and in your default PATH.
For the analysis of somatic mutations Palimpsest requires one mandatory input file -- a mutational catalogue file (vcf) describing somatic mutations in the tumour series. For the analysis of the clonality of somatic mutations (optional), a further two files are required -- a copy number alteration file (cna_data) providing genome-wide absolute copy number estimates, and a minimal sample annotation file (annot_data) indicating gender and tumour purity.
The input files should have the following columns (the header is required, but the order of the columns can change). Example input files are provided with the package.
1]. VCF : somatic mutation data
Sample: Sample identifier. Any alphanumeric string.Type: Mutation type [SNV (e.g. C > A), INS for insertions (e.g. C > CAAA), DEL for deletions (e.g. CTAC > C)]. Although Palimpsest has double base substitution (DBS) extraction and analysis capacaties, DBS mutations are encoded by 2 consecutive lines of SNVs and so their type should remain as such.CHROM: Chromosome. Between chr1 and chr22 or the chrX or chrY ('chr' prefix required).POS: Mutation position. A positive integer.REF: Reference base(s): Each base must be one of A,C,G,T (upper case). Multiple bases are permitted for deletions only, where the value in the POS field refers to the position of the first base in the string.ALT: Alternate base(s): Each base must be one of A,C,G,T (upper case). Multiple bases are permitted for insertions only, where the value in the POS field refers to the position of the first base in the string.Tumor_Varcount: Number of variant bases at the position in the tumour sample.Tumor_Depth: Tumour sample sequencing depth at the position.Normal_Depth: Normal sample sequencing depth at the position.Gene_Name: OPTIONAL column for representing mutated gene name.Driver: OPTIONAL column indicating the driver events to be annotated in tumour history plots.
Optional:
2]. cna_data: copy number alteration data
Sample: Sample identifier. Any alphanumeric string.CHROM: Chromosome. Between chr1 and chr22 or the chrX or chrY ('chr' prefix required).POS_START: Start position of segmented chromosome.POS_END: End position of segmented chromosome.LogR: LogR information.Nmin: Minor allele copy number.Nmaj: Major allele copy number.ntot: Total copy number of segmented chromosome.Ploidy: Tumour ploidy.
3]. annot_data: sample annotation data
Sample: Sample identifier. Any alphanumeric string.Gender: Gender information for patient [M/F].Purity: Tumour purity estimate (Represented as fraction; ranging between 0.01 - 1).
4]. sv_data: structural variant data
Sample: Sample identifier. Any alphanumeric string.Type: Type of structural variant: INV/DEL/DUP/BND.CHROM_1: Chromosome of the first breakpoint. Between chr1 and chr22 or the chrX or chrY ('chr' prefix required).POS_1: Position of the first breakpoint. A positive integer.CHROM_2: Chromosome of the second breakpoint. Between chr1 and chr22 or the chrX or chrY ('chr' prefix required).POS_2: Position of the second breakpoint. A positive integer.Tumor_Varcount: Column for variant allele count information.Tumor_Depth: Column for tumour sequencing depth information.Normal_Depth: Column for normal sequencing depth information.Driver: OPTIONAL column indicating the driver events to be annotated in tumour history plots.
- The RUNNING_PALIMPSEST_EXAMPLE folder contains an example dataset and an R script of a typical Palimpsest analysis using this data. Please try!
- Introduction to Palimpsest provides a comprehensive example of the Palimpsest workflow with detailed explanations of each function.
- Please refer to the following paper for extensive description of the statistical methods used in the package: Letouzé, E., Shinde, J., et al. (2017) Mutational signatures reveal the dynamic interplay of risk factors and cellular processes during liver tumorigenesis. Nat. Commun.
Shinde, J. et al. (2018) Palimpsest: an R package for studying mutational and structural variant signatures along clonal evolution in cancer. Bioinformatics.
Figure 1. (A) Workflow illustrating a typical analysis with Palimpsest. Taking as input somatic mutations, copy-number alterations (CNAs) and structural variants, the package classifies variants as clonal and subclonal, extracts mutational and structural variant signatures separately in early clonal and late subclonal events, and estimates the probability of each alteration being due to each process. The timing of chromosome duplications is also estimated from the ratio of duplicated/non-duplicated mutations to reconstruct the complete natural history of the tumour. (B) Example of output representing, for one tumour, the number of clonal and subclonal mutations, their distribution per mutation signature, the driver alterations (colored according to the most likely causal mutational process) and CNA timing.
Copyright (C) 2019 Benedict Monteiro & Jayendra Shinde
Palimpsest is a free software: you can redistribute it and/or modify it under the terms of the GNU General Public License as published by the Free Software Foundation, either version 3 of the License, or (at your option) any later version.
This program is distributed in the hope that it will be useful, but WITHOUT ANY WARRANTY; without even the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU General Public License for more details.
You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses/.