Best correlation between the expression of a given gene and the rest of genes in scRNAseq data

Question

[mmalumbresm]

I am analyzing single cell RNA seq data with scanpy and I would like to find a method to find which gene(s) display the best correlation in their expression profile when compared to a gene of my choice. Any suggestion? Python-based would be great but any suggestion welcome! Thanks for your help!