This is a test version of Biostars. For the public version, visit https://www.biostars.org.
retrieving attributes of a list generated by read.fasta()

Hi All, I am new to R and I struggle with the following: I generated a list with a MusTr <- read.fasta(file = "Mus_musculus.GRCm38.cdna.all.fa", as.string = TRUE) function, which looks like this

$ENSMUST00000177564.1
[1] "atcggagggatacgag"
attr(,"name")
[1] "ENSMUST00000177564.1"
attr(,"Annot")
[1] ">ENSMUST00000177564.1 cdna chromosome:GRCm38:14:54122226:54122241:1 gene:ENSMUSG00000096176.1 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene gene_symbol:Trdd2 description:T cell receptor delta diversity 2 [Source:MGI Symbol;Acc:MGI:4439546]"
attr(,"class")
[1] "SeqFastadna"

$ENSMUST00000196221.1
[1] "atggcatat"
attr(,"name")
[1] "ENSMUST00000196221.1"
attr(,"Annot")
[1] ">ENSMUST00000196221.1 cdna chromosome:GRCm38:14:54113468:54113476:1 gene:ENSMUSG00000096749.2 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene gene_symbol:Trdd1 description:T cell receptor delta diversity 1 [Source:MGI Symbol;Acc:MGI:4439547]"
attr(,"class")
[1] "SeqFastadna"

$ENSMUST00000179664.1
[1] "atggcatatca"
attr(,"name")
[1] "ENSMUST00000179664.1"
attr(,"Annot")
[1] ">ENSMUST00000179664.1 cdna chromosome:GRCm38:14:54113468:54113478:1 gene:ENSMUSG00000096749.2 gene_biotype:TR_D_gene transcript_biotype:processed_transcript gene_symbol:Trdd1 description:T cell receptor delta diversity 1 [Source:MGI Symbol;Acc:MGI:4439547]"
attr(,"class")
[1] "SeqFastadna"

$ENSMUST00000178537.1
[1] "gggacagggggc"
attr(,"name")
[1] "ENSMUST00000178537.1"
attr(,"Annot")
[1] ">ENSMUST00000178537.1 cdna chromosome:GRCm38:6:41533201:41533212:1 gene:ENSMUSG00000095668.1 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene gene_symbol:Trbd1 description:T cell receptor beta, D region 1 [Source:MGI Symbol;Acc:MGI:4439571]"
attr(,"class")
[1] "SeqFastadna"

I would like to retrieve the attributes of each element (e.g. as a vector). For a single element of the list attr works (here to retrieve "Annot" attribute):

attr(MusTr$ENSMUST00000196221.1, "Annot", exact = FALSE)
[1] ">ENSMUST00000196221.1 cdna chromosome:GRCm38:14:54113468:54113476:1 gene:ENSMUSG00000096749.2 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene gene_symbol:Trdd1 description:T cell receptor delta diversity 1 [Source:MGI Symbol;Acc:MGI:4439547]"

How do I achieve the same for multiple/all elements of the list? Thanks in advance for any suggestions.

rna-seq r
• 4,449 views • link
updated by cpad0112 cpad0112 • written by lech.kaczmarczyk • lech.kaczmarczyk

The Ensembl perl API is the tool to use for this kind of job. Any particular reason you need to use R for this ? If you have to, check the biomaRt package or the mygene package. They would give you access to some annotations but are not as flexible or comprehensive as the Ensembl API.

• 0 views • link
written by Jean-Karim Heriche Jean-Karim Heriche

I would (potentially) find R helpful in manipulating FASTA files, through converting them to data frames, and then easy exporting with write.fasta. Otherwise, no particular reason. Thanks for your suggestion.

• 0 views • link
written by lech.kaczmarczyk • lech.kaczmarczyk

1 answer

little bit late. Try the following in R: I guess you are using seqinR for reading fasta file.

$ lapply(MusTr, function(x) attr(x, "Annot", exact = FALSE)) or $ getAnnot(MusTr)

• 0 views • link
written by cpad0112 cpad0112

Log in to answer this question.

More posts like this
  • Editing header of a fasta file
    written by Melisa • Melisa
    Hello everybody, I've been using sed but for simple steps and now I can't do this: I have this header: >ENSP00000451042.1 pep chromosome:GRCh38:14:22438547:22438554:1 gene:ENSG00000223997.1 transcript:ENST00000415118.1 …
  • Help with finding introns
    written by elisheva elisheva
    Hello everyone!! I have a question that I guess the answer is simple: I have a fasta file of human cDNA that looks like this: …
  • Extract Information from fasta File
    written by saamar.rajput saamar.rajput
    I have big fasta file with information like this, >ENST00000448914.1 cdna chromosome:GRCh38:14:22449113:22449125:1 gene:ENSG00000228985.1 gene_biotype:TR_D_gene transcript_biotype:TR_D_gene gene_symbol:TRDD3 description:T cell receptor delta diversity 3 [Source:HGNC Symbol;Acc:HGNC:12256] ACTGGGGGATACG …