Batch Effects in Single-cell RNA-sequencing Data Are Corrected by Matching Mutual Nearest Neighbors

Overview

Journal Nat Biotechnol

Specialty Biotechnology

Date 2018 Apr 3

PMID 29608177

Citations 977

Authors

Laleh Haghverdi

Aaron T L Lun

Michael D Morgan

John C Marioni

Affiliations

Soon will be listed here.

Abstract

Large-scale single-cell RNA sequencing (scRNA-seq) data sets that are produced in different laboratories and at different times contain batch effects that may compromise the integration and interpretation of the data. Existing scRNA-seq analysis methods incorrectly assume that the composition of cell populations is either known or identical across batches. We present a strategy for batch correction based on the detection of mutual nearest neighbors (MNNs) in the high-dimensional expression space. Our approach does not rely on predefined or equal population compositions across batches; instead, it requires only that a subset of the population be shared between batches. We demonstrate the superiority of our approach compared with existing methods by using both simulated and real scRNA-seq data sets. Using multiple droplet-based scRNA-seq data sets, we demonstrate that our MNN batch-effect-correction method can be scaled to large numbers of cells.

Citing Articles

A single-cell transcriptomic atlas reveals the cell differentiation trajectory and the response to virus invasion in swelling clove of garlic.

Gao S, Li F, Zeng Z, He Q, Mostafa H, Zhang S Hortic Res. 2025; 12(4):uhae365.

PMID: 40070403 PMC: 11894531. DOI: 10.1093/hr/uhae365.

An model for cardiac organoid production: The combined role of geometrical confinement and substrate stiffness.

Santoro R, Piacentini L, Vavassori C, Benzoni P, Colombo G, Banfi C Mater Today Bio. 2025; 31:101566.

PMID: 40061214 PMC: 11889630. DOI: 10.1016/j.mtbio.2025.101566.

Mosquito Cell Atlas: A single-nucleus transcriptomic atlas of the adult mosquito.

Goldman O, DeFoe A, Qi Y, Jiao Y, Weng S, Houri-Zeevi L bioRxiv. 2025; .

PMID: 40060408 PMC: 11888250. DOI: 10.1101/2025.02.25.639765.

MAEST: accurately spatial domain detection in spatial transcriptomics with graph masked autoencoder.

Zhu P, Shu H, Wang Y, Wang X, Zhao Y, Hu J Brief Bioinform. 2025; 26(2).

PMID: 40052440 PMC: 11886571. DOI: 10.1093/bib/bbaf086.

Focal adhesion in the tumour metastasis: from molecular mechanisms to therapeutic targets.

Liu Z, Zhang X, Ben T, Li M, Jin Y, Wang T Biomark Res. 2025; 13(1):38.

PMID: 40045379 PMC: 11884212. DOI: 10.1186/s40364-025-00745-7.

References

Segerstolpe A, Palasantza A, Eliasson P, Andersson E, Andreasson A, Sun X . Single-Cell Transcriptome Profiling of Human Pancreatic Islets in Health and Type 2 Diabetes. Cell Metab. 2016; 24(4):593-607. PMC: 5069352. DOI: 10.1016/j.cmet.2016.08.020. View

Tung P, Blischak J, Hsiao C, Knowles D, Burnett J, Pritchard J . Batch effects and the effective design of single-cell gene expression studies. Sci Rep. 2017; 7:39921. PMC: 5206706. DOI: 10.1038/srep39921. View

Bendall S, Davis K, Amir E, Tadmor M, Simonds E, Chen T . Single-cell trajectory detection uncovers progression and regulatory coordination in human B cell development. Cell. 2014; 157(3):714-25. PMC: 4045247. DOI: 10.1016/j.cell.2014.04.005. View

Angerer P, Haghverdi L, Buttner M, Theis F, Marr C, Buettner F . destiny: diffusion maps for large-scale single-cell data in R. Bioinformatics. 2015; 32(8):1241-3. DOI: 10.1093/bioinformatics/btv715. View

Scialdone A, Tanaka Y, Jawaid W, Moignard V, Wilson N, Macaulay I . Resolving early mesoderm diversification through single-cell expression profiling. Nature. 2016; 535(7611):289-293. PMC: 4947525. DOI: 10.1038/nature18633. View

