These details are for use by our clinical analysts for variant classification and comprehensive data visualization as well as external professionals in the field. As one of St. Jude’s goals is to save kids worldwide, providing sufficient tools to clinicians and experts with no or little access to costly software for variant classification and visualization, would help them to better understand their patient’s cancer and potentially provide them with better care and treatment.
Q: Are there limits on the size of VCF files?
Uploaded files must not exceed 2 MB. If an uploaded file is larger than 2 MB, the cancer predisposition gene list filter will be automatically enabled unless you are using a custom gene list. This reduces the processing burden on the system by removing variants outside of targeted genes. You can archive the file in a .gz or .bgz if needed.
Q. What genome versions are supported?
A. Pecan PIE will accept variants mapped to either GRCh37-lite/hg19 or GRCh38. GRCh38 variants are automatically lifted over to 37, as the system uses 37 internally. The liftover process is able to compensate for strand and reference/variant allele swaps which can occur. A native hg38 version is in development, but is not yet available.
Pecan PIE only works for human data.
Q. What genes are on the curated gene list?
A. The list consists of disease-related genes, both cancer and non-cancer, see the file_download Excel spreadsheet for details. Filtering the source variants to a target list of genes reduces the processing burden on the system.
When browsing the results the view may be filtered to disease sub-categories of interest.
You can also specify your own custom list of genes to process when submitting your VCF file (see the advanced options panel).
Q. Why is the classification column blank in my results?
Q. This column displays the classification assigned by the St. Jude Germline Committee reviewers. If a variant was not classified by this committee before, this field will be blank.
Pecan PIE provides classifications from the Medal Ceremony pipeline, which may assign variants gold, silver, or bronze medals. An “Unknown” medal may be assigned for non-disease-predisposition genes, variants present in the ExAC (ex-TCGA) database at an allele frequency > 0.1%, or variants without functional annotations (which includes most silent variants).
Q. What do the medals mean?
A. The medal column is a rough indicator of the likelihood of the variant being clinically significant as predicted by the medal ceremony software. Variants with gold medals are most likely to be significant, and those with no medal are least likely. More details can be found in the Analysis of Results <results> section.
Q. Why are some of my variants missing?
A. Currently only coding and splice-related variants in disease-related genes make it to the medaling process. Intergenic, intronic, and UTR variants are excluded, as are those in non-coding transcripts.
Q. Why does the ExAC allele frequency shown differ from the ExAC portal?
A. The reported ExAC frequency may differ for several reasons:
- PIE uses the TCGA-subtracted distribution of ExAC rather than the main distribution.
- PIE reports the primary allele frequencies in the ExAC database, specifically the AC, AN, and AF fields from the VCF distribution. The ExAC portal appears to use the “adjusted” frequencies which may be different.
Q. Is Pecan PIE free?
A. Pecan PIE is free for non-commercial use. St. Jude covers the cost of running the pipeline and hosting. DNANexus accounts are required to keep track of your jobs in the cloud so that you can retrieve and manage from multiple locations. Accounts also make it possible to alert you of job completion via email.
Frequently asked questions
If you have any questions not covered here, feel free to reach out on our contact form.
Q: Which files are supported?
PIE works with variants in VCF format:
- Uploaded files must be compliant with the VCF specification.
- VCF files may be either uncompressed, or compressed with
bgzip
only.bgzip
is part of the htslib/tabix packages (see below).
Improperly formatted VCF files will not work with PIE. Some common problems include:
- Missing header line
- Missing required columns
- Files were compressed by gzip, zip, or any method other than the required bgzip
To verify compatibility of your VCF you can try one of these methods:
- Compressing your VCF with bgzip and indexing it with tabix, both programs from the HTSlib package (some systems also use the earlier, pre-HTSlib “tabix” package). This process will only succeed for compliant VCF files, and can help diagnose failures.
- Running “vcf-validator” program from the vcftools package.
While the VCF specification also requires that variants be sorted by chromosome name and position, PIE is now often able to automatically correct sorting issues in uploaded files. PIE requires sorted data in order to query data for targeted genes.