maftools: Error in validateMaf(maf = maf, isTCGA = isTCGA, rdup = removeDuplicatedVariants, : missing required fields from MAF: Chromosome
3
0
Entering edit mode
Ramiro Magno ▴ 100
@ramiro-magno-12376
Last seen 5.5 years ago
CBMR, Faro, Portugal

Hi All,

I have downloaded a Maf file using TCGAbiolinks (version 2.5.9):

maf <- GDCquery_Maf("BRCA", pipelines = "muse", save.csv = TRUE)

Then I tried to load it as a Maf object using the maftools package (version 1.2.30):

BRCA_maf_csv <- "GDCdata/TCGA.BRCA.muse.d9876b23-3e7d-4d7b-bc1b-3b4393cd2afb.DR-7.0.somatic.maf.csv"
tcga_BRCA_maf <- read.maf(maf = BRCA_maf_csv, removeSilent = TRUE, useAll = FALSE)

but I get this unexpected error: 

Error in validateMaf(maf = maf, isTCGA = isTCGA, rdup = removeDuplicatedVariants,  :
  missing required fields from MAF: Chromosome

This is unexpected because

head -n1 TCGA.BRCA.muse.d9876b23-3e7d-4d7b-bc1b-3b4393cd2afb.DR-7.0.somatic.maf.csv

does show that the Chromosome field is present: 

Hugo_Symbol,Entrez_Gene_Id,Center,NCBI_Build,Chromosome,Start_Position,End_Position,Strand,Variant_Classification,Variant_Type,Reference_Allele,Tumor_Seq_Allele1,Tumor_Seq_Allele2,dbSNP_RS,dbSNP_Val_Status,Tumor_Sample_Barcode,Matched_Norm_Sample_Barcode,Match_Norm_Seq_Allele1,Match_Norm_Seq_Allele2,Tumor_Validation_Allele1,Tumor_Validation_Allele2,Match_Norm_Validation_Allele1,Match_Norm_Validation_Allele2,Verification_Status,Validation_Status,Mutation_Status,Sequencing_Phase,Sequence_Source,Validation_Method,Score,BAM_File,Sequencer,Tumor_Sample_UUID,Matched_Norm_Sample_UUID,HGVSc,HGVSp,HGVSp_Short,Transcript_ID,Exon_Number,t_depth,t_ref_count,t_alt_count,n_depth,n_ref_count,n_alt_count,all_effects,Allele,Gene,Feature,Feature_type,One_Consequence,Consequence,cDNA_position,CDS_position,Protein_position,Amino_acids,Codons,Existing_variation,ALLELE_NUM,DISTANCE,TRANSCRIPT_STRAND,SYMBOL,SYMBOL_SOURCE,HGNC_ID,BIOTYPE,CANONICAL,CCDS,ENSP,SWISSPROT,TREMBL,UNIPARC,RefSeq,SIFT,PolyPhen,EXON,INTRON,DOMAINS,GMAF,AFR_MAF,AMR_MAF,ASN_MAF,EAS_MAF,EUR_MAF,SAS_MAF,AA_MAF,EA_MAF,CLIN_SIG,SOMATIC,PUBMED,MOTIF_NAME,MOTIF_POS,HIGH_INF_POS,MOTIF_SCORE_CHANGE,IMPACT,PICK,VARIANT_CLASS,TSL,HGVS_OFFSET,PHENO,MINIMISED,ExAC_AF,ExAC_AF_Adj,ExAC_AF_AFR,ExAC_AF_AMR,ExAC_AF_EAS,ExAC_AF_FIN,ExAC_AF_NFE,ExAC_AF_OTH,ExAC_AF_SAS,GENE_PHENO,FILTER,CONTEXT,src_vcf_id,tumor_bam_uuid,normal_bam_uuid,case_id,GDC_FILTER,COSMIC,MC3_Overlap,GDC_Validation_Status
software error maftools tcgabiolinks • 3.6k views
ADD COMMENT
1
Entering edit mode
Ramiro Magno ▴ 100
@ramiro-magno-12376
Last seen 5.5 years ago
CBMR, Faro, Portugal

My apologies...

read.maf expects a tab-delimited file and the downloaded MAF file with TCGAbiolinks is CSV.

ADD COMMENT
2
Entering edit mode
snijesh ▴ 200
@snijesh-20358
Last seen 8 months ago
India

You can try: using read.delim or fread.

my_maf = data.table::fread("path/to/maf/sample.maf")
my_maf = read.maf(maf = my_maf)
ADD COMMENT
0
Entering edit mode
@tiagochst-7121
Last seen 18 months ago
Miami, US

 Instead of the file you can pass the object.

maf <- GDCquery_Maf("BRCA", pipelines = "muse", save.csv = TRUE)
tcga_BRCA_maf <- read.maf(maf = maf, removeSilent = TRUE, useAll = FALSE)

ADD COMMENT

Login before adding your answer.

Traffic: 458 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6