import scanpy as sc genes = ['brca1', 'brca2'] sc.tl.pca(adata, svd_solver='arpack') sc.pl.pca(adata, color=genes)