2.0.2

This is a bugfix release. It fixes the following problem(s):

• There was a bug in the aggregate report creation. When parsing the all_epitopes file as input to the report creation any N/A amino acid changes would be parsed as null values. This release fixes this issue so that these amino acid change values are correctly parsed as string values.
• This release removes support for Python 3.5 since that Python version has reached end of life and is no longer supported.

2.0.1

This is a bugfix release. It fixes the following problem(s):

• NetMHCstabpan and NetCons have moved to a new server resulting in no results being returned from the old server URL. This results in empty filtered.tsv report files when either the --netmhc-stab or --net-chop-method option were enabled. This release fixes our usage of these tools to use the new server URL.

1.5.14

This is a bugfix release. It fixes the following problem(s):

• NetMHCstabpan and NetCons have moved to a new server resulting in no results being returned from the old server URL. This results in empty filtered.tsv report files when either the --netmhc-stab or --net-chop-method option were enabled. This release fixes our usage of these tools to use the new server URL.

2.0.0

This version adds the following features, outlined below. Please note that pVACtools 2.0 is not backwards-compatible and certain changes will break old workflows.

Breaking changes

• pVACtools now supports variable epitope lengths for class II prediction algorithms. The previous option --epitope-length (-e) no longer exists. It has been replaced with --class-i-epitope-length (-e1) and --class-ii-epitope-length (-e2) for class I and class II epitope lengths, respectively. The defaults are [8, 9, 10, 11] and [12, 13, 14, 15, 16, 17, 18], respectively.
• The --peptide-sequence-length option has been removed. The peptide sequence length is now determined by the epitope length(s) to determine the flanking sequence length before and after the mutation.
• pVACtools no longer depends on conda. pVACtools remains compatible with Python 3.5 and above but users may chose any environment manager to set up an appropriate Python environment.
• When using standalone IEDB, pVACtools is now only compatible with IEDB 3.1 and above. Please see Installation for instructions on installing the latest IEDB version.
• pVACseq is no longer dependent on annotations with the VEP Downstream plugin. This dependency has been replaced with the VEP Frameshift plugin. This requires changes to your existing VEP installation in order to install the Frameshift plugin. Existing VCFs that were previously annotated to work with pVACtools 1.5 and below will no longer work with version 2.0 and above and will need to be reannotated. Please see our documentation on Annotating your VCF with VEP for more information.
• The filtered.condensed.tsv report has been removed and replaced with the all_epitopes.aggregated.tsv report. We believe that this new report will provide a more useful summary of your results. Please see the Output Files sections of each tool for more information on this new report.

New features

• pVACtools now provides binding affinity percentile rank information, in addition to the raw ic50 binding affinity values. Users may filter on the percentile rank by using the new --percentile-threshold argument.
• Users now have the option of calculating the reference proteome similarity of their filtered epitopes. For this, the peptide sequence for the remaining variants is mapped to the reference proteome using BLAST. Variants where this yields a hit to a reference proteome are marked accordingly and a .reference_matches file provides more information about the matches. This option can be enabled using the --run-reference-proteome-similarity option.
• Users may now use the options all, all_class_i, or all_class_ii instead of specific prediction algorithms in order to run all prediction algorithms, all class I prediction algorithms, or all class II prediction algorithms, respectively.
• For successful pVACvector runs, we now output a _results.dna.fa file with the most likely nucleic acid sequence for the predicted vector.

Minor Updates

• When running pVACseq with a proximal variants VCF we would previously assume that your ran VEP with the --pick option and only process the first transcript annotation for a variant. With this update we will now associate the correct transcript for a proximal variant with the matching transcript of the main somatic variant of interest.
• The pvacseq generate_protein_fasta command now allows users to provide a proximal variants VCF using the --phased-proximal-variants-vcf option.
• The pvacseq generate_protein_fasta command now supports multi-sample VCFs. Users may use the --sample-name to provide the sample name of the sample they wish to process.
• pVACseq and pVACfuse would previously error out if the intermediate TSV parsed from the input was empty. In 2.0 the tool will no longer error out but exit with an appropriate message.
• pVACvector would previously error out when no valid path was found. In 2.0 pVACvector will not longer error out but exit with an appropriate message.
• We now set consistent file permissions on all output files.
• We’ve updated our license to BSD 3-Cause Clear. Please note that the individual licenses of our dependent tools remain in place. These can be viewed by on the Tools Used By pVACtools page.

1.5.13

This is a hotfix release. It fixes the following issues:

• MHCflurry would previously be called once per peptide sequence. Because a large overhead for MHCflurry is the creation of the model and the model gets created every time MHCflurry is called, this would cause very long runtimes when using this prediction algorithm. This version updates the processing to pre-calculate all the epitopes of an intermediate fasta file and make only one call to MHCflurry for all of them.
• This version was updated to set consistent file permissions for all of the output files created by pVACseq, pVACfuse, pVACbind, and pVACvector.

1.5.12

This is a hotfix release. It fixes the following issues:

• Vaxrank was pinned to an older version because newer versions made backwards-incompatible changes to some code that pVACseq was using. However this was causing installation issues since this older version of the vaxrank package has an indirect dependency on an old version of pysam. We updated the usage of this module so that we could use the latest version of vaxrank and, thus, newer versions of pysam.
• This version adds error handling for when a normal sample name is provided to pVACseq but the input VCF is a single-sample VCF.

1.5.11

This is a hotfix release. It fixes the following issues:

• The standalone pvacbind top_score_filter command would throw an error because it wasn’t set up correctly. This has now been fixed.
• The standalone pvacfuse generate_protein_fasta would fail when run with the --input-tsv option because it wasn’t able to associate TSV entries with the fasta entries correctly. Using this option will now correctly limit the output to only entries from the input TSV file.
• In certain situation the trimming of problematic peptides in pVACvector would not work correctly. This issue has now been addressed.

1.5.10

This is a hotfix release. It fixes the following issues:

• A stray character at the end of one file was causing a syntax error in Python 3.8. The character has been removed. pVACtools should now be 3.8 compatible although some dependencies might not be compatible yet.

1.5.9

This is a hotfix release. It fixes the following issues:

• Some variant consequences supported by pVACseq would not actually result in a amino acid change (e.g., inframe_insertion&incomplete_terminal_codon_variant). These types of variants were not being filtered out correctly and would cause an error when trying to create the peptide fasta sequences. This issue has been fixed and these variants are now being filtered out upstream.
• Running pVACseq on a non-human VCF would cause an error in the top score filter since the transcript identifiers would not match the expected format. This has been fixed and the top score filter now supports non-human transcripts.
• There was an error in setting up the standalone pvacfuse generate_protein_fasta command which would result in an error when trying to run this command. This has now been fixed.

1.5.8

This is a hotfix release. It fixes the following issues:

• The pvacseq run, pvacfuse run, and pvacbind run commands would previously error out when running with both MHC class I and MHC class II algorithms but one or the other would not produce an all_eptiopes.tsv file. This version fixes this bug.
• MHCflurry version 1.6.0 and higher changed the output file headers. This would cause errors when trying to parse these output files. pVACtools now supports both the old and the new headers.
• AGFusion updated their output file naming convention in newer versions and is now outputting .exons.csv files instead of .exons.txt files. pVACfuse is now able to process either version.