ComBat-seq: batch effect adjustment for RNA-seq count data

Abstract

The benefit of integrating batches of genomic data to increase statistical power is often hindered by batch effects, or unwanted variation in data caused by differences in technical factors across batches. It is therefore critical to effectively address batch effects in genomic data to overcome these challenges. Many existing methods for batch effects adjustment assume the data follow a continuous, bell-shaped Gaussian distribution. However in RNA-seq studies the data are typically skewed, over-dispersed counts, so this assumption is not appropriate and may lead to erroneous results. Negative binomial regression models have been used previously to better capture the properties of counts. We developed a batch correction method, ComBat-seq, using a negative binomial regression model that retains the integer nature of count data in RNA-seq studies, making the batch adjusted data compatible with common differential expression software packages that require integer counts. We show in realistic simulations that the ComBat-seq adjusted data results in better statistical power and control of false positives in differential expression compared to data adjusted by the other available methods. We further demonstrated in a real data example that ComBat-seq successfully removes batch effects and recovers the biological signal in the data.

Cited by papers we hold 11

• CCDC138 overexpression predicts poor prognosis and highlights ciliopathy-linked mechanisms in uterine corpus endometrial carcinoma · 2025
  Not checked
• TLR4+group 2 innate lymphoid cells contribute to persistent type 2 immunity in airway diseases · 2025
  Not checked
• Mitochondrial integrated stress response controls lung epithelial cell fate · 2023
  No notice
• Exposure of iPSC-derived human microglia to brain substrates enables the generation and manipulation of diverse transcriptional states in vitro · 2023
  No notice
• Genetic variation in the immunoglobulin heavy chain locus shapes the human antibody repertoire · 2023
  No notice
• pyComBat, a Python tool for batch effects correction in high-throughput molecular data using empirical Bayes methods · 2023
  No notice
• Machine learning for multi-omics data integration in cancer · 2022
  Not checked
• Monocyte Transcriptional Responses to Mycobacterium tuberculosis Associate with Resistance to Tuberculin Skin Test and Interferon Gamma Release Assay Conversion · 2022
  Not checked
• Longitudinal single-cell RNA-seq analysis reveals stress-promoted chemoresistance in metastatic ovarian cancer · 2022
  No notice
• Flimma: a federated and privacy-aware tool for differential gene expression analysis · 2021
  Not checked
• Cobolt: integrative analysis of multimodal single-cell sequencing data · 2021
  No notice

Reads

A person counts once a day per paper. We keep no cookie and no address: a reader is a hash that changes every day. Any month under five reads shows no number.