Tag: VEP

An FGFR2 mutation as the potential cause of a new phenotype including early-onset osteoporosis and bone fractures: a case report | BMC Medical Genomics

Anamnesis vitae A 13 year old male born was as result of the VII pregnancy, from unrelated parents. Other pregnancies resulted in: I-II silent miscarriage in the second trimester; III – female, born in 2003 (III-3 Fig. 1) that has the following phenotypic features: genu valgum, hip dysplasia, combined thoracolumbar scoliosis,…

Continue Reading An FGFR2 mutation as the potential cause of a new phenotype including early-onset osteoporosis and bone fractures: a case report | BMC Medical Genomics

Annotate variants with ensembl rest api

Annotate variants with ensembl rest api 0 I have a variant file (.vcf.gz), and I want to annotate this file using the Ensembl Rest API, particularly the Vep Rest API. I am new to this variant annotation; however, I have seen a couple of codes from the Ensembl page on…

Continue Reading Annotate variants with ensembl rest api

Genomics England hiring PhD Bioinformatics Intern in London, England, United Kingdom

Company DescriptionGenomics England partners with the NHS to provide whole genome sequencing diagnostics. We also equip researchers to find the causes of disease and develop new treatments – with patients and participants at the heart of it all. Our mission is to continue refining, scaling, and evolving our ability to…

Continue Reading Genomics England hiring PhD Bioinformatics Intern in London, England, United Kingdom

PhD Bioinformatics Intern Job in Greater London, Pharmaceuticals & Life Sciences Career, Intern/Graduate Jobs in Genomics England

Company Description Genomics England partners with the NHS to provide whole genome sequencing diagnostics. We also equip researchers to find the causes of disease and develop new treatments – with patients and participants at the heart of it all. Our mission is to continue refining, scaling, and evolving our…

Continue Reading PhD Bioinformatics Intern Job in Greater London, Pharmaceuticals & Life Sciences Career, Intern/Graduate Jobs in Genomics England

Does GNOMAD use all LOFTEE LoF filters?

Does GNOMAD use all LOFTEE LoF filters? 0 Hi all, I have some lof variants and I want to know if they have already been detected in GNOMAD (I am essentially curious if my lof variants are novel). In order to make this comparison to GNOMAD I have run my…

Continue Reading Does GNOMAD use all LOFTEE LoF filters?

Intradermal DNA vaccine delivery using vacuum-controlled, needle-free electroporation

Intradermal delivery of DNA vaccines via electroporation (ID-EP) has shown clinical promise, but the use of needle electrodes is typically required to achieve consistent results. Here, delivery of a DNA vaccine targeting the Middle East Respiratory Syndrome Coronavirus (MERS-CoV) is achieved using noninvasive intradermal vacuum-EP (ID-VEP), which functions by pulling…

Continue Reading Intradermal DNA vaccine delivery using vacuum-controlled, needle-free electroporation

Loftee no splice site annotations

Loftee no splice site annotations 1 Hello! I am using Loftee in my VEP pipeline and after some fights with my code everything works now, but the splice site annotations…meaning that I dont get them. There is no error at all, but my vcf files do not contain a single…

Continue Reading Loftee no splice site annotations

IJMS | Free Full-Text | Down-Regulation of Rice Glutelin by CRISPR-Cas9 Gene Editing Decreases Carbohydrate Content and Grain Weight and Modulates Synthesis of Seed Storage Proteins during Seed Maturation

3.1. Atypical Morphology of Protein Bodies in the Endosperm of Seeds from Glutelin Gene-Edited Transgenic Rice ER-derived PB-I vacuoles begin to form in response to the aggregation of β-/γ-zeins (maize) and 10-kDa prolamins (rice) and then expand by co-aggregation of α-/δ-zeins (maize) and 13-/16-kDa prolamins (rice) [3,22,23]. Saito’s group reported…

Continue Reading IJMS | Free Full-Text | Down-Regulation of Rice Glutelin by CRISPR-Cas9 Gene Editing Decreases Carbohydrate Content and Grain Weight and Modulates Synthesis of Seed Storage Proteins during Seed Maturation

[maftools]Too many multi_hit and missense mutation

[maftools]Too many multi_hit and missense mutation 0 Describe the issue When using maftools to plot mutational summary data, I encountered some issues: I use WES data to generate a filtered VCF file, and then utilize VEP for annotation to obtain an MAF file. The MAF file contains an excessive number…

Continue Reading [maftools]Too many multi_hit and missense mutation

AlphaMissense_hg38.tsv.gz tabix indexing issue VEP

AlphaMissense_hg38.tsv.gz tabix indexing issue VEP 1 Hi! I’m experiencing the issue with the tabix indexing of AlphaMissense_hg38.tsv.gz file. I’ve installed alphamissense plugin in VEP and have downloaded the requested files and launch the tabix command before use it. The error tells me that it can not be indexed directly because…

Continue Reading AlphaMissense_hg38.tsv.gz tabix indexing issue VEP

Appropriate genome reference for converting TCGA VCF files to MAF

