Hello everyone,

I am currently in an internship to see differential expression between 2 treatments. I managed to find the DE genes and I would like to annotate it but I am working on a "non model" species (absent from geneid in goana). So I have my a dataframe with GO terms associated with genes ("annotation") and one DGELRT including 8 genes wich are DE ("tr"). I extracted DE genes to associate it with the GO terms as follows :

tr <- glmTreat(fit, coef = 4, lfc=log2(1.5))

top_genes <- topTags(qlf)[1:8,]

top_genes <- rownames(top_genes)
top_genes_couleur_id <- as.matrix(annotation[annotation$V1 %in% top_genes_couleur,])

But I have to admit that even by browing tutorials with several packages, I can't annotate the DE genes ....

PS : some genes have several GO terms

Thank you for reading, and have a great day !