Variant Effect Predictor Input form
When you reach the VEP web interface, you will be presented with a form to enter your data and alter various options.
- First select the correct species for your data. Ensembl hosts many vertebrate genomes; genomes for plants, protists and fungi can be found at Ensembl Genomes.
- You can optionally choose a name for the data you upload - this can make it easier for you to identify jobs and files that you have uploaded to the VEP at a later point.
You have three options for uploading your data:
- File upload - click the "Choose file" button and locate the file on your system
- Paste file - simply copy and paste the contents of your file into the large text box
- File URL - point the VEP to a file hosted on a publically accessible address. This can be either a http:// or ftp:// address.
The format of your data is automatically detected; see the examples or the input format documentation.
- For pasted data you can get an instant preview of the results of your first variant by clicking the button that appears when you paste your data. This quickly shows you the consequence type, the IDs of any overlapping variants, genes, transcripts and regulatory features, as well as SIFT and PolyPhen predictions. To see the full results set submit your job as normal.
For some species you can select which transcript database to use. The default is to use Ensembl transcripts, which offer the most rich annotation through VEP.
GENCODE Basic is a subset of the GENCODE gene set, and is intended to provide a simplified, high-quality subset of the GENCODE transcript annotations that will be useful to the majority of users. GENCODE Basic includes all genes in the GENCODE gene set, with a representative subset of the transcripts (splice variants).
You can also select to use RefSeq transcripts from the otherfeatures database; note though that these transcripts are simply aligned to the reference genome and the database is missing much of the annotation found when using the main Ensembl database (e.g. protein domains, CCDS identifiers).
- Gene symbol
Add the Ensembl protein identifer (ENSP). Equivalent to --protein.
Variants and frequency data
Find co-located known variants - report known variants from the Ensembl Variation database that overlap with your input. A list of variant sources imported can be viewed here. Note that this feature is only available for species with an Ensembl Variation database. Equivalent to --check_existing.
VEP will by default compares the alleles of your input variant to that of the existing variant; VEP will only report the existing variant ID if none of the alleles in your input variant are novel.
For example, if your input variant has alleles A/G, and the existing variant has alleles A/T, then the existing variant will not be reported. If instead your input variant has alleles A/T, then the existing variant will be reported.
To disable this allele matching, select the option "Yes but don't compare alleles" for the option "Find co-located known variants".
For known variants VEP can also provide PubMed IDs of publications citing the variant (equivalent to --pubmed).
- Frequency data for co-located variants
VEP can also report allele frequency (AF) data for existing variants from several major genotyping projects, the 1000 Genomes Project, the NHLBI-ESP and gnomAD; this only applies when you have selected human as your species.
- 1000 Genomes global - the combined phase 3 population (i.e. all individuals from all populations). Equivalent to --af
- 1000 Genomes continental - the four continent-level populations - AFR (African), AMR (American), ASN (Asian) and EUR (European). Equivalent to --af_1kg
- ESP - AA (African American) and EA (European American) populations. Equivalent to --af_esp
- gnomAD - combined, AFR, AMR, ASJ, EAS, FIN, NFE, OTH, SAS populations. Equivalent to --af_gnomad
- PubMed IDs for citations of co-located variants
Report the PubMed IDs of any publications that cite the co-located variant(s).
- Include flagged variants
Variants flagged as failed by the Ensembl Variation quality control.
- Transcript biotype
- Exon and intron numbers
Report the exon or intron number that a variant falls in as NUMBER / TOTAL, i.e. exon 2/5 means the variant falls in the 2nd of 5 exons in the transcript. Equivalent to --numbers.
- Transcript support level
- Identify canonical transcripts
- Upstream/Downstream distance (bp)
Change the distance to assign the upstream and downstream consequences. Equivalent to --distance.
- miRNA structure
Determines where in the secondary structure of a miRNA a variant falls. Equivalent to the VEP plugin miRNA.
- Protein domains
- Get regulatory region consequences
In addition to predicting consequences with overlapping transcripts, VEP can find overlaps with known regulatory regions as determined in the Ensembl Regulatory build.
Using this option, VEP will also report if a variant falls in a transcription factor binding motif, and give a score that reflects whether the altered motif sequence is more or less similar to the consensus.
Get regulatory consequences is equivalent to --regulatory.
- SIFT predictions
SIFT predicts whether an amino acid substitution affects protein function based on sequence homology and the physical properties of amino acids. Only available in popular species. For both SIFT and PolyPhen VEP can report either a score between 0 and 1, a prediction in words, or both. Equivalent to --sift.
- PolyPhen predictions
PolyPhen is a tool which predicts possible impact of an amino acid substitution on the structure and function of a human protein using straightforward physical and comparative considerations. Equivalent to --polyphen.
Provides a rank of genic intolerance and consequent susceptibility to disease based on the ratio of Loss-of-function (LoF)to synonymous mutations. Equivalent to the VEP plugin LoFtool.
Looks up the BLOSUM 62 substitution matrix score for the reference and alternative amino acids predicted for a missense mutation. Equivalent to the VEP plugin Blosum62.
- Ancestral allele
- By frequency
Filter variants by minor allele frequency (MAF). Two options are provided:
- Exclude common variants
Filter out variants that are co-located with an existing variant that has a frequency greater than 0.01 (1%) in the 1000 Genomes global population. Equivalent to --filter_common in the VEP script.
- Advanced filtering
Enabling this option allows you to specify a population and frequency to compare to, as well whether matching variants should be included or excluded from the results.
- Exclude common variants
- Return results for variants in coding regions only
Exclude variants that don't fall in a coding region of a transcript. Equivalent to --coding_only.
- Restrict results
For many variants VEP will report multiple consequence types - typically this is because the variant overlaps more than one transcript. For each of these options VEP uses consequence ranks that are subjectively determined by Ensembl. This table gives all of the consquence types predicted by Ensembl, ordered by rank. Note that enabling one of these options not only loses potentially relevant data, but in some cases may be scientifically misleading. Options:
- Show one selected consequence
Pick one consequence type across all those predicted for the variant; the output will include transcript- or feature-specific information. Consequences are chosen by the canonical, biotype status and length of the transcript, along with the ranking of the consequence type according to this table. This is the best method to use if you are interested only in one consequence per variant. Equivalent to --pick.
- Show one selected consequence per gene
Pick one consequence type for each gene using the same criteria as above. Note that if a variant overlaps more than one gene, output for each gene will be reported. Equivalent to --per_gene.
- Show only list of consequences per variant
Give a comma-separated list of all observed consequence types for each variant. No transcript-specific or gene-specific output will be given. Equivalent to --summary.
- Show most severe per variant
Only the most severe of all observed consequence types is reported for each variant. No transcript-specific or gene-specific output will be given. Equivalent to --most_severe.
- Show one selected consequence
Once you have clicked "Run", your input will be checked and submitted to the VEP as a job. All jobs associated with your session or account are shown in the "Recent Tickets" table. You may submit multiple jobs simultaneously.
The "Jobs" column of the table shows the current status of the job.
- Queued - your job is waiting to be submitted to the system
- Running - your job is currently running
- Done - your job is finished - click the [View results] link to be taken to the results page
- Failed - there is a problem with your job - click the magnifying glass icon to see more details
You may delete a job by clicking the trash can icon . If you are logged in to Ensembl, you can save the job by clicking the save icon .
You may also resubmit a job (for example, to re-run with the same data but change some parameters) by clicking the edit icon .
You can see a summary of the options that you selected for your VEP job by clicking on the magnifying glass icon .