Option 1: View demo data
Please wait a couple seconds after clicking and you should be redirected to the Visualize and Explore tab.
Option 2: Upload your own data Files
(Required) Please upload the aggregate report file. Note that this will be the data displayed in the main table in the Explore tab.
(Required) Please upload the corresponding metrics file for the main file that you have chosen.
(Optional) If you would like, you can upload an additional aggregate report file generated with either Class I or Class II results to supplement your main table. (E.g. if you uploaded Class I data as the main table, you can upload your Class II report here as supplemental data)
(Optional) Additionally, you can upload a gene-of-interest list in a tsv format, where each row is a single gene name. These genes (if in your aggregate report) will be highlighted in the Gene Name column.
Basic Instructions: How to explore your data using pVACview?
Step 1: Upload your own data / Load demo data
You can either choose to explore a demo dataset that we have prepared from the HCC1395 cell line, or choose to upload your own datasets.
If you are uploading your own datasets, the two required inputs are output files you obtain after running the pVACseq pipeline. The aggregated tsv file is a list of all predicted epitopes and their binding affinity scores with additional variant information and the metrics json file contains additional transcript and peptide level information.
You have the option of uploading an additional file to supplement the data you are exploring. This includes: additional class I or II information and a gene-of-interest tsv file.
Step 2: Exploring your data
To explore the different aspects of your neoantigen candidates, you will need to navigate to the Aggregate Report of Best Candidate by Variant on the visualize and explore tab. For detailed variant, transcript and peptide information for each candidate listed, you will need to click on the Investigate button for the specific row of interest. This will prompt both the transcript and peptide table to reload with the matching information.
By hovering over each column header, you will be able to see a brief description of the corresponding column and for more details, you can click on the tooltip located at the top right of the aggregate report table.
After investigating each candidate, you can label the candidate using the dropdown menu located at the second to last column of the table. Choices include:
Accept, Reject or Review.
Step 3: Exporting your data
When you have either finished ranking your neoantigen candidates or need to pause and would like to save your current evaluations, you can export the current main aggregate report using the export page.
Navigate to the export tab, and you will be able to name your file prior to downloading in either tsv or excel format. The excel format is user-friendly for downstream visualization and manipulation. However, if you plan on to continuing editing the aggregate report and would like to load it back in pVACview with the previous evaluations preloaded, you will need to download the file in a tsv format. This serves as a way to save your progress as your evaluations are cleared upon closing or refreshing the pVACview app.
Main table full column descriptions
If using pVACview with pVACtools output, the user is required to provide at least the following two files:
all_epitopes.aggregated.tsv
all_epitopes.aggregated.metrics.json
The
all_epitopes.aggregated.tsv
file is an aggregated version of the all_epitopes TSV.
It presents the best-scoring (lowest binding affinity) epitope for each variant, along with
additional binding affinity, expression, and coverage information for that epitope.
It also gives information about the total number of well-scoring epitopes for each variant,
the number of transcripts covered by those epitopes, and the HLA alleles that those
epitopes are well-binding to. Here, a well-binding or well-scoring epitope is any epitope that has a stronger
binding affinity than the
aggregate_inclusion_binding_threshold
described below. The report then bins variants into
tiers that offer suggestions about the suitability of variants for use in vaccines.
The
all_epitopes.aggregated.metrics.json
complements the
all_epitopes_aggregated.tsv
and is required for the tool's proper functioning.
Column Names : Description
ID
:
A unique identifier for the variant
HLA Alleles
:
For each HLA allele in the run, the number of this variant’s
epitopes that bound well to the HLA allele (with
lowest
or
median
mutant binding affinity <
aggregate_inclusion_binding_threshold
)
Gene
:
The Ensembl gene name of the affected gene
AA Change
:
The amino acid change for the mutation
Num Passing Transcripts
:
The number of transcripts
for this mutation that resulted in at least one well-binding peptide (
lowest
or
median
mutant binding affinity <
aggregate_inclusion_binding_threshold
)
Best Peptide
:
The best-binding mutant epitope sequence (lowest binding affinity)
prioritizing epitope sequences that resulted from a protein_coding transcript with a TSL below the
maximum transcript support level and having no problematic positions.
Best Transcript
:
Transcript corresponding to the best peptide with the lowest TSL and shortest length.
TSL
:
Transcript support level of the best peptide
Pos
:
The one-based position of the start of the mutation within the epitope sequence.
0
if the start of the mutation is before the epitope (as can occur downstream of frameshift mutations)
Prob Pos
:
If you specify problematic amino acids when running pVACseq, the number of problematic peptides within the best peptide.
Num Passing Peptides
:
The number of unique well-binding peptides for this mutation.
IC50 MT
:
Lowest
or
Median
ic50 binding affinity of
the best-binding mutant epitope across all prediction algorithms used.
IC50 WT
:
Lowest
or
Median
ic50 binding affinity of
the corresponding wildtype epitope across all prediction algorithms used.
%ile MT
:
Lowest
or
Median
binding affinity percentile rank
of the best-binding mutant epitope across all prediction algorithms used (those that provide percentile output)
%ile WT
:
Lowest
or
Median
binding affinity percentile rank of the
corresponding wildtype epitope across all prediction algorithms used (those that provide percentile output)
RNA Expr
:
Gene expression value for the annotated gene containing the variant.
RNA VAF
:
Tumor RNA variant allele frequency (VAF) at this position.
Allele Expr
:
RNA Expr * RNA VAF
RNA Depth
:
Tumor RNA depth at this position.
DNA VAF
:
Tumor DNA variant allele frequency (VAF) at this position.
Tier
:
A tier suggesting the suitability of variants for use in vaccines.
Evaluation
:
Column to store the evaluation of each variant when evaluating the run in pVACview.
Can be
Accept,
Reject
or
Review
.
How is the Tiering column determined / How are the Tiers assigned?
Tier : Criteria
Pass
:
(MT binding < binding threshold) AND allele expr filter pass AND vaf clonal filter pass
AND tsl filter pass AND anchor residue filter pass
Anchor
:
(MT binding < binding threshold) AND allele expr filter pass AND vaf clonal filter pass
AND tsl filter pass AND anchor residue filter fail
Subclonal
:
(MT binding < binding threshold) AND allele expr filter pass AND vaf clonal filter fail
AND tsl filter pass AND anchor residue filter pass
LowExpr
:
(MT binding < binding threshold) AND low expression criteria met AND allele expr filter pass
AND vaf clonal filter pass AND tsl filter pass AND anchor residue filter pass
Poor
:
Best peptide for current variant FAILS in two or more categories
NoExpr
:
((gene expr == 0) OR (RNA VAF == 0)) AND low expression criteria not met
Here we list out the exact criteria for passing each respective filter:
Allele Expr Filter:
(allele expr >= allele expr cutoff) OR (rna_vaf == 'NA') OR (gene_expr == 'NA')
VAF Clonal Filter:
(dna vaf < vaf subclonal) OR (dna_vaf == 'NA')
TSL Filter:
(TSL != 'NA') AND (TSL < maximum transcript support level)
Anchor Residue Filter:
1.
(Mutation(s) is at anchor(s)) AND
((WT binding < binding threshold) OR (WT percentile < percentile threshold))
OR
2.
Mutation(s) not or not entirely at anchor(s)
Low Expression Criteria:
(allele expr > 0) OR ((gene expr == 0) AND (RNA Depth > RNA Coverage Cutoff) AND (RNA VAF > RNA vaf cutoff))
Note that if a percentile threshold has been provided, then the
%ile MT
column is also required to be lower than
the given threshold to qualify for tiers, including Pass, Anchor, Subclonal and LowExpr.
Option 1: View NeoFox demo data
Please wait a couple seconds after clicking for the data to load.
Option 2: Upload your own neofox data files
(Required) Please upload your neofox output file. This file should be a table generated by NeoFox with the suffix “_neoantigen_candidates_annotated.tsv“
NeoFox (NEOantigen Feature toolbOX)
NeoFox (NEOantigen Feature toolbOX) is a python package that annotates a given set of neoantigen candidate sequences with relevant neoantigen features.
The tool covers neoepitope prediction by MHC binding and ligand prediction, similarity/foreignness of a neoepitope candidate sequence, combinatorial features and machine learning approaches by running a wide range of published toolsets on the given input data. For more detailed information on the specific neoantigen-related algorithms and how to generate your own NeoFox results, please refer to the link below:
Option 1: View Demo data
After clicking the "Load demo data" button, select your desired grouping and sorting parameters in the "Choose How to Visualize Data" panel and click "Visualize". Please wait a couple seconds for the data to load.
Option 2: View NeoPredPipe demo data
After clicking the "Load demo data" button, select your desired grouping and sorting parameters in the "Choose How to Visualize Data" panel and click "Visualize". Please wait a couple seconds for the data to load.
Option 3: View antigen.garnish demo data
After clicking the "Load demo data" button, select your desired grouping and sorting parameters in the "Choose How to Visualize Data" panel and click "Visualize". Please wait a couple seconds for the data to load.
Option 4: Upload your own custom data files
(Required) Please upload your TSV file.
Choose How to Visualize Data
Group peptides together by a certain feature. For example, grouping by variant would allow user to explore all proposed peptides for one variant at a time.
Order peptides by a certain feature. For example, ordering peptides by binding scores to find the best binders.
Choose what features you would like to consider for each group of peptides.
Example Neoantigen Prediction Pipelines
Vaxrank: A computational tool for designing personalized cancer vaccines
Therapeutic vaccines targeting mutant tumor antigens (“neoantigens”) are an increasingly popular form of personalized cancer immunotherapy. Vaxrank is a computational tool for selecting neoantigen vaccine peptides from tumor mutations, tumor RNA data, and patient HLA type. Vaxrank is freely available at www.github.com/openvax/vaxrank under the Apache 2.0 open source license and can also be installed from the Python Package Index.
NeoPredPipe: high-throughput neoantigen prediction and recognition potential pipeline
NeoPredPipe (Neoantigen Prediction Pipeline) is offered as a contiguous means of predicting putative neoantigens and their corresponding recognition potentials for both single and multi-region tumor samples. This tool allows a user to process neoantigens predicted from single- or multi-region vcf files using ANNOVAR and netMHCpan.
antigen.garnish.2: Tumor neoantigen prediction
Human and mouse ensemble tumor neoantigen prediction from SNVs and complex variants. Immunogenicity filtering based on the Tumor Neoantigen Selection Alliance (TESLA).
Bug reports or feature requests can be submitted on the
pVACtools Github page.
You may also contact us by email at
help@pvactools.org
.