Lawlor N, George J, Bolisetty M, Kursawe R, Sun L, Sivakamasundari V . Single-cell transcriptomes identify human islet cell signatures and reveal cell-type-specific expression changes in type 2 diabetes. Genome Res. 2016; 27(2):208-222. PMC: 5287227. DOI: 10.1101/gr.212720.116. View

Risso D, Ngai J, Speed T, Dudoit S . Normalization of RNA-seq data using factor analysis of control genes or samples. Nat Biotechnol. 2014; 32(9):896-902. PMC: 4404308. DOI: 10.1038/nbt.2931. View

Jaitin D, Kenigsberg E, Keren-Shaul H, Elefant N, Paul F, Zaretsky I . Massively parallel single-cell RNA-seq for marker-free decomposition of tissues into cell types. Science. 2014; 343(6172):776-9. PMC: 4412462. DOI: 10.1126/science.1247651. View

Dobin A, Davis C, Schlesinger F, Drenkow J, Zaleski C, Jha S . STAR: ultrafast universal RNA-seq aligner. Bioinformatics. 2012; 29(1):15-21. PMC: 3530905. DOI: 10.1093/bioinformatics/bts635. View

10.

Grun D, Muraro M, Boisset J, Wiebrands K, Lyubimova A, Dharmadhikari G . De Novo Prediction of Stem Cell Identity using Single-Cell Transcriptome Data. Cell Stem Cell. 2016; 19(2):266-277. PMC: 4985539. DOI: 10.1016/j.stem.2016.05.010. View

11.

Macosko E, Basu A, Satija R, Nemesh J, Shekhar K, Goldman M . Highly Parallel Genome-wide Expression Profiling of Individual Cells Using Nanoliter Droplets. Cell. 2015; 161(5):1202-1214. PMC: 4481139. DOI: 10.1016/j.cell.2015.05.002. View

12.

Brandani G, Schor M, MacPhee C, Grubmuller H, Zachariae U, Marenduzzo D . Quantifying disorder through conditional entropy: an application to fluid mixing. PLoS One. 2013; 8(6):e65617. PMC: 3677935. DOI: 10.1371/journal.pone.0065617. View

13.

Picelli S, Bjorklund A, Faridani O, Sagasser S, Winberg G, Sandberg R . Smart-seq2 for sensitive full-length transcriptome profiling in single cells. Nat Methods. 2013; 10(11):1096-8. DOI: 10.1038/nmeth.2639. View

14.

Nestorowa S, Hamey F, Pijuan Sala B, Diamanti E, Shepherd M, Laurenti E . A single-cell resolution map of mouse hematopoietic stem and progenitor cell differentiation. Blood. 2016; 128(8):e20-31. PMC: 5305050. DOI: 10.1182/blood-2016-05-716480. View

15.

Xu C, Su Z . Identification of cell types from single-cell transcriptomes using a novel clustering method. Bioinformatics. 2015; 31(12):1974-80. PMC: 6280782. DOI: 10.1093/bioinformatics/btv088. View

16.

Leek J . svaseq: removing batch effects and other unwanted noise from sequencing data. Nucleic Acids Res. 2014; 42(21). PMC: 4245966. DOI: 10.1093/nar/gku864. View

17.

Ritchie M, Phipson B, Wu D, Hu Y, Law C, Shi W . limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res. 2015; 43(7):e47. PMC: 4402510. DOI: 10.1093/nar/gkv007. View

18.

Johnson W, Li C, Rabinovic A . Adjusting batch effects in microarray expression data using empirical Bayes methods. Biostatistics. 2006; 8(1):118-27. DOI: 10.1093/biostatistics/kxj037. View

19.

Zheng G, Terry J, Belgrader P, Ryvkin P, Bent Z, Wilson R . Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017; 8:14049. PMC: 5241818. DOI: 10.1038/ncomms14049. View

20.

Brennecke P, Anders S, Kim J, Kolodziejczyk A, Zhang X, Proserpio V . Accounting for technical noise in single-cell RNA-seq experiments. Nat Methods. 2013; 10(11):1093-5. DOI: 10.1038/nmeth.2645. View