Appropriate genome reference for converting TCGA VCF files to MAF 0 I have a directory of MAF files obtained from TCGA and I want to convert it to VCF format. Reference: GRCh38.d1.vd1 Reference Sequence Source: gdc.cancer.gov/about-data/gdc-data-processing/gdc-reference-files maf2vcf.pl –input-maf maf/* –output-dir VCF -ref-fasta /home/melchua/.vep/homo_sapiens/GRCh38/GRCh38.d1.vd1.fa.tar.gz Traceback: Use of uninitialized value $lines in…

Continue Reading Appropriate genome reference for converting TCGA VCF files to MAF

Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

Identification of CH cases from WGS in ISL and UKB We used WGS from 45,510 Icelanders and 130,709 British ancestry participants from the UKB17,18. Average sequencing depth was 33× for UKB and 38× for ISL. Participants with prior diagnoses of hematological disorders or grossly abnormal hematology measurements on entry were…

Continue Reading Genetics and epidemiology of mutational barcode-defined clonal hematopoiesis

In VEP annotation, how is the codon field interpreted?

In VEP annotation, how is the codon field interpreted? 0 After annotating with VEP a VCF file, we obtain different fields. One of them is called Codons which represents the affected codon in the transcript of the gene. Below is a screenshot of Insertions from a sample: HGVSp_Short RefSeq Codons…

Continue Reading In VEP annotation, how is the codon field interpreted?

vcf – VEP annotation INFO field Ensembl IDs and locations

I have a vcf file that I annoteted with VEP, for human data. I have run VEP to annotate my files with some additional parameters (as shown below in the ##VEP-command-line). However, my output is rather strange (mainly the INFO column). ##VEP=”v108″ time=”2023-04-27 15:13:08″ cache=”workflow/resources/variants/cache_vep/homo_sapiens/108_GRCh38″ ensembl-funcgen=108.56bb136 ensembl-variation=108.a885ada ensembl-io=108.58d13c1 ensembl=108.d8a9c80 1000genomes=”phase3″…

Continue Reading vcf – VEP annotation INFO field Ensembl IDs and locations

How to make 4-column annotation file

How to make 4-column annotation file 1 Hi everyone, I’ve been following regenie to perform analyses. In its 2nd step a 4-column annotation file is needed. The format is 1:55039839:T:C PCSK9 LoF 1:55039842:G:A PCSK9 missense According to snpEFF documentation, its output is in vcf format, which I think might be…

Continue Reading How to make 4-column annotation file

Locally annotating SNP IDs and Gene names of called variants

Locally annotating SNP IDs and Gene names of called variants 0 I have GWAS results after variant calling. The VCF file only had CHR (1:22) and POS (12345678 etc) information but the ID column has all “.”, namely no rsIDs in it. After GWAS analysis I have a list of…

Continue Reading Locally annotating SNP IDs and Gene names of called variants

Issues with FilterRNAMutationsNoPoN in RNA-Mutect Workflow for Tumour Mutational Burden Estimation

Issues with FilterRNAMutationsNoPoN in RNA-Mutect Workflow for Tumour Mutational Burden Estimation 0 Hello community, I am currently working on estimating tumour mutational burden using the RNA-Mutect-WMN pipeline. The process requires an output from RNA-Mutect, which in turn needs the output of Mutect2. I have proceeded with the following steps: Acquired…

Continue Reading Issues with FilterRNAMutationsNoPoN in RNA-Mutect Workflow for Tumour Mutational Burden Estimation

multiple alternate alleles for a single variant in R

Ensembl’s POST “vep/:species/region”: multiple alternate alleles for a single variant in R 1 I am working with this POST request from Ensembl. Let us suppose, as a starting example, that we have the following variant (in the form of a R code snippet): { “variants” : [“21 26960070 26960070 G…

Continue Reading multiple alternate alleles for a single variant in R

Eukaryotic Genome Annotation | Genome Annotation Pipeline

PASAPASA was originally developed at The Institute for Genomic Research in 2002 as an effort to automatically improve gene structures in Arabidopsis thaliana. Since then, it has been applied to numerous Eukaryotic genome annotation projects including Rice, Aspergillus species, Plasmodium falciparum, Schistosoma mansoni, Aedes aegypti, mouse, human, among others. Functions of PASA…

Continue Reading Eukaryotic Genome Annotation | Genome Annotation Pipeline

NGS one-liner to call variants

Tutorial:NGS one-liner to call variants 0 This is a tutorial about creating a pipeline for sequence analysis in a single line. It is made for capture/amplicon short read sequencing in mind for human DNA and tested with reference exome sequencing data described here. I share the process and debuging steps…

Continue Reading NGS one-liner to call variants

NGS oneliner

Tutorial:NGS oneliner 0 This is a tutorial about creating a pipeline for sequence analysis in a single line.I share the process and debuging steps gone through while putting it together.Source is available at: github.com/barslmn/ngsoneliner/I couldn’t make a longer post, complete version of this post: omics.sbs/blog/NGSoneliner/NGSoneliner.html Pipeline # fastp –in1 “$R1″…

Continue Reading NGS oneliner

VEP ensemble docker and plugins

VEP ensemble docker and plugins 0 Hi, I have used VEP in local using docker and cache and ran this command sudo docker run \ -v /mnt/dodl_drive/sarek_cc/vep:/data ensemblorg/ensembl-vep \ vep -i input_sample.vcf \ –cache \ –output_file output_sample.vcf \ –everything which ran successfully and the output VCF also is as expected….

Continue Reading VEP ensemble docker and plugins

Simulation of undiagnosed patients with novel genetic conditions

Simulated patient initialization We simulate patients for each of the 2134 diseases in Orphanet20 (orphadata.org, accessed October 29, 2019) that do not correspond to a group of clinically heterogeneous disorders (i.e., Orphanet’s “Category” classification31), have at least one associated phenotype, and have at least one causal gene. For Orphanet diseases…

Continue Reading Simulation of undiagnosed patients with novel genetic conditions

GenomicStar hiring Bioinformatics Data Analyst in Chennai, Tamil Nadu, India

Position: Bioinformatics Data Analyst Company: GenomicStar – A SaaS based Precision Medicine Platform Job Purpose: The Bioinformatics Data Analyst will play a critical role in the development, validation, and implementation of data analysis pipelines for various genetic test panels, incorporating state-of-the-art AI techniques. This role is integral to our commitment…

Continue Reading GenomicStar hiring Bioinformatics Data Analyst in Chennai, Tamil Nadu, India

AlphaMissense Plugin VEP

AlphaMissense Plugin VEP 0 I’ve installed alphamissense plugin in VEP, but I can’t use it. I’ve downloaded the requested files and launch the tabix command before use it. Then I’ve launched the command but I got this error: WARNING: Failed to instantiate plugin AlphaMissense: ERROR: No file specified Try using…

Continue Reading AlphaMissense Plugin VEP

Issue VEP installation MacOS

Issue VEP installation MacOS 0 Hi, I’m trying to install VEP on Mac. I’ve tried on the Anaconda Navigator, but I couldn’t install. I also tried through the terminal, but also in this way I can’t install. The last error that I’ve got is: User cram/cram_io.c:61:10: fatal error: ‘lzma.h’ file…

Continue Reading Issue VEP installation MacOS

The Biostar Herald for Tuesday, September 26, 2023

The Biostar Herald publishes user submitted links of bioinformatics relevance. It aims to provide a summary of interesting and relevant information you may have missed. You too can submit links here. This edition of the Herald was brought to you by contribution from Istvan Albert, Shred, and was edited by…

Continue Reading The Biostar Herald for Tuesday, September 26, 2023

Filtering VCF to divide with equal sizes

Filtering VCF to divide with equal sizes 1 Hello everyone! I have a very large VCF file (>400gb), and I want to divide it to use with VEP. VEP recommends separating the vcf, so I generated a list of contigs, based on the header, with 3^7 bases for each chromosome….

Continue Reading Filtering VCF to divide with equal sizes

AlphaMissense Classifies Mutation Pathogenicity | Inside Precision Medicine

Credit: iStock/agsandrew A new, freely available AI catalog has classified the potential effects of millions of missense genetic mutations, which could help establish the cause of diseases such as cystic fibrosis, sickle-cell anemia, and cancer. The AlphaMissense resource from Google DeepMind categorized 89% of all 71 million possible missense variants….

Continue Reading AlphaMissense Classifies Mutation Pathogenicity | Inside Precision Medicine

Human Genomics Team Leader at European Molecular Biology Laboratory (EMBL)

About the team/job Recent advances in technology and analysis have revealed significant opportunities to further our understanding of the diversity and impact of human genetics (such as the first full length telomere-to-telomere human genome and the first release of the human pangenome). As a result, we are excited to launch…

Continue Reading Human Genomics Team Leader at European Molecular Biology Laboratory (EMBL)

Recommended ways to merge multiple filesets with disjoint samples

Hi Chris, Thanks for the quick reply. I am getting used to eccentricities and the mentality of “right tool for right job”. BCF is highly optimized but the file floats which ends up taking entire space. My 2TiB dataset ended up taking 4.3 TiB in bcf file format while the…

Continue Reading Recommended ways to merge multiple filesets with disjoint samples

sarek: Introduction

Introduction nf-core/sarek is a workflow designed to detect variants on whole genome or targeted sequencing data. Initially designed for Human, and Mouse, it can work on any species with a reference genome. Sarek can also handle tumour / normal pairs and could include additional relapses. The pipeline is built using…

Continue Reading sarek: Introduction

Finding Unique values on specific INFO field of the VCF file (dbNSFP, vep annotated multisample VCF)

Hello everyone! I searched the forum but coundn`t find a question that is like mine I have a multisample annotated VCF file (with dbNSFP plugin) to which I have filtered using VEP, like so: /scratch/ensembl-vep-109/ensembl-vep/filter_vep –force_overwrite –input_file {1} –output_file /home/filtering/2/2_{1/.}.vcf –only_matched –filter “(clinvar_clnsig is Pathogenic) or (clinvar_clnsig is Likely_pathogenic) After…

Continue Reading Finding Unique values on specific INFO field of the VCF file (dbNSFP, vep annotated multisample VCF)

Finding Unique values on specific INFO field of the VCF file

Finding Unique values on specific INFO field of the VCF file 0 Hello everyone! I searched the forum but coundn`t find a question that is like mine I have a VCF file to which I have filtered using VEP, like so: /scratch/ensembl-vep-109/ensembl-vep/filter_vep –force_overwrite –input_file {1} –output_file /home/filtering/2/2_{1/.}.vcf –only_matched –filter “(clinvar_clnsig…

Continue Reading Finding Unique values on specific INFO field of the VCF file

Ann.vcf format

Ann.vcf format 0 Hello everyone! I used a pipeline named Sarek to run my WES data. The pipeline used two programs to annotate VEP and SnpEff and give to me an ann.vcf file. I don’t know how to use or open this file in excel for do the filters and…

Continue Reading Ann.vcf format

Installing ensemble VEP by running INSTALL.pl takes too long.

Installing ensemble VEP by running INSTALL.pl takes too long. 0 How long does it usually take to install ensemble VEP with INSTALL.pl? I am running INSTALL.pl and I think it’s taking too long. I am trying to install ensemble VEP with this code: git clone github.com/Ensembl/ensembl-vep.git cd ensembl-vep perl INSTALL.pl…

Continue Reading Installing ensemble VEP by running INSTALL.pl takes too long.

Can I annotate a 32GB vcf file containing 300 samples at once using ANNOVAR?

Can I annotate a 32GB vcf file containing 300 samples at once using ANNOVAR? 1 I have a 32 GB VCF file containing 300 samples. I also have access to a cluster computing system. Do you recommend annotating this large file all at once, or splitting it to two files…

Continue Reading Can I annotate a 32GB vcf file containing 300 samples at once using ANNOVAR?

Confusion about transcript ablation

I’m analyzing the WES data of a patient, after calling variants by GATK, I use Ensembl Variant Effect Predictor (VEP) to annotate my vcf file. Here is one record from the output file: #Uploaded_variation Location Allele Gene Feature Feature_type Consequence cDNA_position CDS_position Protein_position Amino_acids Codons Existing_variation Extra chr11_64341844_GTTGTGGTCTGAGGTCTTGGGCCATCAGTGATGTCACAACCAGATGGCCCAAGACCCCAGACCACAACCCCATGTCTGGT/- chr11:64341844-64341923- ENSG00000278359…

Continue Reading Confusion about transcript ablation

Offline Usage of the Go Plugin

Offline Usage of the Go Plugin 1 I’m using Veb in offline mode. I want to add the GO terms. From what I’ve read, it seems like I need the ‘–database’ connection. But isn’t there another way to run the GO plugin offline? ensembl veb go-terms • 229 views •…

Continue Reading Offline Usage of the Go Plugin

Is there any method for using WES data from Illumina Hiseq into ABSOLUTE?

Is there any method for using WES data from Illumina Hiseq into ABSOLUTE? 2 Hey namhaesly, You need to create segmentation files for your data. You can use VarScan2 (the copynumber) function, or any other program that computes copy number information. My tool of choice is ExomeCNV though I don’t…

Continue Reading Is there any method for using WES data from Illumina Hiseq into ABSOLUTE?

Genome-wide prediction of disease variant effects with a deep protein language model

This study did not require any ethical approval. ESM1b In this study, we have leveraged and expanded the use of ESM1b, a protein language model developed by MetaAI20. The code and pretrained parameters for ESM1b (and other ESM models) were taken from the model’s official GitHub repository at github.com/facebookresearch/esm. Throughout…

Continue Reading Genome-wide prediction of disease variant effects with a deep protein language model

VEP run with –af_gnomadg refuses to provide allele frequencies for extremely rare variants

VEP run with –af_gnomadg refuses to provide allele frequencies for extremely rare variants 0 Here are some examples. Now, I will grant you these are all extremely rare but not unknown. HGVSg gnomADg_AF actual chr20:g.63400308G>A 1.97088e-05 chr20:g.63400728C>A 1.31465e-05 chr20:g.63401127G>A 6.57263e-06 chr20:g.63402766G>A 2.6276e-05 chr20:g.63402967G>A 6.57082e-06 chr20:g.63403148G>A 1.97114e-05 chr20:g.63403351C>T 6.56927e-06 chr20:g.63403638G>A 1.31524e-05…

Continue Reading VEP run with –af_gnomadg refuses to provide allele frequencies for extremely rare variants

Problems when trying to use ensembl-vep

Problems when trying to use ensembl-vep 0 Hi. I’ve been using ensembl-vep through Docker before and managed to use it, but now I’m encountering difficulties.I’m running the commands as instructed in the ensembl-vep page.After running the commands docker pull ensemblorg/ensembl-vep and mkdir $HOME/vep_data I run docker run -t -i -v…

Continue Reading Problems when trying to use ensembl-vep

Biomart doesnt work in R for big input data. How to run it in Python ?

Biomart doesnt work in R for big input data. How to run it in Python ? 0 I am trying to use Biomart for a list of variants (with rs ids) to retrieve the consequence_types for each variant but as because my file is too big (80621 entries) i am…

Continue Reading Biomart doesnt work in R for big input data. How to run it in Python ?

Are there any ways to delete multiple jobs at Ensembl VEP?

Are there any ways to delete multiple jobs at Ensembl VEP? 0 I have looked at the VEP documentation regarding deleting the job. “You may delete a job by clicking the trash can icon”. This I have to delete one by one and I have lots of jobs to delete…

Continue Reading Are there any ways to delete multiple jobs at Ensembl VEP?

What does all this mean

Hi, I am trying to understand these Manta somatic variant VCFs that I have annotated with VEP. chr10 11974892 MantaBND:6118:0:1:0:0:0:0 G ]chr1:41850074]G . PASS SVTYPE=BND;MATEID=MantaBND:6118:0:1:0:0:0:1;CIPOS=0,2;HOMLEN=2;HOMSEQ=AG;SOMATIC;SOMATICSCORE=54;BND_DEPTH=62;MATE_BND_DEPTH=0;CHR2=chr1;ENDPOSSV=41850074;SVLEN=0;RO=-+;TRA;LENSV=-1;RGENUPS=ENSCAFG00000045430___None,-9947,+;RGENISEC=ENSCAFG00000000435___BEST3,0,-;RGENDNS=ENSCAFG00000000436___RAB3IP,59491,+;LGENUPS=ENSCAFG00000000405___PLEKHG1,-146330,+;LGENISEC=ENSCAFG00000000413___MTHFD1L,0,+;LGENDNS=ENSCAFG00000000419___AKAP12,173336,+;CSQ=]chr1:41850074]G|feature_truncation|HIGH|MTHFD1L|MTHFD1L_3|Transcript|XM_038526323.1|mRNA||||||||||||1||||Canis_familiaris.UU_Cfam_GSD_1.0.gtf.gz|,]chr1:41850074]G|feature_truncation&intron_variant|HIGH|BEST3|BEST3_2|Transcript|XM_038549568.1|mRNA||4/5||||||||||-1||||Canis_familiaris.UU_Cfam_GSD_1.0.gtf.gz| PR:SR:RDISTDISC1:RDISTDISC2:RCDIS1:RCDIS2:RUT1:RUT2:DRNOISE1:DRNOISE2:CLMQ1:CLMQ2:EVALRCNOISE 7,0:25,0:0:0:0:0:1:1:0,0,0,0,0,0:0,0,0,0,0,0:0,0,186:0,0,2:1,1 9,0:29,16:0:0:0:0:1:1:0,0,0,0,3,3:0,0,0,0,15,15:0,0,487:0,0,41:0,0 I have consulted the manual github.com/Illumina/manta/blob/master/docs/userGuide/README.md#vcf-info-fields but cannot seem to find information regarding these fields; “RGENISEC” ,”RGENDNS”, “LGENUPS”,”LGENISEC”. Additionally I see “RGENUPS”…

Continue Reading What does all this mean

Utilizing PON-P2 for vcf annotation

Utilizing PON-P2 for vcf annotation 0 I have done a lot of investigation into both the VEP plugin and the standalone api script provided for PON-P2, and I am beginning to suspect that the server is no longer functional. Has anybody used it recently or is able to confirm the…

Continue Reading Utilizing PON-P2 for vcf annotation

DECIPHER database for rare diseases joins EMBL-EBI

Closer collaborations with EMBL-EBI data resources set to further develop and grow this unique rare disease platform DECIPHER database for rare diseases joins EMBL-EBI. Credit: Karen Arnott/EMBL-EBI The DECIPHER website and database, which are used globally by researchers and clinicians to interpret and share phenotypic and genotypic data from rare…

Continue Reading DECIPHER database for rare diseases joins EMBL-EBI

SnpSift empty result

SnpSift empty result 0 Hello, I need some help with SnpSift and fitler option. I am trying to run different filters but in every case the output is empty. No warnings, of fail run. all the works are completed and the .log file does not show errors. However, the .vcf…

Continue Reading SnpSift empty result

Generation of polyphen2 & SFT scores for structural variants.

VEP: Generation of polyphen2 & SFT scores for structural variants. 0 Hello, has anybody generated polyphen 2 scores for structural variants reported by Manta using VEP. I have ran VEP but do not have any scores reported and there is no error in my out put report. Of note I…

Continue Reading Generation of polyphen2 & SFT scores for structural variants.

Why does Ensembl VEP provide HGVSg for some variants and not others?

Why does Ensembl VEP provide HGVSg for some variants and not others? 0 docker run \ -u 1000:1000 \ -v $HOME/vep_data:/data \ -v $HOME/vep_data:/.vep \ -v $PWD:/opt/foo ensemblorg/ensembl-vep \ vep \ -i /opt/foo/{input} \ -o /opt/foo/{output} \ –cache \ –force_overwrite \ –fork {threads} \ –format vcf \ –buffer_size 5000 \…

Continue Reading Why does Ensembl VEP provide HGVSg for some variants and not others?

VCF file CSQ flag

VCF file CSQ flag 0 Hello, I am new to exome sequencing. I ran GATK4 and have a g.VCF file. I used VEP to annotate the gVCF file. I have a couple of questions. I noticed not all the rows(variants) have the CSQ flag. Is that normal? Or I missed…

Continue Reading VCF file CSQ flag

variant allelic fraction

variant allelic fraction 0 Hello, I have a question, how can I find the variant allelic fraction from VEP file? thanks in advanced vep • 289 views • link updated 2 hours ago by Ram 39k • written 1 day ago by Jom • 0 Login before adding your answer….

Continue Reading variant allelic fraction

VEP tab/vcf – Different output

VEP tab/vcf – Different output 1 Hello, I am using VEP to annotate a VCF file. I am generating tab and vcf output. However, when looking at the tab output, I find mutations, that are not present in the vcf output (and also not in the input vcf file). Those…

Continue Reading VEP tab/vcf – Different output

Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine

Study design and participants The GCKD study is an ongoing prospective observational study that enrolled 5,217 adult persons with CKD between 2010 and 2012. Patients regularly seen by nephrologists with eGFR between 30 and 60 ml min−1 per 1.73 m2 or eGFR >60 ml min−1 per 1.73 m2 with UACR > 300 mg per g (or urinary protein/creatinine ratio > 500 mg…

Continue Reading Genetic studies of paired metabolomes reveal enzymatic and transport processes at the interface of plasma and urine

What does canonical transcript mean in the context of VEP

What does canonical transcript mean in the context of VEP 1 Just looking at this documentation: useast.ensembl.org/info/docs/tools/vep/script/vep_options.html and seeing: –canonical Adds a flag indicating if the transcript is the canonical transcript for the gene. Not used by default Does anyone know what this means? Does that mean the fully functional…

Continue Reading What does canonical transcript mean in the context of VEP

Why aren’t authors specific about patch version?

I am finding difficulty finding the exact assembly version (e.g. patch version) of GRCh38 used for major databases. For instance, gnomad says “GRCh38”. But the only information on the version, for v3.1 comes from here, which says it “uses an updated version of Variant Effect Predictor (VEP) based on the…

Continue Reading Why aren’t authors specific about patch version?

problems with MAF for MutSigCV (vcf2maf)

I am trying to run MutSIgCV and got stuck with this error: MutSigCV allsamples.md.tc.ir.br.pr.ug.dbsnp.vep.maf \ “$anno”exome_full192.coverage.txt \ “$anno”gene.covariates.txt \ my_results \ “$anno”mutation_type_dictionary_file.txt \ “$anno”chr_files_hg19 ====================================== MutSigCV v1.4 (c) Mike Lawrence and Gaddy Getz Broad Institute of MIT and Harvard ====================================== MutSigCV: PREPROCESS ——————– Loading mutation_file… Error using MutSigCV>MutSig_preprocess (line 246)…

Continue Reading problems with MAF for MutSigCV (vcf2maf)

Retrieving allele-specific information for a variant using VEP annotation

Hi all! I annotated a VCF file using VEP and noticed that it reports several variant IDs to each input variant. For example, this is an excerpt of one of the variant lines (I removed the annotation info that is not relevant to the question): 12 25398284 . C A…

Continue Reading Retrieving allele-specific information for a variant using VEP annotation

The impact of rare protein coding genetic variation on adult cognitive function

The UKB is approved by the North West Multi-centre Research Ethics Committee (www.ukbiobank.ac.uk/learn-more-about-uk-biobank/about-us/ethics). The current study was conducted under UKB application no. 26041. The data in the UKB were collected after written informed consent was obtained from all participants. The Human Research Committee of the MGB approved the Biobank research…

Continue Reading The impact of rare protein coding genetic variation on adult cognitive function

Splitting of VCF file of CSQ field in the INFO column to tabular format.

VCF file will be having seven fixed columns and INFO column. Chromosome, position, ID, ref, alt, qual, filter, and INFO column. This INFO column will be having the variant related information. In the INFO column CSQ field will be having multiple fields – 82 fields fixed with the delimeter “|”…

Continue Reading Splitting of VCF file of CSQ field in the INFO column to tabular format.

vep – How to normalise indel of different transcript direction?

I am trying to reproduce an analysis result, and there is a requirement regarding the format of the result: the ‘3’ rule ‘: all mutations should be expressed in a position close to the end of the gene transcription direction (3’ end). The result I obtained is(annotate from VEP) Chromosome…

Continue Reading vep – How to normalise indel of different transcript direction?

VEP/ CADD error – ERROR: Assembly is GRCh38 but CADD file does not contain GRCh38 in header.

Dear Biostars, I am having a confusing issue with my CADD plugin. This is confusing because when I run VEP for my whole trio – all the plugins work fine. However when I try to run CADD for individual – pivoted files – it no longer does and I get…

Continue Reading VEP/ CADD error – ERROR: Assembly is GRCh38 but CADD file does not contain GRCh38 in header.

Answer: How do you go about grabbing proper genotyping information from dbSNP using entr

You can check out these services: * myvariant.info/v1/api#/ * api.ncbi.nlm.nih.gov/variation/v0/ * rest.ensembl.org/#VEP As far as I know they all return JSON. You can check out a tool like jq or a package in your fav language to parse JSON. **Example:** curl -s ‘rest.ensembl.org/vep/human/id/rs2656176?’ -H ‘Content-type:application/json’ | jq ‘.[].colocated_variants[].frequencies’ { “A”:…

Continue Reading Answer: How do you go about grabbing proper genotyping information from dbSNP using entr

Senior Bioinformatics Scientist @ Natera

We are seeking a bioinformatician with significant oncology or immunology experience to join a multidisciplinary team developing leading-edge genomics analysis tools to understand the immune system’s response to cancer. This highly motivated, detail-oriented individual would join a clinical genomics analysis group and will be responsible for developing and applying bioinformatics…

Continue Reading Senior Bioinformatics Scientist @ Natera

Biobank-scale inference of ancestral recombination graphs enables genealogical analysis of complex traits

ARG-Needle and ASMC-clust algorithms We introduce two algorithms to construct the ARG of a set of samples, called ARG-Needle and ASMC-clust. Both approaches leverage output from the ASMC algorithm11, which takes as input a pair of genotyping array or sequencing samples and outputs a posterior distribution of the TMRCA across…

Continue Reading Biobank-scale inference of ancestral recombination graphs enables genealogical analysis of complex traits

Issues installing PhyloCSF for VEP: opam install ocamlbuild

I am trying to install PhyloCSF: github.com/mlin/PhyloCSF Since it is needed for new VEP99. In order to do that I needed to install ocaml, bubblewrap, opam. Second bullet point: opam install batteries ocaml+twt gsl Did not work, but I found that this one works: opam install ocaml-twt As can be…

Continue Reading Issues installing PhyloCSF for VEP: opam install ocamlbuild

How to install VEP plugins in conda environment

How to install VEP plugins in conda environment 1 I installed VEP through miniconda and use VEP through my conda environment. I would like to install VEP plugins (ex, CADD), but I do not have the INSTALL.pl script VEP instructs to use to install plugins (I think this is because…

Continue Reading How to install VEP plugins in conda environment

Haploview Error

Haploview Error 0 I would like to request for help for the following. When I am giving input of ped and info file generated from plink into haploview. I am getting an error. “File does not any zero valid indiviuals” I am thankful in advance plink Haploview vep • 36…

Continue Reading Haploview Error

adding transcripts ID to the REVEL score

adding transcripts ID to the REVEL score 0 I am annotating variants for clinical diagnosis. One of the data required by my team is REVEL. We work with specific transcripts so when clinicians do variant interpretation, they need to check if our preferred transcript is among the transcripts considered. REVEL…

Continue Reading adding transcripts ID to the REVEL score

VEP tool

VEP tool 0 I am using VEP tool for variant annotation, I ran the VEP command : ./vep -i sample.vcf -o out.txt –offline but it shows the following error : ——————– EXCEPTION ——————– MSG: ERROR: Cache directory /home/shweta/.vep/homo_sapiens not found STACK Bio::EnsEMBL::VEP::CacheDir::dir /home/shweta/ensembl-vep/modules/Bio/EnsEMBL/VEP/CacheDir.pm:305 STACK Bio::EnsEMBL::VEP::CacheDir::init /home/shweta/ensembl-vep/modules/Bio/EnsEMBL/VEP/CacheDir.pm:219 STACK Bio::EnsEMBL::VEP::CacheDir::new /home/shweta/ensembl-vep/modules/Bio/EnsEMBL/VEP/CacheDir.pm:111 STACK…

Continue Reading VEP tool

VEP-like tool for sequence ontology and HGVS annotation of VCF files

Mehari is a software package for annotating VCF files with variant effect/consequence. The program uses hgvs-rs for projecting genomic variants to transcripts and proteins and thus has high prediction quality. Other popular tools offering variant effect/consequence prediction include: Mehari offers predictions that aim to mirror VariantValidator, the gold standard for…

Continue Reading VEP-like tool for sequence ontology and HGVS annotation of VCF files

ACMG assignment tools

Hi, I’ve been asked to semi automate (as far as full automation of the process seems to be impossible) the ACMG annotations of WGS data for our clinicians. What are the options here? I’ve checked several: varsome (paid full access, closed code) InterVar (seems nice, but works only with annovar,…

Continue Reading ACMG assignment tools

What does the decimal ending of a transcript id mean?

What does the decimal ending of a transcript id mean? 1 When working with alternative transcripts, I see a decimal ending to the transcript id. What does the decimal value add to the transcript meaning? For example, a mutation, BRAF p.V600M is annotated by VEP with: csn=”ENST00000479537.1:c.82G>A_p.Val28Met” hgvsc=”ENST00000479537.1:c.82G>A” transcript_id=”ENST00000479537″ &…

Continue Reading What does the decimal ending of a transcript id mean?

VCF file generated by Gridss not recognised by VEP

VCF file generated by Gridss not recognised by VEP 0 Hello, Currently, I am engaged in analyzing genomic Next-Generation Sequencing (NGS) data as a part of my professional activities. Specifically, I am attempting to integrate the Gridss tool into my workflow for the purpose of variant calling (indels). Unfortunately, I…

Continue Reading VCF file generated by Gridss not recognised by VEP

Easy way to split VCF file by chromosome

Easy way to split VCF file by chromosome 2 Hi, Im trying to submit a job on the TOPMED/Michigan imputation server, but it returns an error saying that I need to split my VCF by chromosome Is there an easy way to do this? Will bcftools help? sequence TOPMED chromosome…

Continue Reading Easy way to split VCF file by chromosome

how to seperate VEP INFO column into seperate columns

I have a vcf files like below: #CHROM POS ID REF ALT QUAL FILTER INFO FORMAT treatmentSample chr1 857100 . C T 1756.06 PASS AC=2;AF=1;AN=2;DP=60;ExcessHet=3.0103;FS=0;MLEAC=2;MLEAF=1;MQ=60;QD=29.27;SOR=1.812;CSQ=chr1:857100|T|SNV|ENSG00000228794|ENST00000445118|LINC01128||1|MODIFIER|non_coding_transcript_exon_variant||||5/5|||||||||||||||||| GT:AD:DP:GQ:PL 1/1:0,60:60:99:1770,180,0 Does anyone know how to seperate INFO columns into different columns? And also how to separate treatmentSample column following the FORMAT ORDER? I…

Continue Reading how to seperate VEP INFO column into seperate columns

Fasta file and GTF file for STAR alignment

Fasta file and GTF file for STAR alignment 3 Hello there, The top-level fasta file will include chromsomes, regions not assembled into chromosomes and N padded haplotype/patch regions. See more here: ftp.ensembl.org/pub/release-92/fasta/mus_musculus/dna/README. If you are only looking for reference genome assembly chromosome level sequences then use the primary_assembly.fa file. The…

Continue Reading Fasta file and GTF file for STAR alignment

Variant effect predictor in ictalurus punctatus

Variant effect predictor in ictalurus punctatus 0 Hello everyone I performed a VEP analysis today using a catfish reference genome and my VCFs are from the species Pimmelodus Yuma. The problem is that the predictors did not give me any results in the missense variants. Is this normal? Here are…

Continue Reading Variant effect predictor in ictalurus punctatus

Bump hg38 VEP version to 94

Ensembl Variant Effect Predictor (VEP) VEP determines the effect of your variants (SNPs, insertions, deletions, CNVs or structural variants) on genes, transcripts, and protein sequence, as well as …… Read more > The GA4GH Variation Representation Specification (4) Context precision: VRS does not lump together or link variation contexts; a variation…

Continue Reading Bump hg38 VEP version to 94

Availability of information on genes in Gnomad VCF data

Availability of information on genes in Gnomad VCF data 1 Hi , Im new to gnomad and genetics in general and i was wondering does the gnomad genome data that is downlaoded in the vcf format on variants contains information of what is the nearest gene and is the genomic…

Continue Reading Availability of information on genes in Gnomad VCF data

BioPerl 1.7.2-GCCcore-8.2.0-Perl-5.28.1 (for BlueBEAR, BEARCloud VMs, and CaStLeS VMs)

BioPerl 1.7.2-GCCcore-8.2.0-Perl-5.28.1 There is a newer version of BioPerl Bioperl is the product of a community effort to produce Perl code which is useful in biology. Examples include Sequence objects, Alignment objects and database searching objects. Accessing BioPerl 1.7.2-GCCcore-8.2.0-Perl-5.28.1 To load the module for BioPerl 1.7.2-GCCcore-8.2.0-Perl-5.28.1 please use this command…

Continue Reading BioPerl 1.7.2-GCCcore-8.2.0-Perl-5.28.1 (for BlueBEAR, BEARCloud VMs, and CaStLeS VMs)

Installing ensembl-vep

Installing ensembl-vep 1 Hello all , I tried installing ensembl-Vep on mu ubuntu using following commands git clone github.com/Ensembl/ensembl-vep cd ensembl-vep perl INSTALL.pl I got the error as :- ERROR: DBI module not found. VEP requires the DBI perl module to function I tried installing DBI etc but the error…

Continue Reading Installing ensembl-vep

hgmd variants from Ensembl

hgmd variants from Ensembl 0 I know that Ensembl has hgmd variants amongst other annotations as given here: hgmd variants How can I download these variants from Ensembl? I downloaded VEP because it uses multiple tracks including HGMD: VEP Is it possible to use VEP to extract just the HGMD…

Continue Reading hgmd variants from Ensembl

Unexprected Ensembl-vep results

Unexprected Ensembl-vep results 0 Hi.I got a VCF from an individual that shows symptoms of a known disease with known mutations. I run it with Ensembl-vep, expecting to find some of those mutations in the results, yet, all the consequences in the results are “intergenic-variant”.The command I used was: –cache…

Continue Reading Unexprected Ensembl-vep results

python – Matching two files(vcf to maf) using a dictionaries, and appending the contents

annotation_file ##INFO=<ID=ClinVar_CLNSIG,Number=.,xxx ##INFO=<ID=ClinVar_CLNREVSTAT,Number=.,yyy ##INFO=<ID=ClinVar_CLNDN,Number=.zzz #CHROM POS ID REF ALT QUAL FILTER INFO chr1 10145 . AAC A 101.83 . AC=2;AF=0.067;AN=30;aaa chr1 10146 . AC A 98.25 . AC=2;AF=0.083;AN=24;bbb chr1 10146 . AC * 79.25 . AC=2;AF=0.083;AN=24;ccc chr1 10439 . AC A 81.33 . AC=1;AF=0.008333;AN=120;ddd chr1 10450 . T G 53.09…

Continue Reading python – Matching two files(vcf to maf) using a dictionaries, and appending the contents

08 compare visualization results of different annotation software

stay In the first two sections , We compared the differences vcf Use of annotation software , And convert the demerit recorded after the annotation into maf File format , because snpeff The comment result cannot be converted to maf, So we will compare later ANNOVAR、VEP、GATK Funcatator The results of…

Continue Reading 08 compare visualization results of different annotation software

Annotating with CADD, gnomad, Clinvar & dbNSFP on UKB RAP – Feature Requests

dint May 9, 2022, 1:33pm #1 i’m just wondering if you can specify cadd, gnomad, clinvar and dbNSFP options when annotating with hail on dxjupyterlab_spark_cluster o the UKB RAP? From the hail website, the following command can be used on your matrix file to annotate with these features: db =…

Continue Reading Annotating with CADD, gnomad, Clinvar & dbNSFP on UKB RAP – Feature Requests

Bioinformatics Analyst II – Remote in Danville, PA for Geisinger

Details Posted: 22-Apr-22 Location: Danville, Pennsylvania Type: Full Time Salary: Open Categories: Operations Job Summary Primary accountability is to leverage the organization’s data assets exome sequencing data (>180,000 individuals) from MyCode Community Health Initiative to improve quality, efficiency and generate knowledge specifically in the field of bioinformatics within health research….

Continue Reading Bioinformatics Analyst II – Remote in Danville, PA for Geisinger

wrong number of fields ?

Error occurence after merging files with bcftools: wrong number of fields ? 1 I have multiple vcf of CASES and CONTROLS variations annotated by VEP, SNPEff, SnpSift. first pair vcf -> only variations| CASES and CONTROLS second pair vcf -> variations + SnpEff | CASES and CONTROLS third pair vcf->…

Continue Reading wrong number of fields ?

split gtex genotype data by chromosomes.

Hello, I used and edited the command line to use –vcf to import vcf file. I used these commands: for chr in $(seq 1 22); do      plink –vcf /dbGAP/GTEx_Analysis_2017-06-05_v8_WholeExomeSeq_979Indiv_VEP_annot.vcf.gz            –chr $chr            –recode            –out…

Continue Reading split gtex genotype data by chromosomes.

In Vivo Methods to Assess Retinal Ganglion Cell and Optic Nerve Function and Structure in Large Animals

Here we demostrate several in vivo tests (flash visual evoked potential, pattern electroretinogram and optic coherence tomography) in goat and rhesus macaque to understand the structure and function of the optic nerve and its neurons. Hi, everyone. This is Yikui Zhang from Wenzhou Eye Hospital. The upper nerve lacks Exxons…

Continue Reading In Vivo Methods to Assess Retinal Ganglion Cell and Optic Nerve Function and Structure in Large Animals

Ensembl VEP gnomAD annotated allele frequencies different from gnomAD browser

I’ve annotated some variants using VEP, and was looking at the minor allele frequencies. Some of the variants had very different MAFs in the annotation than I expected (I expected MAF < 1%, whereas some annotated MAFs were >50%). I looked up the same variants on the gnomAD v3 browser,…

Continue Reading Ensembl VEP gnomAD annotated allele frequencies different from gnomAD browser

variant – Where should you put you cache for ensembl-vep using conda

I’ve installed vep in conda like so: conda install ensembl-vep=105.0-0 And then I installed the human cache like so: vep_install -a cf -s homo_sapiens -y GRCh38 -c /mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/refs/vep –CONVERT But when I try and run vep I get an error: vep –dir_cache /mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/refs/vep -i /mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/cohort.norm_recalibrated.vcf -o /mnt/gpfs/live/rd01__/ritd-ag-project-rd018o-mdflo13/data/test/manual/results/variants/cohort.norm_recalibrated_vep.vcf Am I doing…

Continue Reading variant – Where should you put you cache for ensembl-vep using conda

linux – How to fix Perl from anaconda not installing bioperl? Bailing out the installation for BioPerl

vep -i examples/homo_sapiens_GRCh38.vcf –database Can’t locate Bio/PrimarySeqI.pm in @INC (you may need to install the Bio::PrimarySeqI module) (@INC contains: /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0/modules /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0 /home/youssef/anaconda3/envs/ngs1/lib/site_perl/5.26.2/x86_64-linux-thread-multi /home/youssef/anaconda3/envs/ngs1/lib/site_perl/5.26.2 /home/youssef/anaconda3/envs/ngs1/lib/5.26.2/x86_64-linux-thread-multi /home/youssef/anaconda3/envs/ngs1/lib/5.26.2 .) at /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0/Bio/EnsEMBL/Slice.pm line 75. BEGIN failed–compilation aborted at /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0/Bio/EnsEMBL/Slice.pm line 75. Compilation failed in require at /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0/Bio/EnsEMBL/Feature.pm line 84. BEGIN failed–compilation aborted at /home/youssef/anaconda3/envs/ngs1/share/ensembl-vep-88.9-0/Bio/EnsEMBL/Feature.pm…

Continue Reading linux – How to fix Perl from anaconda not installing bioperl? Bailing out the installation for BioPerl

Fail to install Ensembl VEP

Fail to install Ensembl VEP 2 Hello I have been encounter some issues on installing the Ensembl VEP. I am using Ubuntu 20.04 system, with the newest perl v5.34. perl -MBio::Root::Version -e ‘print $Bio::Root::Version::VERSION,”n”‘ 1.006924 However, when I install the VEP, it failed. Hello! This installer is configured to install…

Continue Reading Fail to install Ensembl VEP

VEP issue: ERROR: Cache assembly version (GRCh37) and database or selected assembly version (GRCh38) do not match

Describe the issue VEP give errors even my query and reference has same assembly version Command :$: ./vep -i examples/homo_sapiens_GRCh37.vcf –cache –refseq cache reference details while running install.pl ? 458 NB: Remember to use –refseq when running the VEP with this cache! downloading ftp.ensembl.org/pub/release-104/variation/indexed_vep_cache/homo_sapiens_refseq_vep_104_GRCh37.tar.gz unpacking homo_sapiens_refseq_vep_104_GRCh37.tar.gz converting cache, this may…

Continue Reading VEP issue: ERROR: Cache assembly version (GRCh37) and database or selected assembly version (GRCh38) do not match

Installing Bio::DB::Big::File on Ubuntu (Building and testing…)

Installing Bio::DB::Big::File on Ubuntu (Building and testing…) 0 Hi everyone, I’m trying to install VEP, I already downloaded all the recommended PERL libraries/dependencies/modules But I still have a problem with Bio::DB::Big::File and Bio::DB::Big::AutoSQL They asked for lib big wig. I’ve found three ways to install this (in my understanding): sudo…

Continue Reading Installing Bio::DB::Big::File on Ubuntu (Building and testing…)

Failed to instantiate plugin dbNSFP in VEP

Failed to instantiate plugin dbNSFP in VEP 0 Hi Team, My VEP (version 105, installed by perl INSTALL.pl) works well. But I face some problems to use dbNSFP plugin (also installed by perl INSTALL.pl) with VEP tool. My dbNSFP version 4.2a was installed by the following code without any warning…

Continue Reading Failed to instantiate plugin dbNSFP in VEP