Best correlation between a given gene and the rest of genes

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?
Thanks for your help!