Hi everyone,

I'm currently working with a dataset that contains Ensembl gene IDs and Ensembl transcript IDs, and I need to map these to RefSeq mRNA IDs and RefSeq Gene IDs. Here's a snippet of my dataset:

> head(isoform.ensembl)

         gene_stable  transcript_stable
1 ENSMPUG00000000002 ENSMPUT00000000002
2 ENSMPUG00000000003 ENSMPUT00000000003
3 ENSMPUG00000000004 ENSMPUT00000000004
4 ENSMPUG00000000005 ENSMPUT00000000005
5 ENSMPUG00000000006 ENSMPUT00000000006
6 ENSMPUG00000000007 ENSMPUT00000000007

I’ve tried using the getBM() function with Ensembl gene and transcript IDs as filters and retrieving the corresponding RefSeq mRNA and RefSeq transcript. However, I keep getting empty results, even for human data. For example, I focus on ferret first:

# Connect to the Ensembl BioMart database for ferret
ensembl.mu <- biomaRt::useMart("ensembl", dataset = "mpfuro_gene_ensembl")
#
getBM(
    attributes = c("ensembl_gene_id", "ensembl_transcript_id", 
                   "refseq_mrna", "refseq_peptide"),
    filters = "ensembl_transcript_id",
    values = isoform.ensembl$transcript_stable,
    mart = ensembl.mu
)

The code works, however, I just got empty results:

ensembl_gene_id ensembl_transcript_id refseq_mrna refseq_gene
1   ENSMPUG00000000002    ENSMPUT00000000002            
2   ENSMPUG00000000003    ENSMPUT00000000003            
3   ENSMPUG00000000004    ENSMPUT00000000004            
4   ENSMPUG00000000005    ENSMPUT00000000005            
5   ENSMPUG00000000006    ENSMPUT00000000006            
6   ENSMPUG00000000007    ENSMPUT00000000007            
7   ENSMPUG00000000008    ENSMPUT00000000008

Even If I tried using humans, I still got empty columns.

I know that there is a lot of different post related to the same idea, however, not of them worked so far.

Please, any comment to solve this issue will be really appreciated!