Entering edit mode
Chelsea • 0@8759ce1a
Last seen 4 weeks ago
I am using the DESeq pipeline to create a heatmap of RNA transcripts deferentially expressed from 3 treatments, I have created a heatmap with the geneIDs but I am struggling to find a way to add gene names or alternatively produce a list of gene names that are deferentially expressed. I have a GFF3 file and a cds file for the genome I am using. I have tried the GenomicFeatures package but I cant access the vignettes for this package. Is there an updated package?
hse <- makeTranscriptDbFromGFF( "Red5_v1.0.gff3.gz", format="gff" ) Error in makeTranscriptDbFromGFF("Red5_v1.0.gff3.gz", format = "gff") : could not find function "makeTranscriptDbFromGFF"
Hi, if it doesn't have to be GenomicFeature package, then there are plenty of other options. I'm usually using biomaRt R package.
In this case, rownames of my results were ensembl_gene_id, but you can of course adjust it to the values you're having. Then if you need you can replace the rownames with your external_gene_name and delete the additional column, that was created after merging.
Of course, it doesn't fit if your target genome is not in the package, and it's the reason you're creating custom database.