Backport PR #3844 on branch 1.11.x (docs: fix score_genes citation) (#3850)

Co-authored-by: MLubetzki <[email protected]>

Author: meeseeksmachine

docs/references.bib

@@ -869,6 +869,17 @@ @article{Tarashansky2019
   month     = {sep},
 }
 
+@article{Tirosh2016,
+  author    = {Tirosh, Itay and Izar, Benjamin and Prakadan, Sanjay M and Wadsworth, Marc H and Treacy, Daniel and Trombetta, John J and Rotem, Asaf and Rodman, Christopher and Lian, Christine and Murphy, George and others},
+  title     = {Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq},
+  journal   = {Science},
+  volume    = {352},
+  number    = {6282},
+  pages     = {189--196},
+  year      = {2016},
+  publisher = {American Association for the Advancement of Science},
+}
+
 @misc{Traag2017,
   author    = {Traag, Vincent},
   doi       = {10.5281/ZENODO.35117},

src/scanpy/tools/_score_genes.py

@@ -69,13 +69,14 @@ def score_genes(  # noqa: PLR0913
     use_raw: bool | None = None,
     layer: str | None = None,
 ) -> AnnData | None:
-    """Score a set of genes :cite:p:`Satija2015`.
+    """Score a set of genes :cite:p:`Tirosh2016`.
 
     The score is the average expression of a set of genes after subtraction by
     the average expression of a reference set of genes. The reference set is
     randomly sampled from the `gene_pool` for each binned expression value.
 
-    This reproduces the approach in Seurat :cite:p:`Satija2015` and has been implemented
+    This reproduces the approach in Seurat :cite:p:`Tirosh2016` ("MITF and AXL expression
+    programs and cell scores" in materials and methods) and has been implemented
     for Scanpy by Davide Cittaro.
 
     Parameters