Marker genes for each significant cluster were found using the function sc.tl.rank_genes_groups with default parameters.