HESTData: provide util to map ensemble ID to gene name

[konst-int-i]

This PR

• Provides util to map genes encoded using the ensemble ID (common for several Spatial Transcriptomics samples) to the gene name
• Addresses Mapping gene names to Ensembl IDs #70

To dos

• Might want to make this bi-directional to have ensembleID_to_gene and gene_to_ensembleID
• There is some redundancy with get_gene_db() which is not needed when using the use_cache when we query the biomart annotations. Are we good to remove this function @pauldoucet as I don't see it used anywhere?

Run instructions

from hest import iter_hest, ensembleID_to_gene

# three samples with ensemblIDs as var_names
id_list = ['SPA118', 'SPA117', 'SPA116']

for st in iter_hest('/home/iain/kh/ssd/hest_data/', id_list=id_list):
    
    print(any([var_name.startswith("ENSG") for var_name in st.adata.var_names]))
    print(st.adata.var_names[:5])
    
    st_updated = ensembleID_to_gene(st)
    
    print(any([var_name.startswith("ENSG") for var_name in st_updated.adata.var_names]))
    print(st_updated.adata.var_names[:5])

Expected output:

True
Index(['ENSG00000000003', 'ENSG00000000005', 'ENSG00000000419',
       'ENSG00000000457', 'ENSG00000000460'],
      dtype='object')
False
Index(['TSPAN6', 'TNMD', 'DPM1', 'SCYL3', 'FIRRM'], dtype='object', name='gene_name')
True
Index(['ENSG00000000003', 'ENSG00000000005', 'ENSG00000000419',
       'ENSG00000000457', 'ENSG00000000460'],
      dtype='object')
False
Index(['TSPAN6', 'TNMD', 'DPM1', 'SCYL3', 'FIRRM'], dtype='object', name='gene_name')
True
Index(['ENSG00000000003', 'ENSG00000000005', 'ENSG00000000419',
       'ENSG00000000457', 'ENSG00000000460'],
      dtype='object')
False
Index(['TSPAN6', 'TNMD', 'DPM1', 'SCYL3', 'FIRRM'], dtype='object', name='gene_name')

[pauldoucet]

Hi @konst-int-i,
Have you tried on TENX24? It gives me 0 valid genes, weird

[konst-int-i]

Hi @konst-int-i, Have you tried on TENX24? It gives me 0 valid genes, weird

Yes, that's because TENX24 doesn't contain any ensemble IDs and currently the default behavior was to invalidate genes without a mapping. Changed it to keep genes without mapping in 4e369f9

[pauldoucet]

looks good to me, thanks for the addition!