Question: Which keytype from match the annotations in the ensembl gtf file?
gravatar for alejandro.colaneri
2.8 years ago by
United States
alejandro.colaneri10 wrote:

Hi, I downloaded the following gtf file from ensembl: Mus_musculus.GRCm38.91.gtf

Using DESeq2 I got a res object that have ENSEMBL symbols on it. I want to use AnnotationDbi to further annotate my res objet with gene symbol and gene names. However I am not sure which kind of keytypes are in my res object. I was expecting either ENSEMBL OR ENSEMBLTRANS however when I look at the length of different object I found this:

length(keys(, keytype = "ENSEMBL"))
[1] 24420
length(keys(, keytype = "ENSEMBLTRANS"))
[1] 19245
[1] 31379     6

There are more keys in my res object than in any of the ENSEMBL related keys from the object. How that is possible? Please find a brief summary of what I did:

I used the ensembl gtf file "Mus_musculus.GRCm38.91.gtf" it to create index and to run alignments with STAR

echo "STAR --runThreadN 12 --genomeDir /genome --sjdbGTFfile /Mus_musculus.GRCm38.91.gtf --readFilesIn /RNAseq/140311_UNC14-SN744_0402_AC3LB8ACX/${i}_R1.fastq.gz /${i}_R2.fastq.gz --outFileNamePrefix /star_mapping/aligments_forSalmon/${i} --outSAMtype BAM Unsorted --readFilesCommand zcat" >> $i.sub

Then I created a summarize experiment object using genomicAligments

se <- summarizeOverlaps(exonsByGene, BamFileList(bamFiles), mode="Union", singleEnd=FALSE, ignore.strand=FALSE, fragments=FALSE)

Later on I run DESeq2 to perform differential expression analysis and I got the res objetc Now I want to add genesymbol, For that I am extracting the keys in my result table with myKeys=rownames(res), however I m having problem to find out which kind of keytypes those are

Hope someone can help me Thanks ALe rna-seq deseq2 • 1.9k views
ADD COMMENTlink written 2.8 years ago by alejandro.colaneri10
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 2.3.0
Traffic: 1084 users visited in the last hour