Abstract:ABSTRACT Functional genomics techniques, such as transposon insertion sequencing and RNA-sequencing, are key to studying relative differences in bacterial mutant fitness or gene expression under selective conditions. However, certain stress conditions, mutations, or antibiotics can directly interfere with DNA synthesis, resulting in systematic changes in local DNA copy numbers along the chromosome. This can lead to artifacts in sequencing-based functional genomics data when comparing antibiotic treatment to an unstressed control. Further, relative differences in gene-wise read counts may result from alterations in chromosomal replication dynamics, rather than selection or direct gene regulation. We term this artifact “chromosomal location bias” and implement a principled statistical approach to correct it by calculating local normalization factors along the chromosome. These normalization factors are then directly incorporated into statistical analyses using standard RNA-sequencing analysis methods without modifying the read counts themselves, preserving important information about the mean-variance relationship in the data. We illustrate the utility of this approach by generating and analyzing a ciprofloxacin-treated transposon insertion sequencing data set in Escherichia coli as a case study. We show that ciprofloxacin treatment generates chromosomal location bias in the resulting data, and we further demonstrate that failing to correct for this bias leads to false predictions of mutant drug sensitivity as measured by minimum inhibitory concentrations. We have developed an R package and user-friendly graphical Shiny application, ChromoCorrect, that detects and corrects for chromosomal bias in read count data, enabling the application of functional genomics technologies to the study of antibiotic stress. IMPORTANCE Altered gene dosage due to changes in DNA replication has been observed under a variety of stresses with a variety of experimental techniques. However, the implications of changes in gene dosage for sequencing-based functional genomics assays are rarely considered. We present a statistically principled approach to correcting for the effect of changes in gene dosage, enabling testing for differences in the fitness effects or regulation of individual genes in the presence of confounding differences in DNA copy number. We show that failing to correct for these effects can lead to incorrect predictions of resistance phenotype when applying functional genomics assays to investigate antibiotic stress, and we provide a user-friendly application to detect and correct for changes in DNA copy number.

Cgcorrect: a Method to Correct for Confounding Cell-Cell Variation Due to Cell Growth in Single-Cell Transcriptomics

Parameterizing Cell-to-cell Regulatory Heterogeneities Via Stochastic Transcriptional Profiles

Assessment of batch-correction methods for scRNA-seq data with a new test metric

A test metric for assessing single-cell RNA-seq batch correction

Data Exploration, Quality Control and Testing in Single-Cell qPCR-Based Gene Expression Experiments

Multi-Experiment Nonlinear Mixed Effect Modeling of Single-Cell Translation Kinetics after Transfection

Evaluating batch correction methods for image-based cell profiling

A method to correct for local alterations in DNA copy number that bias functional genomics assays applied to antibiotic-treated bacteria

A benchmark of computational methods for correcting biases of established and unknown origin in CRISPR-Cas9 screening data

iCellR: Combined Coverage Correction and Principal Component Alignment for Batch Alignment in Single-Cell Sequencing Analysis

Single-Cell Transcriptome Profiling Simulation Reveals the Impact of Sequencing Parameters and Algorithms on Clustering

Single Cell Genomics: when Stochasticity Meets Precision

Differences in molecular sampling and data processing explain variation among single-cell and single-nucleus RNA-seq experiments

Gene count normalization in single-cell imaging-based spatially resolved transcriptomics

: estimating the cell composition of heterogeneous tissue with varying cell sizes using gene expression

GTestimate: Improving relative gene expression estimation in scRNA-seq using the Good-Turing estimator

Normalizing need not be the norm: count-based math for analyzing single-cell data

Cell Cycle and Cell Size Dependent Gene Expression Reveals Distinct Subpopulations at Single-Cell Level

Accurate Estimation of Cell-Type Composition from Gene Expression Data

scCDC: a computational method for gene-specific contamination detection and correction in single-cell and single-nucleus RNA-seq data

Challenges and opportunities to computationally deconvolve heterogeneous tissue with varying cell sizes using single cell RNA-sequencing datasets