MainINDELClass2: Calculate Neoantigen Candidates on INDELs for MHC Class2

Description Usage Arguments Value

View source: R/MainINDELClass2.R

Description

Calculate Neoantigen Candidates on INDELs for MHC Class2

Usage

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
MainINDELClass2(input_file, hla_file, file_name_in_hla_table = input_file,
  refflat_file = paste(hmdir, "lib/refFlat.txt", sep = "/"),
  refmrna_file = paste(hmdir, "lib/refMrna.fa", sep = "/"), hmdir = getwd(),
  job_id = "NO_job_id", export_dir = paste("result", file_name_in_hla_table,
  job_id, "INDEL", sep = "."), rnaexp_file = NA, rnabam_file = NA,
  cnv_file = NA, purity = 1, netMHCIIpan_dir = paste(hmdir,
  "lib/netMHCIIpan-3.1/netMHCIIpan", sep = "/"), refdna_file = NA,
  samtools_dir = NA, bcftools_dir = NA, chr_column = NA,
  mutation_start_column = NA, mutation_end_column = NA,
  mutation_ref_column = NA, mutation_alt_column = NA, nm_id_column = NA,
  depth_normal_column = NA, depth_tumor_column = NA,
  ambiguous_between_exon = 0, ambiguous_codon = 0, peptide_length = c(15))

Arguments

input_file

(Required) An input vcf file annotated by,

e.g., ANNOVAR (http://annovar.openbioinformatics.org/en/latest/) or other softwares.

See by data(sample_vcf); sample_vcf;

hla_file

(Required) A tab separated file indicating HLA types. The 1st column is input_file name, and the following columns indicate HLA types.

See by data(sample_hla_table_c2); sample_hla_table_c2;

file_name_in_hla_table

If the name (1st column) in HLA table is not the same as input_file, indicate the corresponding name (Default=input_file).

refflat_file

refFlat file to be used in constructing peptide. (Default=paste(hmdir, "lib/refFlat.txt", sep="").

See "https://github.com/hase62/Neoantimon"

refmrna_file

refMrna file to be used in constructing peptide (Default=paste(hmdir, "lib/refMrna.fa", sep="").

See "https://github.com/hase62/Neoantimon"

hmdir

Home directory for the analysis (Default = getwd()).

job_id

Job-Id to be attached in output files (Default = "NO_job_id").

export_dir

The directory will be stored results (Default = "paste("result", file_name_in_hla_table, job_id, sep=".")")

rnaexp_file

A file including RNA expressions (Default=NA). The 1st, 2nd and 3rd columns are "GeneSymbol Chr:Exonstart-Exonend (locus) ExpressionAmount", respectively. The 1st row should be any header.

See by data(sample_rna_exp); sample_rna_exp;

rnabam_file

RNA bam file to calculate variant allele frequency of RNA at each mutation (Default=NA).

cnv_file

A file including copy number variation to calculate cancer cell fraction probability (CCFP) (Default=NA). The format is according to ASCAT output files. The columns are "SNPName Chromosome Position LogR segmentedLogR BAF segmentedBAF CopyNumber MinorAllele RawCopyNumber" The 1st row should be the above header.

See data(sample_copynum); sample_copynum;

purity

Tumor purity or tumor contents ratio required to calculate CCFP (Default=1).

netMHCIIpan_dir

The file directory to netMHCpan (Default="lib/netMHCIIpan-3.1/netMHCpan").

refdna_file

refdna_file information to be used to calculate RNA VAF (Default=NA).

See "https://github.com/hase62/Neoantimon"

samtools_dir

The file directory to samtools_0_x_x (Default="samtools"). It shouled be indicated when you indicate RNA-bam and try to calculate RNA VAF.

bcftools_dir

The file directory to netMHCpan (Default="bcftools"). It shouled be indicated when you indicate RNA-bam and try to calculate RNA VAF. samtools 0_x_x includes bcftools in the directory.

chr_column

The column number describing chromosome number in input_file (Default=NA, but will automatically search "Chr" in header).

mutation_start_column

The column number describing mutation start Position in input_file (Default=NA, but will automatically search "Start" in header) .

mutation_end_column

The column number describing mutation end Position in input_file (Default=NA, but will automatically search "End" in header).

mutation_ref_column

The column number describing mutation Ref in input_file (Default=NA, but will automatically search "Ref" in header).

mutation_alt_column

The column number describing mutation Alt in input_file (Default=NA, but will automatically search "Alt" in header).

nm_id_column

The column number describing NM IDs in input_file such as

"SLCO1C1:NM_001145944:exon7:c.692_693insG:p.L231fs" (Default=NA).

depth_normal_column

The column number describing the read count from normal cells (Default = NA).

depth_tumor_column

The column number describing the read count from tumor cells (Default = NA).

ambiguous_between_exon

The maximum number to permit the differences between Exon-Lengths from refFlat and refMrna (Default=0).

ambiguous_codon

The maximum number to permit the differences between inputfile- and refMrna-oriented translation start/end position (Default=0).

peptide_length

Peptide Length to be generated (Default = 15 in HLA Class2).

Value

void (Calculated Neoantigen Files will be generated as .tsv files.):

HLA: HLA type used to calculate neoantigen.

Pos: The position of a fraction of peptide used to be evaluated from the full-length peptide.

Gene: Gene symbol used to be evaluated in NetMHCpan.

Evaluated_Mutant_Peptide_Core: The core peptide of the mutant peptide to be evaluated in NetMHCpan.

Evaluated_Mutant_Peptide: The mutant peptide to be evaluated.

Mut_IC50: IC50 value for evaluated mutant peptide.

Mut_Rank: Rank value for evaluated mutanat peptide.

Chr: Chromosome Number of the mutation.

NM_ID: NM_ID used to construct peptides from the mutation.

Change: The annotation to be described in .vcf file.

Ref: reference type nucleic acid base.

Alt: alternative type nucleic acid base.

Prob: A probability of reference nucleic acid base described in .vcf file.

Mutation_Prob: A probability of alternative nucleic acid base described in .vcf file.

Exon_Start: The exon start position of the corrsponding NM_ID.

Exon_End: The exon end position of the corrsponding NM_ID.

Mutation_Position: The mutation position of the corrsponding NM_ID.

Total_Depth: The sum depth of the reference and alternative nucleic acid base.

Tumor_Depth: The depth of the alternative nucleic acid base.

Wt_Peptide: The full-length of the wild-type peptide.

Mutant_Peptide: The full-length of the mutant peptide.

Total_RNA: The expression amount of the corresponding RNA.

Tumor_RNA_Ratio: The variant allele frequency of the corresponding RNA.

Tumor_RNA: The modified amount of the corresponding RNA level based on RNA Reads.

Tumor_RNA_based_on_DNA: The modified amount of the corresponding RNA level based on DNA Reads.

MutRatio: The mean value of the cancer cell fraction probability.

MutRatio_Min: The 1% percentile of the cancer cell fraction probability.

MutRatio_Max: The 99% percentile of the cancer cell fraction probability.


hase62/Neoantimon documentation built on May 7, 2018, 6:34 p.m.